Thu Apr 13 19:55:45 2023

tetrahedron_monte_carlo_tests():
  Python version:3.8.10
  Test tetrahedron_monte_carlo().

tetrahedron_monte_carlo_test01():
  Integrate tetrahedron_integrand_03()
  Integration region is the unit tetrahedron.
  Use an increasing number of points P_NUM.

     P_NUM      1/(x**2+y**2+1)

         1        0.131660
         2        0.131248
         4        0.139526
         8        0.135605
        16        0.133513
        32        0.143330
        64        0.136048
       128        0.139543
       256        0.140963
       512        0.140268
      1024        0.140928
      2048        0.141540
      4096        0.141531
      8192        0.141172
     16384        0.141345
     32768        0.141287
     65536        0.141342

tetrahedron_monte_carlo_test02():
  Integrate tetrahedron_unit_integrand_user()
  Integration region is over a general tetrahedron.
  Use an increasing number of points P_NUM.

  Tetrahedron vertices:
[[1. 2. 3.]
 [4. 1. 2.]
 [2. 4. 4.]
 [3. 2. 5.]]

     P_NUM

         1      0.208991
         2      0.203066
         4      0.210918
         8      0.196011
        16      0.228833
        32      0.224907
        64      0.218999
       128      0.219930
       256      0.215671
       512      0.216788
      1024      0.217698
      2048      0.216515
      4096      0.217791
      8192      0.216975
     16384      0.216920
     32768      0.217319
     65536      0.217262

tetrahedron_monte_carlo_tests():
  Normal end of execution.
Thu Apr 13 19:55:47 2023