1 files changed, 188 insertions, 206 deletions
diff --git a/src/Lagrange-tensor-product/3d.coeffs/3d.cube.order2.smooth0/coeffs-I.compute.c b/src/Lagrange-tensor-product/3d.coeffs/3d.cube.order2.smooth0/coeffs-I.compute.c
index d629ea0..fc1b8a5 100644
--- a/src/Lagrange-tensor-product/3d.coeffs/3d.cube.order2.smooth0/coeffs-I.compute.c
+++ b/src/Lagrange-tensor-product/3d.coeffs/3d.cube.order2.smooth0/coeffs-I.compute.c
@@ -1,245 +1,227 @@
-fp t82;
+fp t81;
 fp t84;
-fp t114;
-fp t83;
-fp t104;
-fp t26;
-fp t106;
-fp t34;
-fp t154;
+fp t119;
 fp t85;
-fp t129;
-fp t105;
-fp t25;
-fp t29;
-fp t153;
-fp t81;
-fp t130;
-fp t107;
-fp t36;
-fp t152;
-fp t108;
-fp t38;
-fp t46;
-fp t151;
-fp t109;
-fp t45;
-fp t150;
-fp t27;
-fp t131;
-fp t110;
-fp t35;
-fp t149;
-fp t42;
-fp t48;
-fp t148;
-fp t147;
-fp t77;
-fp t111;
-fp t51;
-fp t78;
-fp t128;
-fp t146;
-fp t123;
-fp t60;
-fp t127;
-fp t145;
 fp t115;
-fp t64;
-fp t31;
-fp t119;
-fp t144;
-fp t52;
-fp t143;
-fp t124;
+fp t78;
+fp t56;
+fp t29;
 fp t142;
-fp t70;
-fp t33;
+fp t118;
+fp t65;
 fp t141;
+fp t83;
+fp t116;
+fp t54;
+fp t30;
+fp t120;
 fp t63;
-fp t44;
 fp t140;
 fp t139;
-fp t103;
-fp t41;
-fp t68;
+fp t117;
+fp t50;
+fp t32;
 fp t138;
-fp t47;
+fp t69;
 fp t137;
+fp t103;
+fp t43;
 fp t136;
-fp t50;
 fp t135;
-fp t126;
-fp t65;
 fp t134;
+fp t122;
+fp t106;
+fp t44;
 fp t133;
-fp t24;
+fp t121;
+fp t55;
 fp t132;
+fp t105;
+fp t45;
+fp t123;
+fp t52;
+fp t131;
+fp t130;
+fp t104;
+fp t47;
+fp t129;
+fp t60;
+fp t128;
+fp t127;
+fp t126;
 fp t125;
-fp t121;
-fp t80;
+fp t124;
+fp t114;
 fp t79;
-fp t76;
-fp t32;
-fp t122;
-fp t120;
-fp t118;
-fp t117;
-fp t116;
+fp t35;
+fp t82;
 fp t113;
 fp t112;
+fp t111;
+fp t110;
+fp t109;
+fp t108;
+fp t107;
 fp t102;
 fp t101;
 fp t100;
+fp t77;
+fp t21;
+fp t22;
+fp t23;
+fp t24;
 fp t99;
+fp t25;
+fp t80;
+fp t37;
+fp t41;
 fp t98;
-fp t59;
-fp t62;
+fp t27;
+fp t36;
+fp t42;
 fp t97;
+fp t28;
+fp t38;
+fp t39;
 fp t96;
 fp t95;
 fp t94;
 fp t93;
+fp t72;
 fp t92;
 fp t91;
 fp t90;
 fp t89;
+fp t71;
 fp t88;
 fp t87;
 fp t86;
-      t82 = RATIONAL(1.0,8.0);
-      t84 = z*z;
-      t114 = t82*t84;
-      t83 = x*x;
-      t104 = t83*t114;
-      t26 = y*t104;
-      t106 = x*t114;
-      t34 = y*t106;
-      t154 = t26+t34;
-      t85 = y*y;
-      t129 = z*t85;
-      t105 = t82*t129;
-      t25 = t83*t105;
-      t29 = x*t105;
-      t153 = t25+t29;
-      t81 = RATIONAL(-1.0,8.0);
-      t130 = x*t81;
-      t107 = t84*t130;
-      t36 = y*t107;
-      t152 = t26+t36;
-      t108 = t81*t129;
-      t38 = x*t108;
-      t46 = t83*t108;
-      t151 = t38+t46;
-      t109 = y*t81*t83;
-      t45 = t84*t109;
-      t150 = t34+t45;
-      t27 = z*t109;
-      t131 = y*z;
-      t110 = t82*t131;
-      t35 = x*t110;
-      t149 = t27+t35;
-      t42 = t130*t131;
-      t48 = t83*t110;
-      t148 = t42+t48;
-      t147 = t36+t45;
-      t77 = RATIONAL(1.0,4.0);
-      t111 = t77*t83;
-      t51 = z*t111;
-      t78 = RATIONAL(-1.0,4.0);
-      t128 = z*t78;
-      t146 = (t128+t51)*y;
-      t123 = t78*t84;
-      t60 = x*t123;
-      t127 = y*t77;
-      t145 = x*t127+y*t60;
-      t115 = t78*t85;
-      t64 = x*t115;
-      t31 = t84*t64;
-      t119 = t77*t85;
-      t144 = x*t119+t31;
-      t52 = z*t119;
-      t143 = (t128+t52)*x;
-      t124 = t77*t84;
-      t142 = x*t124+t31;
-      t70 = t84*t119;
-      t33 = x*t70;
-      t141 = t60+t33;
-      t63 = t85*t111;
-      t44 = z*t63;
-      t140 = z*t115+t44;
-      t139 = (t64+t77*x)*z;
+fp t40;
+fp t26;
+      t81 = RATIONAL(-1.0,4.0);
+      t84 = y*y;
+      t119 = t81*t84;
+      t85 = x*x;
+      t115 = t85*t84;
+      t78 = RATIONAL(1.0,4.0);
+      t56 = t78*t115;
+      t29 = z*t56;
+      t142 = t29+z*t119;
+      t118 = t81*t85;
+      t65 = z*t118;
+      t141 = t29+t65;
+      t83 = z*z;
+      t116 = t85*t83;
+      t54 = t78*t116;
+      t30 = y*t54;
+      t120 = t81*t83;
+      t63 = y*t120;
+      t140 = t30+t63;
+      t139 = t30+y*t118;
+      t117 = t84*t83;
+      t50 = t78*t117;
+      t32 = x*t50;
+      t138 = t32+x*t120;
+      t69 = x*t119;
+      t137 = t32+t69;
       t103 = t83*t115;
-      t41 = t84*t103;
-      t68 = t84*t111;
-      t138 = t41+t68;
-      t47 = t85*t104;
-      t137 = t85*t106+t47;
-      t136 = t85*t107+t47;
-      t50 = y*t124;
-      t135 = (t78*y+t50)*x;
-      t126 = t78*t83;
-      t65 = y*t126;
-      t134 = (t127+t65)*z;
-      t133 = t64+t33;
-      t24 = z*t103;
-      t132 = t24+t52;
-      t125 = t85+t84;
-      t121 = t83*t84;
-      t80 = RATIONAL(-1.0,2.0);
+      t43 = t81*t103;
+      t136 = t43+t50;
+      t135 = t43+t54;
+      t134 = t43+t56;
+      t122 = t78*t84;
+      t106 = z*t115;
+      t44 = t81*t106;
+      t133 = t44+z*t122;
+      t121 = t78*t85;
+      t55 = z*t121;
+      t132 = t44+t55;
+      t105 = y*t116;
+      t45 = t81*t105;
+      t123 = t78*t83;
+      t52 = y*t123;
+      t131 = t45+t52;
+      t130 = t45+y*t121;
+      t104 = x*t117;
+      t47 = t81*t104;
+      t129 = t47+x*t123;
+      t60 = x*t122;
+      t128 = t47+t60;
+      t127 = x*y;
+      t126 = y*z;
+      t125 = t78*z;
+      t124 = t81*z;
+      t114 = x*t126;
       t79 = RATIONAL(1.0,2.0);
-      t76 = t79*t85;
-      t32 = t76*t121;
-      t122 = t80*t121+t32;
-      t120 = t85*t84;
-      t118 = y*t68+t41;
-      t117 = t84*t65+t41;
-      t116 = t84+t83;
-      t113 = t83+t85;
-      t112 = t80*t85;
-      t102 = 1.0+t120;
-      t101 = t121+1.0;
-      t100 = t83*t85+1.0;
-      t99 = t29+t46+t136;
-      t98 = z*t126+t44+t138;
-      t59 = t84*t112;
-      t62 = t83*t112;
-      t97 = t59+t32+t76+t62;
-      t96 = t79*t83+t62+t122;
-      t95 = t25+t38+t136;
-      t94 = y*t123+t70+t118;
-      t93 = t50+t70+t117;
-      t92 = t59+t79*t84+t122;
-      t91 = t63+t65+t118;
-      t90 = t63+y*t111+t117;
-      t89 = t35+t48+t137;
-      t88 = t27+t42+t137;
-      t87 = t51+t24+t138;
-      t86 = -t120+t125;
-      coeffs_I->coeff_m1_m1_m1 = t99+t148+t150;
-      coeffs_I->coeff_0_m1_m1 = t94+t134+t140;
-      coeffs_I->coeff_p1_m1_m1 = t89+t147+t151;
-      coeffs_I->coeff_m1_0_m1 = t98+t139+t141;
-      coeffs_I->coeff_0_0_m1 = (t100*t80+t113*t79)*z+t92;
-      coeffs_I->coeff_p1_0_m1 = t98+t142+t143;
-      coeffs_I->coeff_m1_p1_m1 = t99+t149+t152;
-      coeffs_I->coeff_0_p1_m1 = t93+t140+t146;
-      coeffs_I->coeff_p1_p1_m1 = t88+t151+t154;
-      coeffs_I->coeff_m1_m1_0 = t91+t133+t145;
-      coeffs_I->coeff_0_m1_0 = (t101*t80+t116*t79)*y+t97;
-      coeffs_I->coeff_p1_m1_0 = t91+t135+t144;
-      coeffs_I->coeff_m1_0_0 = (t102*t80+t125*t79)*x+t96;
+      t35 = t79*t103;
+      t82 = RATIONAL(-1.0,2.0);
+      t113 = t35+t82*t115;
+      t112 = t83+t84;
+      t111 = t83+t85;
+      t110 = t84+t85;
+      t109 = t83*t127;
+      t108 = t84*x*z;
+      t107 = t85*t126;
+      t102 = 1.0+t116;
+      t101 = 1.0+t115;
+      t100 = t117+1.0;
+      t77 = RATIONAL(1.0,8.0);
+      t21 = t77*t103;
+      t22 = t77*t106;
+      t23 = t77*t105;
+      t24 = t77*t107;
+      t99 = t21+t22+t23+t24;
+      t25 = t77*t104;
+      t80 = RATIONAL(-1.0,8.0);
+      t37 = t80*t105;
+      t41 = t80*t109;
+      t98 = t21+t25+t37+t41;
+      t27 = t77*t109;
+      t36 = t80*t106;
+      t42 = t80*t114;
+      t97 = t21+t27+t36+t42;
+      t28 = t77*t114;
+      t38 = t80*t107;
+      t39 = t80*t104;
+      t96 = t21+t28+t38+t39;
+      t95 = t136+(t55+t124)*y;
+      t94 = z*t60+x*t124+t135;
+      t93 = x*t52+t81*t127+t134;
+      t72 = t82*t116;
+      t92 = t79*t85+t72+t113;
+      t91 = t136+(t65+t125)*y;
+      t90 = z*t69+x*t125+t135;
+      t89 = x*t63+t78*t127+t134;
+      t71 = t82*t117;
+      t88 = t72+t71+t35+t79*t83;
+      t87 = t79*t84+t71+t113;
+      t86 = -t115+t110;
+      t40 = t80*t108;
+      t26 = t77*t108;
+      coeffs_I->coeff_m1_m1_m1 = t26+t39+t24+t37+t97;
+      coeffs_I->coeff_0_m1_m1 = t91+t140+t142;
+      coeffs_I->coeff_p1_m1_m1 = t28+t24+t40+t36+t98;
+      coeffs_I->coeff_m1_0_m1 = t90+t138+t141;
+      coeffs_I->coeff_0_0_m1 = (t101*t82+t110*t79)*z+t88;
+      coeffs_I->coeff_p1_0_m1 = t94+t129+t141;
+      coeffs_I->coeff_m1_p1_m1 = t23+t41+t26+t36+t96;
+      coeffs_I->coeff_0_p1_m1 = t95+t131+t142;
+      coeffs_I->coeff_p1_p1_m1 = t38+t40+t25+t23+t97;
+      coeffs_I->coeff_m1_m1_0 = t89+t137+t139;
+      coeffs_I->coeff_0_m1_0 = (t102*t82+t111*t79)*y+t87;
+      coeffs_I->coeff_p1_m1_0 = t93+t128+t139;
+      coeffs_I->coeff_m1_0_0 = (t100*t82+t112*t79)*x+t92;
       coeffs_I->coeff_0_0_0 = RATIONAL(1.0,1.0)+(-1.0+t86)*t83-t86;
-      coeffs_I->coeff_p1_0_0 = (t102*t79+t125*t80)*x+t96;
-      coeffs_I->coeff_m1_p1_0 = t90+t133+t135;
-      coeffs_I->coeff_0_p1_0 = (t101*t79+t116*t80)*y+t97;
-      coeffs_I->coeff_p1_p1_0 = t90+t144+t145;
-      coeffs_I->coeff_m1_m1_p1 = t95+t149+t150;
-      coeffs_I->coeff_0_m1_p1 = t94+t132+t146;
-      coeffs_I->coeff_p1_m1_p1 = t88+t147+t153;
-      coeffs_I->coeff_m1_0_p1 = t87+t141+t143;
-      coeffs_I->coeff_0_0_p1 = (t100*t79+t113*t80)*z+t92;
-      coeffs_I->coeff_p1_0_p1 = t87+t139+t142;
-      coeffs_I->coeff_m1_p1_p1 = t95+t148+t152;
-      coeffs_I->coeff_0_p1_p1 = t93+t132+t134;
-      coeffs_I->coeff_p1_p1_p1 = t89+t153+t154;
+      coeffs_I->coeff_p1_0_0 = (t100*t79+t112*t82)*x+t92;
+      coeffs_I->coeff_m1_p1_0 = t93+t130+t137;
+      coeffs_I->coeff_0_p1_0 = (t102*t79+t111*t82)*y+t87;
+      coeffs_I->coeff_p1_p1_0 = t89+t128+t130;
+      coeffs_I->coeff_m1_m1_p1 = t37+t27+t22+t40+t96;
+      coeffs_I->coeff_0_m1_p1 = t95+t133+t140;
+      coeffs_I->coeff_p1_m1_p1 = t42+t22+t26+t38+t98;
+      coeffs_I->coeff_m1_0_p1 = t94+t132+t138;
+      coeffs_I->coeff_0_0_p1 = (t101*t79+t110*t82)*z+t88;
+      coeffs_I->coeff_p1_0_p1 = t90+t129+t132;
+      coeffs_I->coeff_m1_p1_p1 = t41+t42+t40+t39+t99;
+      coeffs_I->coeff_0_p1_p1 = t91+t131+t133;
+      coeffs_I->coeff_p1_p1_p1 = t25+t26+t27+t28+t99;