46{
50
51 const std::string monte_carlo_pi_program =
52 " var samples[2 * 10^8] := [(rnd_01^2 + rnd_01^2) <= 1]; "
53 " 4 * sum(samples) / samples[]; ";
54
56
57 symbol_table_t symbol_table;
58 symbol_table.add_function("rnd_01",rnd01);
59
60 expression_t expression;
61 expression.register_symbol_table(symbol_table);
62
63 parser_t parser;
64 parser.compile(monte_carlo_pi_program,expression);
65
66 const T approximate_pi = expression.value();
67
68 const T real_pi = T(3.141592653589793238462643383279502);
69
70 printf("pi ~ %20.17f\terror: %20.17f\n",
71 approximate_pi,
72 std::abs(real_pi - approximate_pi));
73}