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

1 files changed, 3 insertions, 3 deletions

diff --git a/src/Lagrange-tensor-product/2d.coeffs/2d.cube.order2.smooth0/coeffs-dxx.compute.c b/src/Lagrange-tensor-product/2d.coeffs/2d.cube.order2.smooth0/coeffs-dxx.compute.c
index d67f268..0599d25 100644
--- a/src/Lagrange-tensor-product/2d.coeffs/2d.cube.order2.smooth0/coeffs-dxx.compute.c
+++ b/src/Lagrange-tensor-product/2d.coeffs/2d.cube.order2.smooth0/coeffs-dxx.compute.c
@@ -7,14 +7,14 @@ fp t57;
       t62 = y*y;
       t61 = RATIONAL(1.0,2.0);
       t60 = t61*t62;
-      t59 = -t62+RATIONAL(1.0,1.0);
+      t59 = RATIONAL(1.0,1.0)-t62;
       t58 = RATIONAL(-1.0,2.0)*y+t60;
-      t57 = t60+t61*y;
+      t57 = t61*y+t60;
       coeffs_dxx->coeff_m1_m1 = t58;
       coeffs_dxx->coeff_0_m1 = y-t62;
       coeffs_dxx->coeff_p1_m1 = t58;
       coeffs_dxx->coeff_m1_0 = t59;
-      coeffs_dxx->coeff_0_0 = RATIONAL(2.0,1.0)*t62+RATIONAL(-2.0,1.0);
+      coeffs_dxx->coeff_0_0 = RATIONAL(-2.0,1.0)+RATIONAL(2.0,1.0)*t62;
       coeffs_dxx->coeff_p1_0 = t59;
       coeffs_dxx->coeff_m1_p1 = t57;
       coeffs_dxx->coeff_0_p1 = -t62-y;