23 March 2023 9:14:28.222 AM triangle_twb_rule_test(): FORTRAN90 version Test triangle_twb_rule(). triangle_unit_quad_test(): Approximate monomial integrals in triangle with TWB rules. Monomial: x^ 0 y^ 0 1 1 0.500000 2 3 0.500000 4 6 0.500000 5 10 0.500000 7 15 0.500000 9 21 0.500000 11 28 0.500000 13 36 0.500000 14 45 0.500000 16 55 0.500000 18 66 0.500000 20 78 0.500000 21 91 0.500000 23 105 0.500000 25 120 0.500000 Exact 0.500000 Monomial: x^ 0 y^ 1 1 1 0.166667 2 3 0.166667 4 6 0.166667 5 10 0.166667 7 15 0.166667 9 21 0.166667 11 28 0.166667 13 36 0.166667 14 45 0.166667 16 55 0.166667 18 66 0.166667 20 78 0.166667 21 91 0.166667 23 105 0.166667 25 120 0.166667 Exact 0.166667 Monomial: x^ 1 y^ 0 1 1 0.166667 2 3 0.166667 4 6 0.166667 5 10 0.166667 7 15 0.166667 9 21 0.166667 11 28 0.166667 13 36 0.166667 14 45 0.166667 16 55 0.166667 18 66 0.166667 20 78 0.166667 21 91 0.166667 23 105 0.166667 25 120 0.166667 Exact 0.166667 Monomial: x^ 0 y^ 2 1 1 0.555556E-01 2 3 0.833333E-01 4 6 0.833333E-01 5 10 0.833333E-01 7 15 0.833333E-01 9 21 0.833333E-01 11 28 0.833333E-01 13 36 0.833333E-01 14 45 0.833333E-01 16 55 0.833333E-01 18 66 0.833333E-01 20 78 0.833333E-01 21 91 0.833333E-01 23 105 0.833333E-01 25 120 0.833333E-01 Exact 0.833333E-01 Monomial: x^ 1 y^ 1 1 1 0.555556E-01 2 3 0.416667E-01 4 6 0.416667E-01 5 10 0.416667E-01 7 15 0.416667E-01 9 21 0.416667E-01 11 28 0.416667E-01 13 36 0.416667E-01 14 45 0.416667E-01 16 55 0.416667E-01 18 66 0.416667E-01 20 78 0.416667E-01 21 91 0.416667E-01 23 105 0.416667E-01 25 120 0.416667E-01 Exact 0.416667E-01 Monomial: x^ 2 y^ 0 1 1 0.555556E-01 2 3 0.833333E-01 4 6 0.833333E-01 5 10 0.833333E-01 7 15 0.833333E-01 9 21 0.833333E-01 11 28 0.833333E-01 13 36 0.833333E-01 14 45 0.833333E-01 16 55 0.833333E-01 18 66 0.833333E-01 20 78 0.833333E-01 21 91 0.833333E-01 23 105 0.833333E-01 25 120 0.833333E-01 Exact 0.833333E-01 Monomial: x^ 0 y^ 3 1 1 0.185185E-01 2 3 0.509259E-01 4 6 0.500000E-01 5 10 0.500000E-01 7 15 0.500000E-01 9 21 0.500000E-01 11 28 0.500000E-01 13 36 0.500000E-01 14 45 0.500000E-01 16 55 0.500000E-01 18 66 0.500000E-01 20 78 0.500000E-01 21 91 0.500000E-01 23 105 0.500000E-01 25 120 0.500000E-01 Exact 0.500000E-01 Monomial: x^ 1 y^ 2 1 1 0.185185E-01 2 3 0.162037E-01 4 6 0.166667E-01 5 10 0.166667E-01 7 15 0.166667E-01 9 21 0.166667E-01 11 28 0.166667E-01 13 36 0.166667E-01 14 45 0.166667E-01 16 55 0.166667E-01 18 66 0.166667E-01 20 78 0.166667E-01 21 91 0.166667E-01 23 105 0.166667E-01 25 120 0.166667E-01 Exact 0.166667E-01 Monomial: x^ 2 y^ 1 1 1 0.185185E-01 2 3 0.162037E-01 4 6 0.166667E-01 5 10 0.166667E-01 7 15 0.166667E-01 9 21 0.166667E-01 11 28 0.166667E-01 13 36 0.166667E-01 14 45 0.166667E-01 16 55 0.166667E-01 18 66 0.166667E-01 20 78 0.166667E-01 21 91 0.166667E-01 23 105 0.166667E-01 25 120 0.166667E-01 Exact 0.166667E-01 Monomial: x^ 3 y^ 0 1 1 0.185185E-01 2 3 0.509259E-01 4 6 0.500000E-01 5 10 0.500000E-01 7 15 0.500000E-01 9 21 0.500000E-01 11 28 0.500000E-01 13 36 0.500000E-01 14 45 0.500000E-01 16 55 0.500000E-01 18 66 0.500000E-01 20 78 0.500000E-01 21 91 0.500000E-01 23 105 0.500000E-01 25 120 0.500000E-01 Exact 0.500000E-01 Monomial: x^ 0 y^ 4 1 1 0.617284E-02 2 3 0.331790E-01 4 6 0.333333E-01 5 10 0.333333E-01 7 15 0.333333E-01 9 21 0.333333E-01 11 28 0.333333E-01 13 36 0.333333E-01 14 45 0.333333E-01 16 55 0.333333E-01 18 66 0.333333E-01 20 78 0.333333E-01 21 91 0.333333E-01 23 105 0.333333E-01 25 120 0.333333E-01 Exact 0.333333E-01 Monomial: x^ 1 y^ 3 1 1 0.617284E-02 2 3 0.887346E-02 4 6 0.833333E-02 5 10 0.833333E-02 7 15 0.833333E-02 9 21 0.833333E-02 11 28 0.833333E-02 13 36 0.833333E-02 14 45 0.833333E-02 16 55 0.833333E-02 18 66 0.833333E-02 20 78 0.833333E-02 21 91 0.833333E-02 23 105 0.833333E-02 25 120 0.833333E-02 Exact 0.833333E-02 Monomial: x^ 2 y^ 2 1 1 0.617284E-02 2 3 0.424383E-02 4 6 0.555556E-02 5 10 0.555556E-02 7 15 0.555556E-02 9 21 0.555556E-02 11 28 0.555556E-02 13 36 0.555556E-02 14 45 0.555556E-02 16 55 0.555556E-02 18 66 0.555556E-02 20 78 0.555556E-02 21 91 0.555556E-02 23 105 0.555556E-02 25 120 0.555556E-02 Exact 0.555556E-02 Monomial: x^ 3 y^ 1 1 1 0.617284E-02 2 3 0.887346E-02 4 6 0.833333E-02 5 10 0.833333E-02 7 15 0.833333E-02 9 21 0.833333E-02 11 28 0.833333E-02 13 36 0.833333E-02 14 45 0.833333E-02 16 55 0.833333E-02 18 66 0.833333E-02 20 78 0.833333E-02 21 91 0.833333E-02 23 105 0.833333E-02 25 120 0.833333E-02 Exact 0.833333E-02 Monomial: x^ 4 y^ 0 1 1 0.617284E-02 2 3 0.331790E-01 4 6 0.333333E-01 5 10 0.333333E-01 7 15 0.333333E-01 9 21 0.333333E-01 11 28 0.333333E-01 13 36 0.333333E-01 14 45 0.333333E-01 16 55 0.333333E-01 18 66 0.333333E-01 20 78 0.333333E-01 21 91 0.333333E-01 23 105 0.333333E-01 25 120 0.333333E-01 Exact 0.333333E-01 Monomial: x^ 0 y^ 5 1 1 0.205761E-02 2 3 0.219907E-01 4 6 0.239351E-01 5 10 0.238095E-01 7 15 0.238095E-01 9 21 0.238095E-01 11 28 0.238095E-01 13 36 0.238095E-01 14 45 0.238095E-01 16 55 0.238095E-01 18 66 0.238095E-01 20 78 0.238095E-01 21 91 0.238095E-01 23 105 0.238095E-01 25 120 0.238095E-01 Exact 0.238095E-01 Monomial: x^ 1 y^ 4 1 1 0.205761E-02 2 3 0.559414E-02 4 6 0.469912E-02 5 10 0.476190E-02 7 15 0.476190E-02 9 21 0.476190E-02 11 28 0.476190E-02 13 36 0.476190E-02 14 45 0.476190E-02 16 55 0.476190E-02 18 66 0.476190E-02 20 78 0.476190E-02 21 91 0.476190E-02 23 105 0.476190E-02 25 120 0.476190E-02 Exact 0.476190E-02 Monomial: x^ 2 y^ 3 1 1 0.205761E-02 2 3 0.173611E-02 4 6 0.239351E-02 5 10 0.238095E-02 7 15 0.238095E-02 9 21 0.238095E-02 11 28 0.238095E-02 13 36 0.238095E-02 14 45 0.238095E-02 16 55 0.238095E-02 18 66 0.238095E-02 20 78 0.238095E-02 21 91 0.238095E-02 23 105 0.238095E-02 25 120 0.238095E-02 Exact 0.238095E-02 Monomial: x^ 3 y^ 2 1 1 0.205761E-02 2 3 0.173611E-02 4 6 0.239351E-02 5 10 0.238095E-02 7 15 0.238095E-02 9 21 0.238095E-02 11 28 0.238095E-02 13 36 0.238095E-02 14 45 0.238095E-02 16 55 0.238095E-02 18 66 0.238095E-02 20 78 0.238095E-02 21 91 0.238095E-02 23 105 0.238095E-02 25 120 0.238095E-02 Exact 0.238095E-02 Monomial: x^ 4 y^ 1 1 1 0.205761E-02 2 3 0.559414E-02 4 6 0.469912E-02 5 10 0.476190E-02 7 15 0.476190E-02 9 21 0.476190E-02 11 28 0.476190E-02 13 36 0.476190E-02 14 45 0.476190E-02 16 55 0.476190E-02 18 66 0.476190E-02 20 78 0.476190E-02 21 91 0.476190E-02 23 105 0.476190E-02 25 120 0.476190E-02 Exact 0.476190E-02 Monomial: x^ 5 y^ 0 1 1 0.205761E-02 2 3 0.219907E-01 4 6 0.239351E-01 5 10 0.238095E-01 7 15 0.238095E-01 9 21 0.238095E-01 11 28 0.238095E-01 13 36 0.238095E-01 14 45 0.238095E-01 16 55 0.238095E-01 18 66 0.238095E-01 20 78 0.238095E-01 21 91 0.238095E-01 23 105 0.238095E-01 25 120 0.238095E-01 Exact 0.238095E-01 triangle_twb_rule_test(): Normal end of execution. 23 March 2023 9:14:28.223 AM