diff options
author | jthorn <jthorn@df1f8a13-aa1d-4dd4-9681-27ded5b42416> | 2003-04-15 17:12:59 +0000 |
---|---|---|
committer | jthorn <jthorn@df1f8a13-aa1d-4dd4-9681-27ded5b42416> | 2003-04-15 17:12:59 +0000 |
commit | 96309eb56861333d6e3111c5a3cb54c9426a1d61 (patch) | |
tree | 4ce9825883ea56e3a4f16bd910974de4a4c48044 /src/GeneralizedPolynomial-Uniform/Lagrange-tensor-product/1d.coeffs | |
parent | 5a8da73fce7752d0e3ee52fe8fb1149009b3305c (diff) |
new machine-generated coefficient files for this interpolator
(with Lagrange now split into tensor-product and maximum-degree
interpolators)
git-svn-id: http://svn.cactuscode.org/arrangements/CactusBase/LocalInterp/trunk@144 df1f8a13-aa1d-4dd4-9681-27ded5b42416
Diffstat (limited to 'src/GeneralizedPolynomial-Uniform/Lagrange-tensor-product/1d.coeffs')
17 files changed, 268 insertions, 0 deletions
diff --git a/src/GeneralizedPolynomial-Uniform/Lagrange-tensor-product/1d.coeffs/1d.cube.order1.smooth0/coeffs-I.compute.c b/src/GeneralizedPolynomial-Uniform/Lagrange-tensor-product/1d.coeffs/1d.cube.order1.smooth0/coeffs-I.compute.c new file mode 100644 index 0000000..c8a5114 --- /dev/null +++ b/src/GeneralizedPolynomial-Uniform/Lagrange-tensor-product/1d.coeffs/1d.cube.order1.smooth0/coeffs-I.compute.c @@ -0,0 +1,2 @@ + coeffs_I->coeff_0 = RATIONAL(1.0,1.0)-x; + coeffs_I->coeff_p1 = x; diff --git a/src/GeneralizedPolynomial-Uniform/Lagrange-tensor-product/1d.coeffs/1d.cube.order1.smooth0/coeffs-dx.compute.c b/src/GeneralizedPolynomial-Uniform/Lagrange-tensor-product/1d.coeffs/1d.cube.order1.smooth0/coeffs-dx.compute.c new file mode 100644 index 0000000..583278a --- /dev/null +++ b/src/GeneralizedPolynomial-Uniform/Lagrange-tensor-product/1d.coeffs/1d.cube.order1.smooth0/coeffs-dx.compute.c @@ -0,0 +1,2 @@ + coeffs_dx->coeff_0 = RATIONAL(-1.0,1.0); + coeffs_dx->coeff_p1 = RATIONAL(1.0,1.0); diff --git a/src/GeneralizedPolynomial-Uniform/Lagrange-tensor-product/1d.coeffs/1d.cube.order2.smooth0/coeffs-I.compute.c b/src/GeneralizedPolynomial-Uniform/Lagrange-tensor-product/1d.coeffs/1d.cube.order2.smooth0/coeffs-I.compute.c new file mode 100644 index 0000000..6e17264 --- /dev/null +++ b/src/GeneralizedPolynomial-Uniform/Lagrange-tensor-product/1d.coeffs/1d.cube.order2.smooth0/coeffs-I.compute.c @@ -0,0 +1,9 @@ +fp t3; +fp t2; +fp t1; + t3 = x*x; + t2 = RATIONAL(1.0,2.0); + t1 = t2*t3; + coeffs_I->coeff_m1 = RATIONAL(-1.0,2.0)*x+t1; + coeffs_I->coeff_0 = RATIONAL(1.0,1.0)-t3; + coeffs_I->coeff_p1 = t2*x+t1; diff --git a/src/GeneralizedPolynomial-Uniform/Lagrange-tensor-product/1d.coeffs/1d.cube.order2.smooth0/coeffs-dx.compute.c b/src/GeneralizedPolynomial-Uniform/Lagrange-tensor-product/1d.coeffs/1d.cube.order2.smooth0/coeffs-dx.compute.c new file mode 100644 index 0000000..5647f75 --- /dev/null +++ b/src/GeneralizedPolynomial-Uniform/Lagrange-tensor-product/1d.coeffs/1d.cube.order2.smooth0/coeffs-dx.compute.c @@ -0,0 +1,3 @@ + coeffs_dx->coeff_m1 = x+RATIONAL(-1.0,2.0); + coeffs_dx->coeff_0 = RATIONAL(-2.0,1.0)*x; + coeffs_dx->coeff_p1 = RATIONAL(1.0,2.0)+x; diff --git a/src/GeneralizedPolynomial-Uniform/Lagrange-tensor-product/1d.coeffs/1d.cube.order2.smooth0/coeffs-dxx.compute.c b/src/GeneralizedPolynomial-Uniform/Lagrange-tensor-product/1d.coeffs/1d.cube.order2.smooth0/coeffs-dxx.compute.c new file mode 100644 index 0000000..19689e7 --- /dev/null +++ b/src/GeneralizedPolynomial-Uniform/Lagrange-tensor-product/1d.coeffs/1d.cube.order2.smooth0/coeffs-dxx.compute.c @@ -0,0 +1,5 @@ +fp t4; + t4 = RATIONAL(1.0,1.0); + coeffs_dxx->coeff_m1 = t4; + coeffs_dxx->coeff_0 = RATIONAL(-2.0,1.0); + coeffs_dxx->coeff_p1 = t4; diff --git a/src/GeneralizedPolynomial-Uniform/Lagrange-tensor-product/1d.coeffs/1d.cube.order3.smooth0/coeffs-I.compute.c b/src/GeneralizedPolynomial-Uniform/Lagrange-tensor-product/1d.coeffs/1d.cube.order3.smooth0/coeffs-I.compute.c new file mode 100644 index 0000000..6ee0b5a --- /dev/null +++ b/src/GeneralizedPolynomial-Uniform/Lagrange-tensor-product/1d.coeffs/1d.cube.order3.smooth0/coeffs-I.compute.c @@ -0,0 +1,16 @@ +fp t10; +fp t9; +fp t8; +fp t7; +fp t6; +fp t5; + t10 = x*x; + t9 = x*t10; + t8 = RATIONAL(-1.0,2.0); + t7 = RATIONAL(1.0,2.0); + t6 = RATIONAL(-1.0,6.0); + t5 = t7*t10; + coeffs_I->coeff_m1 = RATIONAL(-1.0,3.0)*x+t5+t6*t9; + coeffs_I->coeff_0 = RATIONAL(1.0,1.0)+t8*x-t10+t7*t9; + coeffs_I->coeff_p1 = x+t5+t8*t9; + coeffs_I->coeff_p2 = t6*x+RATIONAL(1.0,6.0)*t9; diff --git a/src/GeneralizedPolynomial-Uniform/Lagrange-tensor-product/1d.coeffs/1d.cube.order3.smooth0/coeffs-dx.compute.c b/src/GeneralizedPolynomial-Uniform/Lagrange-tensor-product/1d.coeffs/1d.cube.order3.smooth0/coeffs-dx.compute.c new file mode 100644 index 0000000..1cf3fc6 --- /dev/null +++ b/src/GeneralizedPolynomial-Uniform/Lagrange-tensor-product/1d.coeffs/1d.cube.order3.smooth0/coeffs-dx.compute.c @@ -0,0 +1,8 @@ +fp t13; +fp t12; + t13 = x*x; + t12 = RATIONAL(-1.0,2.0); + coeffs_dx->coeff_m1 = x+t12*t13+RATIONAL(-1.0,3.0); + coeffs_dx->coeff_0 = t12+RATIONAL(-2.0,1.0)*x+RATIONAL(3.0,2.0)*t13; + coeffs_dx->coeff_p1 = x+RATIONAL(1.0,1.0)+RATIONAL(-3.0,2.0)*t13; + coeffs_dx->coeff_p2 = RATIONAL(1.0,2.0)*t13+RATIONAL(-1.0,6.0); diff --git a/src/GeneralizedPolynomial-Uniform/Lagrange-tensor-product/1d.coeffs/1d.cube.order3.smooth0/coeffs-dxx.compute.c b/src/GeneralizedPolynomial-Uniform/Lagrange-tensor-product/1d.coeffs/1d.cube.order3.smooth0/coeffs-dxx.compute.c new file mode 100644 index 0000000..470882c --- /dev/null +++ b/src/GeneralizedPolynomial-Uniform/Lagrange-tensor-product/1d.coeffs/1d.cube.order3.smooth0/coeffs-dxx.compute.c @@ -0,0 +1,6 @@ +fp t14; + t14 = RATIONAL(1.0,1.0); + coeffs_dxx->coeff_m1 = t14-x; + coeffs_dxx->coeff_0 = RATIONAL(-2.0,1.0)+RATIONAL(3.0,1.0)*x; + coeffs_dxx->coeff_p1 = t14+RATIONAL(-3.0,1.0)*x; + coeffs_dxx->coeff_p2 = x; diff --git a/src/GeneralizedPolynomial-Uniform/Lagrange-tensor-product/1d.coeffs/1d.cube.order4.smooth0/coeffs-I.compute.c b/src/GeneralizedPolynomial-Uniform/Lagrange-tensor-product/1d.coeffs/1d.cube.order4.smooth0/coeffs-I.compute.c new file mode 100644 index 0000000..c30526f --- /dev/null +++ b/src/GeneralizedPolynomial-Uniform/Lagrange-tensor-product/1d.coeffs/1d.cube.order4.smooth0/coeffs-I.compute.c @@ -0,0 +1,24 @@ +fp t25; +fp t23; +fp t29; +fp t21; +fp t22; +fp t28; +fp t24; +fp t20; +fp t19; + t25 = x*x; + t23 = t25*t25; + t29 = RATIONAL(1.0,24.0)*t23+RATIONAL(-1.0,24.0)*t25; + t21 = RATIONAL(-1.0,6.0); + t22 = RATIONAL(2.0,3.0); + t28 = t21*t23+t22*t25; + t24 = x*t25; + t20 = RATIONAL(1.0,12.0); + t19 = RATIONAL(-1.0,12.0); + coeffs_I->coeff_m2 = t20*x+t19*t24+t29; + coeffs_I->coeff_m1 = RATIONAL(-2.0,3.0)*x+RATIONAL(1.0,6.0)*t24+t28; + coeffs_I->coeff_0 = RATIONAL(-5.0,4.0)*t25+RATIONAL(1.0,1.0)+RATIONAL(1.0 +,4.0)*t23; + coeffs_I->coeff_p1 = t21*t24+t22*x+t28; + coeffs_I->coeff_p2 = t20*t24+t19*x+t29; diff --git a/src/GeneralizedPolynomial-Uniform/Lagrange-tensor-product/1d.coeffs/1d.cube.order4.smooth0/coeffs-dx.compute.c b/src/GeneralizedPolynomial-Uniform/Lagrange-tensor-product/1d.coeffs/1d.cube.order4.smooth0/coeffs-dx.compute.c new file mode 100644 index 0000000..14b1551 --- /dev/null +++ b/src/GeneralizedPolynomial-Uniform/Lagrange-tensor-product/1d.coeffs/1d.cube.order4.smooth0/coeffs-dx.compute.c @@ -0,0 +1,17 @@ +fp t34; +fp t37; +fp t36; +fp t40; +fp t35; +fp t39; + t34 = RATIONAL(-1.0,12.0); + t37 = x*x; + t36 = x*t37; + t40 = RATIONAL(1.0,6.0)*t36+t34*x; + t35 = RATIONAL(-2.0,3.0); + t39 = t35*t36+RATIONAL(4.0,3.0)*x; + coeffs_dx->coeff_m2 = RATIONAL(-1.0,4.0)*t37+RATIONAL(1.0,12.0)+t40; + coeffs_dx->coeff_m1 = t35+RATIONAL(1.0,2.0)*t37+t39; + coeffs_dx->coeff_0 = t36+RATIONAL(-5.0,2.0)*x; + coeffs_dx->coeff_p1 = RATIONAL(-1.0,2.0)*t37+RATIONAL(2.0,3.0)+t39; + coeffs_dx->coeff_p2 = t34+RATIONAL(1.0,4.0)*t37+t40; diff --git a/src/GeneralizedPolynomial-Uniform/Lagrange-tensor-product/1d.coeffs/1d.cube.order4.smooth0/coeffs-dxx.compute.c b/src/GeneralizedPolynomial-Uniform/Lagrange-tensor-product/1d.coeffs/1d.cube.order4.smooth0/coeffs-dxx.compute.c new file mode 100644 index 0000000..6f78a84 --- /dev/null +++ b/src/GeneralizedPolynomial-Uniform/Lagrange-tensor-product/1d.coeffs/1d.cube.order4.smooth0/coeffs-dxx.compute.c @@ -0,0 +1,13 @@ +fp t46; +fp t48; +fp t44; +fp t47; + t46 = x*x; + t48 = RATIONAL(-2.0,1.0)*t46+RATIONAL(4.0,3.0); + t44 = RATIONAL(1.0,2.0); + t47 = t44*t46+RATIONAL(-1.0,12.0); + coeffs_dxx->coeff_m2 = RATIONAL(-1.0,2.0)*x+t47; + coeffs_dxx->coeff_m1 = x+t48; + coeffs_dxx->coeff_0 = RATIONAL(-5.0,2.0)+RATIONAL(3.0,1.0)*t46; + coeffs_dxx->coeff_p1 = -x+t48; + coeffs_dxx->coeff_p2 = t44*x+t47; diff --git a/src/GeneralizedPolynomial-Uniform/Lagrange-tensor-product/1d.coeffs/1d.cube.order5.smooth0/coeffs-I.compute.c b/src/GeneralizedPolynomial-Uniform/Lagrange-tensor-product/1d.coeffs/1d.cube.order5.smooth0/coeffs-I.compute.c new file mode 100644 index 0000000..2fe57b6 --- /dev/null +++ b/src/GeneralizedPolynomial-Uniform/Lagrange-tensor-product/1d.coeffs/1d.cube.order5.smooth0/coeffs-I.compute.c @@ -0,0 +1,28 @@ +fp t54; +fp t55; +fp t57; +fp t58; +fp t64; +fp t63; +fp t60; +fp t59; +fp t51; + t54 = RATIONAL(1.0,24.0); + t55 = RATIONAL(-1.0,24.0); + t57 = x*x; + t58 = t57*t57; + t64 = t54*t58+t55*t57; + t63 = RATIONAL(-1.0,6.0)*t58+RATIONAL(2.0,3.0)*t57; + t60 = x*t57; + t59 = t60*t57; + t51 = t55*t60; + coeffs_I->coeff_m2 = RATIONAL(1.0,20.0)*x+t51+RATIONAL(-1.0,120.0)*t59+ +t64; + coeffs_I->coeff_m1 = RATIONAL(-1.0,2.0)*x+t51+t54*t59+t63; + coeffs_I->coeff_0 = RATIONAL(1.0,1.0)+RATIONAL(-1.0,3.0)*x+RATIONAL(-5.0, +4.0)*t57+RATIONAL(5.0,12.0)*t60+RATIONAL(1.0,4.0)*t58+RATIONAL(-1.0,12.0)*t59; + coeffs_I->coeff_p1 = x+RATIONAL(-7.0,12.0)*t60+RATIONAL(1.0,12.0)*t59+t63 +; + coeffs_I->coeff_p2 = RATIONAL(-1.0,4.0)*x+RATIONAL(7.0,24.0)*t60+t55*t59+ +t64; + coeffs_I->coeff_p3 = RATIONAL(1.0,30.0)*x+RATIONAL(1.0,120.0)*t59+t51; diff --git a/src/GeneralizedPolynomial-Uniform/Lagrange-tensor-product/1d.coeffs/1d.cube.order5.smooth0/coeffs-dx.compute.c b/src/GeneralizedPolynomial-Uniform/Lagrange-tensor-product/1d.coeffs/1d.cube.order5.smooth0/coeffs-dx.compute.c new file mode 100644 index 0000000..0fc33a0 --- /dev/null +++ b/src/GeneralizedPolynomial-Uniform/Lagrange-tensor-product/1d.coeffs/1d.cube.order5.smooth0/coeffs-dx.compute.c @@ -0,0 +1,21 @@ +fp t71; +fp t70; +fp t76; +fp t75; +fp t72; +fp t67; + t71 = x*x; + t70 = t71*x; + t76 = RATIONAL(-2.0,3.0)*t70+RATIONAL(4.0,3.0)*x; + t75 = RATIONAL(-1.0,12.0)*x+RATIONAL(1.0,6.0)*t70; + t72 = t71*t71; + t67 = RATIONAL(-1.0,8.0)*t71; + coeffs_dx->coeff_m2 = RATIONAL(-1.0,24.0)*t72+RATIONAL(1.0,20.0)+t67+t75; + coeffs_dx->coeff_m1 = t67+RATIONAL(-1.0,2.0)+RATIONAL(5.0,24.0)*t72+t76; + coeffs_dx->coeff_0 = RATIONAL(-5.0,2.0)*x+t70+RATIONAL(-5.0,12.0)*t72+ +RATIONAL(-1.0,3.0)+RATIONAL(5.0,4.0)*t71; + coeffs_dx->coeff_p1 = RATIONAL(5.0,12.0)*t72+RATIONAL(-7.0,4.0)*t71+ +RATIONAL(1.0,1.0)+t76; + coeffs_dx->coeff_p2 = RATIONAL(7.0,8.0)*t71+RATIONAL(-1.0,4.0)+RATIONAL( +-5.0,24.0)*t72+t75; + coeffs_dx->coeff_p3 = RATIONAL(1.0,24.0)*t72+t67+RATIONAL(1.0,30.0); diff --git a/src/GeneralizedPolynomial-Uniform/Lagrange-tensor-product/1d.coeffs/1d.cube.order5.smooth0/coeffs-dxx.compute.c b/src/GeneralizedPolynomial-Uniform/Lagrange-tensor-product/1d.coeffs/1d.cube.order5.smooth0/coeffs-dxx.compute.c new file mode 100644 index 0000000..2382d07 --- /dev/null +++ b/src/GeneralizedPolynomial-Uniform/Lagrange-tensor-product/1d.coeffs/1d.cube.order5.smooth0/coeffs-dxx.compute.c @@ -0,0 +1,17 @@ +fp t83; +fp t86; +fp t85; +fp t82; +fp t78; + t83 = x*x; + t86 = RATIONAL(1.0,2.0)*t83+RATIONAL(-1.0,12.0); + t85 = RATIONAL(-2.0,1.0)*t83+RATIONAL(4.0,3.0); + t82 = x*t83; + t78 = RATIONAL(-1.0,4.0)*x; + coeffs_dxx->coeff_m2 = RATIONAL(-1.0,6.0)*t82+t78+t86; + coeffs_dxx->coeff_m1 = RATIONAL(5.0,6.0)*t82+t78+t85; + coeffs_dxx->coeff_0 = RATIONAL(-5.0,3.0)*t82+RATIONAL(5.0,2.0)*x+RATIONAL +(3.0,1.0)*t83+RATIONAL(-5.0,2.0); + coeffs_dxx->coeff_p1 = RATIONAL(5.0,3.0)*t82+RATIONAL(-7.0,2.0)*x+t85; + coeffs_dxx->coeff_p2 = RATIONAL(-5.0,6.0)*t82+RATIONAL(7.0,4.0)*x+t86; + coeffs_dxx->coeff_p3 = t78+RATIONAL(1.0,6.0)*t82; diff --git a/src/GeneralizedPolynomial-Uniform/Lagrange-tensor-product/1d.coeffs/1d.cube.order6.smooth0/coeffs-I.compute.c b/src/GeneralizedPolynomial-Uniform/Lagrange-tensor-product/1d.coeffs/1d.cube.order6.smooth0/coeffs-I.compute.c new file mode 100644 index 0000000..28b9679 --- /dev/null +++ b/src/GeneralizedPolynomial-Uniform/Lagrange-tensor-product/1d.coeffs/1d.cube.order6.smooth0/coeffs-I.compute.c @@ -0,0 +1,41 @@ +fp t103; +fp t106; +fp t107; +fp t105; +fp t113; +fp t112; +fp t100; +fp t96; +fp t101; +fp t111; +fp t104; +fp t99; +fp t98; +fp t97; + t103 = x*x; + t106 = t103*t103; + t107 = x*t103; + t105 = t107*t107; + t113 = RATIONAL(1.0,720.0)*t105+RATIONAL(1.0,180.0)*t103+RATIONAL(-1.0, +144.0)*t106; + t112 = RATIONAL(-1.0,120.0)*t105+RATIONAL(-3.0,40.0)*t103+RATIONAL(1.0, +12.0)*t106; + t100 = RATIONAL(1.0,48.0); + t96 = RATIONAL(-13.0,48.0); + t101 = RATIONAL(3.0,4.0); + t111 = t100*t105+t96*t106+t101*t103; + t104 = t107*t103; + t99 = RATIONAL(-1.0,60.0); + t98 = RATIONAL(1.0,60.0); + t97 = RATIONAL(-1.0,48.0); + coeffs_I->coeff_m3 = t99*x+t100*t107+RATIONAL(-1.0,240.0)*t104+t113; + coeffs_I->coeff_m2 = RATIONAL(3.0,20.0)*x+RATIONAL(-1.0,6.0)*t107+t98* +t104+t112; + coeffs_I->coeff_m1 = RATIONAL(-3.0,4.0)*x+RATIONAL(13.0,48.0)*t107+t97* +t104+t111; + coeffs_I->coeff_0 = RATIONAL(7.0,18.0)*t106+RATIONAL(-49.0,36.0)*t103+ +RATIONAL(1.0,1.0)+RATIONAL(-1.0,36.0)*t105; + coeffs_I->coeff_p1 = t96*t107+t101*x+t100*t104+t111; + coeffs_I->coeff_p2 = RATIONAL(1.0,6.0)*t107+RATIONAL(-3.0,20.0)*x+t99* +t104+t112; + coeffs_I->coeff_p3 = t97*t107+t98*x+RATIONAL(1.0,240.0)*t104+t113; diff --git a/src/GeneralizedPolynomial-Uniform/Lagrange-tensor-product/1d.coeffs/1d.cube.order6.smooth0/coeffs-dx.compute.c b/src/GeneralizedPolynomial-Uniform/Lagrange-tensor-product/1d.coeffs/1d.cube.order6.smooth0/coeffs-dx.compute.c new file mode 100644 index 0000000..d4d1004 --- /dev/null +++ b/src/GeneralizedPolynomial-Uniform/Lagrange-tensor-product/1d.coeffs/1d.cube.order6.smooth0/coeffs-dx.compute.c @@ -0,0 +1,32 @@ +fp t123; +fp t125; +fp t128; +fp t126; +fp t133; +fp t132; +fp t131; +fp t127; + t123 = RATIONAL(-3.0,20.0); + t125 = x*x; + t128 = x*t125; + t126 = t128*t125; + t133 = t123*x+RATIONAL(-1.0,20.0)*t126+RATIONAL(1.0,3.0)*t128; + t132 = RATIONAL(1.0,120.0)*t126+RATIONAL(1.0,90.0)*x+RATIONAL(-1.0,36.0)* +t128; + t131 = RATIONAL(3.0,2.0)*x+RATIONAL(-13.0,12.0)*t128+RATIONAL(1.0,8.0)* +t126; + t127 = t125*t125; + coeffs_dx->coeff_m3 = RATIONAL(-1.0,48.0)*t127+RATIONAL(-1.0,60.0)+ +RATIONAL(1.0,16.0)*t125+t132; + coeffs_dx->coeff_m2 = RATIONAL(-1.0,2.0)*t125+RATIONAL(3.0,20.0)+RATIONAL +(1.0,12.0)*t127+t133; + coeffs_dx->coeff_m1 = RATIONAL(-5.0,48.0)*t127+RATIONAL(13.0,16.0)*t125+ +RATIONAL(-3.0,4.0)+t131; + coeffs_dx->coeff_0 = RATIONAL(-1.0,6.0)*t126+RATIONAL(14.0,9.0)*t128+ +RATIONAL(-49.0,18.0)*x; + coeffs_dx->coeff_p1 = RATIONAL(5.0,48.0)*t127+RATIONAL(3.0,4.0)+RATIONAL( +-13.0,16.0)*t125+t131; + coeffs_dx->coeff_p2 = RATIONAL(1.0,2.0)*t125+RATIONAL(-1.0,12.0)*t127+ +t123+t133; + coeffs_dx->coeff_p3 = RATIONAL(-1.0,16.0)*t125+RATIONAL(1.0,60.0)+ +RATIONAL(1.0,48.0)*t127+t132; diff --git a/src/GeneralizedPolynomial-Uniform/Lagrange-tensor-product/1d.coeffs/1d.cube.order6.smooth0/coeffs-dxx.compute.c b/src/GeneralizedPolynomial-Uniform/Lagrange-tensor-product/1d.coeffs/1d.cube.order6.smooth0/coeffs-dxx.compute.c new file mode 100644 index 0000000..02aa677 --- /dev/null +++ b/src/GeneralizedPolynomial-Uniform/Lagrange-tensor-product/1d.coeffs/1d.cube.order6.smooth0/coeffs-dxx.compute.c @@ -0,0 +1,24 @@ +fp t144; +fp t145; +fp t150; +fp t149; +fp t141; +fp t148; +fp t143; + t144 = x*x; + t145 = t144*t144; + t150 = RATIONAL(-13.0,4.0)*t144+RATIONAL(3.0,2.0)+RATIONAL(5.0,8.0)*t145; + t149 = t144+RATIONAL(-3.0,20.0)+RATIONAL(-1.0,4.0)*t145; + t141 = RATIONAL(-1.0,12.0); + t148 = RATIONAL(1.0,90.0)+RATIONAL(1.0,24.0)*t145+t141*t144; + t143 = x*t144; + coeffs_dxx->coeff_m3 = RATIONAL(1.0,8.0)*x+t141*t143+t148; + coeffs_dxx->coeff_m2 = -x+RATIONAL(1.0,3.0)*t143+t149; + coeffs_dxx->coeff_m1 = RATIONAL(-5.0,12.0)*t143+RATIONAL(13.0,8.0)*x+t150 +; + coeffs_dxx->coeff_0 = RATIONAL(14.0,3.0)*t144+RATIONAL(-5.0,6.0)*t145+ +RATIONAL(-49.0,18.0); + coeffs_dxx->coeff_p1 = RATIONAL(5.0,12.0)*t143+RATIONAL(-13.0,8.0)*x+t150 +; + coeffs_dxx->coeff_p2 = x+RATIONAL(-1.0,3.0)*t143+t149; + coeffs_dxx->coeff_p3 = RATIONAL(-1.0,8.0)*x+RATIONAL(1.0,12.0)*t143+t148; |