diff options
Diffstat (limited to 'src/GeneralizedPolynomial-Uniform/Lagrange/1d.coeffs')
34 files changed, 268 insertions, 259 deletions
diff --git a/src/GeneralizedPolynomial-Uniform/Lagrange/1d.coeffs/1d.cube.order1.smooth0/coeff-I.compute.c b/src/GeneralizedPolynomial-Uniform/Lagrange/1d.coeffs/1d.cube.order1.smooth0/coeff-I.compute.c deleted file mode 100644 index b8ea11e..0000000 --- a/src/GeneralizedPolynomial-Uniform/Lagrange/1d.coeffs/1d.cube.order1.smooth0/coeff-I.compute.c +++ /dev/null @@ -1,2 +0,0 @@ - coeff_I_0 = RATIONAL(1.0,1.0)-x; - coeff_I_p1 = x; diff --git a/src/GeneralizedPolynomial-Uniform/Lagrange/1d.coeffs/1d.cube.order1.smooth0/coeff-dx.compute.c b/src/GeneralizedPolynomial-Uniform/Lagrange/1d.coeffs/1d.cube.order1.smooth0/coeff-dx.compute.c deleted file mode 100644 index 4f5fb07..0000000 --- a/src/GeneralizedPolynomial-Uniform/Lagrange/1d.coeffs/1d.cube.order1.smooth0/coeff-dx.compute.c +++ /dev/null @@ -1,2 +0,0 @@ - coeff_dx_0 = RATIONAL(-1.0,1.0); - coeff_dx_p1 = RATIONAL(1.0,1.0); diff --git a/src/GeneralizedPolynomial-Uniform/Lagrange/1d.coeffs/1d.cube.order1.smooth0/coeffs-I.compute.c b/src/GeneralizedPolynomial-Uniform/Lagrange/1d.coeffs/1d.cube.order1.smooth0/coeffs-I.compute.c new file mode 100644 index 0000000..c8a5114 --- /dev/null +++ b/src/GeneralizedPolynomial-Uniform/Lagrange/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/1d.coeffs/1d.cube.order1.smooth0/coeffs-dx.compute.c b/src/GeneralizedPolynomial-Uniform/Lagrange/1d.coeffs/1d.cube.order1.smooth0/coeffs-dx.compute.c new file mode 100644 index 0000000..583278a --- /dev/null +++ b/src/GeneralizedPolynomial-Uniform/Lagrange/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/1d.coeffs/1d.cube.order2.smooth0/coeff-I.compute.c b/src/GeneralizedPolynomial-Uniform/Lagrange/1d.coeffs/1d.cube.order2.smooth0/coeff-I.compute.c deleted file mode 100644 index 4d0de20..0000000 --- a/src/GeneralizedPolynomial-Uniform/Lagrange/1d.coeffs/1d.cube.order2.smooth0/coeff-I.compute.c +++ /dev/null @@ -1,9 +0,0 @@ -fp t3, - t2, - t1; - t3 = x*x; - t2 = RATIONAL(1.0,2.0); - t1 = t2*t3; - coeff_I_m1 = RATIONAL(-1.0,2.0)*x+t1; - coeff_I_0 = RATIONAL(1.0,1.0)-t3; - coeff_I_p1 = t2*x+t1; diff --git a/src/GeneralizedPolynomial-Uniform/Lagrange/1d.coeffs/1d.cube.order2.smooth0/coeff-dx.compute.c b/src/GeneralizedPolynomial-Uniform/Lagrange/1d.coeffs/1d.cube.order2.smooth0/coeff-dx.compute.c deleted file mode 100644 index e2532c9..0000000 --- a/src/GeneralizedPolynomial-Uniform/Lagrange/1d.coeffs/1d.cube.order2.smooth0/coeff-dx.compute.c +++ /dev/null @@ -1,3 +0,0 @@ - coeff_dx_m1 = x+RATIONAL(-1.0,2.0); - coeff_dx_0 = RATIONAL(-2.0,1.0)*x; - coeff_dx_p1 = RATIONAL(1.0,2.0)+x; diff --git a/src/GeneralizedPolynomial-Uniform/Lagrange/1d.coeffs/1d.cube.order2.smooth0/coeff-dxx.compute.c b/src/GeneralizedPolynomial-Uniform/Lagrange/1d.coeffs/1d.cube.order2.smooth0/coeff-dxx.compute.c deleted file mode 100644 index 4e6fd95..0000000 --- a/src/GeneralizedPolynomial-Uniform/Lagrange/1d.coeffs/1d.cube.order2.smooth0/coeff-dxx.compute.c +++ /dev/null @@ -1,5 +0,0 @@ -fp t4; - t4 = RATIONAL(1.0,1.0); - coeff_dxx_m1 = t4; - coeff_dxx_0 = RATIONAL(-2.0,1.0); - coeff_dxx_p1 = t4; diff --git a/src/GeneralizedPolynomial-Uniform/Lagrange/1d.coeffs/1d.cube.order2.smooth0/coeffs-I.compute.c b/src/GeneralizedPolynomial-Uniform/Lagrange/1d.coeffs/1d.cube.order2.smooth0/coeffs-I.compute.c new file mode 100644 index 0000000..6e17264 --- /dev/null +++ b/src/GeneralizedPolynomial-Uniform/Lagrange/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/1d.coeffs/1d.cube.order2.smooth0/coeffs-dx.compute.c b/src/GeneralizedPolynomial-Uniform/Lagrange/1d.coeffs/1d.cube.order2.smooth0/coeffs-dx.compute.c new file mode 100644 index 0000000..5647f75 --- /dev/null +++ b/src/GeneralizedPolynomial-Uniform/Lagrange/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/1d.coeffs/1d.cube.order2.smooth0/coeffs-dxx.compute.c b/src/GeneralizedPolynomial-Uniform/Lagrange/1d.coeffs/1d.cube.order2.smooth0/coeffs-dxx.compute.c new file mode 100644 index 0000000..19689e7 --- /dev/null +++ b/src/GeneralizedPolynomial-Uniform/Lagrange/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/1d.coeffs/1d.cube.order3.smooth0/coeff-I.compute.c b/src/GeneralizedPolynomial-Uniform/Lagrange/1d.coeffs/1d.cube.order3.smooth0/coeff-I.compute.c deleted file mode 100644 index 2df6cd1..0000000 --- a/src/GeneralizedPolynomial-Uniform/Lagrange/1d.coeffs/1d.cube.order3.smooth0/coeff-I.compute.c +++ /dev/null @@ -1,16 +0,0 @@ -fp t10, - t9, - t8, - t7, - t6, - 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; - coeff_I_m1 = RATIONAL(-1.0,3.0)*x+t5+t6*t9; - coeff_I_0 = RATIONAL(1.0,1.0)+t8*x-t10+t7*t9; - coeff_I_p1 = x+t5+t8*t9; - coeff_I_p2 = t6*x+RATIONAL(1.0,6.0)*t9; diff --git a/src/GeneralizedPolynomial-Uniform/Lagrange/1d.coeffs/1d.cube.order3.smooth0/coeff-dx.compute.c b/src/GeneralizedPolynomial-Uniform/Lagrange/1d.coeffs/1d.cube.order3.smooth0/coeff-dx.compute.c deleted file mode 100644 index 4306fc4..0000000 --- a/src/GeneralizedPolynomial-Uniform/Lagrange/1d.coeffs/1d.cube.order3.smooth0/coeff-dx.compute.c +++ /dev/null @@ -1,8 +0,0 @@ -fp t13, - t12; - t13 = x*x; - t12 = RATIONAL(-1.0,2.0); - coeff_dx_m1 = x+t12*t13+RATIONAL(-1.0,3.0); - coeff_dx_0 = t12+RATIONAL(-2.0,1.0)*x+RATIONAL(3.0,2.0)*t13; - coeff_dx_p1 = RATIONAL(-3.0,2.0)*t13+RATIONAL(1.0,1.0)+x; - coeff_dx_p2 = RATIONAL(1.0,2.0)*t13+RATIONAL(-1.0,6.0); diff --git a/src/GeneralizedPolynomial-Uniform/Lagrange/1d.coeffs/1d.cube.order3.smooth0/coeff-dxx.compute.c b/src/GeneralizedPolynomial-Uniform/Lagrange/1d.coeffs/1d.cube.order3.smooth0/coeff-dxx.compute.c deleted file mode 100644 index d5e199b..0000000 --- a/src/GeneralizedPolynomial-Uniform/Lagrange/1d.coeffs/1d.cube.order3.smooth0/coeff-dxx.compute.c +++ /dev/null @@ -1,6 +0,0 @@ -fp t14; - t14 = RATIONAL(1.0,1.0); - coeff_dxx_m1 = t14-x; - coeff_dxx_0 = RATIONAL(-2.0,1.0)+RATIONAL(3.0,1.0)*x; - coeff_dxx_p1 = RATIONAL(-3.0,1.0)*x+t14; - coeff_dxx_p2 = x; diff --git a/src/GeneralizedPolynomial-Uniform/Lagrange/1d.coeffs/1d.cube.order3.smooth0/coeffs-I.compute.c b/src/GeneralizedPolynomial-Uniform/Lagrange/1d.coeffs/1d.cube.order3.smooth0/coeffs-I.compute.c new file mode 100644 index 0000000..6ee0b5a --- /dev/null +++ b/src/GeneralizedPolynomial-Uniform/Lagrange/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/1d.coeffs/1d.cube.order3.smooth0/coeffs-dx.compute.c b/src/GeneralizedPolynomial-Uniform/Lagrange/1d.coeffs/1d.cube.order3.smooth0/coeffs-dx.compute.c new file mode 100644 index 0000000..1cf3fc6 --- /dev/null +++ b/src/GeneralizedPolynomial-Uniform/Lagrange/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/1d.coeffs/1d.cube.order3.smooth0/coeffs-dxx.compute.c b/src/GeneralizedPolynomial-Uniform/Lagrange/1d.coeffs/1d.cube.order3.smooth0/coeffs-dxx.compute.c new file mode 100644 index 0000000..470882c --- /dev/null +++ b/src/GeneralizedPolynomial-Uniform/Lagrange/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/1d.coeffs/1d.cube.order4.smooth0/coeff-I.compute.c b/src/GeneralizedPolynomial-Uniform/Lagrange/1d.coeffs/1d.cube.order4.smooth0/coeff-I.compute.c deleted file mode 100644 index 05b6c68..0000000 --- a/src/GeneralizedPolynomial-Uniform/Lagrange/1d.coeffs/1d.cube.order4.smooth0/coeff-I.compute.c +++ /dev/null @@ -1,24 +0,0 @@ -fp t23, - t24, - t29, - t21, - t22, - t28, - t25, - t20, - t19; - t23 = x*x; - t24 = t23*t23; - t29 = RATIONAL(-1.0,24.0)*t23+RATIONAL(1.0,24.0)*t24; - t21 = RATIONAL(2.0,3.0); - t22 = RATIONAL(-1.0,6.0); - t28 = t21*t23+t22*t24; - t25 = x*t23; - t20 = RATIONAL(1.0,12.0); - t19 = RATIONAL(-1.0,12.0); - coeff_I_m2 = t20*x+t19*t25+t29; - coeff_I_m1 = RATIONAL(-2.0,3.0)*x+RATIONAL(1.0,6.0)*t25+t28; - coeff_I_0 = RATIONAL(-5.0,4.0)*t23+RATIONAL(1.0,1.0)+RATIONAL(1.0,4.0)* -t24; - coeff_I_p1 = t22*t25+t21*x+t28; - coeff_I_p2 = t20*t25+t19*x+t29; diff --git a/src/GeneralizedPolynomial-Uniform/Lagrange/1d.coeffs/1d.cube.order4.smooth0/coeff-dx.compute.c b/src/GeneralizedPolynomial-Uniform/Lagrange/1d.coeffs/1d.cube.order4.smooth0/coeff-dx.compute.c deleted file mode 100644 index 5ac1469..0000000 --- a/src/GeneralizedPolynomial-Uniform/Lagrange/1d.coeffs/1d.cube.order4.smooth0/coeff-dx.compute.c +++ /dev/null @@ -1,17 +0,0 @@ -fp t36, - t37, - t34, - t40, - t35, - t39; - t36 = x*x; - t37 = x*t36; - t34 = RATIONAL(-1.0,12.0); - t40 = t34*x+RATIONAL(1.0,6.0)*t37; - t35 = RATIONAL(-2.0,3.0); - t39 = t35*t37+RATIONAL(4.0,3.0)*x; - coeff_dx_m2 = RATIONAL(-1.0,4.0)*t36+RATIONAL(1.0,12.0)+t40; - coeff_dx_m1 = t35+RATIONAL(1.0,2.0)*t36+t39; - coeff_dx_0 = t37+RATIONAL(-5.0,2.0)*x; - coeff_dx_p1 = RATIONAL(2.0,3.0)+RATIONAL(-1.0,2.0)*t36+t39; - coeff_dx_p2 = t34+RATIONAL(1.0,4.0)*t36+t40; diff --git a/src/GeneralizedPolynomial-Uniform/Lagrange/1d.coeffs/1d.cube.order4.smooth0/coeff-dxx.compute.c b/src/GeneralizedPolynomial-Uniform/Lagrange/1d.coeffs/1d.cube.order4.smooth0/coeff-dxx.compute.c deleted file mode 100644 index 56d8567..0000000 --- a/src/GeneralizedPolynomial-Uniform/Lagrange/1d.coeffs/1d.cube.order4.smooth0/coeff-dxx.compute.c +++ /dev/null @@ -1,13 +0,0 @@ -fp t44, - t46, - t48, - t47; - t44 = RATIONAL(1.0,2.0); - t46 = x*x; - t48 = t44*t46+RATIONAL(-1.0,12.0); - t47 = RATIONAL(-2.0,1.0)*t46+RATIONAL(4.0,3.0); - coeff_dxx_m2 = RATIONAL(-1.0,2.0)*x+t48; - coeff_dxx_m1 = x+t47; - coeff_dxx_0 = RATIONAL(-5.0,2.0)+RATIONAL(3.0,1.0)*t46; - coeff_dxx_p1 = -x+t47; - coeff_dxx_p2 = t44*x+t48; diff --git a/src/GeneralizedPolynomial-Uniform/Lagrange/1d.coeffs/1d.cube.order4.smooth0/coeffs-I.compute.c b/src/GeneralizedPolynomial-Uniform/Lagrange/1d.coeffs/1d.cube.order4.smooth0/coeffs-I.compute.c new file mode 100644 index 0000000..c30526f --- /dev/null +++ b/src/GeneralizedPolynomial-Uniform/Lagrange/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/1d.coeffs/1d.cube.order4.smooth0/coeffs-dx.compute.c b/src/GeneralizedPolynomial-Uniform/Lagrange/1d.coeffs/1d.cube.order4.smooth0/coeffs-dx.compute.c new file mode 100644 index 0000000..14b1551 --- /dev/null +++ b/src/GeneralizedPolynomial-Uniform/Lagrange/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/1d.coeffs/1d.cube.order4.smooth0/coeffs-dxx.compute.c b/src/GeneralizedPolynomial-Uniform/Lagrange/1d.coeffs/1d.cube.order4.smooth0/coeffs-dxx.compute.c new file mode 100644 index 0000000..6f78a84 --- /dev/null +++ b/src/GeneralizedPolynomial-Uniform/Lagrange/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/1d.coeffs/1d.cube.order5.smooth0/coeff-I.compute.c b/src/GeneralizedPolynomial-Uniform/Lagrange/1d.coeffs/1d.cube.order5.smooth0/coeff-I.compute.c deleted file mode 100644 index f32ea37..0000000 --- a/src/GeneralizedPolynomial-Uniform/Lagrange/1d.coeffs/1d.cube.order5.smooth0/coeff-I.compute.c +++ /dev/null @@ -1,25 +0,0 @@ -fp t58, - t56, - t64, - t54, - t55, - t63, - t60, - t57, - t50; - t58 = x*x; - t56 = t58*t58; - t64 = RATIONAL(2.0,3.0)*t58+RATIONAL(-1.0,6.0)*t56; - t54 = RATIONAL(1.0,24.0); - t55 = RATIONAL(-1.0,24.0); - t63 = t54*t56+t55*t58; - t60 = t58*x; - t57 = t58*t60; - t50 = t55*t60; - coeff_I_m2 = RATIONAL(1.0,20.0)*x+t50+RATIONAL(-1.0,120.0)*t57+t63; - coeff_I_m1 = RATIONAL(-1.0,2.0)*x+t50+t54*t57+t64; - coeff_I_0 = RATIONAL(1.0,1.0)+RATIONAL(-1.0,3.0)*x+RATIONAL(-5.0,4.0)*t58 -+RATIONAL(5.0,12.0)*t60+RATIONAL(1.0,4.0)*t56+RATIONAL(-1.0,12.0)*t57; - coeff_I_p1 = x+RATIONAL(-7.0,12.0)*t60+RATIONAL(1.0,12.0)*t57+t64; - coeff_I_p2 = RATIONAL(-1.0,4.0)*x+RATIONAL(7.0,24.0)*t60+t55*t57+t63; - coeff_I_p3 = RATIONAL(1.0,30.0)*x+RATIONAL(1.0,120.0)*t57+t50; diff --git a/src/GeneralizedPolynomial-Uniform/Lagrange/1d.coeffs/1d.cube.order5.smooth0/coeff-dx.compute.c b/src/GeneralizedPolynomial-Uniform/Lagrange/1d.coeffs/1d.cube.order5.smooth0/coeff-dx.compute.c deleted file mode 100644 index 8fb00d3..0000000 --- a/src/GeneralizedPolynomial-Uniform/Lagrange/1d.coeffs/1d.cube.order5.smooth0/coeff-dx.compute.c +++ /dev/null @@ -1,21 +0,0 @@ -fp t71, - t72, - t76, - t75, - t70, - t67; - t71 = x*x; - t72 = x*t71; - t76 = RATIONAL(-1.0,12.0)*x+RATIONAL(1.0,6.0)*t72; - t75 = RATIONAL(-2.0,3.0)*t72+RATIONAL(4.0,3.0)*x; - t70 = t71*t71; - t67 = RATIONAL(-1.0,8.0)*t71; - coeff_dx_m2 = t67+RATIONAL(-1.0,24.0)*t70+RATIONAL(1.0,20.0)+t76; - coeff_dx_m1 = RATIONAL(-1.0,2.0)+t67+RATIONAL(5.0,24.0)*t70+t75; - coeff_dx_0 = RATIONAL(5.0,4.0)*t71+t72+RATIONAL(-5.0,12.0)*t70+RATIONAL( --5.0,2.0)*x+RATIONAL(-1.0,3.0); - coeff_dx_p1 = RATIONAL(5.0,12.0)*t70+RATIONAL(-7.0,4.0)*t71+RATIONAL(1.0, -1.0)+t75; - coeff_dx_p2 = RATIONAL(-1.0,4.0)+RATIONAL(7.0,8.0)*t71+RATIONAL(-5.0,24.0 -)*t70+t76; - coeff_dx_p3 = RATIONAL(1.0,24.0)*t70+RATIONAL(1.0,30.0)+t67; diff --git a/src/GeneralizedPolynomial-Uniform/Lagrange/1d.coeffs/1d.cube.order5.smooth0/coeff-dxx.compute.c b/src/GeneralizedPolynomial-Uniform/Lagrange/1d.coeffs/1d.cube.order5.smooth0/coeff-dxx.compute.c deleted file mode 100644 index 4c7a405..0000000 --- a/src/GeneralizedPolynomial-Uniform/Lagrange/1d.coeffs/1d.cube.order5.smooth0/coeff-dxx.compute.c +++ /dev/null @@ -1,17 +0,0 @@ -fp t82, - t86, - t85, - t83, - t79; - t82 = x*x; - t86 = RATIONAL(1.0,2.0)*t82+RATIONAL(-1.0,12.0); - t85 = RATIONAL(-2.0,1.0)*t82+RATIONAL(4.0,3.0); - t83 = t82*x; - t79 = RATIONAL(-1.0,4.0)*x; - coeff_dxx_m2 = t79+RATIONAL(-1.0,6.0)*t83+t86; - coeff_dxx_m1 = t79+RATIONAL(5.0,6.0)*t83+t85; - coeff_dxx_0 = RATIONAL(-5.0,3.0)*t83+RATIONAL(5.0,2.0)*x+RATIONAL(3.0,1.0 -)*t82+RATIONAL(-5.0,2.0); - coeff_dxx_p1 = RATIONAL(5.0,3.0)*t83+RATIONAL(-7.0,2.0)*x+t85; - coeff_dxx_p2 = RATIONAL(7.0,4.0)*x+RATIONAL(-5.0,6.0)*t83+t86; - coeff_dxx_p3 = t79+RATIONAL(1.0,6.0)*t83; diff --git a/src/GeneralizedPolynomial-Uniform/Lagrange/1d.coeffs/1d.cube.order5.smooth0/coeffs-I.compute.c b/src/GeneralizedPolynomial-Uniform/Lagrange/1d.coeffs/1d.cube.order5.smooth0/coeffs-I.compute.c new file mode 100644 index 0000000..2fe57b6 --- /dev/null +++ b/src/GeneralizedPolynomial-Uniform/Lagrange/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/1d.coeffs/1d.cube.order5.smooth0/coeffs-dx.compute.c b/src/GeneralizedPolynomial-Uniform/Lagrange/1d.coeffs/1d.cube.order5.smooth0/coeffs-dx.compute.c new file mode 100644 index 0000000..0fc33a0 --- /dev/null +++ b/src/GeneralizedPolynomial-Uniform/Lagrange/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/1d.coeffs/1d.cube.order5.smooth0/coeffs-dxx.compute.c b/src/GeneralizedPolynomial-Uniform/Lagrange/1d.coeffs/1d.cube.order5.smooth0/coeffs-dxx.compute.c new file mode 100644 index 0000000..2382d07 --- /dev/null +++ b/src/GeneralizedPolynomial-Uniform/Lagrange/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/1d.coeffs/1d.cube.order6.smooth0/coeff-I.compute.c b/src/GeneralizedPolynomial-Uniform/Lagrange/1d.coeffs/1d.cube.order6.smooth0/coeff-I.compute.c deleted file mode 100644 index b189aa6..0000000 --- a/src/GeneralizedPolynomial-Uniform/Lagrange/1d.coeffs/1d.cube.order6.smooth0/coeff-I.compute.c +++ /dev/null @@ -1,38 +0,0 @@ -fp t103, - t107, - t106, - t102, - t113, - t112, - t96, - t101, - t100, - t111, - t105, - t99, - t98, - t97; - t103 = x*x; - t107 = x*t103; - t106 = t107*t107; - t102 = t103*t103; - t113 = RATIONAL(-1.0,144.0)*t102+RATIONAL(1.0,720.0)*t106+RATIONAL(1.0, -180.0)*t103; - t112 = RATIONAL(-1.0,120.0)*t106+RATIONAL(-3.0,40.0)*t103+RATIONAL(1.0, -12.0)*t102; - t96 = RATIONAL(-13.0,48.0); - t101 = RATIONAL(3.0,4.0); - t100 = RATIONAL(1.0,48.0); - t111 = t96*t102+t100*t106+t101*t103; - t105 = t107*t103; - t99 = RATIONAL(-1.0,60.0); - t98 = RATIONAL(1.0,60.0); - t97 = RATIONAL(-1.0,48.0); - coeff_I_m3 = t99*x+t100*t107+RATIONAL(-1.0,240.0)*t105+t113; - coeff_I_m2 = RATIONAL(3.0,20.0)*x+RATIONAL(-1.0,6.0)*t107+t98*t105+t112; - coeff_I_m1 = RATIONAL(-3.0,4.0)*x+RATIONAL(13.0,48.0)*t107+t97*t105+t111; - coeff_I_0 = RATIONAL(7.0,18.0)*t102+RATIONAL(-49.0,36.0)*t103+RATIONAL( -1.0,1.0)+RATIONAL(-1.0,36.0)*t106; - coeff_I_p1 = t96*t107+t101*x+t100*t105+t111; - coeff_I_p2 = RATIONAL(1.0,6.0)*t107+RATIONAL(-3.0,20.0)*x+t99*t105+t112; - coeff_I_p3 = t97*t107+t98*x+RATIONAL(1.0,240.0)*t105+t113; diff --git a/src/GeneralizedPolynomial-Uniform/Lagrange/1d.coeffs/1d.cube.order6.smooth0/coeff-dx.compute.c b/src/GeneralizedPolynomial-Uniform/Lagrange/1d.coeffs/1d.cube.order6.smooth0/coeff-dx.compute.c deleted file mode 100644 index ab1b33f..0000000 --- a/src/GeneralizedPolynomial-Uniform/Lagrange/1d.coeffs/1d.cube.order6.smooth0/coeff-dx.compute.c +++ /dev/null @@ -1,31 +0,0 @@ -fp t123, - t125, - t126, - t129, - t127, - t133, - t132, - t131; - t123 = RATIONAL(-3.0,20.0); - t125 = x*x; - t126 = x*t125; - t129 = t125*t125; - t127 = x*t129; - t133 = RATIONAL(-1.0,20.0)*t127+t123*x+RATIONAL(1.0,3.0)*t126; - t132 = RATIONAL(1.0,90.0)*x+RATIONAL(-1.0,36.0)*t126+RATIONAL(1.0,120.0)* -t127; - t131 = RATIONAL(-13.0,12.0)*t126+RATIONAL(3.0,2.0)*x+RATIONAL(1.0,8.0)* -t127; - coeff_dx_m3 = RATIONAL(-1.0,60.0)+RATIONAL(1.0,16.0)*t125+RATIONAL(-1.0, -48.0)*t129+t132; - coeff_dx_m2 = RATIONAL(3.0,20.0)+RATIONAL(1.0,12.0)*t129+RATIONAL(-1.0, -2.0)*t125+t133; - coeff_dx_m1 = RATIONAL(-3.0,4.0)+RATIONAL(13.0,16.0)*t125+RATIONAL(-5.0, -48.0)*t129+t131; - coeff_dx_0 = RATIONAL(-1.0,6.0)*t127+RATIONAL(14.0,9.0)*t126+RATIONAL( --49.0,18.0)*x; - coeff_dx_p1 = RATIONAL(5.0,48.0)*t129+RATIONAL(3.0,4.0)+RATIONAL(-13.0, -16.0)*t125+t131; - coeff_dx_p2 = t123+RATIONAL(1.0,2.0)*t125+RATIONAL(-1.0,12.0)*t129+t133; - coeff_dx_p3 = RATIONAL(1.0,60.0)+RATIONAL(1.0,48.0)*t129+RATIONAL(-1.0, -16.0)*t125+t132; diff --git a/src/GeneralizedPolynomial-Uniform/Lagrange/1d.coeffs/1d.cube.order6.smooth0/coeff-dxx.compute.c b/src/GeneralizedPolynomial-Uniform/Lagrange/1d.coeffs/1d.cube.order6.smooth0/coeff-dxx.compute.c deleted file mode 100644 index 46979da..0000000 --- a/src/GeneralizedPolynomial-Uniform/Lagrange/1d.coeffs/1d.cube.order6.smooth0/coeff-dxx.compute.c +++ /dev/null @@ -1,22 +0,0 @@ -fp t144, - t143, - t141, - t150, - t149, - t148, - t145; - t144 = x*x; - t143 = t144*t144; - t141 = RATIONAL(-1.0,12.0); - t150 = RATIONAL(1.0,90.0)+RATIONAL(1.0,24.0)*t143+t141*t144; - t149 = RATIONAL(-3.0,20.0)+RATIONAL(-1.0,4.0)*t143+t144; - t148 = RATIONAL(3.0,2.0)+RATIONAL(-13.0,4.0)*t144+RATIONAL(5.0,8.0)*t143; - t145 = x*t144; - coeff_dxx_m3 = RATIONAL(1.0,8.0)*x+t141*t145+t150; - coeff_dxx_m2 = RATIONAL(1.0,3.0)*t145-x+t149; - coeff_dxx_m1 = RATIONAL(-5.0,12.0)*t145+RATIONAL(13.0,8.0)*x+t148; - coeff_dxx_0 = RATIONAL(14.0,3.0)*t144+RATIONAL(-5.0,6.0)*t143+RATIONAL( --49.0,18.0); - coeff_dxx_p1 = RATIONAL(5.0,12.0)*t145+RATIONAL(-13.0,8.0)*x+t148; - coeff_dxx_p2 = RATIONAL(-1.0,3.0)*t145+x+t149; - coeff_dxx_p3 = RATIONAL(-1.0,8.0)*x+RATIONAL(1.0,12.0)*t145+t150; diff --git a/src/GeneralizedPolynomial-Uniform/Lagrange/1d.coeffs/1d.cube.order6.smooth0/coeffs-I.compute.c b/src/GeneralizedPolynomial-Uniform/Lagrange/1d.coeffs/1d.cube.order6.smooth0/coeffs-I.compute.c new file mode 100644 index 0000000..28b9679 --- /dev/null +++ b/src/GeneralizedPolynomial-Uniform/Lagrange/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/1d.coeffs/1d.cube.order6.smooth0/coeffs-dx.compute.c b/src/GeneralizedPolynomial-Uniform/Lagrange/1d.coeffs/1d.cube.order6.smooth0/coeffs-dx.compute.c new file mode 100644 index 0000000..d4d1004 --- /dev/null +++ b/src/GeneralizedPolynomial-Uniform/Lagrange/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/1d.coeffs/1d.cube.order6.smooth0/coeffs-dxx.compute.c b/src/GeneralizedPolynomial-Uniform/Lagrange/1d.coeffs/1d.cube.order6.smooth0/coeffs-dxx.compute.c new file mode 100644 index 0000000..02aa677 --- /dev/null +++ b/src/GeneralizedPolynomial-Uniform/Lagrange/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; |