Tue May 20 21:35:02 2025 disk01_integrals_test(): python version: 3.10.12 numpy version: 1.26.4 Test disk01_integrals() 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 8 0 0.176028 0.171806 0.0042 6 8 0.000974424 0.000958738 1.6e-05 8 0 0.176028 0.171806 0.0042 4 0 0.396806 0.392699 0.0041 6 0 0.250022 0.245437 0.0046 2 8 0.014656 0.0143172 0.00034 0 0 3.14159 3.14159 4.4e-16 8 6 0.000968171 0.000958738 9.4e-06 8 6 0.000968171 0.000958738 9.4e-06 2 2 0.132446 0.1309 0.0015 4 8 0.00313637 0.00306796 6.8e-05 4 2 0.0497094 0.0490874 0.00062 8 6 0.000968171 0.000958738 9.4e-06 6 4 0.00620701 0.00613592 7.1e-05 0 6 0.250647 0.245437 0.0052 0 4 0.399785 0.392699 0.0071 6 0 0.250022 0.245437 0.0046 6 8 0.000974424 0.000958738 1.6e-05 0 0 3.14159 3.14159 4.4e-16 2 0 0.786627 0.785398 0.0012 disk01_sample_test(): disk01_sample samples the unit disk. Sample points in the unit disk. Row: 0 1 Col 0 : -0.00856771 -0.568361 1 : 0.00281232 0.677809 2 : 0.124902 -0.0844892 3 : -0.176306 -0.0952298 4 : 0.36493 -0.54455 5 : 0.528522 -0.186506 6 : -0.223454 0.944505 7 : 0.702058 -0.711261 8 : -0.646663 0.722339 9 : 0.444941 -0.083374 disk01_integrals_test(): Normal end of execution. Tue May 20 21:35:03 2025