26 August 2021 10:37:41.150 AM annulus_rule_test(): FORTRAN90 version Test annulus_rule(). annulus_area_test(): annulus_area() computes the area of an annulus with center = (CX,CY), inner radius R1 and outer radius R2. ( CX CY ) R1 R2 Area (-1.374480, 0.180316) 0.921038 1.025212 0.636951 ( 2.081647,-0.845150) 0.426664 1.071092 3.032252 (-0.444014, 3.118613) 0.744686 1.263893 3.276271 (-2.386573, 4.479252) 0.278195 0.327792 0.094420 (-3.928417,-0.988068) 0.978676 1.693600 6.001932 (-3.386770,-0.319858) 0.297915 0.407274 0.242275 (-0.493218, 0.361936) 0.340924 1.192111 4.099463 (-4.563584,-1.721205) 0.300979 0.352498 0.105766 ( 0.480293,-1.200268) 0.422284 1.321778 4.928443 ( 0.836076,-1.320093) 0.399172 1.367127 5.371172 annulus_rule_compute_test(): Test annulus_rule_compute(). W, X, Y for annulus quadrature: 1 0.545415E-01 0.578382 0.00000 2 0.872665E-01 0.790569 0.00000 3 0.545415E-01 0.956804 0.00000 4 0.545415E-01 0.500894 0.289191 5 0.872665E-01 0.684653 0.395285 6 0.545415E-01 0.828617 0.478402 7 0.545415E-01 0.289191 0.500894 8 0.872665E-01 0.395285 0.684653 9 0.545415E-01 0.478402 0.828617 10 0.545415E-01 0.354157E-16 0.578382 11 0.872665E-01 0.484084E-16 0.790569 12 0.545415E-01 0.585873E-16 0.956804 13 0.545415E-01 -0.289191 0.500894 14 0.872665E-01 -0.395285 0.684653 15 0.545415E-01 -0.478402 0.828617 16 0.545415E-01 -0.500894 0.289191 17 0.872665E-01 -0.684653 0.395285 18 0.545415E-01 -0.828617 0.478402 19 0.545415E-01 -0.578382 0.708314E-16 20 0.872665E-01 -0.790569 0.968168E-16 21 0.545415E-01 -0.956804 0.117175E-15 22 0.545415E-01 -0.500894 -0.289191 23 0.872665E-01 -0.684653 -0.395285 24 0.545415E-01 -0.828617 -0.478402 25 0.545415E-01 -0.289191 -0.500894 26 0.872665E-01 -0.395285 -0.684653 27 0.545415E-01 -0.478402 -0.828617 28 0.545415E-01 -0.106247E-15 -0.578382 29 0.872665E-01 -0.145225E-15 -0.790569 30 0.545415E-01 -0.175762E-15 -0.956804 31 0.545415E-01 0.289191 -0.500894 32 0.872665E-01 0.395285 -0.684653 33 0.545415E-01 0.478402 -0.828617 34 0.545415E-01 0.500894 -0.289191 35 0.872665E-01 0.684653 -0.395285 36 0.545415E-01 0.828617 -0.478402 annulus_rule_monomial_test(): annulus_rule_compute() supplies a quadrature rule for the annulus centered at ( 0.0000, 0.0000) with R1 = 0.0000 R2 = 1.0000 NR NT 1 X^2 Y^2 X^4 X^2Y^2 Y^4 X^6 4 16 3.14159 0.785398 0.785398 0.392699 0.130900 0.392699 0.245437 8 32 3.14159 0.785398 0.785398 0.392699 0.130900 0.392699 0.245437 16 64 3.14159 0.785398 0.785398 0.392699 0.130900 0.392699 0.245437 32 128 3.14159 0.785398 0.785398 0.392699 0.130900 0.392699 0.245437 64 256 3.14159 0.785398 0.785398 0.392699 0.130900 0.392699 0.245437 Exact 3.14159 0.785398 0.785398 0.392699 0.130900 0.392699 0.245437 annulus_rule_monomial_test(): annulus_rule_compute() supplies a quadrature rule for the annulus centered at ( 0.0000, 0.0000) with R1 = 0.5000 R2 = 1.0000 NR NT 1 X^2 Y^2 X^4 X^2Y^2 Y^4 X^6 4 16 2.35619 0.736311 0.736311 0.386563 0.128854 0.386563 0.244478 8 32 2.35619 0.736311 0.736311 0.386563 0.128854 0.386563 0.244478 16 64 2.35619 0.736311 0.736311 0.386563 0.128854 0.386563 0.244478 32 128 2.35619 0.736311 0.736311 0.386563 0.128854 0.386563 0.244478 64 256 2.35619 0.736311 0.736311 0.386563 0.128854 0.386563 0.244478 annulus_rule_monomial_test(): annulus_rule_compute() supplies a quadrature rule for the annulus centered at ( 1.0000, 0.0000) with R1 = 0.0000 R2 = 1.0000 NR NT 1 X^2 Y^2 X^4 X^2Y^2 Y^4 X^6 4 16 3.14159 3.92699 0.785398 8.24668 0.916298 0.392699 21.0585 8 32 3.14159 3.92699 0.785398 8.24668 0.916298 0.392699 21.0585 16 64 3.14159 3.92699 0.785398 8.24668 0.916298 0.392699 21.0585 32 128 3.14159 3.92699 0.785398 8.24668 0.916298 0.392699 21.0585 64 256 3.14159 3.92699 0.785398 8.24668 0.916298 0.392699 21.0585 annulus_rule_test(): Normal end of execution. 26 August 2021 10:37:41.156 AM