generate coeff-*.store.c files with non-I coefficients multiplied by "factor"

--> needed to handle 1/dx factors for derivatives


git-svn-id: http://svn.cactuscode.org/arrangements/CactusBase/LocalInterp/trunk@43 df1f8a13-aa1d-4dd4-9681-27ded5b42416

19 files changed, 170 insertions, 170 deletions

diff --git a/src/GeneralizedPolynomial-Uniform/1d.coeffs/1d.cube.order4.smooth0/coeff-I.compute.c b/src/GeneralizedPolynomial-Uniform/1d.coeffs/1d.cube.order4.smooth0/coeff-I.compute.c
index 318dbb1..861d59f 100644
--- a/src/GeneralizedPolynomial-Uniform/1d.coeffs/1d.cube.order4.smooth0/coeff-I.compute.c
+++ b/src/GeneralizedPolynomial-Uniform/1d.coeffs/1d.cube.order4.smooth0/coeff-I.compute.c
@@ -1,24 +1,24 @@
 fp t23,
    t24,
-   t29,
    t21,
    t22,
+   t29,
    t28,
    t25,
    t20,
    t19;
       t23 = x*x;
       t24 = t23*t23;
-      t29 = RATIONAL(1.0,24.0)*t24+RATIONAL(-1.0,24.0)*t23;
-      t21 = RATIONAL(-1.0,6.0);
-      t22 = RATIONAL(2.0,3.0);
-      t28 = t22*t23+t21*t24;
+      t21 = RATIONAL(2.0,3.0);
+      t22 = RATIONAL(-1.0,6.0);
+      t29 = t22*t24+t21*t23;
+      t28 = RATIONAL(-1.0,24.0)*t23+RATIONAL(1.0,24.0)*t24;
       t25 = x*t23;
-      t20 = RATIONAL(-1.0,12.0);
-      t19 = RATIONAL(1.0,12.0);
-      coeff_I_m2 = t19*x+t20*t25+t29;
-      coeff_I_m1 = RATIONAL(-2.0,3.0)*x+RATIONAL(1.0,6.0)*t25+t28;
+      t20 = RATIONAL(1.0,12.0);
+      t19 = RATIONAL(-1.0,12.0);
+      coeff_I_m2 = t20*x+t19*t25+t28;
+      coeff_I_m1 = RATIONAL(-2.0,3.0)*x+RATIONAL(1.0,6.0)*t25+t29;
       coeff_I_0 = RATIONAL(-5.0,4.0)*t23+RATIONAL(1.0,1.0)+RATIONAL(1.0,4.0)*
 t24;
-      coeff_I_p1 = t21*t25+t22*x+t28;
-      coeff_I_p2 = t19*t25+t20*x+t29;
+      coeff_I_p1 = t22*t25+t21*x+t29;
+      coeff_I_p2 = t20*t25+t19*x+t28;
diff --git a/src/GeneralizedPolynomial-Uniform/1d.coeffs/1d.cube.order4.smooth0/coeff-dx.compute.c b/src/GeneralizedPolynomial-Uniform/1d.coeffs/1d.cube.order4.smooth0/coeff-dx.compute.c
index 283d97d..5557616 100644
--- a/src/GeneralizedPolynomial-Uniform/1d.coeffs/1d.cube.order4.smooth0/coeff-dx.compute.c
+++ b/src/GeneralizedPolynomial-Uniform/1d.coeffs/1d.cube.order4.smooth0/coeff-dx.compute.c
@@ -1,17 +1,17 @@
-fp t36,
+fp t35,
+   t36,
    t37,
-   t35,
    t40,
    t34,
    t39;
+      t35 = RATIONAL(-2.0,3.0);
       t36 = x*x;
       t37 = x*t36;
-      t35 = RATIONAL(-2.0,3.0);
-      t40 = RATIONAL(4.0,3.0)*x+t35*t37;
+      t40 = t35*t37+RATIONAL(4.0,3.0)*x;
       t34 = RATIONAL(-1.0,12.0);
-      t39 = RATIONAL(1.0,6.0)*t37+t34*x;
+      t39 = t34*x+RATIONAL(1.0,6.0)*t37;
       coeff_dx_m2 = RATIONAL(-1.0,4.0)*t36+RATIONAL(1.0,12.0)+t39;
       coeff_dx_m1 = t35+RATIONAL(1.0,2.0)*t36+t40;
-      coeff_dx_0 = t37+RATIONAL(-5.0,2.0)*x;
-      coeff_dx_p1 = RATIONAL(-1.0,2.0)*t36+RATIONAL(2.0,3.0)+t40;
-      coeff_dx_p2 = t34+RATIONAL(1.0,4.0)*t36+t39;
+      coeff_dx_0 = RATIONAL(-5.0,2.0)*x+t37;
+      coeff_dx_p1 = RATIONAL(2.0,3.0)+RATIONAL(-1.0,2.0)*t36+t40;
+      coeff_dx_p2 = RATIONAL(1.0,4.0)*t36+t34+t39;
diff --git a/src/GeneralizedPolynomial-Uniform/1d.coeffs/1d.cube.order5.smooth0/coeff-I.compute.c b/src/GeneralizedPolynomial-Uniform/1d.coeffs/1d.cube.order5.smooth0/coeff-I.compute.c
index 3b005c4..8385b6e 100644
--- a/src/GeneralizedPolynomial-Uniform/1d.coeffs/1d.cube.order5.smooth0/coeff-I.compute.c
+++ b/src/GeneralizedPolynomial-Uniform/1d.coeffs/1d.cube.order5.smooth0/coeff-I.compute.c
@@ -1,25 +1,25 @@
 fp t56,
    t61,
-   t64,
-   t55,
    t54,
+   t55,
+   t64,
    t63,
    t58,
    t57,
    t51;
       t56 = x*x;
       t61 = t56*t56;
-      t64 = RATIONAL(2.0,3.0)*t56+RATIONAL(-1.0,6.0)*t61;
-      t55 = RATIONAL(-1.0,24.0);
       t54 = RATIONAL(1.0,24.0);
-      t63 = t55*t56+t54*t61;
+      t55 = RATIONAL(-1.0,24.0);
+      t64 = t55*t56+t54*t61;
+      t63 = RATIONAL(2.0,3.0)*t56+RATIONAL(-1.0,6.0)*t61;
       t58 = t61*x;
-      t57 = x*t56;
+      t57 = t56*x;
       t51 = t55*t57;
-      coeff_I_m2 = RATIONAL(1.0,20.0)*x+t51+RATIONAL(-1.0,120.0)*t58+t63;
-      coeff_I_m1 = RATIONAL(-1.0,2.0)*x+t51+t54*t58+t64;
+      coeff_I_m2 = RATIONAL(1.0,20.0)*x+t51+RATIONAL(-1.0,120.0)*t58+t64;
+      coeff_I_m1 = RATIONAL(-1.0,2.0)*x+t51+t54*t58+t63;
       coeff_I_0 = RATIONAL(1.0,1.0)+RATIONAL(-1.0,3.0)*x+RATIONAL(-5.0,4.0)*t56
 +RATIONAL(5.0,12.0)*t57+RATIONAL(1.0,4.0)*t61+RATIONAL(-1.0,12.0)*t58;
-      coeff_I_p1 = x+RATIONAL(-7.0,12.0)*t57+RATIONAL(1.0,12.0)*t58+t64;
-      coeff_I_p2 = RATIONAL(-1.0,4.0)*x+RATIONAL(7.0,24.0)*t57+t55*t58+t63;
+      coeff_I_p1 = x+RATIONAL(-7.0,12.0)*t57+RATIONAL(1.0,12.0)*t58+t63;
+      coeff_I_p2 = RATIONAL(-1.0,4.0)*x+RATIONAL(7.0,24.0)*t57+t55*t58+t64;
       coeff_I_p3 = RATIONAL(1.0,30.0)*x+RATIONAL(1.0,120.0)*t58+t51;
diff --git a/src/GeneralizedPolynomial-Uniform/1d.coeffs/1d.cube.order5.smooth0/coeff-dx.compute.c b/src/GeneralizedPolynomial-Uniform/1d.coeffs/1d.cube.order5.smooth0/coeff-dx.compute.c
index 166a52f..a816996 100644
--- a/src/GeneralizedPolynomial-Uniform/1d.coeffs/1d.cube.order5.smooth0/coeff-dx.compute.c
+++ b/src/GeneralizedPolynomial-Uniform/1d.coeffs/1d.cube.order5.smooth0/coeff-dx.compute.c
@@ -6,16 +6,16 @@ fp t70,
    t68;
       t70 = x*x;
       t71 = x*t70;
-      t76 = RATIONAL(-1.0,12.0)*x+RATIONAL(1.0,6.0)*t71;
+      t76 = RATIONAL(1.0,6.0)*t71+RATIONAL(-1.0,12.0)*x;
       t75 = RATIONAL(4.0,3.0)*x+RATIONAL(-2.0,3.0)*t71;
       t72 = t70*t70;
       t68 = RATIONAL(-1.0,8.0)*t70;
-      coeff_dx_m2 = RATIONAL(1.0,20.0)+RATIONAL(-1.0,24.0)*t72+t68+t76;
+      coeff_dx_m2 = RATIONAL(-1.0,24.0)*t72+t68+RATIONAL(1.0,20.0)+t76;
       coeff_dx_m1 = t68+RATIONAL(-1.0,2.0)+RATIONAL(5.0,24.0)*t72+t75;
-      coeff_dx_0 = RATIONAL(-1.0,3.0)+RATIONAL(-5.0,2.0)*x+t71+RATIONAL(5.0,4.0
-)*t70+RATIONAL(-5.0,12.0)*t72;
-      coeff_dx_p1 = RATIONAL(-7.0,4.0)*t70+RATIONAL(1.0,1.0)+RATIONAL(5.0,12.0)
-*t72+t75;
-      coeff_dx_p2 = RATIONAL(-1.0,4.0)+RATIONAL(7.0,8.0)*t70+RATIONAL(-5.0,24.0
-)*t72+t76;
-      coeff_dx_p3 = RATIONAL(1.0,24.0)*t72+RATIONAL(1.0,30.0)+t68;
+      coeff_dx_0 = RATIONAL(-5.0,12.0)*t72+t71+RATIONAL(5.0,4.0)*t70+RATIONAL(
+-5.0,2.0)*x+RATIONAL(-1.0,3.0);
+      coeff_dx_p1 = RATIONAL(5.0,12.0)*t72+RATIONAL(1.0,1.0)+RATIONAL(-7.0,4.0)
+*t70+t75;
+      coeff_dx_p2 = RATIONAL(-5.0,24.0)*t72+RATIONAL(7.0,8.0)*t70+RATIONAL(-1.0
+,4.0)+t76;
+      coeff_dx_p3 = RATIONAL(1.0,30.0)+t68+RATIONAL(1.0,24.0)*t72;
diff --git a/src/GeneralizedPolynomial-Uniform/1d.coeffs/1d.cube.order5.smooth0/coeff-dxx.compute.c b/src/GeneralizedPolynomial-Uniform/1d.coeffs/1d.cube.order5.smooth0/coeff-dxx.compute.c
index dc55053..f9af108 100644
--- a/src/GeneralizedPolynomial-Uniform/1d.coeffs/1d.cube.order5.smooth0/coeff-dxx.compute.c
+++ b/src/GeneralizedPolynomial-Uniform/1d.coeffs/1d.cube.order5.smooth0/coeff-dxx.compute.c
@@ -4,14 +4,14 @@ fp t82,
    t83,
    t77;
       t82 = x*x;
-      t86 = RATIONAL(-2.0,1.0)*t82+RATIONAL(4.0,3.0);
-      t85 = RATIONAL(1.0,2.0)*t82+RATIONAL(-1.0,12.0);
-      t83 = t82*x;
+      t86 = RATIONAL(1.0,2.0)*t82+RATIONAL(-1.0,12.0);
+      t85 = RATIONAL(-2.0,1.0)*t82+RATIONAL(4.0,3.0);
+      t83 = x*t82;
       t77 = RATIONAL(-1.0,4.0)*x;
-      coeff_dxx_m2 = RATIONAL(-1.0,6.0)*t83+t77+t85;
-      coeff_dxx_m1 = t77+RATIONAL(5.0,6.0)*t83+t86;
-      coeff_dxx_0 = RATIONAL(-5.0,3.0)*t83+RATIONAL(5.0,2.0)*x+RATIONAL(3.0,1.0
-)*t82+RATIONAL(-5.0,2.0);
-      coeff_dxx_p1 = RATIONAL(5.0,3.0)*t83+RATIONAL(-7.0,2.0)*x+t86;
-      coeff_dxx_p2 = RATIONAL(7.0,4.0)*x+RATIONAL(-5.0,6.0)*t83+t85;
-      coeff_dxx_p3 = RATIONAL(1.0,6.0)*t83+t77;
+      coeff_dxx_m2 = t77+RATIONAL(-1.0,6.0)*t83+t86;
+      coeff_dxx_m1 = t77+RATIONAL(5.0,6.0)*t83+t85;
+      coeff_dxx_0 = RATIONAL(5.0,2.0)*x+RATIONAL(3.0,1.0)*t82+RATIONAL(-5.0,3.0
+)*t83+RATIONAL(-5.0,2.0);
+      coeff_dxx_p1 = RATIONAL(-7.0,2.0)*x+RATIONAL(5.0,3.0)*t83+t85;
+      coeff_dxx_p2 = RATIONAL(-5.0,6.0)*t83+RATIONAL(7.0,4.0)*x+t86;
+      coeff_dxx_p3 = t77+RATIONAL(1.0,6.0)*t83;
diff --git a/src/GeneralizedPolynomial-Uniform/1d.coeffs/1d.cube.order6.smooth0/coeff-I.compute.c b/src/GeneralizedPolynomial-Uniform/1d.coeffs/1d.cube.order6.smooth0/coeff-I.compute.c
index d6cc566..1f7b51b 100644
--- a/src/GeneralizedPolynomial-Uniform/1d.coeffs/1d.cube.order6.smooth0/coeff-I.compute.c
+++ b/src/GeneralizedPolynomial-Uniform/1d.coeffs/1d.cube.order6.smooth0/coeff-I.compute.c
@@ -1,38 +1,38 @@
-fp t104,
-   t108,
+fp t96,
+   t100,
+   t101,
+   t104,
    t107,
-   t103,
+   t102,
+   t106,
    t113,
-   t101,
-   t96,
-   t97,
    t112,
    t111,
-   t102,
-   t100,
+   t103,
    t99,
-   t98;
-      t104 = x*x;
-      t108 = t104*t104;
-      t107 = x*t104;
-      t103 = t107*t107;
-      t113 = RATIONAL(-1.0,120.0)*t103+RATIONAL(1.0,12.0)*t108+RATIONAL(-3.0,
-40.0)*t104;
-      t101 = RATIONAL(3.0,4.0);
+   t98,
+   t97;
       t96 = RATIONAL(-13.0,48.0);
-      t97 = RATIONAL(1.0,48.0);
-      t112 = t96*t108+t101*t104+t97*t103;
-      t111 = RATIONAL(-1.0,144.0)*t108+RATIONAL(1.0,180.0)*t104+RATIONAL(1.0,
-720.0)*t103;
-      t102 = x*t108;
-      t100 = RATIONAL(1.0,60.0);
-      t99 = RATIONAL(-1.0,48.0);
-      t98 = RATIONAL(-1.0,60.0);
-      coeff_I_m3 = t98*x+t97*t107+RATIONAL(-1.0,240.0)*t102+t111;
-      coeff_I_m2 = RATIONAL(3.0,20.0)*x+RATIONAL(-1.0,6.0)*t107+t100*t102+t113;
-      coeff_I_m1 = RATIONAL(-3.0,4.0)*x+RATIONAL(13.0,48.0)*t107+t99*t102+t112;
-      coeff_I_0 = RATIONAL(7.0,18.0)*t108+RATIONAL(-49.0,36.0)*t104+RATIONAL(
-1.0,1.0)+RATIONAL(-1.0,36.0)*t103;
-      coeff_I_p1 = t96*t107+t101*x+t97*t102+t112;
-      coeff_I_p2 = RATIONAL(1.0,6.0)*t107+RATIONAL(-3.0,20.0)*x+t98*t102+t113;
-      coeff_I_p3 = t99*t107+t100*x+RATIONAL(1.0,240.0)*t102+t111;
+      t100 = RATIONAL(1.0,48.0);
+      t101 = RATIONAL(3.0,4.0);
+      t104 = x*x;
+      t107 = t104*x;
+      t102 = t107*t107;
+      t106 = t104*t104;
+      t113 = t100*t102+t101*t104+t96*t106;
+      t112 = RATIONAL(1.0,720.0)*t102+RATIONAL(1.0,180.0)*t104+RATIONAL(-1.0,
+144.0)*t106;
+      t111 = RATIONAL(1.0,12.0)*t106+RATIONAL(-3.0,40.0)*t104+RATIONAL(-1.0,
+120.0)*t102;
+      t103 = t104*t107;
+      t99 = RATIONAL(-1.0,60.0);
+      t98 = RATIONAL(1.0,60.0);
+      t97 = RATIONAL(-1.0,48.0);
+      coeff_I_m3 = t99*x+t100*t107+RATIONAL(-1.0,240.0)*t103+t112;
+      coeff_I_m2 = RATIONAL(3.0,20.0)*x+RATIONAL(-1.0,6.0)*t107+t98*t103+t111;
+      coeff_I_m1 = RATIONAL(-3.0,4.0)*x+RATIONAL(13.0,48.0)*t107+t97*t103+t113;
+      coeff_I_0 = RATIONAL(7.0,18.0)*t106+RATIONAL(-49.0,36.0)*t104+RATIONAL(
+1.0,1.0)+RATIONAL(-1.0,36.0)*t102;
+      coeff_I_p1 = t96*t107+t101*x+t100*t103+t113;
+      coeff_I_p2 = RATIONAL(1.0,6.0)*t107+RATIONAL(-3.0,20.0)*x+t99*t103+t111;
+      coeff_I_p3 = t97*t107+t98*x+RATIONAL(1.0,240.0)*t103+t112;
diff --git a/src/GeneralizedPolynomial-Uniform/1d.coeffs/1d.cube.order6.smooth0/coeff-dx.compute.c b/src/GeneralizedPolynomial-Uniform/1d.coeffs/1d.cube.order6.smooth0/coeff-dx.compute.c
index 252c658..738bb32 100644
--- a/src/GeneralizedPolynomial-Uniform/1d.coeffs/1d.cube.order6.smooth0/coeff-dx.compute.c
+++ b/src/GeneralizedPolynomial-Uniform/1d.coeffs/1d.cube.order6.smooth0/coeff-dx.compute.c
@@ -1,31 +1,31 @@
 fp t125,
-   t126,
-   t129,
+   t128,
    t124,
-   t123,
    t133,
+   t123,
    t132,
-   t131;
+   t131,
+   t127;
       t125 = x*x;
-      t126 = x*t125;
-      t129 = t125*t125;
-      t124 = x*t129;
+      t128 = x*t125;
+      t124 = t128*t125;
+      t133 = RATIONAL(-1.0,36.0)*t128+RATIONAL(1.0,90.0)*x+RATIONAL(1.0,120.0)*
+t124;
       t123 = RATIONAL(-3.0,20.0);
-      t133 = RATIONAL(-1.0,20.0)*t124+t123*x+RATIONAL(1.0,3.0)*t126;
-      t132 = RATIONAL(-13.0,12.0)*t126+RATIONAL(1.0,8.0)*t124+RATIONAL(3.0,2.0)
-*x;
-      t131 = RATIONAL(1.0,90.0)*x+RATIONAL(-1.0,36.0)*t126+RATIONAL(1.0,120.0)*
+      t132 = RATIONAL(-1.0,20.0)*t124+t123*x+RATIONAL(1.0,3.0)*t128;
+      t131 = RATIONAL(-13.0,12.0)*t128+RATIONAL(3.0,2.0)*x+RATIONAL(1.0,8.0)*
 t124;
-      coeff_dx_m3 = RATIONAL(1.0,16.0)*t125+RATIONAL(-1.0,48.0)*t129+RATIONAL(
--1.0,60.0)+t131;
-      coeff_dx_m2 = RATIONAL(1.0,12.0)*t129+RATIONAL(3.0,20.0)+RATIONAL(-1.0,
-2.0)*t125+t133;
-      coeff_dx_m1 = RATIONAL(-3.0,4.0)+RATIONAL(-5.0,48.0)*t129+RATIONAL(13.0,
-16.0)*t125+t132;
-      coeff_dx_0 = RATIONAL(-1.0,6.0)*t124+RATIONAL(-49.0,18.0)*x+RATIONAL(14.0
-,9.0)*t126;
-      coeff_dx_p1 = RATIONAL(-13.0,16.0)*t125+RATIONAL(5.0,48.0)*t129+RATIONAL(
-3.0,4.0)+t132;
-      coeff_dx_p2 = t123+RATIONAL(-1.0,12.0)*t129+RATIONAL(1.0,2.0)*t125+t133;
-      coeff_dx_p3 = RATIONAL(1.0,60.0)+RATIONAL(-1.0,16.0)*t125+RATIONAL(1.0,
-48.0)*t129+t131;
+      t127 = t125*t125;
+      coeff_dx_m3 = RATIONAL(1.0,16.0)*t125+RATIONAL(-1.0,48.0)*t127+RATIONAL(
+-1.0,60.0)+t133;
+      coeff_dx_m2 = RATIONAL(-1.0,2.0)*t125+RATIONAL(3.0,20.0)+RATIONAL(1.0,
+12.0)*t127+t132;
+      coeff_dx_m1 = RATIONAL(-3.0,4.0)+RATIONAL(13.0,16.0)*t125+RATIONAL(-5.0,
+48.0)*t127+t131;
+      coeff_dx_0 = RATIONAL(-1.0,6.0)*t124+RATIONAL(14.0,9.0)*t128+RATIONAL(
+-49.0,18.0)*x;
+      coeff_dx_p1 = RATIONAL(-13.0,16.0)*t125+RATIONAL(3.0,4.0)+RATIONAL(5.0,
+48.0)*t127+t131;
+      coeff_dx_p2 = t123+RATIONAL(-1.0,12.0)*t127+RATIONAL(1.0,2.0)*t125+t132;
+      coeff_dx_p3 = RATIONAL(1.0,48.0)*t127+RATIONAL(-1.0,16.0)*t125+RATIONAL(
+1.0,60.0)+t133;
diff --git a/src/GeneralizedPolynomial-Uniform/1d.coeffs/1d.cube.order6.smooth0/coeff-dxx.compute.c b/src/GeneralizedPolynomial-Uniform/1d.coeffs/1d.cube.order6.smooth0/coeff-dxx.compute.c
index 55b1f9d..e63d0f2 100644
--- a/src/GeneralizedPolynomial-Uniform/1d.coeffs/1d.cube.order6.smooth0/coeff-dxx.compute.c
+++ b/src/GeneralizedPolynomial-Uniform/1d.coeffs/1d.cube.order6.smooth0/coeff-dxx.compute.c
@@ -1,22 +1,22 @@
 fp t143,
    t145,
    t150,
-   t141,
+   t140,
    t149,
    t148,
    t144;
       t143 = x*x;
       t145 = t143*t143;
-      t150 = RATIONAL(-1.0,4.0)*t145+t143+RATIONAL(-3.0,20.0);
-      t141 = RATIONAL(-1.0,12.0);
-      t149 = t141*t143+RATIONAL(1.0,24.0)*t145+RATIONAL(1.0,90.0);
-      t148 = RATIONAL(3.0,2.0)+RATIONAL(-13.0,4.0)*t143+RATIONAL(5.0,8.0)*t145;
-      t144 = x*t143;
-      coeff_dxx_m3 = t141*t144+RATIONAL(1.0,8.0)*x+t149;
-      coeff_dxx_m2 = RATIONAL(1.0,3.0)*t144-x+t150;
-      coeff_dxx_m1 = RATIONAL(-5.0,12.0)*t144+RATIONAL(13.0,8.0)*x+t148;
-      coeff_dxx_0 = RATIONAL(-5.0,6.0)*t145+RATIONAL(14.0,3.0)*t143+RATIONAL(
--49.0,18.0);
-      coeff_dxx_p1 = RATIONAL(5.0,12.0)*t144+RATIONAL(-13.0,8.0)*x+t148;
-      coeff_dxx_p2 = RATIONAL(-1.0,3.0)*t144+x+t150;
-      coeff_dxx_p3 = RATIONAL(-1.0,8.0)*x+RATIONAL(1.0,12.0)*t144+t149;
+      t150 = RATIONAL(-13.0,4.0)*t143+RATIONAL(5.0,8.0)*t145+RATIONAL(3.0,2.0);
+      t140 = RATIONAL(-1.0,12.0);
+      t149 = t140*t143+RATIONAL(1.0,24.0)*t145+RATIONAL(1.0,90.0);
+      t148 = RATIONAL(-3.0,20.0)+RATIONAL(-1.0,4.0)*t145+t143;
+      t144 = t143*x;
+      coeff_dxx_m3 = t140*t144+RATIONAL(1.0,8.0)*x+t149;
+      coeff_dxx_m2 = RATIONAL(1.0,3.0)*t144-x+t148;
+      coeff_dxx_m1 = RATIONAL(-5.0,12.0)*t144+RATIONAL(13.0,8.0)*x+t150;
+      coeff_dxx_0 = RATIONAL(14.0,3.0)*t143+RATIONAL(-49.0,18.0)+RATIONAL(-5.0,
+6.0)*t145;
+      coeff_dxx_p1 = RATIONAL(5.0,12.0)*t144+RATIONAL(-13.0,8.0)*x+t150;
+      coeff_dxx_p2 = x+RATIONAL(-1.0,3.0)*t144+t148;
+      coeff_dxx_p3 = RATIONAL(1.0,12.0)*t144+RATIONAL(-1.0,8.0)*x+t149;
diff --git a/src/GeneralizedPolynomial-Uniform/1d.coeffs/1d.cube.size2/coeff-dx.store.c b/src/GeneralizedPolynomial-Uniform/1d.coeffs/1d.cube.size2/coeff-dx.store.c
index 18f835c..eb49618 100644
--- a/src/GeneralizedPolynomial-Uniform/1d.coeffs/1d.cube.size2/coeff-dx.store.c
+++ b/src/GeneralizedPolynomial-Uniform/1d.coeffs/1d.cube.size2/coeff-dx.store.c
@@ -1,2 +1,2 @@
-COEFF(0) = coeff_dx_0;
-COEFF(1) = coeff_dx_p1;
+COEFF(0) = factor * coeff_dx_0;
+COEFF(1) = factor * coeff_dx_p1;
diff --git a/src/GeneralizedPolynomial-Uniform/1d.coeffs/1d.cube.size3/coeff-dx.store.c b/src/GeneralizedPolynomial-Uniform/1d.coeffs/1d.cube.size3/coeff-dx.store.c
index 1454448..672d338 100644
--- a/src/GeneralizedPolynomial-Uniform/1d.coeffs/1d.cube.size3/coeff-dx.store.c
+++ b/src/GeneralizedPolynomial-Uniform/1d.coeffs/1d.cube.size3/coeff-dx.store.c
@@ -1,3 +1,3 @@
-COEFF(-1) = coeff_dx_m1;
-COEFF(0) = coeff_dx_0;
-COEFF(1) = coeff_dx_p1;
+COEFF(-1) = factor * coeff_dx_m1;
+COEFF(0) = factor * coeff_dx_0;
+COEFF(1) = factor * coeff_dx_p1;
diff --git a/src/GeneralizedPolynomial-Uniform/1d.coeffs/1d.cube.size3/coeff-dxx.store.c b/src/GeneralizedPolynomial-Uniform/1d.coeffs/1d.cube.size3/coeff-dxx.store.c
index be09b3f..72b7d92 100644
--- a/src/GeneralizedPolynomial-Uniform/1d.coeffs/1d.cube.size3/coeff-dxx.store.c
+++ b/src/GeneralizedPolynomial-Uniform/1d.coeffs/1d.cube.size3/coeff-dxx.store.c
@@ -1,3 +1,3 @@
-COEFF(-1) = coeff_dxx_m1;
-COEFF(0) = coeff_dxx_0;
-COEFF(1) = coeff_dxx_p1;
+COEFF(-1) = factor * coeff_dxx_m1;
+COEFF(0) = factor * coeff_dxx_0;
+COEFF(1) = factor * coeff_dxx_p1;
diff --git a/src/GeneralizedPolynomial-Uniform/1d.coeffs/1d.cube.size4/coeff-dx.store.c b/src/GeneralizedPolynomial-Uniform/1d.coeffs/1d.cube.size4/coeff-dx.store.c
index 77e5711..ebf3b47 100644
--- a/src/GeneralizedPolynomial-Uniform/1d.coeffs/1d.cube.size4/coeff-dx.store.c
+++ b/src/GeneralizedPolynomial-Uniform/1d.coeffs/1d.cube.size4/coeff-dx.store.c
@@ -1,4 +1,4 @@
-COEFF(-1) = coeff_dx_m1;
-COEFF(0) = coeff_dx_0;
-COEFF(1) = coeff_dx_p1;
-COEFF(2) = coeff_dx_p2;
+COEFF(-1) = factor * coeff_dx_m1;
+COEFF(0) = factor * coeff_dx_0;
+COEFF(1) = factor * coeff_dx_p1;
+COEFF(2) = factor * coeff_dx_p2;
diff --git a/src/GeneralizedPolynomial-Uniform/1d.coeffs/1d.cube.size4/coeff-dxx.store.c b/src/GeneralizedPolynomial-Uniform/1d.coeffs/1d.cube.size4/coeff-dxx.store.c
index c0a6695..487736f 100644
--- a/src/GeneralizedPolynomial-Uniform/1d.coeffs/1d.cube.size4/coeff-dxx.store.c
+++ b/src/GeneralizedPolynomial-Uniform/1d.coeffs/1d.cube.size4/coeff-dxx.store.c
@@ -1,4 +1,4 @@
-COEFF(-1) = coeff_dxx_m1;
-COEFF(0) = coeff_dxx_0;
-COEFF(1) = coeff_dxx_p1;
-COEFF(2) = coeff_dxx_p2;
+COEFF(-1) = factor * coeff_dxx_m1;
+COEFF(0) = factor * coeff_dxx_0;
+COEFF(1) = factor * coeff_dxx_p1;
+COEFF(2) = factor * coeff_dxx_p2;
diff --git a/src/GeneralizedPolynomial-Uniform/1d.coeffs/1d.cube.size5/coeff-dx.store.c b/src/GeneralizedPolynomial-Uniform/1d.coeffs/1d.cube.size5/coeff-dx.store.c
index e6ae190..1dc93a4 100644
--- a/src/GeneralizedPolynomial-Uniform/1d.coeffs/1d.cube.size5/coeff-dx.store.c
+++ b/src/GeneralizedPolynomial-Uniform/1d.coeffs/1d.cube.size5/coeff-dx.store.c
@@ -1,5 +1,5 @@
-COEFF(-2) = coeff_dx_m2;
-COEFF(-1) = coeff_dx_m1;
-COEFF(0) = coeff_dx_0;
-COEFF(1) = coeff_dx_p1;
-COEFF(2) = coeff_dx_p2;
+COEFF(-2) = factor * coeff_dx_m2;
+COEFF(-1) = factor * coeff_dx_m1;
+COEFF(0) = factor * coeff_dx_0;
+COEFF(1) = factor * coeff_dx_p1;
+COEFF(2) = factor * coeff_dx_p2;
diff --git a/src/GeneralizedPolynomial-Uniform/1d.coeffs/1d.cube.size5/coeff-dxx.store.c b/src/GeneralizedPolynomial-Uniform/1d.coeffs/1d.cube.size5/coeff-dxx.store.c
index 3644a16..2d5fd79 100644
--- a/src/GeneralizedPolynomial-Uniform/1d.coeffs/1d.cube.size5/coeff-dxx.store.c
+++ b/src/GeneralizedPolynomial-Uniform/1d.coeffs/1d.cube.size5/coeff-dxx.store.c
@@ -1,5 +1,5 @@
-COEFF(-2) = coeff_dxx_m2;
-COEFF(-1) = coeff_dxx_m1;
-COEFF(0) = coeff_dxx_0;
-COEFF(1) = coeff_dxx_p1;
-COEFF(2) = coeff_dxx_p2;
+COEFF(-2) = factor * coeff_dxx_m2;
+COEFF(-1) = factor * coeff_dxx_m1;
+COEFF(0) = factor * coeff_dxx_0;
+COEFF(1) = factor * coeff_dxx_p1;
+COEFF(2) = factor * coeff_dxx_p2;
diff --git a/src/GeneralizedPolynomial-Uniform/1d.coeffs/1d.cube.size6/coeff-dx.store.c b/src/GeneralizedPolynomial-Uniform/1d.coeffs/1d.cube.size6/coeff-dx.store.c
index aabec4f..3a03b04 100644
--- a/src/GeneralizedPolynomial-Uniform/1d.coeffs/1d.cube.size6/coeff-dx.store.c
+++ b/src/GeneralizedPolynomial-Uniform/1d.coeffs/1d.cube.size6/coeff-dx.store.c
@@ -1,6 +1,6 @@
-COEFF(-2) = coeff_dx_m2;
-COEFF(-1) = coeff_dx_m1;
-COEFF(0) = coeff_dx_0;
-COEFF(1) = coeff_dx_p1;
-COEFF(2) = coeff_dx_p2;
-COEFF(3) = coeff_dx_p3;
+COEFF(-2) = factor * coeff_dx_m2;
+COEFF(-1) = factor * coeff_dx_m1;
+COEFF(0) = factor * coeff_dx_0;
+COEFF(1) = factor * coeff_dx_p1;
+COEFF(2) = factor * coeff_dx_p2;
+COEFF(3) = factor * coeff_dx_p3;
diff --git a/src/GeneralizedPolynomial-Uniform/1d.coeffs/1d.cube.size6/coeff-dxx.store.c b/src/GeneralizedPolynomial-Uniform/1d.coeffs/1d.cube.size6/coeff-dxx.store.c
index 5263353..d437161 100644
--- a/src/GeneralizedPolynomial-Uniform/1d.coeffs/1d.cube.size6/coeff-dxx.store.c
+++ b/src/GeneralizedPolynomial-Uniform/1d.coeffs/1d.cube.size6/coeff-dxx.store.c
@@ -1,6 +1,6 @@
-COEFF(-2) = coeff_dxx_m2;
-COEFF(-1) = coeff_dxx_m1;
-COEFF(0) = coeff_dxx_0;
-COEFF(1) = coeff_dxx_p1;
-COEFF(2) = coeff_dxx_p2;
-COEFF(3) = coeff_dxx_p3;
+COEFF(-2) = factor * coeff_dxx_m2;
+COEFF(-1) = factor * coeff_dxx_m1;
+COEFF(0) = factor * coeff_dxx_0;
+COEFF(1) = factor * coeff_dxx_p1;
+COEFF(2) = factor * coeff_dxx_p2;
+COEFF(3) = factor * coeff_dxx_p3;
diff --git a/src/GeneralizedPolynomial-Uniform/1d.coeffs/1d.cube.size7/coeff-dx.store.c b/src/GeneralizedPolynomial-Uniform/1d.coeffs/1d.cube.size7/coeff-dx.store.c
index ddbef48..7b76907 100644
--- a/src/GeneralizedPolynomial-Uniform/1d.coeffs/1d.cube.size7/coeff-dx.store.c
+++ b/src/GeneralizedPolynomial-Uniform/1d.coeffs/1d.cube.size7/coeff-dx.store.c
@@ -1,7 +1,7 @@
-COEFF(-3) = coeff_dx_m3;
-COEFF(-2) = coeff_dx_m2;
-COEFF(-1) = coeff_dx_m1;
-COEFF(0) = coeff_dx_0;
-COEFF(1) = coeff_dx_p1;
-COEFF(2) = coeff_dx_p2;
-COEFF(3) = coeff_dx_p3;
+COEFF(-3) = factor * coeff_dx_m3;
+COEFF(-2) = factor * coeff_dx_m2;
+COEFF(-1) = factor * coeff_dx_m1;
+COEFF(0) = factor * coeff_dx_0;
+COEFF(1) = factor * coeff_dx_p1;
+COEFF(2) = factor * coeff_dx_p2;
+COEFF(3) = factor * coeff_dx_p3;
diff --git a/src/GeneralizedPolynomial-Uniform/1d.coeffs/1d.cube.size7/coeff-dxx.store.c b/src/GeneralizedPolynomial-Uniform/1d.coeffs/1d.cube.size7/coeff-dxx.store.c
index 857bbfb..b65bc24 100644
--- a/src/GeneralizedPolynomial-Uniform/1d.coeffs/1d.cube.size7/coeff-dxx.store.c
+++ b/src/GeneralizedPolynomial-Uniform/1d.coeffs/1d.cube.size7/coeff-dxx.store.c
@@ -1,7 +1,7 @@
-COEFF(-3) = coeff_dxx_m3;
-COEFF(-2) = coeff_dxx_m2;
-COEFF(-1) = coeff_dxx_m1;
-COEFF(0) = coeff_dxx_0;
-COEFF(1) = coeff_dxx_p1;
-COEFF(2) = coeff_dxx_p2;
-COEFF(3) = coeff_dxx_p3;
+COEFF(-3) = factor * coeff_dxx_m3;
+COEFF(-2) = factor * coeff_dxx_m2;
+COEFF(-1) = factor * coeff_dxx_m1;
+COEFF(0) = factor * coeff_dxx_0;
+COEFF(1) = factor * coeff_dxx_p1;
+COEFF(2) = factor * coeff_dxx_p2;
+COEFF(3) = factor * coeff_dxx_p3;