21 January 2020 11:07:17 AM
TETRAHEDRON_FELIPPA_RULE_TEST
C version
Test the TETRAHEDRON_FELIPPA_RULE library.
TETRAHEDRON_UNIT_MONOMIAL_TEST
For the unit tetrahedron,
TETRAHEDRON_UNIT_MONOMIAL returns the exact value of the
integral of X^ALPHA Y^BETA Z^GAMMA
Volume = 0.166667
ALPHA BETA GAMMA INTEGRAL
0 0 0 0.166667
0 0 1 0.0416667
0 0 2 0.0166667
0 0 3 0.00833333
0 0 4 0.0047619
0 1 0 0.0416667
0 1 1 0.00833333
0 1 2 0.00277778
0 1 3 0.00119048
0 2 0 0.0166667
0 2 1 0.00277778
0 2 2 0.000793651
0 3 0 0.00833333
0 3 1 0.00119048
0 4 0 0.0047619
1 0 0 0.0416667
1 0 1 0.00833333
1 0 2 0.00277778
1 0 3 0.00119048
1 1 0 0.00833333
1 1 1 0.00138889
1 1 2 0.000396825
1 2 0 0.00277778
1 2 1 0.000396825
1 3 0 0.00119048
2 0 0 0.0166667
2 0 1 0.00277778
2 0 2 0.000793651
2 1 0 0.00277778
2 1 1 0.000396825
2 2 0 0.000793651
3 0 0 0.00833333
3 0 1 0.00119048
3 1 0 0.00119048
4 0 0 0.0047619
TETRAHEDRON_UNIT_QUAD_TEST
For the unit tetrahedron,
we approximate monomial integrals with:
TETRAHEDRON_UNIT_O01,
TETRAHEDRON_UNIT_O04,
TETRAHEDRON_UNIT_O08,
TETRAHEDRON_UNIT_O08b,
TETRAHEDRON_UNIT_O14,
TETRAHEDRON_UNIT_O14b,
TETRAHEDRON_UNIT_O15,
TETRAHEDRON_UNIT_O15b,
TETRAHEDRON_UNIT_O24,
Monomial exponents: 0 0 0
1 0.166667
4 0.166667
8 0.166667
8 0.166667
14 0.166667
14 0.166667
15 0.166667
15 0.166667
24 0.166667
Exact 0.166667
Monomial exponents: 1 0 0
1 0.0416667
4 0.0416667
8 0.0416667
8 0.0416667
14 0.0416667
14 0.0416667
15 0.0416667
15 0.0416667
24 0.0416667
Exact 0.0416667
Monomial exponents: 0 1 0
1 0.0416667
4 0.0416667
8 0.0416667
8 0.0416667
14 0.0416667
14 0.0416667
15 0.0416667
15 0.0416667
24 0.0416667
Exact 0.0416667
Monomial exponents: 0 0 1
1 0.0416667
4 0.0416667
8 0.0416667
8 0.0416667
14 0.0416667
14 0.0416667
15 0.0416667
15 0.0416667
24 0.0416667
Exact 0.0416667
Monomial exponents: 2 0 0
1 0.0104167
4 0.0166667
8 0.0166667
8 0.0166667
14 0.0166667
14 0.0166667
15 0.0166667
15 0.0166667
24 0.0166667
Exact 0.0166667
Monomial exponents: 1 1 0
1 0.0104167
4 0.00833333
8 0.00833333
8 0.00833333
14 0.00833333
14 0.00833333
15 0.00833333
15 0.00833333
24 0.00833333
Exact 0.00833333
Monomial exponents: 0 2 0
1 0.0104167
4 0.0166667
8 0.0166667
8 0.0166667
14 0.0166667
14 0.0166667
15 0.0166667
15 0.0166667
24 0.0166667
Exact 0.0166667
Monomial exponents: 1 0 1
1 0.0104167
4 0.00833333
8 0.00833333
8 0.00833333
14 0.00833333
14 0.00833333
15 0.00833333
15 0.00833333
24 0.00833333
Exact 0.00833333
Monomial exponents: 0 1 1
1 0.0104167
4 0.00833333
8 0.00833333
8 0.00833333
14 0.00833333
14 0.00833333
15 0.00833333
15 0.00833333
24 0.00833333
Exact 0.00833333
Monomial exponents: 0 0 2
1 0.0104167
4 0.0166667
8 0.0166667
8 0.0166667
14 0.0166667
14 0.0166667
15 0.0166667
15 0.0166667
24 0.0166667
Exact 0.0166667
Monomial exponents: 3 0 0
1 0.00260417
4 0.00868921
8 0.00833333
8 0.00833333
14 0.00833333
14 0.00833333
15 0.00833333
15 0.00833333
24 0.00833333
Exact 0.00833333
Monomial exponents: 2 1 0
1 0.00260417
4 0.00265915
8 0.00277778
8 0.00277778
14 0.00277778
14 0.00277778
15 0.00277778
15 0.00277778
24 0.00277778
Exact 0.00277778
Monomial exponents: 1 2 0
1 0.00260417
4 0.00265915
8 0.00277778
8 0.00277778
14 0.00277778
14 0.00277778
15 0.00277778
15 0.00277778
24 0.00277778
Exact 0.00277778
Monomial exponents: 0 3 0
1 0.00260417
4 0.00868921
8 0.00833333
8 0.00833333
14 0.00833333
14 0.00833333
15 0.00833333
15 0.00833333
24 0.00833333
Exact 0.00833333
Monomial exponents: 2 0 1
1 0.00260417
4 0.00265915
8 0.00277778
8 0.00277778
14 0.00277778
14 0.00277778
15 0.00277778
15 0.00277778
24 0.00277778
Exact 0.00277778
Monomial exponents: 1 1 1
1 0.00260417
4 0.00150751
8 0.00138889
8 0.00138889
14 0.00138889
14 0.00138889
15 0.00138889
15 0.00138889
24 0.00138889
Exact 0.00138889
Monomial exponents: 0 2 1
1 0.00260417
4 0.00265915
8 0.00277778
8 0.00277778
14 0.00277778
14 0.00277778
15 0.00277778
15 0.00277778
24 0.00277778
Exact 0.00277778
Monomial exponents: 1 0 2
1 0.00260417
4 0.00265915
8 0.00277778
8 0.00277778
14 0.00277778
14 0.00277778
15 0.00277778
15 0.00277778
24 0.00277778
Exact 0.00277778
Monomial exponents: 0 1 2
1 0.00260417
4 0.00265915
8 0.00277778
8 0.00277778
14 0.00277778
14 0.00277778
15 0.00277778
15 0.00277778
24 0.00277778
Exact 0.00277778
Monomial exponents: 0 0 3
1 0.00260417
4 0.00868921
8 0.00833333
8 0.00833333
14 0.00833333
14 0.00833333
15 0.00833333
15 0