diff options
Diffstat (limited to 'src/GeneralizedPolynomial-Uniform/Lagrange-tensor-product/1d.coeffs/1d.cube.order3.smooth0')
3 files changed, 30 insertions, 0 deletions
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; |