diff options
author | jthorn <jthorn@0f49ee68-0e4f-0410-9b9c-b2c123ded7ef> | 2004-05-15 14:11:49 +0000 |
---|---|---|
committer | jthorn <jthorn@0f49ee68-0e4f-0410-9b9c-b2c123ded7ef> | 2004-05-15 14:11:49 +0000 |
commit | 527f138c08f3ca615b124860aab14e30353212f9 (patch) | |
tree | 39089cfcd43afecacbe8334c3df6a4c115359f0c /src/Lagrange-tensor-product/2d.coeffs/2d.cube.order2.smooth0/coeffs-dxy.compute.c | |
parent | 214f8eae27c0c862f2b2ca17b683540a0c3d3662 (diff) |
fix a nasty bug which caused the 2D order=3 Lagrange interpolators
to be only first accurate :( :(
(it was a typo in the maple script generating the coefficients)
(I'm still investigating why my interpolator tests didn't spot this
before)
git-svn-id: http://svn.aei.mpg.de/numrel/AEIThorns/AEILocalInterp/trunk@18 0f49ee68-0e4f-0410-9b9c-b2c123ded7ef
Diffstat (limited to 'src/Lagrange-tensor-product/2d.coeffs/2d.cube.order2.smooth0/coeffs-dxy.compute.c')
-rw-r--r-- | src/Lagrange-tensor-product/2d.coeffs/2d.cube.order2.smooth0/coeffs-dxy.compute.c | 26 |
1 files changed, 13 insertions, 13 deletions
diff --git a/src/Lagrange-tensor-product/2d.coeffs/2d.cube.order2.smooth0/coeffs-dxy.compute.c b/src/Lagrange-tensor-product/2d.coeffs/2d.cube.order2.smooth0/coeffs-dxy.compute.c index 0fda87b..aea818d 100644 --- a/src/Lagrange-tensor-product/2d.coeffs/2d.cube.order2.smooth0/coeffs-dxy.compute.c +++ b/src/Lagrange-tensor-product/2d.coeffs/2d.cube.order2.smooth0/coeffs-dxy.compute.c @@ -5,25 +5,25 @@ fp t69; fp t73; fp t71; fp t70; -fp t65; -fp t64; +fp t67; +fp t66; fp t63; t72 = x*y; t68 = RATIONAL(-1.0,2.0); - t74 = t68*x+t72; + t74 = t68*y+t72; t69 = RATIONAL(1.0,2.0); - t73 = t69*x+t72; + t73 = t69*y+t72; t71 = RATIONAL(-1.0,4.0); t70 = RATIONAL(1.0,4.0); - t65 = t68*y; - t64 = t69*y; + t67 = t68*x; + t66 = t69*x; t63 = RATIONAL(-2.0,1.0)*t72; - coeffs_dxy->coeff_m1_m1 = t65+t70+t74; - coeffs_dxy->coeff_0_m1 = t63+x; - coeffs_dxy->coeff_p1_m1 = t71+t64+t74; - coeffs_dxy->coeff_m1_0 = t63+y; + coeffs_dxy->coeff_m1_m1 = t67+t70+t74; + coeffs_dxy->coeff_0_m1 = x+t63; + coeffs_dxy->coeff_p1_m1 = t71+t67+t73; + coeffs_dxy->coeff_m1_0 = y+t63; coeffs_dxy->coeff_0_0 = RATIONAL(4.0,1.0)*t72; - coeffs_dxy->coeff_p1_0 = -y+t63; - coeffs_dxy->coeff_m1_p1 = t65+t71+t73; + coeffs_dxy->coeff_p1_0 = t63-y; + coeffs_dxy->coeff_m1_p1 = t66+t71+t74; coeffs_dxy->coeff_0_p1 = t63-x; - coeffs_dxy->coeff_p1_p1 = t70+t64+t73; + coeffs_dxy->coeff_p1_p1 = t70+t66+t73; |