diff options
| author | jthorn <jthorn@df1f8a13-aa1d-4dd4-9681-27ded5b42416> | 2002-09-02 10:36:32 +0000 |
|---|---|---|
| committer | jthorn <jthorn@df1f8a13-aa1d-4dd4-9681-27ded5b42416> | 2002-09-02 10:36:32 +0000 |
| commit | 80df27de110b6e739fda09ba4d9eae9fdbefc342 (patch) | |
| tree | 3c17be6653cf3bd67a40fdf79f0821ff1afdd333 /src/GeneralizedPolynomial-Uniform/Lagrange/3d.coeffs/3d.cube.order2.smooth0/coeffs-dx.compute.c | |
| parent | 2248cabada859337c9945616de9a3d612d8f2199 (diff) | |
This is a major reworking of this interpolator.
The interpolation code in template.c is now split up into subfunctions
(in common/, Lagrange/, and Hermite/), so it should now compile in
finite cpu/memory even with optimization.
git-svn-id: http://svn.cactuscode.org/arrangements/CactusBase/LocalInterp/trunk@115 df1f8a13-aa1d-4dd4-9681-27ded5b42416
Diffstat (limited to 'src/GeneralizedPolynomial-Uniform/Lagrange/3d.coeffs/3d.cube.order2.smooth0/coeffs-dx.compute.c')
| -rw-r--r-- | src/GeneralizedPolynomial-Uniform/Lagrange/3d.coeffs/3d.cube.order2.smooth0/coeffs-dx.compute.c | 65 |
1 files changed, 65 insertions, 0 deletions
diff --git a/src/GeneralizedPolynomial-Uniform/Lagrange/3d.coeffs/3d.cube.order2.smooth0/coeffs-dx.compute.c b/src/GeneralizedPolynomial-Uniform/Lagrange/3d.coeffs/3d.cube.order2.smooth0/coeffs-dx.compute.c new file mode 100644 index 0000000..9f29978 --- /dev/null +++ b/src/GeneralizedPolynomial-Uniform/Lagrange/3d.coeffs/3d.cube.order2.smooth0/coeffs-dx.compute.c @@ -0,0 +1,65 @@ +fp t91; +fp t95; +fp t104; +fp t94; +fp t86; +fp t103; +fp t93; +fp t88; +fp t102; +fp t89; +fp t101; +fp t92; +fp t100; +fp t99; +fp t87; +fp t98; +fp t97; +fp t96; +fp t90; + t91 = RATIONAL(1.0,9.0)*x; + t95 = RATIONAL(1.0,18.0); + t104 = t91+t95; + t94 = RATIONAL(1.0,12.0); + t86 = t94*y; + t103 = t91+t86; + t93 = RATIONAL(-1.0,12.0); + t88 = t93*y; + t102 = t91+t88; + t89 = t94*z; + t101 = t91+t89; + t92 = RATIONAL(-1.0,18.0); + t100 = t92+t91; + t99 = t86+t101; + t87 = t93*z; + t98 = t87+t103; + t97 = t87+t102; + t96 = t88+t101; + t90 = RATIONAL(-2.0,9.0)*x; + coeffs_dx->coeff_m1_m1_m1 = t92+t99; + coeffs_dx->coeff_0_m1_m1 = t90; + coeffs_dx->coeff_p1_m1_m1 = t95+t97; + coeffs_dx->coeff_m1_0_m1 = t89+t100; + coeffs_dx->coeff_0_0_m1 = t90; + coeffs_dx->coeff_p1_0_m1 = t87+t104; + coeffs_dx->coeff_m1_p1_m1 = t92+t96; + coeffs_dx->coeff_0_p1_m1 = t90; + coeffs_dx->coeff_p1_p1_m1 = t95+t98; + coeffs_dx->coeff_m1_m1_0 = t86+t100; + coeffs_dx->coeff_0_m1_0 = t90; + coeffs_dx->coeff_p1_m1_0 = t95+t102; + coeffs_dx->coeff_m1_0_0 = t100; + coeffs_dx->coeff_0_0_0 = t90; + coeffs_dx->coeff_p1_0_0 = t104; + coeffs_dx->coeff_m1_p1_0 = t88+t100; + coeffs_dx->coeff_0_p1_0 = t90; + coeffs_dx->coeff_p1_p1_0 = t95+t103; + coeffs_dx->coeff_m1_m1_p1 = t92+t98; + coeffs_dx->coeff_0_m1_p1 = t90; + coeffs_dx->coeff_p1_m1_p1 = t95+t96; + coeffs_dx->coeff_m1_0_p1 = t87+t100; + coeffs_dx->coeff_0_0_p1 = t90; + coeffs_dx->coeff_p1_0_p1 = t95+t101; + coeffs_dx->coeff_m1_p1_p1 = t92+t97; + coeffs_dx->coeff_0_p1_p1 = t90; + coeffs_dx->coeff_p1_p1_p1 = t95+t99; |