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.order2.smooth0/coeff-dx.compute.c | |
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.order2.smooth0/coeff-dx.compute.c')
-rw-r--r-- | src/GeneralizedPolynomial-Uniform/Lagrange/3d.coeffs/3d.cube.order2.smooth0/coeff-dx.compute.c | 67 |
1 files changed, 67 insertions, 0 deletions
diff --git a/src/GeneralizedPolynomial-Uniform/Lagrange/3d.coeffs/3d.cube.order2.smooth0/coeff-dx.compute.c b/src/GeneralizedPolynomial-Uniform/Lagrange/3d.coeffs/3d.cube.order2.smooth0/coeff-dx.compute.c new file mode 100644 index 0000000..758757f --- /dev/null +++ b/src/GeneralizedPolynomial-Uniform/Lagrange/3d.coeffs/3d.cube.order2.smooth0/coeff-dx.compute.c @@ -0,0 +1,67 @@ +fp t85, + t90, + t100, + t88, + t99, + t89, + t81, + t98, + t84, + t97, + t87, + t83, + t96, + t82, + t95, + t94, + t93, + t92, + t91, + t86; + t85 = RATIONAL(1.0,9.0)*x; + t90 = RATIONAL(-1.0,18.0); + t100 = t85+t90; + t88 = RATIONAL(1.0,18.0); + t99 = t85+t88; + t89 = RATIONAL(1.0,12.0); + t81 = t89*y; + t98 = t81+t85; + t84 = t89*z; + t97 = t84+t85; + t87 = RATIONAL(-1.0,12.0); + t83 = t87*y; + t96 = t83+t100; + t82 = t87*z; + t95 = t82+t100; + t94 = t83+t99; + t93 = t88+t98; + t92 = t88+t97; + t91 = t81+t97; + t86 = RATIONAL(-2.0,9.0)*x; + coeff_dx_m1_m1_m1 = t90+t91; + coeff_dx_0_m1_m1 = t86; + coeff_dx_p1_m1_m1 = t82+t94; + coeff_dx_m1_0_m1 = t90+t97; + coeff_dx_0_0_m1 = t86; + coeff_dx_p1_0_m1 = t82+t99; + coeff_dx_m1_p1_m1 = t84+t96; + coeff_dx_0_p1_m1 = t86; + coeff_dx_p1_p1_m1 = t82+t93; + coeff_dx_m1_m1_0 = t90+t98; + coeff_dx_0_m1_0 = t86; + coeff_dx_p1_m1_0 = t94; + coeff_dx_m1_0_0 = t100; + coeff_dx_0_0_0 = t86; + coeff_dx_p1_0_0 = t99; + coeff_dx_m1_p1_0 = t96; + coeff_dx_0_p1_0 = t86; + coeff_dx_p1_p1_0 = t93; + coeff_dx_m1_m1_p1 = t81+t95; + coeff_dx_0_m1_p1 = t86; + coeff_dx_p1_m1_p1 = t83+t92; + coeff_dx_m1_0_p1 = t95; + coeff_dx_0_0_p1 = t86; + coeff_dx_p1_0_p1 = t92; + coeff_dx_m1_p1_p1 = t83+t95; + coeff_dx_0_p1_p1 = t86; + coeff_dx_p1_p1_p1 = t88+t91; |