23 March 2023 4:30:45.743 PM TRIANGLE_NCC_RULE_TEST(): FORTRAN90 version Test TRIANGLE_NCC_RULE(). TEST01 TRIANGLE_NCC_RULE_NUM returns the number of rules; TRIANGLE_NCC_DEGREE returns the degree of a rule; TRIANGLE_NCC_ORDER_NUM returns the order of a rule. Number of available rules = 9 Rule Degree Order 1 0 1 2 1 3 3 2 3 4 3 10 5 4 12 6 5 21 7 6 25 8 7 36 9 8 42 TEST02 TRIANGLE_NCC_RULE returns the points and weights of an NCC rule for the triangle. In this test, we simply check that the weights sum to 1. Number of available rules = 9 Rule Sum of weights 1 1.000000000000000 2 1.000000000000000 3 1.000000000000000 4 1.000000000000000 5 1.000000000000000 6 1.000000000000000 7 1.000000000000000 8 0.9999999999999997 9 1.000000000000000 TEST03 TRIANGLE_NCC_RULE returns the points and weights of an NCC rule for the triangle. In this test, we simply check that, for each quadrature point, the barycentric coordinates sum to 1. Rule Suborder Sum of coordinates 1 1 1.000000000000000 2 1 1.000000000000000 3 1 1.000000000000000 4 3 1.000000000000000 1.000000000000000 1.000000000000000 5 3 1.000000000000000 1.000000000000000 1.000000000000000 6 5 1.000000000000000 1.000000000000000 1.000000000000000 1.000000000000000 1.000000000000000 7 6 1.000000000000000 1.000000000000000 0.9999999999999999 1.000000000000000 0.9999999999999999 1.000000000000000 8 8 1.000000000000000 1.000000000000000 1.000000000000000 1.000000000000000 1.000000000000000 1.000000000000000 1.000000000000000 0.9999999999999999 9 9 1.000000000000000 1.000000000000000 1.000000000000000 1.000000000000000 1.000000000000000 1.000000000000000 1.000000000000000 1.000000000000000 1.000000000000000 TEST04 TRIANGLE_NCC_RULE returns the points and weights of an NCC rule for the unit triangle. This routine uses those rules to estimate the integral of monomomials in the unit triangle. Integrate 2.00000 * X ^ 0 * Y ^ 0 Rule QUAD ERROR 1 1.00000 0.00000000 2 1.00000 0.00000000 3 1.00000 0.00000000 4 1.00000 0.00000000 5 1.00000 0.00000000 6 1.00000 0.00000000 7 1.00000 0.00000000 8 1.00000 0.00000000 9 1.00000 0.00000000 Integrate 6.00000 * X ^ 0 * Y ^ 1 Rule QUAD ERROR 1 1.00000 0.00000000 2 1.00000 0.00000000 3 1.00000 0.00000000 4 1.00000 0.00000000 5 1.00000 0.00000000 6 1.00000 0.00000000 7 1.00000 0.00000000 8 1.00000 0.00000000 9 1.00000 0.00000000 Integrate 12.0000 * X ^ 0 * Y ^ 2 Rule QUAD ERROR 1 0.666667 0.33333333 2 2.00000 1.00000000 3 1.00000 0.00000000 4 1.00000 0.00000000 5 1.00000 0.00000000 6 1.00000 0.00000000 7 1.00000 0.00000000 8 1.00000 0.00000000 9 1.00000 0.00000000 Integrate 20.0000 * X ^ 0 * Y ^ 3 Rule QUAD ERROR 1 0.370370 0.62962963 2 3.33333 2.33333333 3 0.833333 0.16666667 4 1.00000 0.00000000 5 1.00000 0.00000000 6 1.00000 0.00000000 7 1.00000 0.00000000 8 1.00000 0.00000000 9 1.00000 0.00000000 Integrate 30.0000 * X ^ 0 * Y ^ 4 Rule QUAD ERROR 1 0.185185 0.81481481 2 5.00000 4.00000000 3 0.625000 0.37500000 4 1.05556 0.05555556 5 1.00000 0.00000000 6 1.00000 0.00000000 7 1.00000 0.00000000 8 1.00000 0.00000000 9 1.00000 0.00000000 Integrate 42.0000 * X ^ 0 * Y ^ 5 Rule QUAD ERROR 1 0.864198E-01 0.91358025 2 7.00000 6.00000000 3 0.437500 0.56250000 4 1.16667 0.16666667 5 0.984375 0.01562500 6 1.00000 0.00000000 7 1.00000 0.00000000 8 1.00000 0.00000000 9 1.00000 0.00000000 Integrate 56.0000 * X ^ 0 * Y ^ 6 Rule QUAD ERROR 1 0.384088E-01 0.96159122 2 9.33333 8.33333333 3 0.291667 0.70833333 4 1.32510 0.32510288 5 0.947917 0.05208333 6 1.00587 0.00586667 7 1.00000 0.00000000 8 1.00000 0.00000000 9 1.00000 0.00000000 Integrate 72.0000 * X ^ 0 * Y ^ 7 Rule QUAD ERROR 1 0.164609E-01 0.98353909 2 12.0000 11.00000000 3 0.187500 0.81250000 4 1.52593 0.52592593 5 0.892969 0.10703125 6 1.02336 0.02336000 7 0.998148 0.00185185 8 1.00000 0.00000000 9 1.00000 0.00000000 Integrate 90.0000 * X ^ 0 * Y ^ 8 Rule QUAD ERROR 1 0.685871E-02 0.99314129 2 15.0000 14.00000000 3 0.117188 0.88281250 4 1.76749 0.76748971 5 0.824707 0.17529297 6 1.05600 0.05600000 7 0.991898 0.00810185 8 1.00071 0.00070974 9 1.00000 0.00000000 Integrate 110.000 * X ^ 0 * Y ^ 9 Rule QUAD ERROR 1 0.279429E-02 0.99720571 2 18.3333 17.33333333 3 0.716146E-01 0.92838542 4 2.04961 1.04961134 5 0.748596 0.25140381 6 1.10528 0.10528000 7 0.979224 0.02077618 8 1.00353 0.00353453 9 0.999765 0.00023524 Integrate 132.000 * X ^ 0 * Y ^ 10 Rule QUAD ERROR 1 0.111772E-02 0.99888228 2 22.0000 21.00000000 3 0.429688E-01 0.95703125 4 2.37235 1.37235178 5 0.669373 0.33062744 6 1.17143 0.17142784 7 0.959105 0.04089506 8 1.01016 0.01016401 9 0.998730 0.00127029 Integrate 6.00000 * X ^ 1 * Y ^ 0 Rule QUAD ERROR 1 1.00000 0.00000000 2 1.00000 0.00000000 3 1.00000 0.00000000 4 1.00000 0.00000000 5 1.00000 0.00000000 6 1.00000 0.00000000 7 1.00000 0.00000000 8 1.00000 0.00000000 9 1.00000 0.00000000 Integrate 24.0000 * X ^ 1 * Y ^ 1 Rule QUAD ERROR 1 1.33333 0.33333333 2 0.00000 1.00000000 3 1.00000 0.00000000 4 1.00000 0.00000000 5 1.00000 0.00000000 6 1.00000 0.00000000 7 1.00000 0.00000000 8 1.00000 0.00000000 9 1.00000 0.00000000 Integrate 60.0000 * X ^ 1 * Y ^ 2 Rule QUAD ERROR 1 1.11111 0.11111111 2 0.00000 1.00000000 3 1.25000 0.25000000 4 1.00000 0.00000000 5 1.00000 0.00000000 6 1.00000 0.00000000 7 1.00000 0.00000000 8 1.00000 0.00000000 9 1.00000 0.00000000 Integrate 120.000 * X ^ 1 * Y ^ 3 Rule QUAD ERROR 1 0.740741 0.25925926 2 0.00000 1.00000000 3 1.25000 0.25000000 4 0.888889 0.11111111 5 1.00000 0.00000000 6 1.00000 0.00000000 7 1.00000 0.00000000 8 1.00000 0.00000000 9 1.00000 0.00000000 Integrate 210.000 * X ^ 1 * Y ^ 4 Rule QUAD ERROR 1 0.432099 0.56790123 2 0.00000 1.00000000 3 1.09375 0.09375000 4 0.777778 0.22222222 5 1.03906 0.03906250 6 1.00000 0.00000000 7 1.00000 0.00000000 8 1.00000 0.00000000 9 1.00000 0.00000000 Integrate 336.000 * X ^ 1 * Y ^ 5 Rule QUAD ERROR 1 0.230453 0.76954733 2 0.00000 1.00000000 3 0.875000 0.12500000 4 0.691358 0.30864198 5 1.09375 0.09375000 6 0.982400 0.01760000 7 1.00000 0.00000000 8 1.00000 0.00000000 9 1.00000 0.00000000 Integrate 504.000 * X ^ 1 * Y ^ 6 Rule QUAD ERROR 1 0.115226 0.88477366 2 0.00000 1.00000000 3 0.656250 0.34375000 4 0.622222 0.37777778 5 1.14023 0.14023437 6 0.944640 0.05536000 7 1.00648 0.00648148 8 1.00000 0.00000000 9 1.00000 0.00000000 Integrate 720.000 * X ^ 1 * Y ^ 7 Rule QUAD ERROR 1 0.548697E-01 0.94513032 2 0.00000 1.00000000 3 0.468750 0.53125000 4 0.559671 0.44032922 5 1.16602 0.16601562 6 0.892800 0.10720000 7 1.02315 0.02314815 8 0.997161 0.00283895 9 1.00000 0.00000000 Integrate 990.000 * X ^ 1 * Y ^ 8 Rule QUAD ERROR 1 0.251486E-01 0.97485139 2 0.00000 1.00000000 3 0.322266 0.67773438 4 0.497942 0.50205761 5 1.16721 0.16720581 6 0.834240 0.16576000 7 1.04893 0.04893261 8 0.987998 0.01200180 9 1.00106 0.00105858 Integrate 1320.00 * X ^ 1 * Y ^ 9 Rule QUAD ERROR 1 0.111772E-01 0.98882284 2 0.00000 1.00000000 3 0.214844 0.78515625 4 0.435909 0.56409084 5 1.14471 0.14471436 6 0.774541 0.22545920 7 1.07982 0.07981824 8 0.970387 0.02961290 9 1.00494 0.00494003 Integrate 12.0000 * X ^ 2 * Y ^ 0 Rule QUAD ERROR 1 0.666667 0.33333333 2 2.00000 1.00000000 3 1.00000 0.00000000 4 1.00000 0.00000000 5 1.00000 0.00000000 6 1.00000 0.00000000 7 1.00000 0.00000000 8 1.00000 0.00000000 9 1.00000 0.00000000 Integrate 60.0000 * X ^ 2 * Y ^ 1 Rule QUAD ERROR 1 1.11111 0.11111111 2 0.00000 1.00000000 3 1.25000 0.25000000 4 1.00000 0.00000000 5 1.00000 0.00000000 6 1.00000 0.00000000 7 1.00000 0.00000000 8 1.00000 0.00000000 9 1.00000 0.00000000 Integrate 180.000 * X ^ 2 * Y ^ 2 Rule QUAD ERROR 1 1.11111 0.11111111 2 0.00000 1.00000000 3 1.87500 0.87500000 4 1.16667 0.16666667 5 1.00000 0.00000000 6 1.00000 0.00000000 7 1.00000 0.00000000 8 1.00000 0.00000000 9 1.00000 0.00000000 Integrate 420.000 * X ^ 2 * Y ^ 3 Rule QUAD ERROR 1 0.864198 0.13580247 2 0.00000 1.00000000 3 2.18750 1.18750000 4 1.16667 0.16666667 5 0.984375 0.01562500 6 1.00000 0.00000000 7 1.00000 0.00000000 8 1.00000 0.00000000 9 1.00000 0.00000000 Integrate 840.000 * X ^ 2 * Y ^ 4 Rule QUAD ERROR 1 0.576132 0.42386831 2 0.00000 1.00000000 3 2.18750 1.18750000 4 1.12346 0.12345679 5 1.05729 0.05729167 6 1.03467 0.03466667 7 1.00000 0.00000000 8 1.00000 0.00000000 9 1.00000 0.00000000 Integrate 1512.00 * X ^ 2 * Y ^ 5 Rule QUAD ERROR 1 0.345679 0.65432099 2 0.00000 1.00000000 3 1.96875 0.96875000 4 1.08889 0.08888889 5 1.20586 0.20585938 6 1.07712 0.07712000 7 0.991667 0.00833333 8 1.00000 0.00000000 9 1.00000 0.00000000 Integrate 2520.00 * X ^ 2 * Y ^ 6 Rule QUAD ERROR 1 0.192044 0.80795610 2 0.00000 1.00000000 3 1.64062 0.64062500 4 1.06584 0.06584362 5 1.38770 0.38769531 6 1.10720 0.10720000 7 0.986883 0.01311728 8 1.00723 0.00722913 9 1.00000 0.00000000 Integrate 3960.00 * X ^ 2 * Y ^ 7 Rule QUAD ERROR 1 0.100594 0.89940558 2 0.00000 1.00000000 3 1.28906 0.28906250 4 1.04115 0.04115226 5 1.56702 0.56701660 6 1.11936 0.11936000 7 1.00154 0.00154321 8 1.02445 0.02444701 9 0.997821 0.00217946 Integrate 5940.00 * X ^ 2 * Y ^ 8 Rule QUAD ERROR 1 0.502972E-01 0.94970279 2 0.00000 1.00000000 3 0.966797 0.03320312 4 1.00343 0.00342936 5 1.72110 0.72109985 6 1.11598 0.11598080 7 1.04469 0.04468879 8 1.04857 0.04857184 9 0.992792 0.00720787 Integrate 20.0000 * X ^ 3 * Y ^ 0 Rule QUAD ERROR 1 0.370370 0.62962963 2 3.33333 2.33333333 3 0.833333 0.16666667 4 1.00000 0.00000000 5 1.00000 0.00000000 6 1.00000 0.00000000 7 1.00000 0.00000000 8 1.00000 0.00000000 9 1.00000 0.00000000 Integrate 120.000 * X ^ 3 * Y ^ 1 Rule QUAD ERROR 1 0.740741 0.25925926 2 0.00000 1.00000000 3 1.25000 0.25000000 4 0.888889 0.11111111 5 1.00000 0.00000000 6 1.00000 0.00000000 7 1.00000 0.00000000 8 1.00000 0.00000000 9 1.00000 0.00000000 Integrate 420.000 * X ^ 3 * Y ^ 2 Rule QUAD ERROR 1 0.864198 0.13580247 2 0.00000 1.00000000 3 2.18750 1.18750000 4 1.16667 0.16666667 5 0.984375 0.01562500 6 1.00000 0.00000000 7 1.00000 0.00000000 8 1.00000 0.00000000 9 1.00000 0.00000000 Integrate 1120.00 * X ^ 3 * Y ^ 3 Rule QUAD ERROR 1 0.768176 0.23182442 2 0.00000 1.00000000 3 2.91667 1.91666667 4 1.26749 0.26748971 5 0.875000 0.12500000 6 0.960000 0.04000000 7 1.00000 0.00000000 8 1.00000 0.00000000 9 1.00000 0.00000000 Integrate 2520.00 * X ^ 3 * Y ^ 4 Rule QUAD ERROR 1 0.576132 0.42386831 2 0.00000 1.00000000 3 3.28125 2.28125000 4 1.29630 0.29629630 5 0.861328 0.13867188 6 0.969600 0.03040000 7 1.00463 0.00462963 8 1.00000 0.00000000 9 1.00000 0.00000000 Integrate 5040.00 * X ^ 3 * Y ^ 5 Rule QUAD ERROR 1 0.384088 0.61591221 2 0.00000 1.00000000 3 3.28125 2.28125000 4 1.32510 0.32510288 5 0.984375 0.01562500 6 1.01760 0.01760000 7 0.981481 0.01851852 8 0.988249 0.01175106 9 1.00000 0.00000000 Integrate 9240.00 * X ^ 3 * Y ^ 6 Rule QUAD ERROR 1 0.234720 0.76527968 2 0.00000 1.00000000 3 3.00781 2.00781250 4 1.37311 0.37311385 5 1.21252 0.21252441 6 1.07712 0.07712000 7 0.931756 0.06824417 8 0.973771 0.02622915 9 1.00269 0.00268809 Integrate 15840.0 * X ^ 3 * Y ^ 7 Rule QUAD ERROR 1 0.134126 0.86587410 2 0.00000 1.00000000 3 2.57812 1.57812500 4 1.42844 0.42844079 5 1.49854 0.49853516 6 1.13203 0.13203200 7 0.879887 0.12011317 8 0.969272 0.03072821 9 1.00305 0.00305176 Integrate 30.0000 * X ^ 4 * Y ^ 0 Rule QUAD ERROR 1 0.185185 0.81481481 2 5.00000 4.00000000 3 0.625000 0.37500000 4 1.05556 0.05555556 5 1.00000 0.00000000 6 1.00000 0.00000000 7 1.00000 0.00000000 8 1.00000 0.00000000 9 1.00000 0.00000000 Integrate 210.000 * X ^ 4 * Y ^ 1 Rule QUAD ERROR 1 0.432099 0.56790123 2 0.00000 1.00000000 3 1.09375 0.09375000 4 0.777778 0.22222222 5 1.03906 0.03906250 6 1.00000 0.00000000 7 1.00000 0.00000000 8 1.00000 0.00000000 9 1.00000 0.00000000 Integrate 840.000 * X ^ 4 * Y ^ 2 Rule QUAD ERROR 1 0.576132 0.42386831 2 0.00000 1.00000000 3 2.18750 1.18750000 4 1.12346 0.12345679 5 1.05729 0.05729167 6 1.03467 0.03466667 7 1.00000 0.00000000 8 1.00000 0.00000000 9 1.00000 0.00000000 Integrate 2520.00 * X ^ 4 * Y ^ 3 Rule QUAD ERROR 1 0.576132 0.42386831 2 0.00000 1.00000000 3 3.28125 2.28125000 4 1.29630 0.29629630 5 0.861328 0.13867188 6 0.969600 0.03040000 7 1.00463 0.00462963 8 1.00000 0.00000000 9 1.00000 0.00000000 Integrate 6300.00 * X ^ 4 * Y ^ 4 Rule QUAD ERROR 1 0.480110 0.51989026 2 0.00000 1.00000000 3 4.10156 3.10156250 4 1.36831 0.36831276 5 0.700684 0.29931641 6 0.939200 0.06080000 7 1.04360 0.04359568 8 1.01401 0.01401202 9 1.00000 0.00000000 Integrate 13860.0 * X ^ 4 * Y ^ 5 Rule QUAD ERROR 1 0.352080 0.64791952 2 0.00000 1.00000000 3 4.51172 3.51171875 4 1.42593 0.42592593 5 0.690857 0.30914307 6 0.960960 0.03904000 7 1.05530 0.05529835 8 1.01239 0.01238571 9 0.998786 0.00121435 Integrate 27720.0 * X ^ 4 * Y ^ 6 Rule QUAD ERROR 1 0.234720 0.76527968 2 0.00000 1.00000000 3 4.51172 3.51171875 4 1.51395 0.51394604 5 0.841248 0.15875244 6 1.01348 0.01347840 7 1.00626 0.00625857 8 0.991878 0.00812162 9 1.00880 0.00880051 Integrate 42.0000 * X ^ 5 * Y ^ 0 Rule QUAD ERROR 1 0.864198E-01 0.91358025 2 7.00000 6.00000000 3 0.437500 0.56250000 4 1.16667 0.16666667 5 0.984375 0.01562500 6 1.00000 0.00000000 7 1.00000 0.00000000 8 1.00000 0.00000000 9 1.00000 0.00000000 Integrate 336.000 * X ^ 5 * Y ^ 1 Rule QUAD ERROR 1 0.230453 0.76954733 2 0.00000 1.00000000 3 0.875000 0.12500000 4 0.691358 0.30864198 5 1.09375 0.09375000 6 0.982400 0.01760000 7 1.00000 0.00000000 8 1.00000 0.00000000 9 1.00000 0.00000000 Integrate 1512.00 * X ^ 5 * Y ^ 2 Rule QUAD ERROR 1 0.345679 0.65432099 2 0.00000 1.00000000 3 1.96875 0.96875000 4 1.08889 0.08888889 5 1.20586 0.20585938 6 1.07712 0.07712000 7 0.991667 0.00833333 8 1.00000 0.00000000 9 1.00000 0.00000000 Integrate 5040.00 * X ^ 5 * Y ^ 3 Rule QUAD ERROR 1 0.384088 0.61591221 2 0.00000 1.00000000 3 3.28125 2.28125000 4 1.32510 0.32510288 5 0.984375 0.01562500 6 1.01760 0.01760000 7 0.981481 0.01851852 8 0.988249 0.01175106 9 1.00000 0.00000000 Integrate 13860.0 * X ^ 5 * Y ^ 4 Rule QUAD ERROR 1 0.352080 0.64791952 2 0.00000 1.00000000 3 4.51172 3.51171875 4 1.42593 0.42592593 5 0.690857 0.30914307 6 0.960960 0.03904000 7 1.05530 0.05529835 8 1.01239 0.01238571 9 0.998786 0.00121435 Integrate 33264.0 * X ^ 5 * Y ^ 5 Rule QUAD ERROR 1 0.281664 0.71833562 2 0.00000 1.00000000 3 5.41406 4.41406250 4 1.47874 0.47873800 5 0.507568 0.49243164 6 0.953948 0.04605184 7 1.14159 0.14158951 8 1.03405 0.03404619 9 0.984924 0.01507568 Integrate 56.0000 * X ^ 6 * Y ^ 0 Rule QUAD ERROR 1 0.384088E-01 0.96159122 2 9.33333 8.33333333 3 0.291667 0.70833333 4 1.32510 0.32510288 5 0.947917 0.05208333 6 1.00587 0.00586667 7 1.00000 0.00000000 8 1.00000 0.00000000 9 1.00000 0.00000000 Integrate 504.000 * X ^ 6 * Y ^ 1 Rule QUAD ERROR 1 0.115226 0.88477366 2 0.00000 1.00000000 3 0.656250 0.34375000 4 0.622222 0.37777778 5 1.14023 0.14023437 6 0.944640 0.05536000 7 1.00648 0.00648148 8 1.00000 0.00000000 9 1.00000 0.00000000 Integrate 2520.00 * X ^ 6 * Y ^ 2 Rule QUAD ERROR 1 0.192044 0.80795610 2 0.00000 1.00000000 3 1.64062 0.64062500 4 1.06584 0.06584362 5 1.38770 0.38769531 6 1.10720 0.10720000 7 0.986883 0.01311728 8 1.00723 0.00722913 9 1.00000 0.00000000 Integrate 9240.00 * X ^ 6 * Y ^ 3 Rule QUAD ERROR 1 0.234720 0.76527968 2 0.00000 1.00000000 3 3.00781 2.00781250 4 1.37311 0.37311385 5 1.21252 0.21252441 6 1.07712 0.07712000 7 0.931756 0.06824417 8 0.973771 0.02622915 9 1.00269 0.00268809 Integrate 27720.0 * X ^ 6 * Y ^ 4 Rule QUAD ERROR 1 0.234720 0.76527968 2 0.00000 1.00000000 3 4.51172 3.51171875 4 1.51395 0.51394604 5 0.841248 0.15875244 6 1.01348 0.01347840 7 1.00626 0.00625857 8 0.991878 0.00812162 9 1.00880 0.00880051 Integrate 72.0000 * X ^ 7 * Y ^ 0 Rule QUAD ERROR 1 0.164609E-01 0.98353909 2 12.0000 11.00000000 3 0.187500 0.81250000 4 1.52593 0.52592593 5 0.892969 0.10703125 6 1.02336 0.02336000 7 0.998148 0.00185185 8 1.00000 0.00000000 9 1.00000 0.00000000 Integrate 720.000 * X ^ 7 * Y ^ 1 Rule QUAD ERROR 1 0.548697E-01 0.94513032 2 0.00000 1.00000000 3 0.468750 0.53125000 4 0.559671 0.44032922 5 1.16602 0.16601562 6 0.892800 0.10720000 7 1.02315 0.02314815 8 0.997161 0.00283895 9 1.00000 0.00000000 Integrate 3960.00 * X ^ 7 * Y ^ 2 Rule QUAD ERROR 1 0.100594 0.89940558 2 0.00000 1.00000000 3 1.28906 0.28906250 4 1.04115 0.04115226 5 1.56702 0.56701660 6 1.11936 0.11936000 7 1.00154 0.00154321 8 1.02445 0.02444701 9 0.997821 0.00217946 Integrate 15840.0 * X ^ 7 * Y ^ 3 Rule QUAD ERROR 1 0.134126 0.86587410 2 0.00000 1.00000000 3 2.57812 1.57812500 4 1.42844 0.42844079 5 1.49854 0.49853516 6 1.13203 0.13203200 7 0.879887 0.12011317 8 0.969272 0.03072821 9 1.00305 0.00305176 Integrate 90.0000 * X ^ 8 * Y ^ 0 Rule QUAD ERROR 1 0.685871E-02 0.99314129 2 15.0000 14.00000000 3 0.117188 0.88281250 4 1.76749 0.76748971 5 0.824707 0.17529297 6 1.05600 0.05600000 7 0.991898 0.00810185 8 1.00071 0.00070974 9 1.00000 0.00000000 Integrate 990.000 * X ^ 8 * Y ^ 1 Rule QUAD ERROR 1 0.251486E-01 0.97485139 2 0.00000 1.00000000 3 0.322266 0.67773438 4 0.497942 0.50205761 5 1.16721 0.16720581 6 0.834240 0.16576000 7 1.04893 0.04893261 8 0.987998 0.01200180 9 1.00106 0.00105858 Integrate 5940.00 * X ^ 8 * Y ^ 2 Rule QUAD ERROR 1 0.502972E-01 0.94970279 2 0.00000 1.00000000 3 0.966797 0.03320312 4 1.00343 0.00342936 5 1.72110 0.72109985 6 1.11598 0.11598080 7 1.04469 0.04468879 8 1.04857 0.04857184 9 0.992792 0.00720787 Integrate 110.000 * X ^ 9 * Y ^ 0 Rule QUAD ERROR 1 0.279429E-02 0.99720571 2 18.3333 17.33333333 3 0.716146E-01 0.92838542 4 2.04961 1.04961134 5 0.748596 0.25140381 6 1.10528 0.10528000 7 0.979224 0.02077618 8 1.00353 0.00353453 9 0.999765 0.00023524 Integrate 1320.00 * X ^ 9 * Y ^ 1 Rule QUAD ERROR 1 0.111772E-01 0.98882284 2 0.00000 1.00000000 3 0.214844 0.78515625 4 0.435909 0.56409084 5 1.14471 0.14471436 6 0.774541 0.22545920 7 1.07982 0.07981824 8 0.970387 0.02961290 9 1.00494 0.00494003 Integrate 132.000 * X ^ 10 * Y ^ 0 Rule QUAD ERROR 1 0.111772E-02 0.99888228 2 22.0000 21.00000000 3 0.429688E-01 0.95703125 4 2.37235 1.37235178 5 0.669373 0.33062744 6 1.17143 0.17142784 7 0.959105 0.04089506 8 1.01016 0.01016401 9 0.998730 0.00127029 TEST05 REFERENCE_TO_PHYSICAL_T3 transforms a rule on the unit (reference) triangle to a rule on an arbitrary (physical) triangle. The reference triangle: 1 0.00000 0.00000 2 1.00000 0.00000 3 0.00000 1.00000 Rule 3 for reference triangle with area = 0.500000 X Y W 1 0.500000 0.500000 0.333333 2 0.500000 0.00000 0.333333 3 0.00000 0.500000 0.333333 The physical triangle: 1 1.00000 2.00000 2 1.00000 1.00000 3 3.00000 2.00000 Rule 3 for physical triangle with area = 1.00000 X Y W 1 2.00000 1.50000 0.333333 2 1.00000 1.50000 0.333333 3 2.00000 2.00000 0.333333 TRIANGLE_NCC_RULE_TEST(): Normal end of execution. 23 March 2023 4:30:45.745 PM