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, 4 insertions, 4 deletions

diff --git a/src/Lagrange-tensor-product/2d.coeffs/2d.cube.order2.smooth0/coeffs-dyy.compute.c b/src/Lagrange-tensor-product/2d.coeffs/2d.cube.order2.smooth0/coeffs-dyy.compute.c
index 9a941cd..dfbebbe 100644
--- a/src/Lagrange-tensor-product/2d.coeffs/2d.cube.order2.smooth0/coeffs-dyy.compute.c
+++ b/src/Lagrange-tensor-product/2d.coeffs/2d.cube.order2.smooth0/coeffs-dyy.compute.c
@@ -7,15 +7,15 @@ fp t75;
       t80 = x*x;
       t79 = RATIONAL(1.0,2.0);
       t78 = t79*t80;
-      t77 = RATIONAL(1.0,1.0)-t80;
+      t77 = -t80+RATIONAL(1.0,1.0);
       t76 = t78+RATIONAL(-1.0,2.0)*x;
-      t75 = t79*x+t78;
+      t75 = t78+t79*x;
       coeffs_dyy->coeff_m1_m1 = t76;
       coeffs_dyy->coeff_0_m1 = t77;
       coeffs_dyy->coeff_p1_m1 = t75;
       coeffs_dyy->coeff_m1_0 = -t80+x;
-      coeffs_dyy->coeff_0_0 = RATIONAL(2.0,1.0)*t80+RATIONAL(-2.0,1.0);
-      coeffs_dyy->coeff_p1_0 = -x-t80;
+      coeffs_dyy->coeff_0_0 = RATIONAL(-2.0,1.0)+RATIONAL(2.0,1.0)*t80;
+      coeffs_dyy->coeff_p1_0 = -t80-x;
       coeffs_dyy->coeff_m1_p1 = t76;
       coeffs_dyy->coeff_0_p1 = t77;
       coeffs_dyy->coeff_p1_p1 = t75;