This commit reorganizes the Maple scripts which generate the interpolation

coefficients, and the coefficient files themselves, to properly support
multiple interpolation operators, and adds some skeleton support for
Hermite interpolation.  I will add full support for Hermite in a future
checkin.


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

6 files changed, 882 insertions, 0 deletions

diff --git a/src/GeneralizedPolynomial-Uniform/Lagrange/2d.coeffs/2d.cube.order4.smooth0/coeff-I.compute.c b/src/GeneralizedPolynomial-Uniform/Lagrange/2d.coeffs/2d.cube.order4.smooth0/coeff-I.compute.c
new file mode 100644
index 0000000..49655f0
--- /dev/null
+++ b/src/GeneralizedPolynomial-Uniform/Lagrange/2d.coeffs/2d.cube.order4.smooth0/coeff-I.compute.c
@@ -0,0 +1,239 @@
+fp t413,
+   t418,
+   t466,
+   t403,
+   t417,
+   t358,
+   t414,
+   t373,
+   t438,
+   t444,
+   t465,
+   t405,
+   t364,
+   t367,
+   t439,
+   t464,
+   t436,
+   t463,
+   t410,
+   t359,
+   t442,
+   t409,
+   t318,
+   t392,
+   t462,
+   t411,
+   t360,
+   t441,
+   t400,
+   t298,
+   t461,
+   t366,
+   t443,
+   t395,
+   t401,
+   t460,
+   t440,
+   t296,
+   t399,
+   t459,
+   t372,
+   t458,
+   t333,
+   t457,
+   t396,
+   t412,
+   t327,
+   t456,
+   t407,
+   t331,
+   t391,
+   t455,
+   t402,
+   t313,
+   t397,
+   t454,
+   t325,
+   t453,
+   t324,
+   t452,
+   t311,
+   t451,
+   t450,
+   t387,
+   t398,
+   t449,
+   t415,
+   t371,
+   t416,
+   t357,
+   t383,
+   t382,
+   t448,
+   t317,
+   t356,
+   t447,
+   t362,
+   t446,
+   t445,
+   t437,
+   t435,
+   t393,
+   t404,
+   t369,
+   t376,
+   t434,
+   t365,
+   t384,
+   t386,
+   t433,
+   t432,
+   t431,
+   t430,
+   t429,
+   t428,
+   t427,
+   t406,
+   t426,
+   t425,
+   t408,
+   t424,
+   t423,
+   t394,
+   t389,
+   t388,
+   t385;
+      t413 = x*x;
+      t418 = y*y;
+      t466 = t418+t413;
+      t403 = RATIONAL(-1.0,60.0);
+      t417 = y*t418;
+      t358 = t403*t417;
+      t414 = t413*x;
+      t373 = t403*t414;
+      t438 = y*t373+x*t358;
+      t444 = x*y;
+      t465 = t438+RATIONAL(37.0,300.0)*t444;
+      t405 = RATIONAL(1.0,60.0);
+      t364 = t405*t414;
+      t367 = t405*t417;
+      t439 = y*t364+x*t367;
+      t464 = t439+RATIONAL(-37.0,300.0)*t444;
+      t436 = t413*t418;
+      t463 = RATIONAL(-99.0,1225.0)+RATIONAL(-1.0,98.0)*t436;
+      t410 = RATIONAL(-1.0,30.0);
+      t359 = t410*t417;
+      t442 = y*t413;
+      t409 = RATIONAL(1.0,70.0);
+      t318 = t409*t442;
+      t392 = RATIONAL(11.0,105.0);
+      t462 = t359+t318+t392*y;
+      t411 = RATIONAL(1.0,30.0);
+      t360 = t411*t417;
+      t441 = y*t414;
+      t400 = RATIONAL(-1.0,120.0);
+      t298 = RATIONAL(-21.0,200.0)*t444;
+      t461 = x*t360+t298+t400*t441;
+      t366 = t410*t414;
+      t443 = x*t418;
+      t395 = RATIONAL(31.0,210.0);
+      t401 = RATIONAL(-1.0,140.0);
+      t460 = t395*x+t366+t401*t443;
+      t440 = x*t417;
+      t296 = RATIONAL(21.0,200.0)*t444;
+      t399 = RATIONAL(1.0,120.0);
+      t459 = t296+y*t366+t399*t440;
+      t372 = t411*t414;
+      t458 = t400*t440+y*t372+t298;
+      t333 = t409*t443;
+      t457 = t366+t333+t392*x;
+      t396 = RATIONAL(-11.0,105.0);
+      t412 = RATIONAL(-1.0,70.0);
+      t327 = t412*t443;
+      t456 = t396*x+t372+t327;
+      t407 = RATIONAL(1.0,35.0);
+      t331 = t407*t442;
+      t391 = RATIONAL(-31.0,420.0);
+      t455 = t331+t391*y+t367;
+      t402 = RATIONAL(-1.0,35.0);
+      t313 = t402*t442;
+      t397 = RATIONAL(31.0,420.0);
+      t454 = t358+t397*y+t313;
+      t325 = t402*t443;
+      t453 = t397*x+t325+t373;
+      t324 = t407*t443;
+      t452 = t364+t324+t391*x;
+      t311 = t412*t442;
+      t451 = t360+t311+t396*y;
+      t450 = x*t359+t296+t399*t441;
+      t387 = RATIONAL(-31.0,210.0);
+      t398 = RATIONAL(1.0,140.0);
+      t449 = t372+t387*x+t398*t443;
+      t415 = t413*t413;
+      t371 = t410*t415;
+      t416 = t418*t418;
+      t357 = t399*t416;
+      t383 = RATIONAL(71.0,5880.0);
+      t382 = RATIONAL(113.0,735.0);
+      t448 = t371+t382*t413+t383*t418+t357+t463;
+      t317 = RATIONAL(1.0,49.0)*t436;
+      t356 = t399*t415;
+      t447 = t317+RATIONAL(51.0,1225.0)+t356+t357+t466*RATIONAL(-289.0,5880.0);
+      t362 = t410*t416;
+      t446 = t382*t418+t356+t383*t413+t362+t463;
+      t445 = t362+RATIONAL(1.0,196.0)*t436+t371+RATIONAL(-24.0,1225.0)+t466*
+RATIONAL(181.0,1470.0);
+      t437 = RATIONAL(-1.0,49.0)*t436+RATIONAL(96.0,1225.0);
+      t435 = RATIONAL(1.0,98.0)*t436+RATIONAL(246.0,1225.0);
+      t393 = RATIONAL(-41.0,196.0);
+      t404 = RATIONAL(1.0,20.0);
+      t369 = t404*t415;
+      t376 = RATIONAL(191.0,5880.0);
+      t434 = t376*t418+t369+t357+t393*t413+t437;
+      t365 = t404*t416;
+      t384 = RATIONAL(-53.0,196.0);
+      t386 = RATIONAL(83.0,735.0);
+      t433 = t365+t371+t384*t418+t386*t413+t435;
+      t432 = t365+t376*t413+t356+t393*t418+t437;
+      t431 = t369+t384*t413+t362+t386*t418+t435;
+      t430 = RATIONAL(-11.0,150.0)*t444+t439+t447;
+      t429 = t359+t395*y+t401*t442+t445;
+      t428 = t360+t387*y+t398*t442+t445;
+      t427 = RATIONAL(11.0,150.0)*t444+t438+t447;
+      t406 = RATIONAL(1.0,84.0);
+      t426 = t367+t311+t406*y+t448;
+      t425 = t364+t406*x+t327+t446;
+      t408 = RATIONAL(-1.0,84.0);
+      t424 = t358+t318+t408*y+t448;
+      t423 = t333+t408*x+t373+t446;
+      t394 = RATIONAL(-17.0,420.0);
+      t389 = RATIONAL(17.0,105.0);
+      t388 = RATIONAL(17.0,420.0);
+      t385 = RATIONAL(-17.0,105.0);
+      coeff_I_m2_m2 = t430+t453+t454;
+      coeff_I_m1_m2 = t424+t456+t459;
+      coeff_I_0_m2 = t331+t394*y+t358+t434;
+      coeff_I_p1_m2 = t424+t457+t458;
+      coeff_I_p2_m2 = t427+t452+t454;
+      coeff_I_m2_m1 = t423+t450+t451;
+      coeff_I_m1_m1 = t428+t449+t465;
+      coeff_I_0_m1 = t360+t385*y+t318+t431;
+      coeff_I_p1_m1 = t428+t460+t464;
+      coeff_I_p2_m1 = t425+t451+t461;
+      coeff_I_m2_0 = t373+t394*x+t324+t432;
+      coeff_I_m1_0 = t385*x+t333+t372+t433;
+      coeff_I_0_0 = t317+t365+t369+RATIONAL(541.0,1225.0)+t466*RATIONAL(-57.0,
+196.0);
+      coeff_I_p1_0 = t389*x+t327+t366+t433;
+      coeff_I_p2_0 = t325+t364+t388*x+t432;
+      coeff_I_m2_p1 = t423+t461+t462;
+      coeff_I_m1_p1 = t429+t449+t464;
+      coeff_I_0_p1 = t359+t389*y+t311+t431;
+      coeff_I_p1_p1 = t429+t460+t465;
+      coeff_I_p2_p1 = t425+t450+t462;
+      coeff_I_m2_p2 = t427+t453+t455;
+      coeff_I_m1_p2 = t426+t456+t458;
+      coeff_I_0_p2 = t313+t388*y+t367+t434;
+      coeff_I_p1_p2 = t426+t457+t459;
+      coeff_I_p2_p2 = t430+t452+t455;
diff --git a/src/GeneralizedPolynomial-Uniform/Lagrange/2d.coeffs/2d.cube.order4.smooth0/coeff-dx.compute.c b/src/GeneralizedPolynomial-Uniform/Lagrange/2d.coeffs/2d.cube.order4.smooth0/coeff-dx.compute.c
new file mode 100644
index 0000000..83ba345
--- /dev/null
+++ b/src/GeneralizedPolynomial-Uniform/Lagrange/2d.coeffs/2d.cube.order4.smooth0/coeff-dx.compute.c
@@ -0,0 +1,147 @@
+fp t531,
+   t534,
+   t535,
+   t513,
+   t570,
+   t508,
+   t569,
+   t537,
+   t536,
+   t502,
+   t552,
+   t568,
+   t512,
+   t555,
+   t567,
+   t556,
+   t566,
+   t509,
+   t533,
+   t510,
+   t565,
+   t530,
+   t506,
+   t507,
+   t564,
+   t563,
+   t471,
+   t562,
+   t561,
+   t523,
+   t511,
+   t560,
+   t528,
+   t504,
+   t559,
+   t558,
+   t557,
+   t478,
+   t554,
+   t553,
+   t516,
+   t476,
+   t551,
+   t483,
+   t550,
+   t549,
+   t548,
+   t547,
+   t546,
+   t545,
+   t544,
+   t543,
+   t542,
+   t541,
+   t540,
+   t484,
+   t477,
+   t475,
+   t473,
+   t469,
+   t467;
+      t531 = RATIONAL(1.0,30.0);
+      t534 = x*x;
+      t535 = t534*x;
+      t513 = t531*t535;
+      t570 = t513+RATIONAL(71.0,2940.0)*x;
+      t508 = RATIONAL(-2.0,15.0)*t535;
+      t569 = t508+RATIONAL(226.0,735.0)*x;
+      t537 = y*y;
+      t536 = y*t537;
+      t502 = RATIONAL(1.0,20.0)*t534;
+      t552 = RATIONAL(1.0,60.0)*t536+y*t502;
+      t568 = t552+RATIONAL(-37.0,300.0)*y;
+      t512 = RATIONAL(-1.0,20.0)*t534;
+      t555 = y*t512+RATIONAL(-1.0,60.0)*t536;
+      t567 = t555+RATIONAL(37.0,300.0)*y;
+      t556 = x*t537;
+      t566 = t508+RATIONAL(181.0,735.0)*x+RATIONAL(1.0,98.0)*t556;
+      t509 = RATIONAL(1.0,10.0)*t534;
+      t533 = RATIONAL(-1.0,70.0);
+      t510 = t533*t537;
+      t565 = t509+t510+RATIONAL(-11.0,105.0);
+      t530 = RATIONAL(1.0,70.0);
+      t506 = t530*t537;
+      t507 = RATIONAL(-1.0,10.0)*t534;
+      t564 = t506+t507+RATIONAL(11.0,105.0);
+      t563 = t510+RATIONAL(1.0,84.0)+t502;
+      t471 = RATIONAL(2.0,49.0)*t556;
+      t562 = t513+t471+RATIONAL(-289.0,2940.0)*x;
+      t561 = t506+t512+RATIONAL(-1.0,84.0);
+      t523 = RATIONAL(-1.0,35.0);
+      t511 = t523*t537;
+      t560 = t511+t512+RATIONAL(31.0,420.0);
+      t528 = RATIONAL(1.0,35.0);
+      t504 = t528*t537;
+      t559 = RATIONAL(-31.0,420.0)+t502+t504;
+      t558 = x*y;
+      t557 = y*t534;
+      t478 = RATIONAL(-1.0,49.0)*t556;
+      t554 = RATIONAL(-21.0,200.0)*y+t478;
+      t553 = RATIONAL(21.0,200.0)*y+t478;
+      t516 = RATIONAL(1.0,5.0)*t535;
+      t476 = RATIONAL(-2.0,49.0)*t556;
+      t551 = t516+t476+RATIONAL(-41.0,98.0)*x;
+      t483 = RATIONAL(1.0,49.0)*t556;
+      t550 = t516+t483+RATIONAL(-53.0,98.0)*x;
+      t549 = t508+RATIONAL(166.0,735.0)*x+t483;
+      t548 = t513+t476+RATIONAL(191.0,2940.0)*x;
+      t547 = t509+RATIONAL(-31.0,210.0)+RATIONAL(1.0,140.0)*t537+t566;
+      t546 = RATIONAL(11.0,150.0)*y+t555+t562;
+      t545 = y*t509+RATIONAL(-1.0,120.0)*t536+t554+t569;
+      t544 = t507+RATIONAL(31.0,210.0)+RATIONAL(-1.0,140.0)*t537+t566;
+      t543 = y*t507+RATIONAL(1.0,120.0)*t536+t553+t569;
+      t542 = t531*t536+RATIONAL(-1.0,40.0)*t557+t554+t570;
+      t541 = RATIONAL(-11.0,150.0)*y+t552+t562;
+      t540 = RATIONAL(1.0,40.0)*t557+RATIONAL(-1.0,30.0)*t536+t553+t570;
+      t484 = t530*t558;
+      t477 = t528*t558;
+      t475 = RATIONAL(2.0,35.0)*t558;
+      t473 = t533*t558;
+      t469 = RATIONAL(-2.0,35.0)*t558;
+      t467 = t523*t558;
+      coeff_dx_m2_m2 = t469+t541+t560;
+      coeff_dx_m1_m2 = t477+t543+t565;
+      coeff_dx_0_m2 = t475+t551;
+      coeff_dx_p1_m2 = t477+t545+t564;
+      coeff_dx_p2_m2 = t469+t546+t559;
+      coeff_dx_m2_m1 = t467+t540+t561;
+      coeff_dx_m1_m1 = t484+t547+t567;
+      coeff_dx_0_m1 = t477+t550;
+      coeff_dx_p1_m1 = t484+t544+t568;
+      coeff_dx_p2_m1 = t467+t542+t563;
+      coeff_dx_m2_0 = t504+t512+RATIONAL(-17.0,420.0)+t548;
+      coeff_dx_m1_0 = RATIONAL(-17.0,105.0)+t509+t506+t549;
+      coeff_dx_0_0 = RATIONAL(-57.0,98.0)*x+t516+t471;
+      coeff_dx_p1_0 = t510+RATIONAL(17.0,105.0)+t507+t549;
+      coeff_dx_p2_0 = t511+RATIONAL(17.0,420.0)+t502+t548;
+      coeff_dx_m2_p1 = t477+t542+t561;
+      coeff_dx_m1_p1 = t473+t547+t568;
+      coeff_dx_0_p1 = t467+t550;
+      coeff_dx_p1_p1 = t473+t544+t567;
+      coeff_dx_p2_p1 = t477+t540+t563;
+      coeff_dx_m2_p2 = t475+t546+t560;
+      coeff_dx_m1_p2 = t467+t545+t565;
+      coeff_dx_0_p2 = t469+t551;
+      coeff_dx_p1_p2 = t467+t543+t564;
+      coeff_dx_p2_p2 = t475+t541+t559;
diff --git a/src/GeneralizedPolynomial-Uniform/Lagrange/2d.coeffs/2d.cube.order4.smooth0/coeff-dxx.compute.c b/src/GeneralizedPolynomial-Uniform/Lagrange/2d.coeffs/2d.cube.order4.smooth0/coeff-dxx.compute.c
new file mode 100644
index 0000000..385d8a7
--- /dev/null
+++ b/src/GeneralizedPolynomial-Uniform/Lagrange/2d.coeffs/2d.cube.order4.smooth0/coeff-dxx.compute.c
@@ -0,0 +1,111 @@
+fp t712,
+   t692,
+   t707,
+   t711,
+   t693,
+   t729,
+   t687,
+   t728,
+   t697,
+   t727,
+   t726,
+   t725,
+   t686,
+   t694,
+   t695,
+   t698,
+   t689,
+   t724,
+   t681,
+   t723,
+   t696,
+   t722,
+   t721,
+   t688,
+   t720,
+   t691,
+   t719,
+   t685,
+   t718,
+   t717,
+   t716,
+   t690,
+   t715,
+   t714,
+   t713,
+   t684,
+   t683,
+   t680,
+   t679,
+   t678,
+   t677,
+   t676,
+   t675;
+      t712 = y*y;
+      t692 = RATIONAL(2.0,49.0)*t712;
+      t707 = RATIONAL(1.0,10.0);
+      t711 = x*x;
+      t693 = t707*t711;
+      t729 = t692+t693+RATIONAL(-289.0,2940.0);
+      t687 = RATIONAL(-1.0,49.0)*t712;
+      t728 = t693+RATIONAL(71.0,2940.0)+t687;
+      t697 = RATIONAL(-2.0,5.0)*t711;
+      t727 = t697+RATIONAL(226.0,735.0)+t687;
+      t726 = t697+RATIONAL(181.0,735.0)+RATIONAL(1.0,98.0)*t712;
+      t725 = x*y;
+      t686 = t707*x;
+      t694 = RATIONAL(-1.0,10.0)*x;
+      t695 = RATIONAL(-1.0,5.0)*x;
+      t698 = RATIONAL(1.0,5.0)*x;
+      t689 = RATIONAL(-2.0,49.0)*t712;
+      t724 = t693+RATIONAL(191.0,2940.0)+t689;
+      t681 = RATIONAL(1.0,49.0)*t712;
+      t723 = t697+RATIONAL(166.0,735.0)+t681;
+      t696 = RATIONAL(3.0,5.0)*t711;
+      t722 = t696+RATIONAL(-53.0,98.0)+t681;
+      t721 = t696+RATIONAL(-41.0,98.0)+t689;
+      t688 = RATIONAL(2.0,35.0)*y;
+      t720 = t688+t729;
+      t691 = RATIONAL(-2.0,35.0)*y;
+      t719 = t691+t729;
+      t685 = RATIONAL(-1.0,35.0)*y;
+      t718 = t685+t727;
+      t717 = t686+t728;
+      t716 = t695+t726;
+      t690 = RATIONAL(1.0,35.0)*y;
+      t715 = t690+t727;
+      t714 = t698+t726;
+      t713 = t694+t728;
+      t684 = RATIONAL(1.0,70.0)*y;
+      t683 = RATIONAL(-1.0,70.0)*y;
+      t680 = y*t695;
+      t679 = y*t698;
+      t678 = RATIONAL(-1.0,20.0)*t725;
+      t677 = y*t686;
+      t676 = y*t694;
+      t675 = RATIONAL(1.0,20.0)*t725;
+      coeff_dxx_m2_m2 = t677+t694+t719;
+      coeff_dxx_m1_m2 = t680+t698+t715;
+      coeff_dxx_0_m2 = t688+t721;
+      coeff_dxx_p1_m2 = t679+t695+t715;
+      coeff_dxx_p2_m2 = t686+t676+t719;
+      coeff_dxx_m2_m1 = t685+t675+t713;
+      coeff_dxx_m1_m1 = t676+t684+t714;
+      coeff_dxx_0_m1 = t690+t722;
+      coeff_dxx_p1_m1 = t684+t677+t716;
+      coeff_dxx_p2_m1 = t685+t678+t717;
+      coeff_dxx_m2_0 = t694+t724;
+      coeff_dxx_m1_0 = t698+t723;
+      coeff_dxx_0_0 = t692+RATIONAL(-57.0,98.0)+t696;
+      coeff_dxx_p1_0 = t695+t723;
+      coeff_dxx_p2_0 = t686+t724;
+      coeff_dxx_m2_p1 = t690+t678+t713;
+      coeff_dxx_m1_p1 = t677+t683+t714;
+      coeff_dxx_0_p1 = t685+t722;
+      coeff_dxx_p1_p1 = t683+t676+t716;
+      coeff_dxx_p2_p1 = t690+t675+t717;
+      coeff_dxx_m2_p2 = t676+t694+t720;
+      coeff_dxx_m1_p2 = t698+t679+t718;
+      coeff_dxx_0_p2 = t691+t721;
+      coeff_dxx_p1_p2 = t695+t680+t718;
+      coeff_dxx_p2_p2 = t686+t677+t720;
diff --git a/src/GeneralizedPolynomial-Uniform/Lagrange/2d.coeffs/2d.cube.order4.smooth0/coeff-dxy.compute.c b/src/GeneralizedPolynomial-Uniform/Lagrange/2d.coeffs/2d.cube.order4.smooth0/coeff-dxy.compute.c
new file mode 100644
index 0000000..9388906
--- /dev/null
+++ b/src/GeneralizedPolynomial-Uniform/Lagrange/2d.coeffs/2d.cube.order4.smooth0/coeff-dxy.compute.c
@@ -0,0 +1,119 @@
+fp t778,
+   t777,
+   t796,
+   t766,
+   t735,
+   t772,
+   t755,
+   t795,
+   t751,
+   t748,
+   t794,
+   t793,
+   t792,
+   t791,
+   t733,
+   t790,
+   t789,
+   t788,
+   t787,
+   t765,
+   t770,
+   t786,
+   t734,
+   t785,
+   t767,
+   t771,
+   t784,
+   t783,
+   t782,
+   t730,
+   t781,
+   t780,
+   t779,
+   t776,
+   t775,
+   t773,
+   t769,
+   t756,
+   t754,
+   t750,
+   t749,
+   t746,
+   t745,
+   t743,
+   t738,
+   t732,
+   t731;
+      t778 = y*y;
+      t777 = x*x;
+      t796 = t778+t777;
+      t766 = RATIONAL(-1.0,35.0);
+      t735 = t766*x;
+      t772 = RATIONAL(1.0,35.0);
+      t755 = t772*y;
+      t795 = t735+t755;
+      t751 = t766*y;
+      t748 = t772*x;
+      t794 = t751+t748;
+      t793 = t735+t751;
+      t792 = t755+t748;
+      t791 = x*y;
+      t733 = RATIONAL(-2.0,49.0)*t791;
+      t790 = t733+RATIONAL(-21.0,200.0);
+      t789 = t733+RATIONAL(21.0,200.0);
+      t788 = t796*RATIONAL(-1.0,20.0);
+      t787 = t796*RATIONAL(1.0,20.0);
+      t765 = RATIONAL(-1.0,40.0);
+      t770 = RATIONAL(1.0,10.0);
+      t786 = t765*t777+t770*t778+t790;
+      t734 = RATIONAL(4.0,49.0)*t791;
+      t785 = t734+RATIONAL(11.0,150.0)+t788;
+      t767 = RATIONAL(1.0,40.0);
+      t771 = RATIONAL(-1.0,10.0);
+      t784 = t771*t778+t767*t777+t789;
+      t783 = t771*t777+t767*t778+t789;
+      t782 = t734+RATIONAL(-11.0,150.0)+t787;
+      t730 = RATIONAL(1.0,49.0)*t791;
+      t781 = RATIONAL(37.0,300.0)+t730+t788;
+      t780 = RATIONAL(-37.0,300.0)+t730+t787;
+      t779 = t765*t778+t770*t777+t790;
+      t776 = RATIONAL(-2.0,35.0);
+      t775 = RATIONAL(-1.0,70.0);
+      t773 = RATIONAL(1.0,70.0);
+      t769 = RATIONAL(2.0,35.0);
+      t756 = t776*y;
+      t754 = t769*y;
+      t750 = t776*x;
+      t749 = t775*x;
+      t746 = t773*y;
+      t745 = t775*y;
+      t743 = t769*x;
+      t738 = t773*x;
+      t732 = RATIONAL(-4.0,49.0)*t791;
+      t731 = RATIONAL(2.0,49.0)*t791;
+      coeff_dxy_m2_m2 = t756+t750+t782;
+      coeff_dxy_m1_m2 = t783+t794;
+      coeff_dxy_0_m2 = t743+t732;
+      coeff_dxy_p1_m2 = t779+t792;
+      coeff_dxy_p2_m2 = t750+t754+t785;
+      coeff_dxy_m2_m1 = t784+t795;
+      coeff_dxy_m1_m1 = t738+t746+t781;
+      coeff_dxy_0_m1 = t748+t731;
+      coeff_dxy_p1_m1 = t745+t738+t780;
+      coeff_dxy_p2_m1 = t786+t793;
+      coeff_dxy_m2_0 = t754+t732;
+      coeff_dxy_m1_0 = t755+t731;
+      coeff_dxy_0_0 = t734;
+      coeff_dxy_p1_0 = t751+t731;
+      coeff_dxy_p2_0 = t732+t756;
+      coeff_dxy_m2_p1 = t786+t792;
+      coeff_dxy_m1_p1 = t749+t746+t780;
+      coeff_dxy_0_p1 = t735+t731;
+      coeff_dxy_p1_p1 = t749+t745+t781;
+      coeff_dxy_p2_p1 = t784+t794;
+      coeff_dxy_m2_p2 = t743+t756+t785;
+      coeff_dxy_m1_p2 = t779+t793;
+      coeff_dxy_0_p2 = t750+t732;
+      coeff_dxy_p1_p2 = t783+t795;
+      coeff_dxy_p2_p2 = t754+t743+t782;
diff --git a/src/GeneralizedPolynomial-Uniform/Lagrange/2d.coeffs/2d.cube.order4.smooth0/coeff-dy.compute.c b/src/GeneralizedPolynomial-Uniform/Lagrange/2d.coeffs/2d.cube.order4.smooth0/coeff-dy.compute.c
new file mode 100644
index 0000000..03c68c3
--- /dev/null
+++ b/src/GeneralizedPolynomial-Uniform/Lagrange/2d.coeffs/2d.cube.order4.smooth0/coeff-dy.compute.c
@@ -0,0 +1,151 @@
+fp t641,
+   t611,
+   t638,
+   t639,
+   t657,
+   t674,
+   t608,
+   t656,
+   t673,
+   t661,
+   t589,
+   t672,
+   t671,
+   t635,
+   t640,
+   t609,
+   t670,
+   t574,
+   t669,
+   t614,
+   t598,
+   t668,
+   t632,
+   t615,
+   t667,
+   t634,
+   t619,
+   t666,
+   t617,
+   t596,
+   t665,
+   t627,
+   t606,
+   t664,
+   t637,
+   t616,
+   t663,
+   t662,
+   t660,
+   t618,
+   t659,
+   t658,
+   t586,
+   t655,
+   t578,
+   t654,
+   t620,
+   t653,
+   t652,
+   t651,
+   t650,
+   t649,
+   t648,
+   t647,
+   t646,
+   t645,
+   t644,
+   t588,
+   t583,
+   t582,
+   t579,
+   t576,
+   t572;
+      t641 = y*y;
+      t611 = RATIONAL(1.0,20.0)*t641;
+      t638 = x*x;
+      t639 = t638*x;
+      t657 = RATIONAL(1.0,60.0)*t639+x*t611;
+      t674 = t657+RATIONAL(-37.0,300.0)*x;
+      t608 = RATIONAL(-1.0,20.0)*t641;
+      t656 = RATIONAL(-1.0,60.0)*t639+x*t608;
+      t673 = t656+RATIONAL(37.0,300.0)*x;
+      t661 = y*t638;
+      t589 = RATIONAL(-1.0,49.0)*t661;
+      t672 = t589+RATIONAL(226.0,735.0)*y;
+      t671 = RATIONAL(1.0,98.0)*t661+RATIONAL(181.0,735.0)*y;
+      t635 = RATIONAL(1.0,30.0);
+      t640 = y*t641;
+      t609 = t635*t640;
+      t670 = t589+t609+RATIONAL(71.0,2940.0)*y;
+      t574 = RATIONAL(2.0,49.0)*t661;
+      t669 = RATIONAL(-289.0,2940.0)*y+t574+t609;
+      t614 = RATIONAL(-1.0,10.0)*t641;
+      t598 = RATIONAL(21.0,200.0)*x;
+      t668 = x*t614+RATIONAL(1.0,120.0)*t639+t598;
+      t632 = RATIONAL(1.0,35.0);
+      t615 = t632*t638;
+      t667 = t611+RATIONAL(-31.0,420.0)+t615;
+      t634 = RATIONAL(1.0,70.0);
+      t619 = t634*t638;
+      t666 = t619+RATIONAL(-1.0,84.0)+t608;
+      t617 = RATIONAL(1.0,10.0)*t641;
+      t596 = RATIONAL(-21.0,200.0)*x;
+      t665 = RATIONAL(-1.0,120.0)*t639+x*t617+t596;
+      t627 = RATIONAL(-1.0,35.0);
+      t606 = t627*t638;
+      t664 = RATIONAL(31.0,420.0)+t606+t608;
+      t637 = RATIONAL(-1.0,70.0);
+      t616 = t637*t638;
+      t663 = t611+t616+RATIONAL(1.0,84.0);
+      t662 = x*y;
+      t660 = x*t641;
+      t618 = RATIONAL(-2.0,15.0)*t640;
+      t659 = t617+t618;
+      t658 = t618+t614;
+      t586 = RATIONAL(-2.0,49.0)*t661;
+      t655 = t586+t609+RATIONAL(191.0,2940.0)*y;
+      t578 = RATIONAL(1.0,49.0)*t661;
+      t654 = t618+RATIONAL(166.0,735.0)*y+t578;
+      t620 = RATIONAL(1.0,5.0)*t640;
+      t653 = RATIONAL(-53.0,98.0)*y+t620+t578;
+      t652 = t620+RATIONAL(-41.0,98.0)*y+t586;
+      t651 = RATIONAL(1.0,140.0)*t638+RATIONAL(-31.0,210.0)+t659+t671;
+      t650 = t619+RATIONAL(11.0,105.0)+t658+t672;
+      t649 = RATIONAL(-11.0,150.0)*x+t657+t669;
+      t648 = RATIONAL(11.0,150.0)*x+t656+t669;
+      t647 = t616+RATIONAL(-11.0,105.0)+t659+t672;
+      t646 = t635*t639+RATIONAL(-1.0,40.0)*t660+t596+t670;
+      t645 = RATIONAL(31.0,210.0)+RATIONAL(-1.0,140.0)*t638+t658+t671;
+      t644 = RATIONAL(-1.0,30.0)*t639+RATIONAL(1.0,40.0)*t660+t598+t670;
+      t588 = t627*t662;
+      t583 = t632*t662;
+      t582 = RATIONAL(2.0,35.0)*t662;
+      t579 = t634*t662;
+      t576 = t637*t662;
+      t572 = RATIONAL(-2.0,35.0)*t662;
+      coeff_dy_m2_m2 = t572+t649+t664;
+      coeff_dy_m1_m2 = t588+t644+t666;
+      coeff_dy_0_m2 = RATIONAL(-17.0,420.0)+t615+t608+t655;
+      coeff_dy_p1_m2 = t583+t646+t666;
+      coeff_dy_p2_m2 = t582+t648+t664;
+      coeff_dy_m2_m1 = t583+t647+t668;
+      coeff_dy_m1_m1 = t579+t651+t673;
+      coeff_dy_0_m1 = RATIONAL(-17.0,105.0)+t617+t619+t654;
+      coeff_dy_p1_m1 = t576+t651+t674;
+      coeff_dy_p2_m1 = t588+t647+t665;
+      coeff_dy_m2_0 = t582+t652;
+      coeff_dy_m1_0 = t583+t653;
+      coeff_dy_0_0 = t574+RATIONAL(-57.0,98.0)*y+t620;
+      coeff_dy_p1_0 = t588+t653;
+      coeff_dy_p2_0 = t572+t652;
+      coeff_dy_m2_p1 = t583+t650+t665;
+      coeff_dy_m1_p1 = t579+t645+t674;
+      coeff_dy_0_p1 = RATIONAL(17.0,105.0)+t616+t614+t654;
+      coeff_dy_p1_p1 = t576+t645+t673;
+      coeff_dy_p2_p1 = t588+t650+t668;
+      coeff_dy_m2_p2 = t572+t648+t667;
+      coeff_dy_m1_p2 = t588+t646+t663;
+      coeff_dy_0_p2 = t611+RATIONAL(17.0,420.0)+t606+t655;
+      coeff_dy_p1_p2 = t583+t644+t663;
+      coeff_dy_p2_p2 = t582+t649+t667;
diff --git a/src/GeneralizedPolynomial-Uniform/Lagrange/2d.coeffs/2d.cube.order4.smooth0/coeff-dyy.compute.c b/src/GeneralizedPolynomial-Uniform/Lagrange/2d.coeffs/2d.cube.order4.smooth0/coeff-dyy.compute.c
new file mode 100644
index 0000000..cc0ccef
--- /dev/null
+++ b/src/GeneralizedPolynomial-Uniform/Lagrange/2d.coeffs/2d.cube.order4.smooth0/coeff-dyy.compute.c
@@ -0,0 +1,115 @@
+fp t833,
+   t808,
+   t853,
+   t807,
+   t852,
+   t834,
+   t818,
+   t851,
+   t850,
+   t849,
+   t829,
+   t815,
+   t817,
+   t806,
+   t820,
+   t804,
+   t848,
+   t847,
+   t812,
+   t846,
+   t819,
+   t845,
+   t803,
+   t844,
+   t843,
+   t798,
+   t842,
+   t841,
+   t840,
+   t839,
+   t838,
+   t837,
+   t800,
+   t836,
+   t835,
+   t816,
+   t814,
+   t811,
+   t810,
+   t809,
+   t805,
+   t802,
+   t801,
+   t799,
+   t797;
+      t833 = x*x;
+      t808 = RATIONAL(-1.0,49.0)*t833;
+      t853 = t808+RATIONAL(71.0,2940.0);
+      t807 = RATIONAL(2.0,49.0)*t833;
+      t852 = t807+RATIONAL(-289.0,2940.0);
+      t834 = y*y;
+      t818 = RATIONAL(-2.0,5.0)*t834;
+      t851 = t818+RATIONAL(181.0,735.0)+RATIONAL(1.0,98.0)*t833;
+      t850 = t818+t808+RATIONAL(226.0,735.0);
+      t849 = x*y;
+      t829 = RATIONAL(1.0,10.0);
+      t815 = t829*y;
+      t817 = RATIONAL(-1.0,5.0)*y;
+      t806 = RATIONAL(-1.0,10.0)*y;
+      t820 = RATIONAL(1.0,5.0)*y;
+      t804 = t829*t834;
+      t848 = t815+t804;
+      t847 = t806+t804;
+      t812 = RATIONAL(1.0,49.0)*t833;
+      t846 = t818+RATIONAL(166.0,735.0)+t812;
+      t819 = RATIONAL(3.0,5.0)*t834;
+      t845 = t819+RATIONAL(-53.0,98.0)+t812;
+      t803 = RATIONAL(-2.0,49.0)*t833;
+      t844 = t803+RATIONAL(191.0,2940.0)+t804;
+      t843 = t803+t819+RATIONAL(-41.0,98.0);
+      t798 = x*t815;
+      t842 = t798+t851;
+      t841 = t848+t853;
+      t840 = t847+t853;
+      t839 = t847+t852;
+      t838 = t820+t850;
+      t837 = t817+t850;
+      t800 = x*t806;
+      t836 = t800+t851;
+      t835 = t848+t852;
+      t816 = RATIONAL(2.0,35.0)*x;
+      t814 = RATIONAL(1.0,70.0)*x;
+      t811 = RATIONAL(-1.0,70.0)*x;
+      t810 = RATIONAL(1.0,35.0)*x;
+      t809 = RATIONAL(-1.0,35.0)*x;
+      t805 = RATIONAL(-2.0,35.0)*x;
+      t802 = x*t817;
+      t801 = x*t820;
+      t799 = RATIONAL(-1.0,20.0)*t849;
+      t797 = RATIONAL(1.0,20.0)*t849;
+      coeff_dyy_m2_m2 = t798+t805+t839;
+      coeff_dyy_m1_m2 = t797+t809+t840;
+      coeff_dyy_0_m2 = t806+t844;
+      coeff_dyy_p1_m2 = t810+t799+t840;
+      coeff_dyy_p2_m2 = t800+t816+t839;
+      coeff_dyy_m2_m1 = t810+t802+t838;
+      coeff_dyy_m1_m1 = t820+t814+t836;
+      coeff_dyy_0_m1 = t820+t846;
+      coeff_dyy_p1_m1 = t811+t820+t842;
+      coeff_dyy_p2_m1 = t809+t801+t838;
+      coeff_dyy_m2_0 = t816+t843;
+      coeff_dyy_m1_0 = t810+t845;
+      coeff_dyy_0_0 = t819+RATIONAL(-57.0,98.0)+t807;
+      coeff_dyy_p1_0 = t809+t845;
+      coeff_dyy_p2_0 = t805+t843;
+      coeff_dyy_m2_p1 = t810+t801+t837;
+      coeff_dyy_m1_p1 = t817+t814+t842;
+      coeff_dyy_0_p1 = t817+t846;
+      coeff_dyy_p1_p1 = t811+t817+t836;
+      coeff_dyy_p2_p1 = t809+t802+t837;
+      coeff_dyy_m2_p2 = t805+t800+t835;
+      coeff_dyy_m1_p2 = t809+t799+t841;
+      coeff_dyy_0_p2 = t815+t844;
+      coeff_dyy_p1_p2 = t797+t810+t841;
+      coeff_dyy_p2_p2 = t816+t798+t835;