27 August 2014 8:07:42.792 PM PYRAMID_FELIPPA_RULE_PRB FORTRAN77 version Test the PYRAMID_FELIPPA_RULE library. PYRAMID_UNIT_MONOMIAL_TEST For the unit pyramid, PYRAMID_UNIT_MONOMIAL returns the exact value of the integral of X^ALPHA Y^BETA Z^GAMMA Volume = 1.33333 ALPHA BETA GAMMA INTEGRAL 0 0 0 1.33333 0 0 1 0.333333 0 0 2 0.133333 0 0 3 0.666667E-01 0 0 4 0.380952E-01 0 1 0 0.00000 0 1 1 0.00000 0 1 2 0.00000 0 1 3 0.00000 0 2 0 0.266667 0 2 1 0.444444E-01 0 2 2 0.126984E-01 0 3 0 0.00000 0 3 1 0.00000 0 4 0 0.114286 1 0 0 0.00000 1 0 1 0.00000 1 0 2 0.00000 1 0 3 0.00000 1 1 0 0.00000 1 1 1 0.00000 1 1 2 0.00000 1 2 0 0.00000 1 2 1 0.00000 1 3 0 0.00000 2 0 0 0.266667 2 0 1 0.444444E-01 2 0 2 0.126984E-01 2 1 0 0.00000 2 1 1 0.00000 2 2 0 0.634921E-01 3 0 0 0.00000 3 0 1 0.00000 3 1 0 0.00000 4 0 0 0.114286 PYRAMID_UNIT_QUAD_TEST For the unit pyramid, we approximate monomial integrals with: PYRAMID_UNIT_O01, PYRAMID_UNIT_O05, PYRAMID_UNIT_O06, PYRAMID_UNIT_O08, PYRAMID_UNIT_O08b, PYRAMID_UNIT_O09, PYRAMID_UNIT_O13, PYRAMID_UNIT_O18, PYRAMID_UNIT_O27, PYRAMID_UNIT_O48. Monomial exponents: 0 0 0 1 1.33333 5 1.33333 6 1.33333 8 1.33333 8 1.45172 9 1.33333 13 1.33333 18 1.33333 27 1.33333 48 1.33333 Exact 1.33333 Monomial exponents: 0 0 1 1 0.333333 5 0.333333 6 0.333333 8 0.333333 8 0.302654 9 0.333333 13 0.333333 18 0.333333 27 0.333333 48 0.333333 Exact 0.333333 Monomial exponents: 2 0 0 1 0.00000 5 0.266667 6 0.266667 8 0.266667 8 0.317135 9 0.266667 13 0.266667 18 0.266667 27 0.266667 48 0.266667 Exact 0.266667 Monomial exponents: 0 2 0 1 0.00000 5 0.266667 6 0.266667 8 0.266667 8 0.317135 9 0.266667 13 0.266667 18 0.266667 27 0.266667 48 0.266667 Exact 0.266667 Monomial exponents: 0 0 2 1 0.833333E-01 5 0.133333 6 0.133333 8 0.133333 8 0.106154 9 0.133333 13 0.133333 18 0.133333 27 0.133333 48 0.133333 Exact 0.133333 Monomial exponents: 2 0 1 1 0.00000 5 0.444444E-01 6 0.444444E-01 8 0.444444E-01 8 0.469792E-01 9 0.444444E-01 13 0.444444E-01 18 0.444444E-01 27 0.444444E-01 48 0.444444E-01 Exact 0.444444E-01 Monomial exponents: 0 2 1 1 0.00000 5 0.444444E-01 6 0.444444E-01 8 0.444444E-01 8 0.469792E-01 9 0.444444E-01 13 0.444444E-01 18 0.444444E-01 27 0.444444E-01 48 0.444444E-01 Exact 0.444444E-01 Monomial exponents: 0 0 3 1 0.208333E-01 5 0.766667E-01 6 0.773148E-01 8 0.666667E-01 8 0.507213E-01 9 0.669312E-01 13 0.666504E-01 18 0.666667E-01 27 0.666667E-01 48 0.666667E-01 Exact 0.666667E-01 Monomial exponents: 4 0 0 1 0.00000 5 0.632099E-01 6 0.634921E-01 8 0.632099E-01 8 0.777846E-01 9 0.634921E-01 13 0.664669E-01 18 0.113778 27 0.114286 48 0.114286 Exact 0.114286 Monomial exponents: 2 2 0 1 0.00000 5 0.632099E-01 6 0.634921E-01 8 0.632099E-01 8 0.777846E-01 9 0.634921E-01 13 0.634921E-01 18 0.632099E-01 27 0.634921E-01 48 0.634921E-01 Exact 0.634921E-01 Monomial exponents: 0 4 0 1 0.00000 5 0.632099E-01 6 0.634921E-01 8 0.632099E-01 8 0.777846E-01 9 0.634921E-01 13 0.664669E-01 18 0.113778 27 0.114286 48 0.114286 Exact 0.114286 Monomial exponents: 2 0 2 1 0.00000 5 0.740741E-02 6 0.740741E-02 8 0.118519E-01 8 0.104830E-01 9 0.126984E-01 13 0.126984E-01 18 0.118519E-01 27 0.126984E-01 48 0.126984E-01 Exact 0.126984E-01 Monomial exponents: 0 2 2 1 0.00000 5 0.740741E-02 6 0.740741E-02 8 0.118519E-01 8 0.104830E-01 9 0.126984E-01 13 0.126984E-01 18 0.118519E-01 27 0.126984E-01 48 0.126984E-01 Exact 0.126984E-01 Monomial exponents: 0 0 4 1 0.520833E-02 5 0.508889E-01 6 0.523148E-01 8 0.355556E-01 8 0.267514E-01 9 0.392038E-01 13 0.379673E-01 18 0.355556E-01 27 0.380952E-01 48 0.380952E-01 Exact 0.380952E-01 Monomial exponents: 4 0 1 1 0.00000 5 0.105350E-01 6 0.105820E-01 8 0.855967E-02 8 0.995664E-02 9 0.793651E-02 13 0.894583E-02 18 0.154074E-01 27 0.142857E-01 48 0.142857E-01 Exact 0.142857E-01 Monomial exponents: 2 2 1 1 0.00000 5 0.105350E-01 6 0.105820E-01 8 0.855967E-02 8 0.995664E-02 9 0.793651E-02 13 0.793651E-02 18 0.855967E-02 27 0.793651E-02 48 0.793651E-02 Exact 0.793651E-02 Monomial exponents: 0 4 1 1 0.00000 5 0.105350E-01 6 0.105820E-01 8 0.855967E-02 8 0.995664E-02 9 0.793651E-02 13 0.894583E-02 18 0.154074E-01 27 0.142857E-01 48 0.142857E-01 Exact 0.142857E-01 Monomial exponents: 2 0 3 1 0.00000 5 0.123457E-02 6 0.123457E-02 8 0.493827E-02 8 0.389044E-02 9 0.466742E-02 13 0.476190E-02 18 0.493827E-02 27 0.476190E-02 48 0.476190E-02 Exact 0.476190E-02 Monomial exponents: 0 2 3 1 0.00000 5 0.123457E-02 6 0.123457E-02 8 0.493827E-02 8 0.389044E-02 9 0.466742E-02 13 0.476190E-02 18 0.493827E-02 27 0.476190E-02 48 0.476190E-02 Exact 0.476190E-02 Monomial exponents: 0 0 5 1 0.130208E-02 5 0.351593E-01 6 0.371881E-01 8 0.192593E-01 8 0.144544E-01 9 0.261261E-01 13 0.238700E-01 18 0.192593E-01 27 0.238095E-01 48 0.238095E-01 Exact 0.238095E-01 PYRAMID_FELIPPA_RULE_PRB Normal end of execution. 27 August 2014 8:07:42.795 PM