1 files changed, 151 insertions, 133 deletions
diff --git a/src/Lagrange-tensor-product/3d.coeffs/3d.cube.order2.smooth0/coeffs-dxz.compute.c b/src/Lagrange-tensor-product/3d.coeffs/3d.cube.order2.smooth0/coeffs-dxz.compute.c
index 66fe899..48bfb8a 100644
--- a/src/Lagrange-tensor-product/3d.coeffs/3d.cube.order2.smooth0/coeffs-dxz.compute.c
+++ b/src/Lagrange-tensor-product/3d.coeffs/3d.cube.order2.smooth0/coeffs-dxz.compute.c
@@ -1,152 +1,170 @@
-fp t596;
-fp t602;
-fp t597;
-fp t618;
-fp t565;
-fp t629;
-fp t593;
-fp t587;
+fp t591;
+fp t598;
+fp t612;
+fp t567;
+fp t558;
+fp t589;
+fp t630;
 fp t592;
+fp t629;
+fp t581;
+fp t566;
+fp t559;
+fp t590;
+fp t580;
+fp t561;
 fp t628;
-fp t617;
-fp t576;
+fp t593;
+fp t585;
+fp t569;
 fp t627;
-fp t583;
+fp t594;
+fp t586;
+fp t574;
+fp t560;
 fp t626;
-fp t624;
-fp t616;
-fp t601;
 fp t625;
+fp t563;
+fp t579;
+fp t624;
+fp t584;
 fp t623;
+fp t565;
 fp t622;
+fp t587;
 fp t621;
-fp t620;
 fp t571;
-fp t615;
-fp t562;
+fp t620;
 fp t619;
-fp t581;
-fp t567;
-fp t614;
-fp t613;
-fp t600;
-fp t612;
 fp t611;
-fp t610;
-fp t609;
-fp t598;
-fp t589;
-fp t590;
-fp t564;
-fp t561;
-fp t608;
-fp t599;
-fp t591;
 fp t573;
-fp t569;
-fp t607;
+fp t618;
+fp t575;
+fp t582;
+fp t617;
+fp t597;
+fp t613;
+fp t609;
+fp t616;
+fp t615;
+fp t614;
+fp t596;
 fp t588;
+fp t557;
 fp t577;
-fp t566;
+fp t610;
+fp t608;
+fp t607;
 fp t606;
-fp t578;
-fp t580;
-fp t563;
 fp t605;
-fp t572;
+fp t562;
+fp t564;
 fp t604;
-fp t570;
+fp t572;
 fp t603;
-fp t579;
-fp t575;
-fp t574;
-      t596 = RATIONAL(1.0,2.0);
-      t602 = y*y;
-      t597 = RATIONAL(-1.0,2.0);
-      t618 = t602*t597;
-      t565 = x*t618;
-      t629 = t565+t596*x;
+fp t570;
+fp t602;
+fp t601;
+fp t600;
+fp t599;
+fp t576;
+fp t568;
+      t591 = RATIONAL(1.0,2.0);
+      t598 = y*y;
+      t612 = t591*t598;
+      t567 = x*t612;
+      t558 = z*t567;
+      t589 = RATIONAL(1.0,8.0);
+      t630 = t558+t589*t598;
+      t592 = RATIONAL(-1.0,8.0);
+      t629 = t558+t592*t598;
+      t581 = t591*z;
+      t566 = y*t581;
+      t559 = x*t566;
+      t590 = RATIONAL(1.0,4.0);
+      t580 = t590*t598;
+      t561 = z*t580;
+      t628 = t559+t561;
       t593 = RATIONAL(-1.0,4.0);
-      t587 = t593*t602;
-      t592 = RATIONAL(1.0,4.0);
-      t628 = t587+t592;
-      t617 = t602*t596;
-      t576 = x*t617;
-      t627 = t576+t597*x;
-      t583 = t592*t602;
-      t626 = t593+t583;
-      t624 = t602*x;
-      t616 = z*t624;
-      t601 = x*z;
-      t625 = -t616+t601;
-      t623 = y*t593;
-      t622 = y*t592;
-      t621 = y*t597;
-      t620 = y*t596;
-      t571 = RATIONAL(-2.0,1.0)*t601;
-      t615 = RATIONAL(2.0,1.0)*t601;
-      t562 = t602*t615;
-      t619 = t571+t562;
-      t581 = z*t620;
-      t567 = z*t621;
-      t614 = (1.0-t602)*z;
-      t613 = x-t624;
-      t600 = y*t601;
-      t612 = -t600-t616;
-      t611 = t600-t616;
-      t610 = -t616+x*t620;
-      t609 = -t616+t567;
-      t598 = RATIONAL(1.0,8.0);
-      t589 = t598*y;
-      t590 = t598*t602;
-      t564 = z*t576;
-      t561 = x*t581;
-      t608 = t589+t590+t564+t561;
-      t599 = RATIONAL(-1.0,8.0);
-      t591 = t599*y;
-      t573 = z*t583;
-      t569 = x*t623;
-      t607 = t591+t573+t569+t564;
-      t588 = t599*t602;
-      t577 = z*t587;
-      t566 = x*t583;
-      t606 = t588+t577+t564+t566;
-      t578 = x*t622;
-      t580 = x*t587;
-      t563 = x*t567;
-      t605 = t564+t578+t580+t563;
-      t572 = z*t617;
-      t604 = t572+t597*z+t625;
-      t570 = z*t618;
-      t603 = t570+t596*z+t625;
-      t579 = z*t622;
-      t575 = x*t621;
-      t574 = z*t623;
-      coeffs_dxz->coeff_m1_m1_m1 = t579+t591+t590+t577+t605;
-      coeffs_dxz->coeff_0_m1_m1 = t575+t576+t611;
-      coeffs_dxz->coeff_p1_m1_m1 = t589+t574+t588+t573+t605;
-      coeffs_dxz->coeff_m1_0_m1 = t604+t627+t628;
-      coeffs_dxz->coeff_0_0_m1 = t613+t619;
-      coeffs_dxz->coeff_p1_0_m1 = t603+t626+t627;
-      coeffs_dxz->coeff_m1_p1_m1 = t580+t574+t569+t577+t608;
-      coeffs_dxz->coeff_0_p1_m1 = -t600+t576+t610;
-      coeffs_dxz->coeff_p1_p1_m1 = t588+t580+t561+t579+t607;
-      coeffs_dxz->coeff_m1_m1_0 = t572+t600+t609;
-      coeffs_dxz->coeff_0_m1_0 = t562+y*t571;
-      coeffs_dxz->coeff_p1_m1_0 = t581+t570+t611;
-      coeffs_dxz->coeff_m1_0_0 = t614+t619;
-      coeffs_dxz->coeff_0_0_0 = (RATIONAL(-4.0,1.0)*t602+RATIONAL(4.0,1.0))*
-t601;
-      coeffs_dxz->coeff_p1_0_0 = -t614+t619;
-      coeffs_dxz->coeff_m1_p1_0 = t581+t572+t612;
-      coeffs_dxz->coeff_0_p1_0 = t562+y*t615;
-      coeffs_dxz->coeff_p1_p1_0 = t570-t600+t609;
-      coeffs_dxz->coeff_m1_m1_p1 = t579+t569+t563+t589+t606;
-      coeffs_dxz->coeff_0_m1_p1 = t565+t600+t610;
-      coeffs_dxz->coeff_p1_m1_p1 = t566+t574+t563+t590+t607;
-      coeffs_dxz->coeff_m1_0_p1 = t604+t626+t629;
-      coeffs_dxz->coeff_0_0_p1 = -t613+t619;
-      coeffs_dxz->coeff_p1_0_p1 = t603+t628+t629;
-      coeffs_dxz->coeff_m1_p1_p1 = t574+t591+t561+t578+t606;
-      coeffs_dxz->coeff_0_p1_p1 = t575+t565+t612;
-      coeffs_dxz->coeff_p1_p1_p1 = t579+t566+t578+t573+t608;
+      t585 = t593*t598;
+      t569 = z*t585;
+      t627 = t559+t569;
+      t594 = RATIONAL(-1.0,2.0);
+      t586 = t594*z;
+      t574 = y*t586;
+      t560 = x*t574;
+      t626 = t560+t561;
+      t625 = t560+t569;
+      t563 = x*t580;
+      t579 = t589*y;
+      t624 = t563+t579;
+      t584 = t592*y;
+      t623 = t563+t584;
+      t565 = z*t612;
+      t622 = t565+t586;
+      t587 = t594*x;
+      t621 = t567+t587;
+      t571 = x*t585;
+      t620 = t571+t579;
+      t619 = t571+t584;
+      t611 = t594*t598;
+      t573 = z*t611;
+      t618 = t573+t581;
+      t575 = x*t611;
+      t582 = t591*x;
+      t617 = t575+t582;
+      t597 = x*z;
+      t613 = x*t598;
+      t609 = z*t613;
+      t616 = t597-t609;
+      t615 = t590*y;
+      t614 = t593*y;
+      t596 = y*t597;
+      t588 = RATIONAL(2.0,1.0);
+      t557 = t588*t609;
+      t577 = RATIONAL(-2.0,1.0)*t597;
+      t610 = t557+t577;
+      t608 = x-t613;
+      t607 = (1.0-t598)*z;
+      t606 = -t596-t609;
+      t605 = t596-t609;
+      t562 = z*t615;
+      t564 = x*t615;
+      t604 = t562+t564+t630;
+      t572 = x*t614;
+      t603 = t562+t572+t629;
+      t570 = z*t614;
+      t602 = t564+t570+t629;
+      t601 = t570+t572+t630;
+      t600 = t580+t593+t616;
+      t599 = t585+t590+t616;
+      t576 = y*t587;
+      t568 = y*t582;
+      coeffs_dxz->coeff_m1_m1_m1 = t604+t619+t625;
+      coeffs_dxz->coeff_0_m1_m1 = t567+t576+t605;
+      coeffs_dxz->coeff_p1_m1_m1 = t602+t620+t626;
+      coeffs_dxz->coeff_m1_0_m1 = t599+t621+t622;
+      coeffs_dxz->coeff_0_0_m1 = t608+t610;
+      coeffs_dxz->coeff_p1_0_m1 = t600+t618+t621;
+      coeffs_dxz->coeff_m1_p1_m1 = t601+t620+t627;
+      coeffs_dxz->coeff_0_p1_m1 = t567+t568+t606;
+      coeffs_dxz->coeff_p1_p1_m1 = t603+t619+t628;
+      coeffs_dxz->coeff_m1_m1_0 = t565+t574+t605;
+      coeffs_dxz->coeff_0_m1_0 = y*t577+t557;
+      coeffs_dxz->coeff_p1_m1_0 = t573+t566+t605;
+      coeffs_dxz->coeff_m1_0_0 = t607+t610;
+      coeffs_dxz->coeff_0_0_0 = (RATIONAL(4.0,1.0)+RATIONAL(-4.0,1.0)*t598)*
+t597;
+      coeffs_dxz->coeff_p1_0_0 = -t607+t610;
+      coeffs_dxz->coeff_m1_p1_0 = t566+t565+t606;
+      coeffs_dxz->coeff_0_p1_0 = t557+t588*t596;
+      coeffs_dxz->coeff_p1_p1_0 = t573+t574+t606;
+      coeffs_dxz->coeff_m1_m1_p1 = t603+t624+t625;
+      coeffs_dxz->coeff_0_m1_p1 = t568+t575+t605;
+      coeffs_dxz->coeff_p1_m1_p1 = t601+t623+t626;
+      coeffs_dxz->coeff_m1_0_p1 = t600+t617+t622;
+      coeffs_dxz->coeff_0_0_p1 = -t608+t610;
+      coeffs_dxz->coeff_p1_0_p1 = t599+t617+t618;
+      coeffs_dxz->coeff_m1_p1_p1 = t602+t623+t627;
+      coeffs_dxz->coeff_0_p1_p1 = t575+t576+t606;
+      coeffs_dxz->coeff_p1_p1_p1 = t604+t624+t628;