13 August 2019 10:18:54 AM
WEDGE_FELIPPA_RULE_TEST
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 1
3 2 6 1
-3 2 6 1
6 3 18 1
-6 2 12 1
7 3 21 1
12 4 48 1
Exact 1
Monomial exponents: 1 0 0
1 1 1 0.333333
3 2 6 0.333333
-3 2 6 0.333333
6 3 18 0.333333
-6 2 12 0.333333
7 3 21 0.333333
12 4 48 0.333333
Exact 0.333333
Monomial exponents: 0 1 0
1 1 1 0.333333
3 2 6 0.333333
-3 2 6 0.333333
6 3 18 0.333333
-6 2 12 0.333333
7 3 21 0.333333
12 4 48 0.333333
Exact 0.333333
Monomial exponents: 2 0 0
1 1 1 0.111111
3 2 6 0.166667
-3 2 6 0.166667
6 3 18 0.166667
-6 2 12 0.166667
7 3 21 0.166667
12 4 48 0.166667
Exact 0.166667
Monomial exponents: 1 1 0
1 1 1 0.111111
3 2 6 0.0833333
-3 2 6 0.0833333
6 3 18 0.0833333
-6 2 12 0.0833333
7 3 21 0.0833333
12 4 48 0.0833333
Exact 0.0833333
Monomial exponents: 0 2 0
1 1 1 0.111111
3 2 6 0.166667
-3 2 6 0.166667
6 3 18 0.166667
-6 2 12 0.166667
7 3 21 0.166667
12 4 48 0.166667
Exact 0.166667
Monomial exponents: 0 0 2
1 1 1 0
3 2 6 0.333333
-3 2 6 0.333333
6 3 18 0.333333
-6 2 12 0.333333
7 3 21 0.333333
12 4 48 0.333333
Exact 0.333333
Monomial exponents: 3 0 0
1 1 1 0.037037
3 2 6 0.101852
-3 2 6 0.0833333
6 3 18 0.1
-6 2 12 0.1
7 3 21 0.1
12 4 48 0.1
Exact 0.1
Monomial exponents: 2 1 0
1 1 1 0.037037
3 2 6 0.0324074
-3 2 6 0.0416667
6 3 18 0.0333333
-6 2 12 0.0333333
7 3 21 0.0333333
12 4 48 0.0333333
Exact 0.0333333
Monomial exponents: 1 2 0
1 1 1 0.037037
3 2 6 0.0324074
-3 2 6 0.0416667
6 3 18 0.0333333
-6 2 12 0.0333333
7 3 21 0.0333333
12 4 48 0.0333333
Exact 0.0333333
Monomial exponents: 0 3 0
1 1 1 0.037037
3 2 6 0.101852
-3 2 6 0.0833333
6 3 18 0.1
-6 2 12 0.1
7 3 21 0.1
12 4 48 0.1
Exact 0.1
Monomial exponents: 1 0 2
1 1 1 0
3 2 6 0.111111
-3 2 6 0.111111
6 3 18 0.111111
-6 2 12 0.111111
7 3 21 0.111111
12 4 48 0.111111
Exact 0.111111
Monomial exponents: 0 1 2
1 1 1 0
3 2 6 0.111111
-3 2 6 0.111111
6 3 18 0.111111
-6 2 12 0.111111
7 3 21 0.111111
12 4 48 0.111111
Exact 0.111111
Monomial exponents: 4 0 0
1 1 1 0.0123457
3 2 6 0.066358
-3 2 6 0.0416667
6 3 18 0.0666667
-6 2 12 0.0638889
7 3 21 0.0666667
12 4 48 0.0666667
Exact 0.0666667
Monomial exponents: 3 1 0
1 1 1 0.0123457
3 2 6 0.0177469
-3 2 6 0.0208333
6 3 18 0.0166667
-6 2 12 0.0180556
7 3 21 0.0166667
12 4 48 0.0166667
Exact 0.0166667
Monomial exponents: 2 2 0
1 1 1 0.0123457
3 2 6 0.00848765
-3 2 6 0.0208333
6 3 18 0.0111111
-6 2 12 0.00972222
7 3 21 0.0111111
12 4 48 0.0111111
Exact 0.0111111
Monomial exponents: 1 3 0
1 1 1 0.0123457
3 2 6 0.0177469
-3