25 August 2014 9:23:55.304 AM TRIANGLE_FELIPPA_RULE_PRB FORTRAN90 version Test the TRIANGLE_FELIPPA_RULE library. TRIANGLE_UNIT_MONOMIAL_TEST For the unit triangle, TRIANGLE_UNIT_MONOMIAL returns the exact value of the integral of X^ALPHA Y^BETA Volume = 0.500000 ALPHA BETA INTEGRAL 0 0 0.500000 0 1 0.166667 0 2 0.833333E-01 0 3 0.500000E-01 0 4 0.333333E-01 1 0 0.166667 1 1 0.416667E-01 1 2 0.166667E-01 1 3 0.833333E-02 2 0 0.833333E-01 2 1 0.166667E-01 2 2 0.555556E-02 3 0 0.500000E-01 3 1 0.833333E-02 4 0 0.333333E-01 TRIANGLE_UNIT_QUAD_TEST For the unit triangle, we approximate monomial integrals with: TRIANGLE_UNIT_O01, TRIANGLE_UNIT_O03, TRIANGLE_UNIT_O03b, TRIANGLE_UNIT_O06, TRIANGLE_UNIT_O06b, TRIANGLE_UNIT_O07, TRIANGLE_UNIT_O012, Exponents: 0 0 Order Error ----- ----- 1 0.500000 3 0.500000 3 0.500000 6 0.500000 6 0.500000 7 0.500000 12 0.500000 Exact 0.500000 Exponents: 1 0 Order Error ----- ----- 1 0.166667 3 0.166667 3 0.166667 6 0.166667 6 0.166667 7 0.166667 12 0.166667 Exact 0.166667 Exponents: 0 1 Order Error ----- ----- 1 0.166667 3 0.166667 3 0.166667 6 0.166667 6 0.166667 7 0.166667 12 0.166667 Exact 0.166667 Exponents: 2 0 Order Error ----- ----- 1 0.555556E-01 3 0.833333E-01 3 0.833333E-01 6 0.833333E-01 6 0.833333E-01 7 0.833333E-01 12 0.833333E-01 Exact 0.833333E-01 Exponents: 1 1 Order Error ----- ----- 1 0.555556E-01 3 0.416667E-01 3 0.416667E-01 6 0.416667E-01 6 0.416667E-01 7 0.416667E-01 12 0.416667E-01 Exact 0.416667E-01 Exponents: 0 2 Order Error ----- ----- 1 0.555556E-01 3 0.833333E-01 3 0.833333E-01 6 0.833333E-01 6 0.833333E-01 7 0.833333E-01 12 0.833333E-01 Exact 0.833333E-01 Exponents: 3 0 Order Error ----- ----- 1 0.185185E-01 3 0.509259E-01 3 0.416667E-01 6 0.500000E-01 6 0.500000E-01 7 0.500000E-01 12 0.500000E-01 Exact 0.500000E-01 Exponents: 2 1 Order Error ----- ----- 1 0.185185E-01 3 0.162037E-01 3 0.208333E-01 6 0.166667E-01 6 0.166667E-01 7 0.166667E-01 12 0.166667E-01 Exact 0.166667E-01 Exponents: 1 2 Order Error ----- ----- 1 0.185185E-01 3 0.162037E-01 3 0.208333E-01 6 0.166667E-01 6 0.166667E-01 7 0.166667E-01 12 0.166667E-01 Exact 0.166667E-01 Exponents: 0 3 Order Error ----- ----- 1 0.185185E-01 3 0.509259E-01 3 0.416667E-01 6 0.500000E-01 6 0.500000E-01 7 0.500000E-01 12 0.500000E-01 Exact 0.500000E-01 Exponents: 4 0 Order Error ----- ----- 1 0.617284E-02 3 0.331790E-01 3 0.208333E-01 6 0.333333E-01 6 0.319444E-01 7 0.333333E-01 12 0.333333E-01 Exact 0.333333E-01 Exponents: 3 1 Order Error ----- ----- 1 0.617284E-02 3 0.887346E-02 3 0.104167E-01 6 0.833333E-02 6 0.902778E-02 7 0.833333E-02 12 0.833333E-02 Exact 0.833333E-02 Exponents: 2 2 Order Error ----- ----- 1 0.617284E-02 3 0.424383E-02 3 0.104167E-01 6 0.555556E-02 6 0.486111E-02 7 0.555556E-02 12 0.555556E-02 Exact 0.555556E-02 Exponents: 1 3 Order Error ----- ----- 1 0.617284E-02 3 0.887346E-02 3 0.104167E-01 6 0.833333E-02 6 0.902778E-02 7 0.833333E-02 12 0.833333E-02 Exact 0.833333E-02 Exponents: 0 4 Order Error ----- ----- 1 0.617284E-02 3 0.331790E-01 3 0.208333E-01 6 0.333333E-01 6 0.319444E-01 7 0.333333E-01 12 0.333333E-01 Exact 0.333333E-01 Exponents: 5 0 Order Error ----- ----- 1 0.205761E-02 3 0.219907E-01 3 0.104167E-01 6 0.239351E-01 6 0.208333E-01 7 0.238095E-01 12 0.238095E-01 Exact 0.238095E-01 Exponents: 4 1 Order Error ----- ----- 1 0.205761E-02 3 0.559414E-02 3 0.520833E-02 6 0.469912E-02 6 0.555556E-02 7 0.476190E-02 12 0.476190E-02 Exact 0.476190E-02 Exponents: 3 2 Order Error ----- ----- 1 0.205761E-02 3 0.173611E-02 3 0.520833E-02 6 0.239351E-02 6 0.208333E-02 7 0.238095E-02 12 0.238095E-02 Exact 0.238095E-02 Exponents: 2 3 Order Error ----- ----- 1 0.205761E-02 3 0.173611E-02 3 0.520833E-02 6 0.239351E-02 6 0.208333E-02 7 0.238095E-02 12 0.238095E-02 Exact 0.238095E-02 Exponents: 1 4 Order Error ----- ----- 1 0.205761E-02 3 0.559414E-02 3 0.520833E-02 6 0.469912E-02 6 0.555556E-02 7 0.476190E-02 12 0.476190E-02 Exact 0.476190E-02 Exponents: 0 5 Order Error ----- ----- 1 0.205761E-02 3 0.219907E-01 3 0.104167E-01 6 0.239351E-01 6 0.208333E-01 7 0.238095E-01 12 0.238095E-01 Exact 0.238095E-01 Exponents: 6 0 Order Error ----- ----- 1 0.685871E-03 3 0.146391E-01 3 0.520833E-02 6 0.180884E-01 6 0.136960E-01 7 0.177753E-01 12 0.178571E-01 Exact 0.178571E-01 Exponents: 5 1 Order Error ----- ----- 1 0.685871E-03 3 0.367584E-02 3 0.260417E-02 6 0.292335E-02 6 0.356867E-02 7 0.301713E-02 12 0.297619E-02 Exact 0.297619E-02 Exponents: 4 2 Order Error ----- ----- 1 0.685871E-03 3 0.975223E-03 3 0.260417E-02 6 0.114099E-02 6 0.113812E-02 7 0.114638E-02 12 0.119048E-02 Exact 0.119048E-02 Exponents: 3 3 Order Error ----- ----- 1 0.685871E-03 3 0.460820E-03 3 0.260417E-02 6 0.949550E-03 6 0.675154E-03 7 0.938524E-03 12 0.892857E-03 Exact 0.892857E-03 Exponents: 2 4 Order Error ----- ----- 1 0.685871E-03 3 0.975223E-03 3 0.260417E-02 6 0.114099E-02 6 0.113812E-02 7 0.114638E-02 12 0.119048E-02 Exact 0.119048E-02 Exponents: 1 5 Order Error ----- ----- 1 0.685871E-03 3 0.367584E-02 3 0.260417E-02 6 0.292335E-02 6 0.356867E-02 7 0.301713E-02 12 0.297619E-02 Exact 0.297619E-02 Exponents: 0 6 Order Error ----- ----- 1 0.685871E-03 3 0.146391E-01 3 0.520833E-02 6 0.180884E-01 6 0.136960E-01 7 0.177753E-01 12 0.178571E-01 Exact 0.178571E-01 Exponents: 7 0 Order Error ----- ----- 1 0.228624E-03 3 0.975580E-02 3 0.260417E-02 6 0.141236E-01 6 0.904064E-02 7 0.136348E-01 12 0.138916E-01 Exact 0.138889E-01 Exponents: 6 1 Order Error ----- ----- 1 0.228624E-03 3 0.244163E-02 3 0.130208E-02 6 0.198238E-02 6 0.232767E-02 7 0.207021E-02 12 0.198275E-02 Exact 0.198413E-02 Exponents: 5 2 Order Error ----- ----- 1 0.228624E-03 3 0.619785E-03 3 0.130208E-02 6 0.583004E-03 6 0.688014E-03 7 0.601687E-03 12 0.661965E-03 Exact 0.661376E-03 Exponents: 4 3 Order Error ----- ----- 1 0.228624E-03 3 0.191115E-03 3 0.130208E-02 6 0.419591E-03 6 0.302212E-03 7 0.418121E-03 12 0.396629E-03 Exact 0.396825E-03 Exponents: 3 4 Order Error ----- ----- 1 0.228624E-03 3 0.191115E-03 3 0.130208E-02 6 0.419591E-03 6 0.302212E-03 7 0.418121E-03 12 0.396629E-03 Exact 0.396825E-03 Exponents: 2 5 Order Error ----- ----- 1 0.228624E-03 3 0.619785E-03 3 0.130208E-02 6 0.583004E-03 6 0.688014E-03 7 0.601687E-03 12 0.661965E-03 Exact 0.661376E-03 Exponents: 1 6 Order Error ----- ----- 1 0.228624E-03 3 0.244163E-02 3 0.130208E-02 6 0.198238E-02 6 0.232767E-02 7 0.207021E-02 12 0.198275E-02 Exact 0.198413E-02 Exponents: 0 7 Order Error ----- ----- 1 0.228624E-03 3 0.975580E-02 3 0.260417E-02 6 0.141236E-01 6 0.904064E-02 7 0.136348E-01 12 0.138916E-01 Exact 0.138889E-01 TRIANGLE_FELIPPA_RULE_PRB Normal end of execution. 25 August 2014 9:23:55.306 AM