20 August 2014 11:55:53 AM WEDGE_FELIPPA_RULE_PRB C++ version Test the WEDGE_FELIPPA_RULE library. TEST01 For the unit wedge, WEDGE_INTEGRAL returns the exact value of the integral of X^ALPHA Y^BETA Z^GAMMA Volume = 1 ALPHA BETA GAMMA INTEGRAL 0 0 0 1 0 0 1 0 0 0 2 0.333333 0 0 3 0 0 0 4 0.2 0 1 0 0.333333 0 1 1 0 0 1 2 0.111111 0 1 3 0 0 2 0 0.166667 0 2 1 0 0 2 2 0.0555556 0 3 0 0.1 0 3 1 0 0 4 0 0.0666667 1 0 0 0.333333 1 0 1 0 1 0 2 0.111111 1 0 3 0 1 1 0 0.0833333 1 1 1 0 1 1 2 0.0277778 1 2 0 0.0333333 1 2 1 0 1 3 0 0.0166667 2 0 0 0.166667 2 0 1 0 2 0 2 0.0555556 2 1 0 0.0333333 2 1 1 0 2 2 0 0.0111111 3 0 0 0.1 3 0 1 0 3 1 0 0.0166667 4 0 0 0.0666667 TEST02 For the unit wedge, we approximate monomial integrals with WEDG_UNIT_RULE. Monomial exponents: 0 0 0 1 1 1 3 2 1 -3 2 1 6 3 1 -6 2 1 7 3 1 12 4 1 Exact 1 Monomial exponents: 1 0 0 1 1 0.333333 3 2 0.333333 -3 2 0.333333 6 3 0.333333 -6 2 0.333333 7 3 0.333333 12 4 0.333333 Exact 0.333333 Monomial exponents: 0 1 0 1 1 0.333333 3 2 0.333333 -3 2 0.333333 6 3 0.333333 -6 2 0.333333 7 3 0.333333 12 4 0.333333 Exact 0.333333 Monomial exponents: 2 0 0 1 1 0.111111 3 2 0.166667 -3 2 0.166667 6 3 0.166667 -6 2 0.166667 7 3 0.166667 12 4 0.166667 Exact 0.166667 Monomial exponents: 1 1 0 1 1 0.111111 3 2 0.0833333 -3 2 0.0833333 6 3 0.0833333 -6 2 0.0833333 7 3 0.0833333 12 4 0.0833333 Exact 0.0833333 Monomial exponents: 0 2 0 1 1 0.111111 3 2 0.166667 -3 2 0.166667 6 3 0.166667 -6 2 0.166667 7 3 0.166667 12 4 0.166667 Exact 0.166667 Monomial exponents: 0 0 2 1 1 0 3 2 0.333333 -3 2 0.333333 6 3 0.333333 -6 2 0.333333 7 3 0.333333 12 4 0.333333 Exact 0.333333 Monomial exponents: 3 0 0 1 1 0.037037 3 2 0.101852 -3 2 0.0833333 6 3 0.1 -6 2 0.1 7 3 0.1 12 4 0.1 Exact 0.1 Monomial exponents: 2 1 0 1 1 0.037037 3 2 0.0324074 -3 2 0.0416667 6 3 0.0333333 -6 2 0.0333333 7 3 0.0333333 12 4 0.0333333 Exact 0.0333333 Monomial exponents: 1 2 0 1 1 0.037037 3 2 0.0324074 -3 2 0.0416667 6 3 0.0333333 -6 2 0.0333333 7 3 0.0333333 12 4 0.0333333 Exact 0.0333333 Monomial exponents: 0 3 0 1 1 0.037037 3 2 0.101852 -3 2 0.0833333 6 3 0.1 -6 2 0.1 7 3 0.1 12 4 0.1 Exact 0.1 Monomial exponents: 1 0 2 1 1 0 3 2 0.111111 -3 2 0.111111 6 3 0.111111 -6 2 0.111111 7 3 0.111111 12 4 0.111111 Exact 0.111111 Monomial exponents: 0 1 2 1 1 0 3 2 0.111111 -3 2 0.111111 6 3 0.111111 -6 2 0.111111 7 3 0.111111 12 4 0.111111 Exact 0.111111 Monomial exponents: 4 0 0 1 1 0.0123457 3 2 0.066358 -3 2 0.0416667 6 3 0.0666667 -6 2 0.0638889 7 3 0.0666667 12 4 0.0666667 Exact 0.0666667 Monomial exponents: 3 1 0 1 1 0.0123457 3 2 0.0177469 -3 2 0.0208333 6 3 0.0166667 -6 2 0.0180556 7 3 0.0166667 12 4 0.0166667 Exact 0.0166667 Monomial exponents: 2 2 0 1 1 0.0123457 3 2 0.00848765 -3 2 0.0208333 6 3 0.0111111 -6 2 0.00972222 7 3 0.0111111 12 4 0.0111111 Exact 0.0111111 Monomial exponents: 1 3 0 1 1 0.0123457 3 2 0.0177469 -3 2 0.0208333 6 3 0.0166667 -6 2 0.0180556 7 3 0.0166667 12 4 0.0166667 Exact 0.0166667 Monomial exponents: 0 4 0 1 1 0.0123457 3 2 0.066358 -3 2 0.0416667 6 3 0.0666667 -6 2 0.0638889 7 3 0.0666667 12 4 0.0666667 Exact 0.0666667 Monomial exponents: 2 0 2 1 1 0 3 2 0.0555556 -3 2 0.0555556 6 3 0.0555556 -6 2 0.0555556 7 3 0.0555556 12 4 0.0555556 Exact 0.0555556 Monomial exponents: 1 1 2 1 1 0 3 2 0.0277778 -3 2 0.0277778 6 3 0.0277778 -6 2 0.0277778 7 3 0.0277778 12 4 0.0277778 Exact 0.0277778 Monomial exponents: 0 2 2 1 1 0 3 2 0.0555556 -3 2 0.0555556 6 3 0.0555556 -6 2 0.0555556 7 3 0.0555556 12 4 0.0555556 Exact 0.0555556 Monomial exponents: 0 0 4 1 1 0 3 2 0.111111 -3 2 0.111111 6 3 0.2 -6 2 0.111111 7 3 0.2 12 4 0.2 Exact 0.2 TEST03 For the unit wedge, write some rules to a file Rule Trig Line Total W_File X_File Order Order Order 0 1 1 1 wedge_felippa_1x1_w.txt wedge_felippa_1x1_x.txt 1 3 2 6 wedge_felippa_3x2_w.txt wedge_felippa_3x2_x.txt 2 -3 2 6 wedge_felippa_3bx2_w.txt wedge_felippa_3bx2_x.txt 3 6 3 18 wedge_felippa_6x3_w.txt wedge_felippa_6x3_x.txt 4 -6 2 12 wedge_felippa_6bx2_w.txt wedge_felippa_6bx2_x.txt 5 7 3 21 wedge_felippa_7x3_w.txt wedge_felippa_7x3_x.txt 6 12 4 48 wedge_felippa_12x4_w.txt wedge_felippa_12x4_x.txt WEDGE_FELIPPA_RULE_PRB Normal end of execution. 20 August 2014 11:55:53 AM