Tue May 20 21:35:04 2025

disk01_monte_carlo_test():
  python version: 3.10.12
  numpy version:  1.26.4
  Test disk01_monte_carlo().

disk01_area_test():
  disk01_area() returns the area of the unit disk.

  disk01_area() = 3.14159

disk01_monomial_integral_test():
  disk01_monomial_integral() computes monomial integrals
  over the interior of the unit disk in 2D.
  Compare with a Monte Carlo value.

  Number of sample points used is 4192

  If any exponent is odd, the integral is zero.
  We will restrict this test to randomly chosen even exponents.

  Ex  Ey     MC-Estimate           Exact      Error

   6   0        0.247631        0.245437      0.0022
   2   0         0.78342        0.785398       0.002
   8   4      0.00307776      0.00306796     9.8e-06
   8   0        0.173702        0.171806      0.0019
   6   8     0.000961245     0.000958738     2.5e-06
   6   2       0.0244692       0.0245437     7.5e-05
   8   0        0.173702        0.171806      0.0019
   8   8     0.000374644     0.000372843     1.8e-06
   8   6      0.00096244     0.000958738     3.7e-06
   0   8        0.172054        0.171806     0.00025
   8   6      0.00096244     0.000958738     3.7e-06
   4   8      0.00306136      0.00306796     6.6e-06
   2   4       0.0481623       0.0490874     0.00093
   2   2        0.128134          0.1309      0.0028
   4   0        0.394464        0.392699      0.0018
   8   4      0.00307776      0.00306796     9.8e-06
   6   6      0.00218451       0.0021914     6.9e-06
   6   8     0.000961245     0.000958738     2.5e-06
   4   0        0.394464        0.392699      0.0018
   6   0        0.247631        0.245437      0.0022

disk01_sample_test():
  disk01_sample() samples the unit disk.

  Sample points in the unit disk.
[[-0.090233    0.75279263  0.19459458 -0.12078011  0.37436995 -0.57839635
   0.92587409 -0.29757552  0.16960368 -0.27430262]
 [-0.1102944   0.01201912  0.12336431  0.91738146  0.71033471 -0.14871476
  -0.19054401 -0.25016607 -0.19655801 -0.61691993]]

disk01_monte_carlo_test():
  Normal end of execution.
Tue May 20 21:35:04 2025