diff options
Diffstat (limited to 'src/GeneralizedPolynomial-Uniform/Hermite/1d.coeffs')
9 files changed, 133 insertions, 0 deletions
diff --git a/src/GeneralizedPolynomial-Uniform/Hermite/1d.coeffs/1d.cube.order2/coeffs-I.compute.c b/src/GeneralizedPolynomial-Uniform/Hermite/1d.coeffs/1d.cube.order2/coeffs-I.compute.c new file mode 100644 index 0000000..b51e6dd --- /dev/null +++ b/src/GeneralizedPolynomial-Uniform/Hermite/1d.coeffs/1d.cube.order2/coeffs-I.compute.c @@ -0,0 +1,13 @@ +fp t4; +fp t3; +fp t2; +fp t1; + t4 = x*x; + t3 = x*t4; + t2 = RATIONAL(1.0,2.0); + t1 = RATIONAL(-1.0,2.0); + coeffs_I->coeff_m1 = t4+(x+t3)*t1; + coeffs_I->coeff_0 = RATIONAL(-5.0,2.0)*t4+RATIONAL(1.0,1.0)+RATIONAL(3.0, +2.0)*t3; + coeffs_I->coeff_p1 = RATIONAL(-3.0,2.0)*t3+t2*x+RATIONAL(2.0,1.0)*t4; + coeffs_I->coeff_p2 = t2*t3+t1*t4; diff --git a/src/GeneralizedPolynomial-Uniform/Hermite/1d.coeffs/1d.cube.order2/coeffs-dx.compute.c b/src/GeneralizedPolynomial-Uniform/Hermite/1d.coeffs/1d.cube.order2/coeffs-dx.compute.c new file mode 100644 index 0000000..d727b4e --- /dev/null +++ b/src/GeneralizedPolynomial-Uniform/Hermite/1d.coeffs/1d.cube.order2/coeffs-dx.compute.c @@ -0,0 +1,8 @@ +fp t6; + t6 = x*x; + coeffs_dx->coeff_m1 = RATIONAL(2.0,1.0)*x+RATIONAL(-3.0,2.0)*t6+RATIONAL( +-1.0,2.0); + coeffs_dx->coeff_0 = RATIONAL(-5.0,1.0)*x+RATIONAL(9.0,2.0)*t6; + coeffs_dx->coeff_p1 = RATIONAL(1.0,2.0)+RATIONAL(4.0,1.0)*x+RATIONAL(-9.0 +,2.0)*t6; + coeffs_dx->coeff_p2 = -x+RATIONAL(3.0,2.0)*t6; diff --git a/src/GeneralizedPolynomial-Uniform/Hermite/1d.coeffs/1d.cube.order2/coeffs-dxx.compute.c b/src/GeneralizedPolynomial-Uniform/Hermite/1d.coeffs/1d.cube.order2/coeffs-dxx.compute.c new file mode 100644 index 0000000..63c1f43 --- /dev/null +++ b/src/GeneralizedPolynomial-Uniform/Hermite/1d.coeffs/1d.cube.order2/coeffs-dxx.compute.c @@ -0,0 +1,4 @@ + coeffs_dxx->coeff_m1 = RATIONAL(2.0,1.0)+RATIONAL(-3.0,1.0)*x; + coeffs_dxx->coeff_0 = RATIONAL(-5.0,1.0)+RATIONAL(9.0,1.0)*x; + coeffs_dxx->coeff_p1 = RATIONAL(-9.0,1.0)*x+RATIONAL(4.0,1.0); + coeffs_dxx->coeff_p2 = RATIONAL(-1.0,1.0)+RATIONAL(3.0,1.0)*x; diff --git a/src/GeneralizedPolynomial-Uniform/Hermite/1d.coeffs/1d.cube.order3/coeffs-I.compute.c b/src/GeneralizedPolynomial-Uniform/Hermite/1d.coeffs/1d.cube.order3/coeffs-I.compute.c new file mode 100644 index 0000000..107d2ed --- /dev/null +++ b/src/GeneralizedPolynomial-Uniform/Hermite/1d.coeffs/1d.cube.order3/coeffs-I.compute.c @@ -0,0 +1,17 @@ +fp t10; +fp t9; +fp t8; +fp t7; + t10 = x*x; + t9 = x*t10; + t8 = RATIONAL(1.0,12.0); + t7 = RATIONAL(-1.0,12.0); + coeffs_I->coeff_m2 = RATIONAL(-1.0,6.0)*t10+(x+t9)*t8; + coeffs_I->coeff_m1 = RATIONAL(-2.0,3.0)*x+RATIONAL(5.0,4.0)*t10+RATIONAL( +-7.0,12.0)*t9; + coeffs_I->coeff_0 = RATIONAL(-7.0,3.0)*t10+RATIONAL(1.0,1.0)+RATIONAL(4.0 +,3.0)*t9; + coeffs_I->coeff_p1 = RATIONAL(-4.0,3.0)*t9+RATIONAL(2.0,3.0)*x+RATIONAL( +5.0,3.0)*t10; + coeffs_I->coeff_p2 = RATIONAL(7.0,12.0)*t9+t7*x+RATIONAL(-1.0,2.0)*t10; + coeffs_I->coeff_p3 = t7*t9+t8*t10; diff --git a/src/GeneralizedPolynomial-Uniform/Hermite/1d.coeffs/1d.cube.order3/coeffs-dx.compute.c b/src/GeneralizedPolynomial-Uniform/Hermite/1d.coeffs/1d.cube.order3/coeffs-dx.compute.c new file mode 100644 index 0000000..22285e9 --- /dev/null +++ b/src/GeneralizedPolynomial-Uniform/Hermite/1d.coeffs/1d.cube.order3/coeffs-dx.compute.c @@ -0,0 +1,11 @@ +fp t12; + t12 = x*x; + coeffs_dx->coeff_m2 = RATIONAL(-1.0,3.0)*x+RATIONAL(1.0,12.0)+RATIONAL( +1.0,4.0)*t12; + coeffs_dx->coeff_m1 = RATIONAL(-2.0,3.0)+RATIONAL(5.0,2.0)*x+RATIONAL( +-7.0,4.0)*t12; + coeffs_dx->coeff_0 = RATIONAL(-14.0,3.0)*x+RATIONAL(4.0,1.0)*t12; + coeffs_dx->coeff_p1 = RATIONAL(-4.0,1.0)*t12+RATIONAL(10.0,3.0)*x+ +RATIONAL(2.0,3.0); + coeffs_dx->coeff_p2 = RATIONAL(7.0,4.0)*t12+RATIONAL(-1.0,12.0)-x; + coeffs_dx->coeff_p3 = RATIONAL(1.0,6.0)*x+RATIONAL(-1.0,4.0)*t12; diff --git a/src/GeneralizedPolynomial-Uniform/Hermite/1d.coeffs/1d.cube.order3/coeffs-dxx.compute.c b/src/GeneralizedPolynomial-Uniform/Hermite/1d.coeffs/1d.cube.order3/coeffs-dxx.compute.c new file mode 100644 index 0000000..399cdbb --- /dev/null +++ b/src/GeneralizedPolynomial-Uniform/Hermite/1d.coeffs/1d.cube.order3/coeffs-dxx.compute.c @@ -0,0 +1,6 @@ + coeffs_dxx->coeff_m2 = RATIONAL(-1.0,3.0)+RATIONAL(1.0,2.0)*x; + coeffs_dxx->coeff_m1 = RATIONAL(-7.0,2.0)*x+RATIONAL(5.0,2.0); + coeffs_dxx->coeff_0 = RATIONAL(-14.0,3.0)+RATIONAL(8.0,1.0)*x; + coeffs_dxx->coeff_p1 = RATIONAL(10.0,3.0)+RATIONAL(-8.0,1.0)*x; + coeffs_dxx->coeff_p2 = RATIONAL(7.0,2.0)*x+RATIONAL(-1.0,1.0); + coeffs_dxx->coeff_p3 = RATIONAL(1.0,6.0)+RATIONAL(-1.0,2.0)*x; diff --git a/src/GeneralizedPolynomial-Uniform/Hermite/1d.coeffs/1d.cube.order4/coeffs-I.compute.c b/src/GeneralizedPolynomial-Uniform/Hermite/1d.coeffs/1d.cube.order4/coeffs-I.compute.c new file mode 100644 index 0000000..f170025 --- /dev/null +++ b/src/GeneralizedPolynomial-Uniform/Hermite/1d.coeffs/1d.cube.order4/coeffs-I.compute.c @@ -0,0 +1,33 @@ +fp t23; +fp t25; +fp t22; +fp t20; +fp t19; +fp t18; +fp t17; +fp t16; +fp t15; +fp t14; +fp t13; + t23 = x*x; + t25 = t23*t23; + t22 = x*t23; + t20 = x*t25; + t19 = RATIONAL(-1.0,24.0); + t18 = RATIONAL(1.0,12.0); + t17 = RATIONAL(-7.0,12.0); + t16 = RATIONAL(5.0,12.0); + t15 = RATIONAL(-1.0,12.0); + t14 = RATIONAL(13.0,12.0); + t13 = t19*t22; + coeffs_I->coeff_m2 = t18*x+RATIONAL(-1.0,8.0)*t23+t13+RATIONAL(1.0,8.0)* +t25+t19*t20; + coeffs_I->coeff_m1 = RATIONAL(-2.0,3.0)*x+t14*t23+t13+t17*t25+RATIONAL( +5.0,24.0)*t20; + coeffs_I->coeff_0 = RATIONAL(-25.0,12.0)*t23+t16*t22+t14*t25+RATIONAL(1.0 +,1.0)+RATIONAL(-5.0,12.0)*t20; + coeffs_I->coeff_p1 = t17*t22+RATIONAL(2.0,3.0)*x+RATIONAL(3.0,2.0)*t23+ +t16*t20-t25; + coeffs_I->coeff_p2 = t15*x+RATIONAL(-11.0,24.0)*t23+RATIONAL(7.0,24.0)* +t22+RATIONAL(-5.0,24.0)*t20+RATIONAL(11.0,24.0)*t25; + coeffs_I->coeff_p3 = t13+RATIONAL(1.0,24.0)*t20+t18*t23+t15*t25; diff --git a/src/GeneralizedPolynomial-Uniform/Hermite/1d.coeffs/1d.cube.order4/coeffs-dx.compute.c b/src/GeneralizedPolynomial-Uniform/Hermite/1d.coeffs/1d.cube.order4/coeffs-dx.compute.c new file mode 100644 index 0000000..bba4a11 --- /dev/null +++ b/src/GeneralizedPolynomial-Uniform/Hermite/1d.coeffs/1d.cube.order4/coeffs-dx.compute.c @@ -0,0 +1,20 @@ +fp t30; +fp t29; +fp t28; +fp t27; + t30 = x*x; + t29 = x*t30; + t28 = t30*t30; + t27 = RATIONAL(-1.0,8.0)*t30; + coeffs_dx->coeff_m2 = RATIONAL(-5.0,24.0)*t28+RATIONAL(1.0,2.0)*t29+ +RATIONAL(1.0,12.0)+t27+RATIONAL(-1.0,4.0)*x; + coeffs_dx->coeff_m1 = RATIONAL(13.0,6.0)*x+RATIONAL(25.0,24.0)*t28+t27+ +RATIONAL(-2.0,3.0)+RATIONAL(-7.0,3.0)*t29; + coeffs_dx->coeff_0 = RATIONAL(5.0,4.0)*t30+RATIONAL(13.0,3.0)*t29+ +RATIONAL(-25.0,6.0)*x+RATIONAL(-25.0,12.0)*t28; + coeffs_dx->coeff_p1 = RATIONAL(2.0,3.0)+RATIONAL(25.0,12.0)*t28+RATIONAL( +-7.0,4.0)*t30+RATIONAL(-4.0,1.0)*t29+RATIONAL(3.0,1.0)*x; + coeffs_dx->coeff_p2 = RATIONAL(-25.0,24.0)*t28+RATIONAL(-1.0,12.0)+ +RATIONAL(11.0,6.0)*t29+RATIONAL(-11.0,12.0)*x+RATIONAL(7.0,8.0)*t30; + coeffs_dx->coeff_p3 = RATIONAL(-1.0,3.0)*t29+RATIONAL(1.0,6.0)*x+RATIONAL +(5.0,24.0)*t28+t27; diff --git a/src/GeneralizedPolynomial-Uniform/Hermite/1d.coeffs/1d.cube.order4/coeffs-dxx.compute.c b/src/GeneralizedPolynomial-Uniform/Hermite/1d.coeffs/1d.cube.order4/coeffs-dxx.compute.c new file mode 100644 index 0000000..cc8dc91 --- /dev/null +++ b/src/GeneralizedPolynomial-Uniform/Hermite/1d.coeffs/1d.cube.order4/coeffs-dxx.compute.c @@ -0,0 +1,21 @@ +fp t37; +fp t36; +fp t35; +fp t34; +fp t33; + t37 = x*x; + t36 = x*t37; + t35 = RATIONAL(-1.0,4.0); + t34 = RATIONAL(-25.0,6.0); + t33 = t35*x; + coeffs_dxx->coeff_m2 = t35+t33+RATIONAL(3.0,2.0)*t37+RATIONAL(-5.0,6.0)* +t36; + coeffs_dxx->coeff_m1 = RATIONAL(13.0,6.0)+t33+RATIONAL(-7.0,1.0)*t37+ +RATIONAL(25.0,6.0)*t36; + coeffs_dxx->coeff_0 = RATIONAL(-25.0,3.0)*t36+RATIONAL(13.0,1.0)*t37+t34+ +RATIONAL(5.0,2.0)*x; + coeffs_dxx->coeff_p1 = RATIONAL(3.0,1.0)+RATIONAL(-7.0,2.0)*x+RATIONAL( +-12.0,1.0)*t37+RATIONAL(25.0,3.0)*t36; + coeffs_dxx->coeff_p2 = RATIONAL(11.0,2.0)*t37+t34*t36+RATIONAL(7.0,4.0)*x ++RATIONAL(-11.0,12.0); + coeffs_dxx->coeff_p3 = RATIONAL(5.0,6.0)*t36+RATIONAL(1.0,6.0)+t33-t37; |