1 files changed, 101 insertions, 74 deletions
diff --git a/src/Lagrange-tensor-product/2d.coeffs/2d.cube.order3.smooth0/coeffs-dxx.compute.c b/src/Lagrange-tensor-product/2d.coeffs/2d.cube.order3.smooth0/coeffs-dxx.compute.c
index 8c53891..fef7036 100644
--- a/src/Lagrange-tensor-product/2d.coeffs/2d.cube.order3.smooth0/coeffs-dxx.compute.c
+++ b/src/Lagrange-tensor-product/2d.coeffs/2d.cube.order3.smooth0/coeffs-dxx.compute.c
@@ -1,74 +1,101 @@
-fp t286;
-fp t287;
-fp t296;
-fp t282;
-fp t276;
-fp t284;
-fp t295;
-fp t285;
-fp t283;
-fp t294;
-fp t274;
-fp t293;
-fp t292;
-fp t291;
-fp t290;
-fp t289;
-fp t270;
-fp t269;
-fp t268;
-fp t266;
-fp t265;
-      t286 = y*y;
-      t287 = t286*y;
-      t296 = x*t287;
-      t282 = RATIONAL(1.0,4.0);
-      t276 = t282*t286;
-      t284 = RATIONAL(1.0,2.0);
-      t295 = t284*y+t276;
-      t285 = RATIONAL(-1.0,2.0);
-      t283 = RATIONAL(-1.0,4.0);
-      t294 = t283*y+t285*t286;
-      t274 = t283*t286;
-      t293 = t285*y+t274;
-      t292 = RATIONAL(3.0,25.0)*t296+RATIONAL(-1.0,200.0);
-      t291 = t284*t286+t282*y;
-      t290 = RATIONAL(-1.0,6.0)*y+t276;
-      t289 = RATIONAL(1.0,6.0)*y+t274;
-      t270 = RATIONAL(1.0,12.0)*y;
-      t269 = RATIONAL(-1.0,12.0)*y;
-      t268 = RATIONAL(-3.0,25.0)*t296;
-      t266 = RATIONAL(3.0,50.0)*t296;
-      t265 = RATIONAL(-3.0,50.0)*t296;
-      coeffs_dxx->coeff_m1_m1 = RATIONAL(-37.0,100.0)*x+RATIONAL(-17.0,150.0)*
-t287+t266+RATIONAL(37.0,200.0)+t290;
-      coeffs_dxx->coeff_0_m1 = RATIONAL(-111.0,200.0)+RATIONAL(13.0,75.0)*t287+
-(RATIONAL(-9.0,50.0)*t287+RATIONAL(111.0,100.0))*x+t289;
-      coeffs_dxx->coeff_p1_m1 = RATIONAL(-1.0,150.0)*t287+RATIONAL(111.0,200.0)
-+(RATIONAL(-111.0,100.0)+RATIONAL(9.0,50.0)*t287)*x+t289;
-      coeffs_dxx->coeff_p2_m1 = RATIONAL(-37.0,200.0)+RATIONAL(-4.0,75.0)*t287+
-t265+RATIONAL(37.0,100.0)*x+t290;
-      coeffs_dxx->coeff_m1_0 = RATIONAL(23.0,100.0)*t287+RATIONAL(133.0,200.0)+
-(RATIONAL(-33.0,100.0)+RATIONAL(1.0,25.0)*t287)*x+t294;
-      coeffs_dxx->coeff_0_0 = RATIONAL(-199.0,200.0)+RATIONAL(-19.0,100.0)*t287
-+RATIONAL(99.0,100.0)*x+t268+t291;
-      coeffs_dxx->coeff_p1_0 = RATIONAL(-99.0,100.0)*x+RATIONAL(-31.0,100.0)*
-t287+t291+t292;
-      coeffs_dxx->coeff_p2_0 = RATIONAL(67.0,200.0)+RATIONAL(27.0,100.0)*t287+(
-RATIONAL(-1.0,25.0)*t287+RATIONAL(33.0,100.0))*x+t294;
-      coeffs_dxx->coeff_m1_p1 = RATIONAL(-13.0,50.0)*t287+RATIONAL(29.0,200.0)+
-(RATIONAL(1.0,50.0)*t287+RATIONAL(-29.0,100.0))*x+t295;
-      coeffs_dxx->coeff_0_p1 = RATIONAL(-87.0,200.0)+t265+RATIONAL(87.0,100.0)*
-x+RATIONAL(7.0,25.0)*t287+t293;
-      coeffs_dxx->coeff_p1_p1 = RATIONAL(-87.0,100.0)*x+t266+RATIONAL(87.0,
-200.0)+RATIONAL(11.0,50.0)*t287+t293;
-      coeffs_dxx->coeff_p2_p1 = RATIONAL(-6.0,25.0)*t287+RATIONAL(-29.0,200.0)+
-(RATIONAL(-1.0,50.0)*t287+RATIONAL(29.0,100.0))*x+t295;
-      coeffs_dxx->coeff_m1_p2 = RATIONAL(43.0,300.0)*t287+RATIONAL(-1.0,100.0)*
-x+t268+t269+RATIONAL(1.0,200.0);
-      coeffs_dxx->coeff_0_p2 = RATIONAL(-79.0,300.0)*t287+RATIONAL(-3.0,200.0)+
-t270+(RATIONAL(3.0,100.0)+RATIONAL(9.0,25.0)*t287)*x;
-      coeffs_dxx->coeff_p1_p2 = t270+RATIONAL(29.0,300.0)*t287+RATIONAL(3.0,
-200.0)+(RATIONAL(-9.0,25.0)*t287+RATIONAL(-3.0,100.0))*x;
-      coeffs_dxx->coeff_p2_p2 = RATIONAL(1.0,100.0)*x+RATIONAL(7.0,300.0)*t287+
-t269+t292;
+fp t374;
+fp t375;
+fp t376;
+fp t388;
+fp t362;
+fp t358;
+fp t369;
+fp t361;
+fp t373;
+fp t356;
+fp t370;
+fp t387;
+fp t386;
+fp t357;
+fp t348;
+fp t363;
+fp t385;
+fp t372;
+fp t355;
+fp t384;
+fp t354;
+fp t345;
+fp t383;
+fp t367;
+fp t382;
+fp t360;
+fp t359;
+fp t381;
+fp t380;
+fp t379;
+fp t378;
+fp t371;
+fp t368;
+fp t366;
+fp t365;
+fp t353;
+fp t352;
+fp t351;
+fp t350;
+fp t346;
+fp t344;
+fp t343;
+      t374 = x*y;
+      t375 = y*y;
+      t376 = t375*y;
+      t388 = x*t376;
+      t362 = RATIONAL(-1.0,2.0);
+      t358 = t362*y;
+      t369 = RATIONAL(-1.0,3.0);
+      t361 = t369*y;
+      t373 = RATIONAL(-1.0,6.0);
+      t356 = t373*y;
+      t370 = RATIONAL(1.0,3.0);
+      t387 = t370*y;
+      t386 = x*t375;
+      t357 = t362*t376;
+      t348 = x*t357;
+      t363 = RATIONAL(1.0,2.0);
+      t385 = t363*t374+t348;
+      t372 = RATIONAL(1.0,6.0);
+      t355 = t372*t376;
+      t384 = t355+t356;
+      t354 = t363*t376;
+      t345 = x*t354;
+      t383 = t345+x*t358;
+      t367 = RATIONAL(3.0,2.0);
+      t382 = t367*t386-t375;
+      t360 = t373*t376;
+      t359 = t363*t375;
+      t381 = y+t357+t359;
+      t380 = t386-x;
+      t379 = t359+t361+t360;
+      t378 = t354+t358+RATIONAL(1.0,1.0)-t375;
+      t371 = RATIONAL(3.0,1.0);
+      t368 = RATIONAL(-3.0,1.0);
+      t366 = RATIONAL(-2.0,1.0);
+      t365 = RATIONAL(-3.0,2.0);
+      t353 = t367*t388;
+      t352 = t365*t388;
+      t351 = x*t355;
+      t350 = t362*t386;
+      t346 = t365*t386;
+      t344 = x*t360;
+      t343 = x*t359;
+      coeffs_dxx->coeff_m1_m1 = t350+t351+x*t387+t379;
+      coeffs_dxx->coeff_0_m1 = -t374+t370*t376+t348+RATIONAL(2.0,3.0)*y+t382;
+      coeffs_dxx->coeff_p1_m1 = t345+t374+t346+t379;
+      coeffs_dxx->coeff_p2_m1 = x*t361+t344+t343;
+      coeffs_dxx->coeff_m1_0 = t378+t380+t385;
+      coeffs_dxx->coeff_0_0 = -t376+y+RATIONAL(2.0,1.0)*t375+t366+t353+(t365*y+
+t371+t368*t375)*x;
+      coeffs_dxx->coeff_p1_0 = t352+(t371*t375+t368+t367*y)*x+t378;
+      coeffs_dxx->coeff_p2_0 = -t380+t383;
+      coeffs_dxx->coeff_m1_p1 = -t374+t350+t345+t381;
+      coeffs_dxx->coeff_0_p1 = t352+t376+(t371*x+t366)*y+t382;
+      coeffs_dxx->coeff_p1_p1 = t353+t346+t368*t374+t381;
+      coeffs_dxx->coeff_p2_p1 = t374+t343+t348;
+      coeffs_dxx->coeff_m1_p2 = t372*t374+t344+t384;
+      coeffs_dxx->coeff_0_p2 = t369*t376+t387+t383;
+      coeffs_dxx->coeff_p1_p2 = t384+t385;
+      coeffs_dxx->coeff_p2_p2 = t351+x*t356;