diff options
author | jthorn <jthorn@df1f8a13-aa1d-4dd4-9681-27ded5b42416> | 2002-08-18 15:12:49 +0000 |
---|---|---|
committer | jthorn <jthorn@df1f8a13-aa1d-4dd4-9681-27ded5b42416> | 2002-08-18 15:12:49 +0000 |
commit | 15c98b3abee2525b24ca89932256ff2cda63c9bc (patch) | |
tree | fba7c114f0b5e7ffdad3e60254fd0e8dd1e42f97 /src/GeneralizedPolynomial-Uniform/Lagrange/3d.coeffs/3d.cube.order1.smooth0 | |
parent | 7be0a94bec4b4e89c433e5380fce2cf84e5e0681 (diff) |
This commit reorganizes the Maple scripts which generate the interpolation
coefficients, and the coefficient files themselves, to properly support
multiple interpolation operators, and adds some skeleton support for
Hermite interpolation. I will add full support for Hermite in a future
checkin.
git-svn-id: http://svn.cactuscode.org/arrangements/CactusBase/LocalInterp/trunk@81 df1f8a13-aa1d-4dd4-9681-27ded5b42416
Diffstat (limited to 'src/GeneralizedPolynomial-Uniform/Lagrange/3d.coeffs/3d.cube.order1.smooth0')
4 files changed, 68 insertions, 0 deletions
diff --git a/src/GeneralizedPolynomial-Uniform/Lagrange/3d.coeffs/3d.cube.order1.smooth0/coeff-I.compute.c b/src/GeneralizedPolynomial-Uniform/Lagrange/3d.coeffs/3d.cube.order1.smooth0/coeff-I.compute.c new file mode 100644 index 0000000..666c48b --- /dev/null +++ b/src/GeneralizedPolynomial-Uniform/Lagrange/3d.coeffs/3d.cube.order1.smooth0/coeff-I.compute.c @@ -0,0 +1,32 @@ +fp t7, + t2, + t8, + t4, + t12, + t1, + t6, + t11, + t3, + t10, + t5, + t9; + t7 = RATIONAL(1.0,4.0); + t2 = t7*y; + t8 = RATIONAL(-1.0,4.0); + t4 = t8*x; + t12 = t2+t4; + t1 = t7*z; + t6 = t7*x; + t11 = t1+t6; + t3 = t8*z; + t10 = t3+t7; + t5 = t8*y; + t9 = t4+t5; + coeff_I_0_0_0 = RATIONAL(1.0,2.0)+t3+t9; + coeff_I_p1_0_0 = t5+t6+t10; + coeff_I_0_p1_0 = t10+t12; + coeff_I_p1_p1_0 = t2+t3+t6; + coeff_I_0_0_p1 = t7+t1+t9; + coeff_I_p1_0_p1 = t5+t11; + coeff_I_0_p1_p1 = t1+t12; + coeff_I_p1_p1_p1 = t2+t8+t11; diff --git a/src/GeneralizedPolynomial-Uniform/Lagrange/3d.coeffs/3d.cube.order1.smooth0/coeff-dx.compute.c b/src/GeneralizedPolynomial-Uniform/Lagrange/3d.coeffs/3d.cube.order1.smooth0/coeff-dx.compute.c new file mode 100644 index 0000000..fd79e4a --- /dev/null +++ b/src/GeneralizedPolynomial-Uniform/Lagrange/3d.coeffs/3d.cube.order1.smooth0/coeff-dx.compute.c @@ -0,0 +1,12 @@ +fp t14, + t13; + t14 = RATIONAL(-1.0,4.0); + t13 = RATIONAL(1.0,4.0); + coeff_dx_0_0_0 = t14; + coeff_dx_p1_0_0 = t13; + coeff_dx_0_p1_0 = t14; + coeff_dx_p1_p1_0 = t13; + coeff_dx_0_0_p1 = t14; + coeff_dx_p1_0_p1 = t13; + coeff_dx_0_p1_p1 = t14; + coeff_dx_p1_p1_p1 = t13; diff --git a/src/GeneralizedPolynomial-Uniform/Lagrange/3d.coeffs/3d.cube.order1.smooth0/coeff-dy.compute.c b/src/GeneralizedPolynomial-Uniform/Lagrange/3d.coeffs/3d.cube.order1.smooth0/coeff-dy.compute.c new file mode 100644 index 0000000..2638b34 --- /dev/null +++ b/src/GeneralizedPolynomial-Uniform/Lagrange/3d.coeffs/3d.cube.order1.smooth0/coeff-dy.compute.c @@ -0,0 +1,12 @@ +fp t16, + t15; + t16 = RATIONAL(-1.0,4.0); + t15 = RATIONAL(1.0,4.0); + coeff_dy_0_0_0 = t16; + coeff_dy_p1_0_0 = t16; + coeff_dy_0_p1_0 = t15; + coeff_dy_p1_p1_0 = t15; + coeff_dy_0_0_p1 = t16; + coeff_dy_p1_0_p1 = t16; + coeff_dy_0_p1_p1 = t15; + coeff_dy_p1_p1_p1 = t15; diff --git a/src/GeneralizedPolynomial-Uniform/Lagrange/3d.coeffs/3d.cube.order1.smooth0/coeff-dz.compute.c b/src/GeneralizedPolynomial-Uniform/Lagrange/3d.coeffs/3d.cube.order1.smooth0/coeff-dz.compute.c new file mode 100644 index 0000000..62d5bec --- /dev/null +++ b/src/GeneralizedPolynomial-Uniform/Lagrange/3d.coeffs/3d.cube.order1.smooth0/coeff-dz.compute.c @@ -0,0 +1,12 @@ +fp t18, + t17; + t18 = RATIONAL(-1.0,4.0); + t17 = RATIONAL(1.0,4.0); + coeff_dz_0_0_0 = t18; + coeff_dz_p1_0_0 = t18; + coeff_dz_0_p1_0 = t18; + coeff_dz_p1_p1_0 = t18; + coeff_dz_0_0_p1 = t17; + coeff_dz_p1_0_p1 = t17; + coeff_dz_0_p1_p1 = t17; + coeff_dz_p1_p1_p1 = t17; |