1 files changed, 101 insertions, 71 deletions
diff --git a/src/Lagrange-tensor-product/2d.coeffs/2d.cube.order3.smooth0/coeffs-dyy.compute.c b/src/Lagrange-tensor-product/2d.coeffs/2d.cube.order3.smooth0/coeffs-dyy.compute.c
index 1b690c6..a6f79fd 100644
--- a/src/Lagrange-tensor-product/2d.coeffs/2d.cube.order3.smooth0/coeffs-dyy.compute.c
+++ b/src/Lagrange-tensor-product/2d.coeffs/2d.cube.order3.smooth0/coeffs-dyy.compute.c
@@ -1,71 +1,101 @@
-fp t350;
-fp t351;
-fp t357;
-fp t342;
-fp t356;
-fp t343;
-fp t355;
-fp t348;
-fp t354;
-fp t353;
-fp t349;
-fp t345;
-fp t341;
-fp t340;
-fp t335;
-fp t334;
-fp t333;
-fp t332;
-      t350 = x*x;
-      t351 = t350*x;
-      t357 = y*t351;
-      t342 = RATIONAL(-1.0,4.0)*t350;
-      t356 = RATIONAL(7.0,20.0)*x+t342+RATIONAL(9.0,20.0);
-      t343 = RATIONAL(1.0,4.0)*t350;
-      t355 = RATIONAL(1.0,20.0)*x+t343+RATIONAL(-3.0,20.0);
-      t348 = RATIONAL(3.0,20.0);
-      t354 = t348*x+t342+RATIONAL(11.0,20.0);
-      t353 = t343+t348+RATIONAL(-11.0,20.0)*x;
-      t349 = RATIONAL(-3.0,10.0);
-      t345 = RATIONAL(-1.0,50.0);
-      t341 = RATIONAL(-1.0,2.0)*t350;
-      t340 = RATIONAL(1.0,2.0)*t350;
-      t335 = RATIONAL(3.0,50.0)*t357;
-      t334 = RATIONAL(-3.0,50.0)*t357;
-      t333 = RATIONAL(-3.0,25.0)*t357;
-      t332 = RATIONAL(3.0,25.0)*t357;
-      coeffs_dyy->coeff_m1_m1 = t335+(RATIONAL(59.0,125.0)*x+RATIONAL(-12.0,
-125.0)+RATIONAL(-17.0,50.0)*t350)*y+t353;
-      coeffs_dyy->coeff_0_m1 = (RATIONAL(-89.0,125.0)+RATIONAL(21.0,250.0)*x+
-RATIONAL(13.0,25.0)*t350+RATIONAL(-9.0,50.0)*t351)*y+t354;
-      coeffs_dyy->coeff_p1_m1 = (RATIONAL(-73.0,125.0)*x+RATIONAL(-36.0,125.0)+
-RATIONAL(9.0,50.0)*t351+t345*t350)*y+t356;
-      coeffs_dyy->coeff_p2_m1 = t334+(RATIONAL(12.0,125.0)+RATIONAL(-4.0,25.0)*
-t350+RATIONAL(7.0,250.0)*x)*y+t355;
-      coeffs_dyy->coeff_m1_0 = RATIONAL(11.0,10.0)*x+t349+t341+(RATIONAL(243.0,
-500.0)+RATIONAL(1.0,25.0)*t351+RATIONAL(69.0,100.0)*t350+RATIONAL(-851.0,500.0)
-*x)*y;
-      coeffs_dyy->coeff_0_0 = t340+t349*x+t333+RATIONAL(-11.0,10.0)+(RATIONAL(
--57.0,100.0)*t350+RATIONAL(771.0,500.0)+RATIONAL(303.0,500.0)*x)*y;
-      coeffs_dyy->coeff_p1_0 = t332+t340+RATIONAL(-9.0,10.0)+RATIONAL(-7.0,10.0
-)*x+(RATIONAL(729.0,500.0)+RATIONAL(-93.0,100.0)*t350+RATIONAL(447.0,500.0)*x)*
-y;
-      coeffs_dyy->coeff_p2_0 = t341+RATIONAL(3.0,10.0)+RATIONAL(-1.0,10.0)*x+(
-RATIONAL(-243.0,500.0)+RATIONAL(101.0,500.0)*x+RATIONAL(81.0,100.0)*t350+
-RATIONAL(-1.0,25.0)*t351)*y;
-      coeffs_dyy->coeff_m1_p1 = (RATIONAL(1.0,50.0)*t351+RATIONAL(-54.0,125.0)+
-RATIONAL(203.0,125.0)*x+RATIONAL(-39.0,50.0)*t350)*y+t353;
-      coeffs_dyy->coeff_0_p1 = t334+(RATIONAL(-213.0,125.0)+RATIONAL(21.0,25.0)
-*t350+RATIONAL(-93.0,250.0)*x)*y+t354;
-      coeffs_dyy->coeff_p1_p1 = t335+(RATIONAL(-141.0,125.0)*x+RATIONAL(33.0,
-50.0)*t350+RATIONAL(-162.0,125.0))*y+t356;
-      coeffs_dyy->coeff_p2_p1 = (RATIONAL(54.0,125.0)+RATIONAL(-18.0,25.0)*t350
-+RATIONAL(-31.0,250.0)*x+t345*t351)*y+t355;
-      coeffs_dyy->coeff_m1_p2 = t333+(RATIONAL(21.0,500.0)+RATIONAL(-197.0,
-500.0)*x+RATIONAL(43.0,100.0)*t350)*y;
-      coeffs_dyy->coeff_0_p2 = (RATIONAL(9.0,25.0)*t351+RATIONAL(-159.0,500.0)*
-x+RATIONAL(-79.0,100.0)*t350+RATIONAL(437.0,500.0))*y;
-      coeffs_dyy->coeff_p1_p2 = (RATIONAL(29.0,100.0)*t350+RATIONAL(-9.0,25.0)*
-t351+RATIONAL(63.0,500.0)+RATIONAL(409.0,500.0)*x)*y;
-      coeffs_dyy->coeff_p2_p2 = t332+(RATIONAL(7.0,100.0)*t350+RATIONAL(-53.0,
-500.0)*x+RATIONAL(-21.0,500.0))*y;
+fp t493;
+fp t494;
+fp t508;
+fp t481;
+fp t507;
+fp t488;
+fp t477;
+fp t495;
+fp t506;
+fp t489;
+fp t505;
+fp t492;
+fp t480;
+fp t482;
+fp t476;
+fp t473;
+fp t491;
+fp t478;
+fp t475;
+fp t479;
+fp t486;
+fp t504;
+fp t503;
+fp t472;
+fp t502;
+fp t462;
+fp t501;
+fp t500;
+fp t499;
+fp t498;
+fp t497;
+fp t490;
+fp t487;
+fp t485;
+fp t484;
+fp t471;
+fp t469;
+fp t468;
+fp t466;
+fp t465;
+fp t464;
+fp t463;
+      t493 = x*y;
+      t494 = x*x;
+      t508 = t494*y;
+      t481 = RATIONAL(-1.0,2.0);
+      t507 = y*t481;
+      t488 = RATIONAL(-1.0,3.0);
+      t477 = t488*x;
+      t495 = x*t494;
+      t506 = t495*y;
+      t489 = RATIONAL(1.0,3.0);
+      t505 = t489*x;
+      t492 = RATIONAL(-1.0,6.0);
+      t480 = t492*x;
+      t482 = RATIONAL(1.0,2.0);
+      t476 = t482*t494;
+      t473 = t482*t495;
+      t491 = RATIONAL(1.0,6.0);
+      t478 = t491*t495;
+      t475 = t492*t495;
+      t479 = t481*t495;
+      t486 = RATIONAL(3.0,2.0);
+      t504 = -t494+t486*t508;
+      t503 = t480+t478;
+      t472 = y*t473;
+      t502 = x*t507+t472;
+      t462 = y*t479;
+      t501 = t462+t482*t493;
+      t500 = x+t479+t476;
+      t499 = t508-y;
+      t498 = t475+t477+t476;
+      t497 = -t494+RATIONAL(1.0,1.0)+t473+t481*x;
+      t490 = RATIONAL(3.0,1.0);
+      t487 = RATIONAL(-3.0,1.0);
+      t485 = RATIONAL(-2.0,1.0);
+      t484 = RATIONAL(-3.0,2.0);
+      t471 = y*t475;
+      t469 = t486*t506;
+      t468 = t484*t506;
+      t466 = t484*t508;
+      t465 = t494*t507;
+      t464 = y*t478;
+      t463 = y*t476;
+      coeffs_dyy->coeff_m1_m1 = t464+t465+y*t505+t498;
+      coeffs_dyy->coeff_0_m1 = t497+t499+t501;
+      coeffs_dyy->coeff_p1_m1 = t465-t493+t472+t500;
+      coeffs_dyy->coeff_p2_m1 = t491*t493+t471+t503;
+      coeffs_dyy->coeff_m1_0 = RATIONAL(2.0,3.0)*x-t493+t462+t489*t495+t504;
+      coeffs_dyy->coeff_0_0 = x+RATIONAL(2.0,1.0)*t494+t485+t469-t495+(t484*x+
+t490+t487*t494)*y;
+      coeffs_dyy->coeff_p1_0 = t495+t468+(t485+t490*y)*x+t504;
+      coeffs_dyy->coeff_p2_0 = t505+t488*t495+t502;
+      coeffs_dyy->coeff_m1_p1 = t472+t493+t466+t498;
+      coeffs_dyy->coeff_0_p1 = t468+(t490*t494+t487+t486*x)*y+t497;
+      coeffs_dyy->coeff_p1_p1 = t466+t469+t487*t493+t500;
+      coeffs_dyy->coeff_p2_p1 = t501+t503;
+      coeffs_dyy->coeff_m1_p2 = y*t477+t471+t463;
+      coeffs_dyy->coeff_0_p2 = -t499+t502;
+      coeffs_dyy->coeff_p1_p2 = t463+t493+t462;
+      coeffs_dyy->coeff_p2_p2 = y*t480+t464;