From 80df27de110b6e739fda09ba4d9eae9fdbefc342 Mon Sep 17 00:00:00 2001 From: jthorn Date: Mon, 2 Sep 2002 10:36:32 +0000 Subject: This is a major reworking of this interpolator. The interpolation code in template.c is now split up into subfunctions (in common/, Lagrange/, and Hermite/), so it should now compile in finite cpu/memory even with optimization. git-svn-id: http://svn.cactuscode.org/arrangements/CactusBase/LocalInterp/trunk@115 df1f8a13-aa1d-4dd4-9681-27ded5b42416 --- .../2d.cube.order1.smooth0/coeff-I.compute.c | 18 -- .../2d.cube.order1.smooth0/coeff-dx.compute.c | 8 - .../2d.cube.order1.smooth0/coeff-dy.compute.c | 8 - .../2d.cube.order1.smooth0/coeffs-I.compute.c | 18 ++ .../2d.cube.order1.smooth0/coeffs-dx.compute.c | 8 + .../2d.cube.order1.smooth0/coeffs-dy.compute.c | 8 + .../2d.cube.order2.smooth0/coeff-I.compute.c | 53 ----- .../2d.cube.order2.smooth0/coeff-dx.compute.c | 21 -- .../2d.cube.order2.smooth0/coeff-dxx.compute.c | 13 -- .../2d.cube.order2.smooth0/coeff-dxy.compute.c | 15 -- .../2d.cube.order2.smooth0/coeff-dy.compute.c | 21 -- .../2d.cube.order2.smooth0/coeff-dyy.compute.c | 13 -- .../2d.cube.order2.smooth0/coeffs-I.compute.c | 53 +++++ .../2d.cube.order2.smooth0/coeffs-dx.compute.c | 21 ++ .../2d.cube.order2.smooth0/coeffs-dxx.compute.c | 13 ++ .../2d.cube.order2.smooth0/coeffs-dxy.compute.c | 15 ++ .../2d.cube.order2.smooth0/coeffs-dy.compute.c | 21 ++ .../2d.cube.order2.smooth0/coeffs-dyy.compute.c | 13 ++ .../2d.cube.order3.smooth0/coeff-I.compute.c | 196 ----------------- .../2d.cube.order3.smooth0/coeff-dx.compute.c | 95 -------- .../2d.cube.order3.smooth0/coeff-dxx.compute.c | 40 ---- .../2d.cube.order3.smooth0/coeff-dxy.compute.c | 52 ----- .../2d.cube.order3.smooth0/coeff-dy.compute.c | 95 -------- .../2d.cube.order3.smooth0/coeff-dyy.compute.c | 40 ---- .../2d.cube.order3.smooth0/coeffs-I.compute.c | 190 ++++++++++++++++ .../2d.cube.order3.smooth0/coeffs-dx.compute.c | 97 ++++++++ .../2d.cube.order3.smooth0/coeffs-dxx.compute.c | 40 ++++ .../2d.cube.order3.smooth0/coeffs-dxy.compute.c | 52 +++++ .../2d.cube.order3.smooth0/coeffs-dy.compute.c | 97 ++++++++ .../2d.cube.order3.smooth0/coeffs-dyy.compute.c | 40 ++++ .../2d.cube.order4.smooth0/coeff-I.compute.c | 243 --------------------- .../2d.cube.order4.smooth0/coeff-dx.compute.c | 155 ------------- .../2d.cube.order4.smooth0/coeff-dxx.compute.c | 107 --------- .../2d.cube.order4.smooth0/coeff-dxy.compute.c | 119 ---------- .../2d.cube.order4.smooth0/coeff-dy.compute.c | 147 ------------- .../2d.cube.order4.smooth0/coeff-dyy.compute.c | 115 ---------- .../2d.cube.order4.smooth0/coeffs-I.compute.c | 239 ++++++++++++++++++++ .../2d.cube.order4.smooth0/coeffs-dx.compute.c | 155 +++++++++++++ .../2d.cube.order4.smooth0/coeffs-dxx.compute.c | 107 +++++++++ .../2d.cube.order4.smooth0/coeffs-dxy.compute.c | 119 ++++++++++ .../2d.cube.order4.smooth0/coeffs-dy.compute.c | 147 +++++++++++++ .../2d.cube.order4.smooth0/coeffs-dyy.compute.c | 107 +++++++++ 42 files changed, 1560 insertions(+), 1574 deletions(-) delete mode 100644 src/GeneralizedPolynomial-Uniform/Lagrange/2d.coeffs/2d.cube.order1.smooth0/coeff-I.compute.c delete mode 100644 src/GeneralizedPolynomial-Uniform/Lagrange/2d.coeffs/2d.cube.order1.smooth0/coeff-dx.compute.c delete mode 100644 src/GeneralizedPolynomial-Uniform/Lagrange/2d.coeffs/2d.cube.order1.smooth0/coeff-dy.compute.c create mode 100644 src/GeneralizedPolynomial-Uniform/Lagrange/2d.coeffs/2d.cube.order1.smooth0/coeffs-I.compute.c create mode 100644 src/GeneralizedPolynomial-Uniform/Lagrange/2d.coeffs/2d.cube.order1.smooth0/coeffs-dx.compute.c create mode 100644 src/GeneralizedPolynomial-Uniform/Lagrange/2d.coeffs/2d.cube.order1.smooth0/coeffs-dy.compute.c delete mode 100644 src/GeneralizedPolynomial-Uniform/Lagrange/2d.coeffs/2d.cube.order2.smooth0/coeff-I.compute.c delete mode 100644 src/GeneralizedPolynomial-Uniform/Lagrange/2d.coeffs/2d.cube.order2.smooth0/coeff-dx.compute.c delete mode 100644 src/GeneralizedPolynomial-Uniform/Lagrange/2d.coeffs/2d.cube.order2.smooth0/coeff-dxx.compute.c delete mode 100644 src/GeneralizedPolynomial-Uniform/Lagrange/2d.coeffs/2d.cube.order2.smooth0/coeff-dxy.compute.c delete mode 100644 src/GeneralizedPolynomial-Uniform/Lagrange/2d.coeffs/2d.cube.order2.smooth0/coeff-dy.compute.c delete mode 100644 src/GeneralizedPolynomial-Uniform/Lagrange/2d.coeffs/2d.cube.order2.smooth0/coeff-dyy.compute.c create mode 100644 src/GeneralizedPolynomial-Uniform/Lagrange/2d.coeffs/2d.cube.order2.smooth0/coeffs-I.compute.c create mode 100644 src/GeneralizedPolynomial-Uniform/Lagrange/2d.coeffs/2d.cube.order2.smooth0/coeffs-dx.compute.c create mode 100644 src/GeneralizedPolynomial-Uniform/Lagrange/2d.coeffs/2d.cube.order2.smooth0/coeffs-dxx.compute.c create mode 100644 src/GeneralizedPolynomial-Uniform/Lagrange/2d.coeffs/2d.cube.order2.smooth0/coeffs-dxy.compute.c create mode 100644 src/GeneralizedPolynomial-Uniform/Lagrange/2d.coeffs/2d.cube.order2.smooth0/coeffs-dy.compute.c create mode 100644 src/GeneralizedPolynomial-Uniform/Lagrange/2d.coeffs/2d.cube.order2.smooth0/coeffs-dyy.compute.c delete mode 100644 src/GeneralizedPolynomial-Uniform/Lagrange/2d.coeffs/2d.cube.order3.smooth0/coeff-I.compute.c delete mode 100644 src/GeneralizedPolynomial-Uniform/Lagrange/2d.coeffs/2d.cube.order3.smooth0/coeff-dx.compute.c delete mode 100644 src/GeneralizedPolynomial-Uniform/Lagrange/2d.coeffs/2d.cube.order3.smooth0/coeff-dxx.compute.c delete mode 100644 src/GeneralizedPolynomial-Uniform/Lagrange/2d.coeffs/2d.cube.order3.smooth0/coeff-dxy.compute.c delete mode 100644 src/GeneralizedPolynomial-Uniform/Lagrange/2d.coeffs/2d.cube.order3.smooth0/coeff-dy.compute.c delete mode 100644 src/GeneralizedPolynomial-Uniform/Lagrange/2d.coeffs/2d.cube.order3.smooth0/coeff-dyy.compute.c create mode 100644 src/GeneralizedPolynomial-Uniform/Lagrange/2d.coeffs/2d.cube.order3.smooth0/coeffs-I.compute.c create mode 100644 src/GeneralizedPolynomial-Uniform/Lagrange/2d.coeffs/2d.cube.order3.smooth0/coeffs-dx.compute.c create mode 100644 src/GeneralizedPolynomial-Uniform/Lagrange/2d.coeffs/2d.cube.order3.smooth0/coeffs-dxx.compute.c create mode 100644 src/GeneralizedPolynomial-Uniform/Lagrange/2d.coeffs/2d.cube.order3.smooth0/coeffs-dxy.compute.c create mode 100644 src/GeneralizedPolynomial-Uniform/Lagrange/2d.coeffs/2d.cube.order3.smooth0/coeffs-dy.compute.c create mode 100644 src/GeneralizedPolynomial-Uniform/Lagrange/2d.coeffs/2d.cube.order3.smooth0/coeffs-dyy.compute.c delete mode 100644 src/GeneralizedPolynomial-Uniform/Lagrange/2d.coeffs/2d.cube.order4.smooth0/coeff-I.compute.c delete mode 100644 src/GeneralizedPolynomial-Uniform/Lagrange/2d.coeffs/2d.cube.order4.smooth0/coeff-dx.compute.c delete mode 100644 src/GeneralizedPolynomial-Uniform/Lagrange/2d.coeffs/2d.cube.order4.smooth0/coeff-dxx.compute.c delete mode 100644 src/GeneralizedPolynomial-Uniform/Lagrange/2d.coeffs/2d.cube.order4.smooth0/coeff-dxy.compute.c delete mode 100644 src/GeneralizedPolynomial-Uniform/Lagrange/2d.coeffs/2d.cube.order4.smooth0/coeff-dy.compute.c delete mode 100644 src/GeneralizedPolynomial-Uniform/Lagrange/2d.coeffs/2d.cube.order4.smooth0/coeff-dyy.compute.c create mode 100644 src/GeneralizedPolynomial-Uniform/Lagrange/2d.coeffs/2d.cube.order4.smooth0/coeffs-I.compute.c create mode 100644 src/GeneralizedPolynomial-Uniform/Lagrange/2d.coeffs/2d.cube.order4.smooth0/coeffs-dx.compute.c create mode 100644 src/GeneralizedPolynomial-Uniform/Lagrange/2d.coeffs/2d.cube.order4.smooth0/coeffs-dxx.compute.c create mode 100644 src/GeneralizedPolynomial-Uniform/Lagrange/2d.coeffs/2d.cube.order4.smooth0/coeffs-dxy.compute.c create mode 100644 src/GeneralizedPolynomial-Uniform/Lagrange/2d.coeffs/2d.cube.order4.smooth0/coeffs-dy.compute.c create mode 100644 src/GeneralizedPolynomial-Uniform/Lagrange/2d.coeffs/2d.cube.order4.smooth0/coeffs-dyy.compute.c (limited to 'src/GeneralizedPolynomial-Uniform/Lagrange/2d.coeffs') diff --git a/src/GeneralizedPolynomial-Uniform/Lagrange/2d.coeffs/2d.cube.order1.smooth0/coeff-I.compute.c b/src/GeneralizedPolynomial-Uniform/Lagrange/2d.coeffs/2d.cube.order1.smooth0/coeff-I.compute.c deleted file mode 100644 index 221b33b..0000000 --- a/src/GeneralizedPolynomial-Uniform/Lagrange/2d.coeffs/2d.cube.order1.smooth0/coeff-I.compute.c +++ /dev/null @@ -1,18 +0,0 @@ -fp t7, - t6, - t5, - t4, - t3, - t2, - t1; - t7 = RATIONAL(-1.0,2.0); - t6 = RATIONAL(1.0,4.0); - t5 = RATIONAL(1.0,2.0); - t4 = t7*y; - t3 = t7*x; - t2 = t5*x; - t1 = t5*y; - coeff_I_0_0 = RATIONAL(3.0,4.0)+t4+t3; - coeff_I_p1_0 = t6+t4+t2; - coeff_I_0_p1 = t3+t6+t1; - coeff_I_p1_p1 = t1+RATIONAL(-1.0,4.0)+t2; diff --git a/src/GeneralizedPolynomial-Uniform/Lagrange/2d.coeffs/2d.cube.order1.smooth0/coeff-dx.compute.c b/src/GeneralizedPolynomial-Uniform/Lagrange/2d.coeffs/2d.cube.order1.smooth0/coeff-dx.compute.c deleted file mode 100644 index 3d89f01..0000000 --- a/src/GeneralizedPolynomial-Uniform/Lagrange/2d.coeffs/2d.cube.order1.smooth0/coeff-dx.compute.c +++ /dev/null @@ -1,8 +0,0 @@ -fp t9, - t8; - t9 = RATIONAL(-1.0,2.0); - t8 = RATIONAL(1.0,2.0); - coeff_dx_0_0 = t9; - coeff_dx_p1_0 = t8; - coeff_dx_0_p1 = t9; - coeff_dx_p1_p1 = t8; diff --git a/src/GeneralizedPolynomial-Uniform/Lagrange/2d.coeffs/2d.cube.order1.smooth0/coeff-dy.compute.c b/src/GeneralizedPolynomial-Uniform/Lagrange/2d.coeffs/2d.cube.order1.smooth0/coeff-dy.compute.c deleted file mode 100644 index 42819a8..0000000 --- a/src/GeneralizedPolynomial-Uniform/Lagrange/2d.coeffs/2d.cube.order1.smooth0/coeff-dy.compute.c +++ /dev/null @@ -1,8 +0,0 @@ -fp t11, - t10; - t11 = RATIONAL(-1.0,2.0); - t10 = RATIONAL(1.0,2.0); - coeff_dy_0_0 = t11; - coeff_dy_p1_0 = t11; - coeff_dy_0_p1 = t10; - coeff_dy_p1_p1 = t10; diff --git a/src/GeneralizedPolynomial-Uniform/Lagrange/2d.coeffs/2d.cube.order1.smooth0/coeffs-I.compute.c b/src/GeneralizedPolynomial-Uniform/Lagrange/2d.coeffs/2d.cube.order1.smooth0/coeffs-I.compute.c new file mode 100644 index 0000000..cda1593 --- /dev/null +++ b/src/GeneralizedPolynomial-Uniform/Lagrange/2d.coeffs/2d.cube.order1.smooth0/coeffs-I.compute.c @@ -0,0 +1,18 @@ +fp t7; +fp t6; +fp t5; +fp t4; +fp t3; +fp t2; +fp t1; + t7 = RATIONAL(-1.0,2.0); + t6 = RATIONAL(1.0,4.0); + t5 = RATIONAL(1.0,2.0); + t4 = t7*y; + t3 = t7*x; + t2 = t5*x; + t1 = t5*y; + coeffs_I->coeff_0_0 = RATIONAL(3.0,4.0)+t4+t3; + coeffs_I->coeff_p1_0 = t6+t4+t2; + coeffs_I->coeff_0_p1 = t3+t6+t1; + coeffs_I->coeff_p1_p1 = t1+RATIONAL(-1.0,4.0)+t2; diff --git a/src/GeneralizedPolynomial-Uniform/Lagrange/2d.coeffs/2d.cube.order1.smooth0/coeffs-dx.compute.c b/src/GeneralizedPolynomial-Uniform/Lagrange/2d.coeffs/2d.cube.order1.smooth0/coeffs-dx.compute.c new file mode 100644 index 0000000..cf806b2 --- /dev/null +++ b/src/GeneralizedPolynomial-Uniform/Lagrange/2d.coeffs/2d.cube.order1.smooth0/coeffs-dx.compute.c @@ -0,0 +1,8 @@ +fp t9; +fp t8; + t9 = RATIONAL(-1.0,2.0); + t8 = RATIONAL(1.0,2.0); + coeffs_dx->coeff_0_0 = t9; + coeffs_dx->coeff_p1_0 = t8; + coeffs_dx->coeff_0_p1 = t9; + coeffs_dx->coeff_p1_p1 = t8; diff --git a/src/GeneralizedPolynomial-Uniform/Lagrange/2d.coeffs/2d.cube.order1.smooth0/coeffs-dy.compute.c b/src/GeneralizedPolynomial-Uniform/Lagrange/2d.coeffs/2d.cube.order1.smooth0/coeffs-dy.compute.c new file mode 100644 index 0000000..e5bb3a5 --- /dev/null +++ b/src/GeneralizedPolynomial-Uniform/Lagrange/2d.coeffs/2d.cube.order1.smooth0/coeffs-dy.compute.c @@ -0,0 +1,8 @@ +fp t11; +fp t10; + t11 = RATIONAL(-1.0,2.0); + t10 = RATIONAL(1.0,2.0); + coeffs_dy->coeff_0_0 = t11; + coeffs_dy->coeff_p1_0 = t11; + coeffs_dy->coeff_0_p1 = t10; + coeffs_dy->coeff_p1_p1 = t10; diff --git a/src/GeneralizedPolynomial-Uniform/Lagrange/2d.coeffs/2d.cube.order2.smooth0/coeff-I.compute.c b/src/GeneralizedPolynomial-Uniform/Lagrange/2d.coeffs/2d.cube.order2.smooth0/coeff-I.compute.c deleted file mode 100644 index 52b6d48..0000000 --- a/src/GeneralizedPolynomial-Uniform/Lagrange/2d.coeffs/2d.cube.order2.smooth0/coeff-I.compute.c +++ /dev/null @@ -1,53 +0,0 @@ -fp t27, - t26, - t18, - t28, - t21, - t34, - t33, - t23, - t15, - t22, - t32, - t17, - t31, - t14, - t30, - t25, - t20, - t29, - t19, - t16, - t13, - t12; - t27 = x*x; - t26 = RATIONAL(1.0,6.0); - t18 = t26*t27; - t28 = y*y; - t21 = t26*t28; - t34 = t18+t21+RATIONAL(-1.0,9.0); - t33 = x*y; - t23 = RATIONAL(-1.0,3.0); - t15 = t23*t28; - t22 = RATIONAL(2.0,9.0); - t32 = t15+t18+t22; - t17 = t23*t27; - t31 = t17+t21+t22; - t14 = t26*y; - t30 = t14+t34; - t25 = RATIONAL(-1.0,6.0); - t20 = t25*y; - t29 = t20+t34; - t19 = t25*x; - t16 = t26*x; - t13 = RATIONAL(1.0,4.0)*t33; - t12 = RATIONAL(-1.0,4.0)*t33; - coeff_I_m1_m1 = t13+t19+t29; - coeff_I_0_m1 = t20+t31; - coeff_I_p1_m1 = t12+t16+t29; - coeff_I_m1_0 = t19+t32; - coeff_I_0_0 = t15+RATIONAL(5.0,9.0)+t17; - coeff_I_p1_0 = t16+t32; - coeff_I_m1_p1 = t19+t12+t30; - coeff_I_0_p1 = t14+t31; - coeff_I_p1_p1 = t16+t13+t30; diff --git a/src/GeneralizedPolynomial-Uniform/Lagrange/2d.coeffs/2d.cube.order2.smooth0/coeff-dx.compute.c b/src/GeneralizedPolynomial-Uniform/Lagrange/2d.coeffs/2d.cube.order2.smooth0/coeff-dx.compute.c deleted file mode 100644 index 14b1a1c..0000000 --- a/src/GeneralizedPolynomial-Uniform/Lagrange/2d.coeffs/2d.cube.order2.smooth0/coeff-dx.compute.c +++ /dev/null @@ -1,21 +0,0 @@ -fp t37, - t42, - t41, - t40, - t39, - t36; - t37 = RATIONAL(1.0,3.0)*x; - t42 = RATIONAL(1.0,4.0)*y+t37; - t41 = t37+RATIONAL(-1.0,4.0)*y; - t40 = RATIONAL(1.0,6.0); - t39 = RATIONAL(-1.0,6.0); - t36 = RATIONAL(-2.0,3.0)*x; - coeff_dx_m1_m1 = t39+t42; - coeff_dx_0_m1 = t36; - coeff_dx_p1_m1 = t40+t41; - coeff_dx_m1_0 = t39+t37; - coeff_dx_0_0 = t36; - coeff_dx_p1_0 = t40+t37; - coeff_dx_m1_p1 = t39+t41; - coeff_dx_0_p1 = t36; - coeff_dx_p1_p1 = t40+t42; diff --git a/src/GeneralizedPolynomial-Uniform/Lagrange/2d.coeffs/2d.cube.order2.smooth0/coeff-dxx.compute.c b/src/GeneralizedPolynomial-Uniform/Lagrange/2d.coeffs/2d.cube.order2.smooth0/coeff-dxx.compute.c deleted file mode 100644 index 3381569..0000000 --- a/src/GeneralizedPolynomial-Uniform/Lagrange/2d.coeffs/2d.cube.order2.smooth0/coeff-dxx.compute.c +++ /dev/null @@ -1,13 +0,0 @@ -fp t52, - t51; - t52 = RATIONAL(-2.0,3.0); - t51 = RATIONAL(1.0,3.0); - coeff_dxx_m1_m1 = t51; - coeff_dxx_0_m1 = t52; - coeff_dxx_p1_m1 = t51; - coeff_dxx_m1_0 = t51; - coeff_dxx_0_0 = t52; - coeff_dxx_p1_0 = t51; - coeff_dxx_m1_p1 = t51; - coeff_dxx_0_p1 = t52; - coeff_dxx_p1_p1 = t51; diff --git a/src/GeneralizedPolynomial-Uniform/Lagrange/2d.coeffs/2d.cube.order2.smooth0/coeff-dxy.compute.c b/src/GeneralizedPolynomial-Uniform/Lagrange/2d.coeffs/2d.cube.order2.smooth0/coeff-dxy.compute.c deleted file mode 100644 index c795e1d..0000000 --- a/src/GeneralizedPolynomial-Uniform/Lagrange/2d.coeffs/2d.cube.order2.smooth0/coeff-dxy.compute.c +++ /dev/null @@ -1,15 +0,0 @@ -fp t55, - t54, - t53; - t55 = RATIONAL(0.0,1.0); - t54 = RATIONAL(1.0,4.0); - t53 = RATIONAL(-1.0,4.0); - coeff_dxy_m1_m1 = t54; - coeff_dxy_0_m1 = t55; - coeff_dxy_p1_m1 = t53; - coeff_dxy_m1_0 = t55; - coeff_dxy_0_0 = t55; - coeff_dxy_p1_0 = t55; - coeff_dxy_m1_p1 = t53; - coeff_dxy_0_p1 = t55; - coeff_dxy_p1_p1 = t54; diff --git a/src/GeneralizedPolynomial-Uniform/Lagrange/2d.coeffs/2d.cube.order2.smooth0/coeff-dy.compute.c b/src/GeneralizedPolynomial-Uniform/Lagrange/2d.coeffs/2d.cube.order2.smooth0/coeff-dy.compute.c deleted file mode 100644 index 65e74bd..0000000 --- a/src/GeneralizedPolynomial-Uniform/Lagrange/2d.coeffs/2d.cube.order2.smooth0/coeff-dy.compute.c +++ /dev/null @@ -1,21 +0,0 @@ -fp t46, - t50, - t49, - t48, - t47, - t43; - t46 = RATIONAL(1.0,3.0)*y; - t50 = RATIONAL(1.0,4.0)*x+t46; - t49 = RATIONAL(-1.0,4.0)*x+t46; - t48 = RATIONAL(1.0,6.0); - t47 = RATIONAL(-1.0,6.0); - t43 = RATIONAL(-2.0,3.0)*y; - coeff_dy_m1_m1 = t47+t50; - coeff_dy_0_m1 = t47+t46; - coeff_dy_p1_m1 = t47+t49; - coeff_dy_m1_0 = t43; - coeff_dy_0_0 = t43; - coeff_dy_p1_0 = t43; - coeff_dy_m1_p1 = t48+t49; - coeff_dy_0_p1 = t46+t48; - coeff_dy_p1_p1 = t48+t50; diff --git a/src/GeneralizedPolynomial-Uniform/Lagrange/2d.coeffs/2d.cube.order2.smooth0/coeff-dyy.compute.c b/src/GeneralizedPolynomial-Uniform/Lagrange/2d.coeffs/2d.cube.order2.smooth0/coeff-dyy.compute.c deleted file mode 100644 index 5d0b896..0000000 --- a/src/GeneralizedPolynomial-Uniform/Lagrange/2d.coeffs/2d.cube.order2.smooth0/coeff-dyy.compute.c +++ /dev/null @@ -1,13 +0,0 @@ -fp t57, - t56; - t57 = RATIONAL(-2.0,3.0); - t56 = RATIONAL(1.0,3.0); - coeff_dyy_m1_m1 = t56; - coeff_dyy_0_m1 = t56; - coeff_dyy_p1_m1 = t56; - coeff_dyy_m1_0 = t57; - coeff_dyy_0_0 = t57; - coeff_dyy_p1_0 = t57; - coeff_dyy_m1_p1 = t56; - coeff_dyy_0_p1 = t56; - coeff_dyy_p1_p1 = t56; diff --git a/src/GeneralizedPolynomial-Uniform/Lagrange/2d.coeffs/2d.cube.order2.smooth0/coeffs-I.compute.c b/src/GeneralizedPolynomial-Uniform/Lagrange/2d.coeffs/2d.cube.order2.smooth0/coeffs-I.compute.c new file mode 100644 index 0000000..c5e53d7 --- /dev/null +++ b/src/GeneralizedPolynomial-Uniform/Lagrange/2d.coeffs/2d.cube.order2.smooth0/coeffs-I.compute.c @@ -0,0 +1,53 @@ +fp t27; +fp t25; +fp t18; +fp t28; +fp t21; +fp t34; +fp t33; +fp t22; +fp t17; +fp t23; +fp t32; +fp t15; +fp t31; +fp t14; +fp t30; +fp t24; +fp t20; +fp t29; +fp t19; +fp t16; +fp t13; +fp t12; + t27 = x*x; + t25 = RATIONAL(1.0,6.0); + t18 = t25*t27; + t28 = y*y; + t21 = t25*t28; + t34 = t18+t21+RATIONAL(-1.0,9.0); + t33 = x*y; + t22 = RATIONAL(-1.0,3.0); + t17 = t22*t27; + t23 = RATIONAL(2.0,9.0); + t32 = t17+t21+t23; + t15 = t22*t28; + t31 = t15+t18+t23; + t14 = t25*y; + t30 = t14+t34; + t24 = RATIONAL(-1.0,6.0); + t20 = t24*y; + t29 = t20+t34; + t19 = t24*x; + t16 = t25*x; + t13 = RATIONAL(-1.0,4.0)*t33; + t12 = RATIONAL(1.0,4.0)*t33; + coeffs_I->coeff_m1_m1 = t19+t12+t29; + coeffs_I->coeff_0_m1 = t20+t32; + coeffs_I->coeff_p1_m1 = t13+t16+t29; + coeffs_I->coeff_m1_0 = t19+t31; + coeffs_I->coeff_0_0 = t15+t17+RATIONAL(5.0,9.0); + coeffs_I->coeff_p1_0 = t16+t31; + coeffs_I->coeff_m1_p1 = t19+t13+t30; + coeffs_I->coeff_0_p1 = t14+t32; + coeffs_I->coeff_p1_p1 = t16+t12+t30; diff --git a/src/GeneralizedPolynomial-Uniform/Lagrange/2d.coeffs/2d.cube.order2.smooth0/coeffs-dx.compute.c b/src/GeneralizedPolynomial-Uniform/Lagrange/2d.coeffs/2d.cube.order2.smooth0/coeffs-dx.compute.c new file mode 100644 index 0000000..7c930d7 --- /dev/null +++ b/src/GeneralizedPolynomial-Uniform/Lagrange/2d.coeffs/2d.cube.order2.smooth0/coeffs-dx.compute.c @@ -0,0 +1,21 @@ +fp t36; +fp t42; +fp t41; +fp t40; +fp t39; +fp t37; + t36 = RATIONAL(1.0,3.0)*x; + t42 = RATIONAL(1.0,4.0)*y+t36; + t41 = t36+RATIONAL(-1.0,4.0)*y; + t40 = RATIONAL(1.0,6.0); + t39 = RATIONAL(-1.0,6.0); + t37 = RATIONAL(-2.0,3.0)*x; + coeffs_dx->coeff_m1_m1 = t39+t42; + coeffs_dx->coeff_0_m1 = t37; + coeffs_dx->coeff_p1_m1 = t40+t41; + coeffs_dx->coeff_m1_0 = t36+t39; + coeffs_dx->coeff_0_0 = t37; + coeffs_dx->coeff_p1_0 = t36+t40; + coeffs_dx->coeff_m1_p1 = t39+t41; + coeffs_dx->coeff_0_p1 = t37; + coeffs_dx->coeff_p1_p1 = t40+t42; diff --git a/src/GeneralizedPolynomial-Uniform/Lagrange/2d.coeffs/2d.cube.order2.smooth0/coeffs-dxx.compute.c b/src/GeneralizedPolynomial-Uniform/Lagrange/2d.coeffs/2d.cube.order2.smooth0/coeffs-dxx.compute.c new file mode 100644 index 0000000..0831674 --- /dev/null +++ b/src/GeneralizedPolynomial-Uniform/Lagrange/2d.coeffs/2d.cube.order2.smooth0/coeffs-dxx.compute.c @@ -0,0 +1,13 @@ +fp t52; +fp t51; + t52 = RATIONAL(-2.0,3.0); + t51 = RATIONAL(1.0,3.0); + coeffs_dxx->coeff_m1_m1 = t51; + coeffs_dxx->coeff_0_m1 = t52; + coeffs_dxx->coeff_p1_m1 = t51; + coeffs_dxx->coeff_m1_0 = t51; + coeffs_dxx->coeff_0_0 = t52; + coeffs_dxx->coeff_p1_0 = t51; + coeffs_dxx->coeff_m1_p1 = t51; + coeffs_dxx->coeff_0_p1 = t52; + coeffs_dxx->coeff_p1_p1 = t51; diff --git a/src/GeneralizedPolynomial-Uniform/Lagrange/2d.coeffs/2d.cube.order2.smooth0/coeffs-dxy.compute.c b/src/GeneralizedPolynomial-Uniform/Lagrange/2d.coeffs/2d.cube.order2.smooth0/coeffs-dxy.compute.c new file mode 100644 index 0000000..f886764 --- /dev/null +++ b/src/GeneralizedPolynomial-Uniform/Lagrange/2d.coeffs/2d.cube.order2.smooth0/coeffs-dxy.compute.c @@ -0,0 +1,15 @@ +fp t55; +fp t54; +fp t53; + t55 = RATIONAL(0.0,1.0); + t54 = RATIONAL(1.0,4.0); + t53 = RATIONAL(-1.0,4.0); + coeffs_dxy->coeff_m1_m1 = t54; + coeffs_dxy->coeff_0_m1 = t55; + coeffs_dxy->coeff_p1_m1 = t53; + coeffs_dxy->coeff_m1_0 = t55; + coeffs_dxy->coeff_0_0 = t55; + coeffs_dxy->coeff_p1_0 = t55; + coeffs_dxy->coeff_m1_p1 = t53; + coeffs_dxy->coeff_0_p1 = t55; + coeffs_dxy->coeff_p1_p1 = t54; diff --git a/src/GeneralizedPolynomial-Uniform/Lagrange/2d.coeffs/2d.cube.order2.smooth0/coeffs-dy.compute.c b/src/GeneralizedPolynomial-Uniform/Lagrange/2d.coeffs/2d.cube.order2.smooth0/coeffs-dy.compute.c new file mode 100644 index 0000000..f18c6ef --- /dev/null +++ b/src/GeneralizedPolynomial-Uniform/Lagrange/2d.coeffs/2d.cube.order2.smooth0/coeffs-dy.compute.c @@ -0,0 +1,21 @@ +fp t46; +fp t50; +fp t49; +fp t48; +fp t47; +fp t43; + t46 = RATIONAL(1.0,3.0)*y; + t50 = RATIONAL(-1.0,4.0)*x+t46; + t49 = RATIONAL(1.0,4.0)*x+t46; + t48 = RATIONAL(1.0,6.0); + t47 = RATIONAL(-1.0,6.0); + t43 = RATIONAL(-2.0,3.0)*y; + coeffs_dy->coeff_m1_m1 = t47+t49; + coeffs_dy->coeff_0_m1 = t47+t46; + coeffs_dy->coeff_p1_m1 = t47+t50; + coeffs_dy->coeff_m1_0 = t43; + coeffs_dy->coeff_0_0 = t43; + coeffs_dy->coeff_p1_0 = t43; + coeffs_dy->coeff_m1_p1 = t48+t50; + coeffs_dy->coeff_0_p1 = t46+t48; + coeffs_dy->coeff_p1_p1 = t48+t49; diff --git a/src/GeneralizedPolynomial-Uniform/Lagrange/2d.coeffs/2d.cube.order2.smooth0/coeffs-dyy.compute.c b/src/GeneralizedPolynomial-Uniform/Lagrange/2d.coeffs/2d.cube.order2.smooth0/coeffs-dyy.compute.c new file mode 100644 index 0000000..47224b8 --- /dev/null +++ b/src/GeneralizedPolynomial-Uniform/Lagrange/2d.coeffs/2d.cube.order2.smooth0/coeffs-dyy.compute.c @@ -0,0 +1,13 @@ +fp t57; +fp t56; + t57 = RATIONAL(-2.0,3.0); + t56 = RATIONAL(1.0,3.0); + coeffs_dyy->coeff_m1_m1 = t56; + coeffs_dyy->coeff_0_m1 = t56; + coeffs_dyy->coeff_p1_m1 = t56; + coeffs_dyy->coeff_m1_0 = t57; + coeffs_dyy->coeff_0_0 = t57; + coeffs_dyy->coeff_p1_0 = t57; + coeffs_dyy->coeff_m1_p1 = t56; + coeffs_dyy->coeff_0_p1 = t56; + coeffs_dyy->coeff_p1_p1 = t56; diff --git a/src/GeneralizedPolynomial-Uniform/Lagrange/2d.coeffs/2d.cube.order3.smooth0/coeff-I.compute.c b/src/GeneralizedPolynomial-Uniform/Lagrange/2d.coeffs/2d.cube.order3.smooth0/coeff-I.compute.c deleted file mode 100644 index a022b8a..0000000 --- a/src/GeneralizedPolynomial-Uniform/Lagrange/2d.coeffs/2d.cube.order3.smooth0/coeff-I.compute.c +++ /dev/null @@ -1,196 +0,0 @@ -fp t160, - t131, - t159, - t133, - t158, - t117, - t66, - t132, - t128, - t87, - t157, - t129, - t84, - t156, - t124, - t68, - t130, - t116, - t77, - t155, - t126, - t70, - t125, - t83, - t154, - t71, - t153, - t108, - t122, - t69, - t152, - t151, - t63, - t150, - t67, - t86, - t149, - t65, - t148, - t85, - t147, - t79, - t146, - t82, - t145, - t144, - t143, - t142, - t123, - t141, - t101, - t140, - t139, - t138, - t137, - t136, - t127, - t121, - t119, - t118, - t115, - t113, - t112, - t111, - t109, - t106, - t105, - t103, - t102, - t100, - t98, - t97, - t96, - t95, - t94, - t93, - t92, - t91, - t90, - t89, - t88, - t81, - t80, - t78, - t76, - t75, - t74, - t72; - t160 = x*y; - t131 = x*x; - t159 = t131*y; - t133 = y*y; - t158 = x*t133; - t117 = RATIONAL(1.0,40.0); - t66 = t117*t159; - t132 = t133*y; - t128 = RATIONAL(-1.0,8.0); - t87 = t128*t132; - t157 = t66+t87; - t129 = RATIONAL(1.0,8.0); - t84 = t129*t132; - t156 = t66+t84; - t124 = RATIONAL(3.0,40.0); - t68 = t124*t158; - t130 = t131*x; - t116 = RATIONAL(1.0,24.0); - t77 = t116*t130; - t155 = t68+t77; - t126 = RATIONAL(-3.0,40.0); - t70 = t126*t159; - t125 = RATIONAL(-1.0,24.0); - t83 = t125*t132; - t154 = t70+t83; - t71 = t124*t159; - t153 = t71+t83; - t108 = RATIONAL(11.0,80.0); - t122 = RATIONAL(-1.0,40.0); - t69 = t122*t158; - t152 = t69+t108*t133; - t151 = RATIONAL(-1.0,100.0)*t160+RATIONAL(37.0,200.0); - t63 = t122*t159; - t150 = t63+t87; - t67 = t117*t158; - t86 = t129*t130; - t149 = t67+t86; - t65 = t126*t158; - t148 = t65+t77; - t85 = t128*t130; - t147 = t67+t85; - t79 = t116*t132; - t146 = t70+t79; - t82 = t125*t130; - t145 = t68+t82; - t144 = t63+t84; - t143 = t65+t82; - t142 = t71+t79; - t123 = RATIONAL(-1.0,50.0); - t141 = t123+RATIONAL(-13.0,100.0)*t160; - t101 = RATIONAL(-17.0,200.0); - t140 = RATIONAL(-9.0,100.0)*t160+t101; - t139 = RATIONAL(11.0,50.0)+t123*t160; - t138 = RATIONAL(2.0,25.0)*t160+RATIONAL(-2.0,25.0); - t137 = RATIONAL(7.0,100.0)*t160+RATIONAL(13.0,100.0); - t136 = t133+t131; - t127 = RATIONAL(1.0,80.0); - t121 = RATIONAL(9.0,80.0); - t119 = RATIONAL(-1.0,80.0); - t118 = RATIONAL(7.0,80.0); - t115 = RATIONAL(3.0,80.0); - t113 = RATIONAL(-19.0,80.0); - t112 = RATIONAL(-17.0,80.0); - t111 = RATIONAL(13.0,80.0); - t109 = RATIONAL(-23.0,80.0); - t106 = RATIONAL(-63.0,400.0); - t105 = RATIONAL(-49.0,400.0); - t103 = RATIONAL(-31.0,400.0); - t102 = RATIONAL(69.0,400.0); - t100 = RATIONAL(-57.0,400.0); - t98 = RATIONAL(117.0,400.0); - t97 = RATIONAL(43.0,1200.0); - t96 = RATIONAL(103.0,400.0); - t95 = RATIONAL(37.0,1200.0); - t94 = RATIONAL(111.0,400.0); - t93 = RATIONAL(89.0,1200.0); - t92 = RATIONAL(-109.0,1200.0); - t91 = RATIONAL(-223.0,1200.0); - t90 = RATIONAL(-157.0,1200.0); - t89 = RATIONAL(-149.0,1200.0); - t88 = RATIONAL(-131.0,1200.0); - t81 = t121*t133; - t80 = t118*t131; - t78 = t118*t133; - t76 = t121*t131; - t75 = t111*t131; - t74 = t111*t133; - t72 = t108*t131; - coeff_I_m1_m1 = t75+RATIONAL(-23.0,200.0)+t74+t92*y+(t92+RATIONAL(6.0, -25.0)*y)*x+t143+t154; - coeff_I_0_m1 = t109*t131+t91*y+t86+t103*x+t139+t152+t153; - coeff_I_p1_m1 = t81+t94*x+t90*y+t80+t141+t147+t153; - coeff_I_p2_m1 = t78+t93*y+t88*x+t115*t131+t140+t154+t155; - coeff_I_m1_0 = t72+t109*t133+t91*x+t103*y+t139+t144+t145; - coeff_I_0_0 = RATIONAL(93.0,200.0)+t100*y+t136*RATIONAL(-21.0,80.0)+( -RATIONAL(-1.0,25.0)*y+t100)*x+t149+t156; - coeff_I_p1_0 = t69+t76+t98*x+t85+t113*t133+t106*y+t151+t156; - coeff_I_p2_0 = t97*x+t127*t131+t105*y+t112*t133+t137+t144+t148; - coeff_I_m1_p1 = t76+t90*x+t78+t94*y+t141+t145+t157; - coeff_I_0_p1 = t113*t131+t81+t106*x+t98*y+t149+t150+t151; - coeff_I_p1_p1 = t72+t101+t96*y+t85+(t96+RATIONAL(3.0,50.0)*y)*x+t150+t152 -; - coeff_I_p2_p1 = t102*y+t74+t119*t131+t95*x+t138+t148+t157; - coeff_I_m1_p2 = t80+t88*y+t93*x+t115*t133+t140+t142+t143; - coeff_I_0_p2 = t127*t133+t86+t97*y+t105*x+t112*t131+t69+t137+t146; - coeff_I_p1_p2 = t95*y+t75+t119*t133+t102*x+t138+t146+t147; - coeff_I_p2_p2 = t89*y+RATIONAL(7.0,200.0)+t136*RATIONAL(-3.0,80.0)+(t89+ -RATIONAL(-3.0,50.0)*y)*x+t142+t155; diff --git a/src/GeneralizedPolynomial-Uniform/Lagrange/2d.coeffs/2d.cube.order3.smooth0/coeff-dx.compute.c b/src/GeneralizedPolynomial-Uniform/Lagrange/2d.coeffs/2d.cube.order3.smooth0/coeff-dx.compute.c deleted file mode 100644 index 5bfa6f4..0000000 --- a/src/GeneralizedPolynomial-Uniform/Lagrange/2d.coeffs/2d.cube.order3.smooth0/coeff-dx.compute.c +++ /dev/null @@ -1,95 +0,0 @@ -fp t201, - t161, - t200, - t184, - t188, - t174, - t187, - t180, - t199, - t164, - t198, - t183, - t175, - t197, - t162, - t196, - t181, - t195, - t194, - t163, - t193, - t186, - t176, - t182, - t192, - t185, - t172, - t191, - t179, - t190, - t189, - t178, - t177, - t169, - t167, - t166, - t165; - t201 = x*y; - t161 = RATIONAL(1.0,20.0)*t201; - t200 = t161+RATIONAL(9.0,40.0)*x; - t184 = RATIONAL(-1.0,40.0); - t188 = y*y; - t174 = t184*t188; - t187 = x*x; - t180 = RATIONAL(-3.0,8.0)*t187; - t199 = t174+t180; - t164 = RATIONAL(-3.0,20.0)*t201; - t198 = t164+RATIONAL(13.0,40.0)*x; - t183 = RATIONAL(1.0,40.0); - t175 = t183*t188; - t197 = t175+t180; - t162 = RATIONAL(-1.0,20.0)*t201; - t196 = t162+RATIONAL(11.0,40.0)*x; - t181 = RATIONAL(3.0,8.0)*t187; - t195 = t174+t181; - t194 = t175+t181; - t163 = RATIONAL(3.0,20.0)*t201; - t193 = RATIONAL(7.0,40.0)*x+t163; - t186 = RATIONAL(-3.0,40.0); - t176 = t186*t188; - t182 = RATIONAL(-1.0,8.0)*t187; - t192 = t176+t182; - t185 = RATIONAL(3.0,40.0); - t172 = t185*t188; - t191 = t172+t182; - t179 = RATIONAL(1.0,8.0)*t187; - t190 = t172+t179; - t189 = t176+t179; - t178 = RATIONAL(2.0,25.0)*y; - t177 = RATIONAL(-1.0,50.0)*y; - t169 = RATIONAL(-9.0,100.0)*y; - t167 = RATIONAL(-1.0,100.0)*y; - t166 = RATIONAL(7.0,100.0)*y; - t165 = RATIONAL(-13.0,100.0)*y; - coeff_dx_m1_m1 = RATIONAL(-109.0,1200.0)+RATIONAL(6.0,25.0)*y+t192+t198; - coeff_dx_0_m1 = RATIONAL(-31.0,400.0)+t177+RATIONAL(-23.0,40.0)*x+t163+ -t195; - coeff_dx_p1_m1 = RATIONAL(111.0,400.0)+t165+t193+t197; - coeff_dx_p2_m1 = t169+RATIONAL(-131.0,1200.0)+t164+t185*x+t190; - coeff_dx_m1_0 = RATIONAL(-223.0,1200.0)+t177+t191+t196; - coeff_dx_0_0 = RATIONAL(-1.0,25.0)*y+t161+RATIONAL(-21.0,40.0)*x+RATIONAL -(-57.0,400.0)+t194; - coeff_dx_p1_0 = t167+RATIONAL(117.0,400.0)+t199+t200; - coeff_dx_p2_0 = t183*x+t162+RATIONAL(43.0,1200.0)+t166+t189; - coeff_dx_m1_p1 = RATIONAL(-157.0,1200.0)+t165+t191+t200; - coeff_dx_0_p1 = RATIONAL(-19.0,40.0)*x+RATIONAL(-63.0,400.0)+t162+t167+ -t194; - coeff_dx_p1_p1 = RATIONAL(3.0,50.0)*y+RATIONAL(103.0,400.0)+t196+t199; - coeff_dx_p2_p1 = t161+t184*x+RATIONAL(37.0,1200.0)+t178+t189; - coeff_dx_m1_p2 = RATIONAL(89.0,1200.0)+t169+t192+t193; - coeff_dx_0_p2 = RATIONAL(-17.0,40.0)*x+t166+RATIONAL(-49.0,400.0)+t164+ -t195; - coeff_dx_p1_p2 = t178+RATIONAL(69.0,400.0)+t197+t198; - coeff_dx_p2_p2 = t163+RATIONAL(-149.0,1200.0)+t186*x+RATIONAL(-3.0,50.0)* -y+t190; diff --git a/src/GeneralizedPolynomial-Uniform/Lagrange/2d.coeffs/2d.cube.order3.smooth0/coeff-dxx.compute.c b/src/GeneralizedPolynomial-Uniform/Lagrange/2d.coeffs/2d.cube.order3.smooth0/coeff-dxx.compute.c deleted file mode 100644 index 29b5494..0000000 --- a/src/GeneralizedPolynomial-Uniform/Lagrange/2d.coeffs/2d.cube.order3.smooth0/coeff-dxx.compute.c +++ /dev/null @@ -1,40 +0,0 @@ -fp t244, - t258, - t246, - t257, - t245, - t256, - t243, - t255, - t250, - t249, - t248, - t247; - t244 = RATIONAL(-3.0,20.0)*y; - t258 = t244+RATIONAL(13.0,40.0); - t246 = RATIONAL(-1.0,20.0)*y; - t257 = RATIONAL(11.0,40.0)+t246; - t245 = RATIONAL(1.0,20.0)*y; - t256 = RATIONAL(9.0,40.0)+t245; - t243 = RATIONAL(3.0,20.0)*y; - t255 = t243+RATIONAL(7.0,40.0); - t250 = RATIONAL(1.0,4.0)*x; - t249 = RATIONAL(-1.0,4.0)*x; - t248 = RATIONAL(3.0,4.0)*x; - t247 = RATIONAL(-3.0,4.0)*x; - coeff_dxx_m1_m1 = t249+t258; - coeff_dxx_0_m1 = t248+RATIONAL(-23.0,40.0)+t243; - coeff_dxx_p1_m1 = t247+t255; - coeff_dxx_p2_m1 = RATIONAL(3.0,40.0)+t250+t244; - coeff_dxx_m1_0 = t249+t257; - coeff_dxx_0_0 = RATIONAL(-21.0,40.0)+t245+t248; - coeff_dxx_p1_0 = t247+t256; - coeff_dxx_p2_0 = RATIONAL(1.0,40.0)+t246+t250; - coeff_dxx_m1_p1 = t249+t256; - coeff_dxx_0_p1 = t248+RATIONAL(-19.0,40.0)+t246; - coeff_dxx_p1_p1 = t247+t257; - coeff_dxx_p2_p1 = RATIONAL(-1.0,40.0)+t250+t245; - coeff_dxx_m1_p2 = t249+t255; - coeff_dxx_0_p2 = t244+RATIONAL(-17.0,40.0)+t248; - coeff_dxx_p1_p2 = t247+t258; - coeff_dxx_p2_p2 = RATIONAL(-3.0,40.0)+t250+t243; diff --git a/src/GeneralizedPolynomial-Uniform/Lagrange/2d.coeffs/2d.cube.order3.smooth0/coeff-dxy.compute.c b/src/GeneralizedPolynomial-Uniform/Lagrange/2d.coeffs/2d.cube.order3.smooth0/coeff-dxy.compute.c deleted file mode 100644 index 20e2106..0000000 --- a/src/GeneralizedPolynomial-Uniform/Lagrange/2d.coeffs/2d.cube.order3.smooth0/coeff-dxy.compute.c +++ /dev/null @@ -1,52 +0,0 @@ -fp t276, - t275, - t274, - t273, - t272, - t271, - t270, - t269, - t268, - t267, - t266, - t265, - t264, - t263, - t262, - t261, - t260, - t259; - t276 = RATIONAL(-1.0,50.0); - t275 = RATIONAL(2.0,25.0); - t274 = RATIONAL(-3.0,20.0); - t273 = RATIONAL(3.0,20.0); - t272 = RATIONAL(-1.0,20.0); - t271 = RATIONAL(1.0,20.0); - t270 = RATIONAL(-9.0,100.0); - t269 = RATIONAL(-1.0,100.0); - t268 = RATIONAL(7.0,100.0); - t267 = RATIONAL(-13.0,100.0); - t266 = t272*y; - t265 = t271*y; - t264 = t273*x; - t263 = t274*x; - t262 = t271*x; - t261 = t272*x; - t260 = t274*y; - t259 = t273*y; - coeff_dxy_m1_m1 = t260+RATIONAL(6.0,25.0)+t263; - coeff_dxy_0_m1 = t276+t264+t266; - coeff_dxy_p1_m1 = t264+t267+t265; - coeff_dxy_p2_m1 = t263+t270+t259; - coeff_dxy_m1_0 = t276+t259+t261; - coeff_dxy_0_0 = t262+RATIONAL(-1.0,25.0)+t265; - coeff_dxy_p1_0 = t269+t262+t266; - coeff_dxy_p2_0 = t260+t268+t261; - coeff_dxy_m1_p1 = t262+t267+t259; - coeff_dxy_0_p1 = t269+t265+t261; - coeff_dxy_p1_p1 = RATIONAL(3.0,50.0)+t266+t261; - coeff_dxy_p2_p1 = t260+t275+t262; - coeff_dxy_m1_p2 = t264+t260+t270; - coeff_dxy_0_p2 = t268+t263+t266; - coeff_dxy_p1_p2 = t263+t275+t265; - coeff_dxy_p2_p2 = t259+RATIONAL(-3.0,50.0)+t264; diff --git a/src/GeneralizedPolynomial-Uniform/Lagrange/2d.coeffs/2d.cube.order3.smooth0/coeff-dy.compute.c b/src/GeneralizedPolynomial-Uniform/Lagrange/2d.coeffs/2d.cube.order3.smooth0/coeff-dy.compute.c deleted file mode 100644 index c766157..0000000 --- a/src/GeneralizedPolynomial-Uniform/Lagrange/2d.coeffs/2d.cube.order3.smooth0/coeff-dy.compute.c +++ /dev/null @@ -1,95 +0,0 @@ -fp t242, - t225, - t228, - t212, - t229, - t220, - t241, - t203, - t240, - t204, - t239, - t205, - t238, - t224, - t213, - t237, - t202, - t236, - t226, - t215, - t222, - t235, - t227, - t218, - t234, - t221, - t233, - t232, - t223, - t231, - t230, - t219, - t217, - t209, - t208, - t207, - t206; - t242 = x*y; - t225 = RATIONAL(-1.0,40.0); - t228 = x*x; - t212 = t225*t228; - t229 = y*y; - t220 = RATIONAL(-3.0,8.0)*t229; - t241 = t212+t220; - t203 = RATIONAL(1.0,20.0)*t242; - t240 = t203+RATIONAL(9.0,40.0)*y; - t204 = RATIONAL(-1.0,20.0)*t242; - t239 = RATIONAL(11.0,40.0)*y+t204; - t205 = RATIONAL(-3.0,20.0)*t242; - t238 = t205+RATIONAL(13.0,40.0)*y; - t224 = RATIONAL(1.0,40.0); - t213 = t224*t228; - t237 = t213+t220; - t202 = RATIONAL(3.0,20.0)*t242; - t236 = RATIONAL(7.0,40.0)*y+t202; - t226 = RATIONAL(3.0,40.0); - t215 = t226*t228; - t222 = RATIONAL(-1.0,8.0)*t229; - t235 = t215+t222; - t227 = RATIONAL(-3.0,40.0); - t218 = t227*t228; - t234 = t218+t222; - t221 = RATIONAL(3.0,8.0)*t229; - t233 = t212+t221; - t232 = t213+t221; - t223 = RATIONAL(1.0,8.0)*t229; - t231 = t218+t223; - t230 = t215+t223; - t219 = RATIONAL(2.0,25.0)*x; - t217 = RATIONAL(-1.0,50.0)*x; - t209 = RATIONAL(-9.0,100.0)*x; - t208 = RATIONAL(7.0,100.0)*x; - t207 = RATIONAL(-1.0,100.0)*x; - t206 = RATIONAL(-13.0,100.0)*x; - coeff_dy_m1_m1 = RATIONAL(6.0,25.0)*x+RATIONAL(-109.0,1200.0)+t234+t238; - coeff_dy_0_m1 = RATIONAL(-223.0,1200.0)+t217+t235+t239; - coeff_dy_p1_m1 = t206+RATIONAL(-157.0,1200.0)+t235+t240; - coeff_dy_p2_m1 = t209+RATIONAL(89.0,1200.0)+t234+t236; - coeff_dy_m1_0 = RATIONAL(-23.0,40.0)*y+t202+RATIONAL(-31.0,400.0)+t217+ -t233; - coeff_dy_0_0 = RATIONAL(-1.0,25.0)*x+RATIONAL(-57.0,400.0)+t203+RATIONAL( --21.0,40.0)*y+t232; - coeff_dy_p1_0 = RATIONAL(-63.0,400.0)+t207+t204+RATIONAL(-19.0,40.0)*y+ -t232; - coeff_dy_p2_0 = t205+RATIONAL(-17.0,40.0)*y+RATIONAL(-49.0,400.0)+t208+ -t233; - coeff_dy_m1_p1 = RATIONAL(111.0,400.0)+t206+t236+t237; - coeff_dy_0_p1 = RATIONAL(117.0,400.0)+t207+t240+t241; - coeff_dy_p1_p1 = RATIONAL(3.0,50.0)*x+RATIONAL(103.0,400.0)+t239+t241; - coeff_dy_p2_p1 = t219+RATIONAL(69.0,400.0)+t237+t238; - coeff_dy_m1_p2 = t209+RATIONAL(-131.0,1200.0)+t205+t226*y+t230; - coeff_dy_0_p2 = t204+t208+t224*y+RATIONAL(43.0,1200.0)+t231; - coeff_dy_p1_p2 = t225*y+t219+t203+RATIONAL(37.0,1200.0)+t231; - coeff_dy_p2_p2 = RATIONAL(-149.0,1200.0)+t202+RATIONAL(-3.0,50.0)*x+t227* -y+t230; diff --git a/src/GeneralizedPolynomial-Uniform/Lagrange/2d.coeffs/2d.cube.order3.smooth0/coeff-dyy.compute.c b/src/GeneralizedPolynomial-Uniform/Lagrange/2d.coeffs/2d.cube.order3.smooth0/coeff-dyy.compute.c deleted file mode 100644 index b449df9..0000000 --- a/src/GeneralizedPolynomial-Uniform/Lagrange/2d.coeffs/2d.cube.order3.smooth0/coeff-dyy.compute.c +++ /dev/null @@ -1,40 +0,0 @@ -fp t277, - t292, - t279, - t291, - t280, - t290, - t278, - t289, - t284, - t283, - t282, - t281; - t277 = RATIONAL(3.0,20.0)*x; - t292 = t277+RATIONAL(7.0,40.0); - t279 = RATIONAL(-1.0,20.0)*x; - t291 = t279+RATIONAL(11.0,40.0); - t280 = RATIONAL(-3.0,20.0)*x; - t290 = t280+RATIONAL(13.0,40.0); - t278 = RATIONAL(1.0,20.0)*x; - t289 = t278+RATIONAL(9.0,40.0); - t284 = RATIONAL(-1.0,4.0)*y; - t283 = RATIONAL(-3.0,4.0)*y; - t282 = RATIONAL(1.0,4.0)*y; - t281 = RATIONAL(3.0,4.0)*y; - coeff_dyy_m1_m1 = t284+t290; - coeff_dyy_0_m1 = t284+t291; - coeff_dyy_p1_m1 = t284+t289; - coeff_dyy_p2_m1 = t284+t292; - coeff_dyy_m1_0 = RATIONAL(-23.0,40.0)+t281+t277; - coeff_dyy_0_0 = t281+RATIONAL(-21.0,40.0)+t278; - coeff_dyy_p1_0 = RATIONAL(-19.0,40.0)+t279+t281; - coeff_dyy_p2_0 = t281+RATIONAL(-17.0,40.0)+t280; - coeff_dyy_m1_p1 = t283+t292; - coeff_dyy_0_p1 = t283+t289; - coeff_dyy_p1_p1 = t283+t291; - coeff_dyy_p2_p1 = t283+t290; - coeff_dyy_m1_p2 = t282+RATIONAL(3.0,40.0)+t280; - coeff_dyy_0_p2 = t282+RATIONAL(1.0,40.0)+t279; - coeff_dyy_p1_p2 = RATIONAL(-1.0,40.0)+t278+t282; - coeff_dyy_p2_p2 = t277+t282+RATIONAL(-3.0,40.0); diff --git a/src/GeneralizedPolynomial-Uniform/Lagrange/2d.coeffs/2d.cube.order3.smooth0/coeffs-I.compute.c b/src/GeneralizedPolynomial-Uniform/Lagrange/2d.coeffs/2d.cube.order3.smooth0/coeffs-I.compute.c new file mode 100644 index 0000000..ac11344 --- /dev/null +++ b/src/GeneralizedPolynomial-Uniform/Lagrange/2d.coeffs/2d.cube.order3.smooth0/coeffs-I.compute.c @@ -0,0 +1,190 @@ +fp t158; +fp t131; +fp t157; +fp t132; +fp t156; +fp t117; +fp t63; +fp t129; +fp t130; +fp t84; +fp t155; +fp t154; +fp t128; +fp t87; +fp t153; +fp t126; +fp t64; +fp t133; +fp t116; +fp t76; +fp t152; +fp t67; +fp t86; +fp t151; +fp t123; +fp t71; +fp t150; +fp t124; +fp t65; +fp t149; +fp t125; +fp t81; +fp t148; +fp t85; +fp t147; +fp t146; +fp t112; +fp t69; +fp t145; +fp t110; +fp t66; +fp t144; +fp t143; +fp t70; +fp t120; +fp t142; +fp t122; +fp t141; +fp t140; +fp t101; +fp t139; +fp t138; +fp t137; +fp t136; +fp t127; +fp t121; +fp t118; +fp t115; +fp t113; +fp t109; +fp t107; +fp t106; +fp t104; +fp t103; +fp t102; +fp t100; +fp t98; +fp t97; +fp t96; +fp t95; +fp t94; +fp t93; +fp t92; +fp t91; +fp t90; +fp t89; +fp t88; +fp t83; +fp t82; +fp t80; +fp t78; +fp t77; +fp t74; +fp t73; + t158 = x*y; + t131 = x*x; + t157 = t131*y; + t132 = y*y; + t156 = x*t132; + t117 = RATIONAL(1.0,40.0); + t63 = t117*t157; + t129 = RATIONAL(1.0,8.0); + t130 = t132*y; + t84 = t129*t130; + t155 = t63+t84; + t154 = RATIONAL(-1.0,100.0)*t158+RATIONAL(37.0,200.0); + t128 = RATIONAL(-1.0,8.0); + t87 = t128*t130; + t153 = t63+t87; + t126 = RATIONAL(-3.0,40.0); + t64 = t126*t156; + t133 = t131*x; + t116 = RATIONAL(1.0,24.0); + t76 = t116*t133; + t152 = t64+t76; + t67 = t117*t156; + t86 = t129*t133; + t151 = t67+t86; + t123 = RATIONAL(-1.0,40.0); + t71 = t123*t156; + t150 = t71+t86; + t124 = RATIONAL(3.0,40.0); + t65 = t124*t156; + t149 = t65+t76; + t125 = RATIONAL(-1.0,24.0); + t81 = t125*t133; + t148 = t64+t81; + t85 = t128*t133; + t147 = t71+t85; + t146 = t65+t81; + t112 = RATIONAL(13.0,80.0); + t69 = t126*t157; + t145 = t69+t112*t131; + t110 = RATIONAL(11.0,80.0); + t66 = t123*t157; + t144 = t66+t110*t131; + t143 = t67+t85; + t70 = t124*t157; + t120 = RATIONAL(7.0,80.0); + t142 = t70+t120*t131; + t122 = RATIONAL(-1.0,50.0); + t141 = t122+RATIONAL(-13.0,100.0)*t158; + t140 = RATIONAL(2.0,25.0)*t158+RATIONAL(-2.0,25.0); + t101 = RATIONAL(-17.0,200.0); + t139 = RATIONAL(-9.0,100.0)*t158+t101; + t138 = RATIONAL(11.0,50.0)+t122*t158; + t137 = RATIONAL(7.0,100.0)*t158+RATIONAL(13.0,100.0); + t136 = t131+t132; + t127 = RATIONAL(1.0,80.0); + t121 = RATIONAL(9.0,80.0); + t118 = RATIONAL(-1.0,80.0); + t115 = RATIONAL(3.0,80.0); + t113 = RATIONAL(-17.0,80.0); + t109 = RATIONAL(-23.0,80.0); + t107 = RATIONAL(-19.0,80.0); + t106 = RATIONAL(-63.0,400.0); + t104 = RATIONAL(-49.0,400.0); + t103 = RATIONAL(-31.0,400.0); + t102 = RATIONAL(69.0,400.0); + t100 = RATIONAL(-57.0,400.0); + t98 = RATIONAL(117.0,400.0); + t97 = RATIONAL(43.0,1200.0); + t96 = RATIONAL(103.0,400.0); + t95 = RATIONAL(111.0,400.0); + t94 = RATIONAL(37.0,1200.0); + t93 = RATIONAL(89.0,1200.0); + t92 = RATIONAL(-109.0,1200.0); + t91 = RATIONAL(-223.0,1200.0); + t90 = RATIONAL(-157.0,1200.0); + t89 = RATIONAL(-131.0,1200.0); + t88 = RATIONAL(-149.0,1200.0); + t83 = t121*t131; + t82 = t125*t130; + t80 = t121*t132; + t78 = t116*t130; + t77 = t120*t132; + t74 = t112*t132; + t73 = t110*t132; + coeffs_I->coeff_m1_m1 = t74+t82+t92*y+RATIONAL(-23.0,200.0)+(t92+RATIONAL +(6.0,25.0)*y)*x+t145+t148; + coeffs_I->coeff_0_m1 = t70+t91*y+t73+t103*x+t82+t109*t131+t138+t150; + coeffs_I->coeff_p1_m1 = t80+t82+t95*x+t90*y+t141+t142+t143; + coeffs_I->coeff_p2_m1 = t93*y+t89*x+t77+t82+t69+t115*t131+t139+t149; + coeffs_I->coeff_m1_0 = t109*t132+t84+t91*x+t103*y+t138+t144+t146; + coeffs_I->coeff_0_0 = RATIONAL(93.0,200.0)+t100*y+t136*RATIONAL(-21.0, +80.0)+(t100+RATIONAL(-1.0,25.0)*y)*x+t151+t155; + coeffs_I->coeff_p1_0 = t107*t132+t106*y+t83+t98*x+t147+t154+t155; + coeffs_I->coeff_p2_0 = t66+t113*t132+t127*t131+t84+t104*y+t97*x+t137+t152 +; + coeffs_I->coeff_m1_p1 = t95*y+t90*x+t77+t83+t141+t146+t153; + coeffs_I->coeff_0_p1 = t87+t106*x+t66+t80+t98*y+t107*t131+t151+t154; + coeffs_I->coeff_p1_p1 = t87+t73+t101+t96*y+(t96+RATIONAL(3.0,50.0)*y)*x+ +t144+t147; + coeffs_I->coeff_p2_p1 = t118*t131+t102*y+t74+t94*x+t140+t152+t153; + coeffs_I->coeff_m1_p2 = t93*x+t89*y+t115*t132+t78+t139+t142+t148; + coeffs_I->coeff_0_p2 = t97*y+t69+t127*t132+t113*t131+t78+t104*x+t137+t150 +; + coeffs_I->coeff_p1_p2 = t102*x+t94*y+t78+t118*t132+t140+t143+t145; + coeffs_I->coeff_p2_p2 = t70+t88*y+t78+RATIONAL(7.0,200.0)+t136*RATIONAL( +-3.0,80.0)+(t88+RATIONAL(-3.0,50.0)*y)*x+t149; diff --git a/src/GeneralizedPolynomial-Uniform/Lagrange/2d.coeffs/2d.cube.order3.smooth0/coeffs-dx.compute.c b/src/GeneralizedPolynomial-Uniform/Lagrange/2d.coeffs/2d.cube.order3.smooth0/coeffs-dx.compute.c new file mode 100644 index 0000000..0ea2d65 --- /dev/null +++ b/src/GeneralizedPolynomial-Uniform/Lagrange/2d.coeffs/2d.cube.order3.smooth0/coeffs-dx.compute.c @@ -0,0 +1,97 @@ +fp t199; +fp t185; +fp t178; +fp t181; +fp t186; +fp t173; +fp t198; +fp t159; +fp t197; +fp t182; +fp t175; +fp t196; +fp t161; +fp t195; +fp t162; +fp t194; +fp t160; +fp t193; +fp t179; +fp t184; +fp t176; +fp t192; +fp t183; +fp t171; +fp t191; +fp t180; +fp t190; +fp t189; +fp t177; +fp t188; +fp t187; +fp t174; +fp t169; +fp t167; +fp t165; +fp t164; +fp t163; + t199 = x*y; + t185 = x*x; + t178 = RATIONAL(-3.0,8.0)*t185; + t181 = RATIONAL(1.0,40.0); + t186 = y*y; + t173 = t181*t186; + t198 = t178+t173; + t159 = RATIONAL(-3.0,20.0)*t199; + t197 = RATIONAL(13.0,40.0)*x+t159; + t182 = RATIONAL(-1.0,40.0); + t175 = t182*t186; + t196 = t178+t175; + t161 = RATIONAL(-1.0,20.0)*t199; + t195 = t161+RATIONAL(11.0,40.0)*x; + t162 = RATIONAL(3.0,20.0)*t199; + t194 = RATIONAL(7.0,40.0)*x+t162; + t160 = RATIONAL(1.0,20.0)*t199; + t193 = RATIONAL(9.0,40.0)*x+t160; + t179 = RATIONAL(-1.0,8.0)*t185; + t184 = RATIONAL(-3.0,40.0); + t176 = t184*t186; + t192 = t179+t176; + t183 = RATIONAL(3.0,40.0); + t171 = t183*t186; + t191 = t179+t171; + t180 = RATIONAL(3.0,8.0)*t185; + t190 = t180+t175; + t189 = t180+t173; + t177 = RATIONAL(1.0,8.0)*t185; + t188 = t171+t177; + t187 = t176+t177; + t174 = RATIONAL(-1.0,50.0)*y; + t169 = RATIONAL(2.0,25.0)*y; + t167 = RATIONAL(-9.0,100.0)*y; + t165 = RATIONAL(-1.0,100.0)*y; + t164 = RATIONAL(7.0,100.0)*y; + t163 = RATIONAL(-13.0,100.0)*y; + coeffs_dx->coeff_m1_m1 = RATIONAL(6.0,25.0)*y+RATIONAL(-109.0,1200.0)+ +t192+t197; + coeffs_dx->coeff_0_m1 = t162+t174+RATIONAL(-31.0,400.0)+RATIONAL(-23.0, +40.0)*x+t190; + coeffs_dx->coeff_p1_m1 = RATIONAL(111.0,400.0)+t163+t194+t198; + coeffs_dx->coeff_p2_m1 = t183*x+RATIONAL(-131.0,1200.0)+t167+t159+t188; + coeffs_dx->coeff_m1_0 = RATIONAL(-223.0,1200.0)+t174+t191+t195; + coeffs_dx->coeff_0_0 = RATIONAL(-1.0,25.0)*y+t160+RATIONAL(-21.0,40.0)*x+ +RATIONAL(-57.0,400.0)+t189; + coeffs_dx->coeff_p1_0 = t165+RATIONAL(117.0,400.0)+t193+t196; + coeffs_dx->coeff_p2_0 = t161+t181*x+RATIONAL(43.0,1200.0)+t164+t187; + coeffs_dx->coeff_m1_p1 = RATIONAL(-157.0,1200.0)+t163+t191+t193; + coeffs_dx->coeff_0_p1 = RATIONAL(-63.0,400.0)+t165+t161+RATIONAL(-19.0, +40.0)*x+t189; + coeffs_dx->coeff_p1_p1 = RATIONAL(103.0,400.0)+RATIONAL(3.0,50.0)*y+t195+ +t196; + coeffs_dx->coeff_p2_p1 = t160+t169+RATIONAL(37.0,1200.0)+t182*x+t187; + coeffs_dx->coeff_m1_p2 = RATIONAL(89.0,1200.0)+t167+t192+t194; + coeffs_dx->coeff_0_p2 = t159+t164+RATIONAL(-17.0,40.0)*x+RATIONAL(-49.0, +400.0)+t190; + coeffs_dx->coeff_p1_p2 = RATIONAL(69.0,400.0)+t169+t197+t198; + coeffs_dx->coeff_p2_p2 = RATIONAL(-149.0,1200.0)+t162+RATIONAL(-3.0,50.0) +*y+t184*x+t188; diff --git a/src/GeneralizedPolynomial-Uniform/Lagrange/2d.coeffs/2d.cube.order3.smooth0/coeffs-dxx.compute.c b/src/GeneralizedPolynomial-Uniform/Lagrange/2d.coeffs/2d.cube.order3.smooth0/coeffs-dxx.compute.c new file mode 100644 index 0000000..ba139bf --- /dev/null +++ b/src/GeneralizedPolynomial-Uniform/Lagrange/2d.coeffs/2d.cube.order3.smooth0/coeffs-dxx.compute.c @@ -0,0 +1,40 @@ +fp t244; +fp t256; +fp t242; +fp t255; +fp t241; +fp t254; +fp t243; +fp t253; +fp t248; +fp t247; +fp t246; +fp t245; + t244 = RATIONAL(-3.0,20.0)*y; + t256 = t244+RATIONAL(13.0,40.0); + t242 = RATIONAL(-1.0,20.0)*y; + t255 = t242+RATIONAL(11.0,40.0); + t241 = RATIONAL(1.0,20.0)*y; + t254 = t241+RATIONAL(9.0,40.0); + t243 = RATIONAL(3.0,20.0)*y; + t253 = t243+RATIONAL(7.0,40.0); + t248 = RATIONAL(-1.0,4.0)*x; + t247 = RATIONAL(1.0,4.0)*x; + t246 = RATIONAL(-3.0,4.0)*x; + t245 = RATIONAL(3.0,4.0)*x; + coeffs_dxx->coeff_m1_m1 = t248+t256; + coeffs_dxx->coeff_0_m1 = t245+t243+RATIONAL(-23.0,40.0); + coeffs_dxx->coeff_p1_m1 = t246+t253; + coeffs_dxx->coeff_p2_m1 = t244+RATIONAL(3.0,40.0)+t247; + coeffs_dxx->coeff_m1_0 = t248+t255; + coeffs_dxx->coeff_0_0 = RATIONAL(-21.0,40.0)+t241+t245; + coeffs_dxx->coeff_p1_0 = t246+t254; + coeffs_dxx->coeff_p2_0 = t247+t242+RATIONAL(1.0,40.0); + coeffs_dxx->coeff_m1_p1 = t248+t254; + coeffs_dxx->coeff_0_p1 = RATIONAL(-19.0,40.0)+t242+t245; + coeffs_dxx->coeff_p1_p1 = t246+t255; + coeffs_dxx->coeff_p2_p1 = RATIONAL(-1.0,40.0)+t241+t247; + coeffs_dxx->coeff_m1_p2 = t248+t253; + coeffs_dxx->coeff_0_p2 = t244+t245+RATIONAL(-17.0,40.0); + coeffs_dxx->coeff_p1_p2 = t246+t256; + coeffs_dxx->coeff_p2_p2 = RATIONAL(-3.0,40.0)+t247+t243; diff --git a/src/GeneralizedPolynomial-Uniform/Lagrange/2d.coeffs/2d.cube.order3.smooth0/coeffs-dxy.compute.c b/src/GeneralizedPolynomial-Uniform/Lagrange/2d.coeffs/2d.cube.order3.smooth0/coeffs-dxy.compute.c new file mode 100644 index 0000000..e9208a2 --- /dev/null +++ b/src/GeneralizedPolynomial-Uniform/Lagrange/2d.coeffs/2d.cube.order3.smooth0/coeffs-dxy.compute.c @@ -0,0 +1,52 @@ +fp t274; +fp t273; +fp t272; +fp t271; +fp t270; +fp t269; +fp t268; +fp t267; +fp t266; +fp t265; +fp t264; +fp t263; +fp t262; +fp t261; +fp t260; +fp t259; +fp t258; +fp t257; + t274 = RATIONAL(-3.0,20.0); + t273 = RATIONAL(3.0,20.0); + t272 = RATIONAL(2.0,25.0); + t271 = RATIONAL(-1.0,20.0); + t270 = RATIONAL(1.0,20.0); + t269 = RATIONAL(-1.0,50.0); + t268 = RATIONAL(-1.0,100.0); + t267 = RATIONAL(7.0,100.0); + t266 = RATIONAL(-9.0,100.0); + t265 = RATIONAL(-13.0,100.0); + t264 = t274*x; + t263 = t273*x; + t262 = t271*x; + t261 = t274*y; + t260 = t270*x; + t259 = t273*y; + t258 = t271*y; + t257 = t270*y; + coeffs_dxy->coeff_m1_m1 = t264+RATIONAL(6.0,25.0)+t261; + coeffs_dxy->coeff_0_m1 = t263+t269+t258; + coeffs_dxy->coeff_p1_m1 = t257+t265+t263; + coeffs_dxy->coeff_p2_m1 = t259+t264+t266; + coeffs_dxy->coeff_m1_0 = t259+t262+t269; + coeffs_dxy->coeff_0_0 = t260+RATIONAL(-1.0,25.0)+t257; + coeffs_dxy->coeff_p1_0 = t258+t268+t260; + coeffs_dxy->coeff_p2_0 = t261+t262+t267; + coeffs_dxy->coeff_m1_p1 = t259+t265+t260; + coeffs_dxy->coeff_0_p1 = t262+t257+t268; + coeffs_dxy->coeff_p1_p1 = t262+RATIONAL(3.0,50.0)+t258; + coeffs_dxy->coeff_p2_p1 = t261+t260+t272; + coeffs_dxy->coeff_m1_p2 = t266+t261+t263; + coeffs_dxy->coeff_0_p2 = t267+t258+t264; + coeffs_dxy->coeff_p1_p2 = t264+t257+t272; + coeffs_dxy->coeff_p2_p2 = t263+t259+RATIONAL(-3.0,50.0); diff --git a/src/GeneralizedPolynomial-Uniform/Lagrange/2d.coeffs/2d.cube.order3.smooth0/coeffs-dy.compute.c b/src/GeneralizedPolynomial-Uniform/Lagrange/2d.coeffs/2d.cube.order3.smooth0/coeffs-dy.compute.c new file mode 100644 index 0000000..c12d9d1 --- /dev/null +++ b/src/GeneralizedPolynomial-Uniform/Lagrange/2d.coeffs/2d.cube.order3.smooth0/coeffs-dy.compute.c @@ -0,0 +1,97 @@ +fp t240; +fp t227; +fp t219; +fp t226; +fp t222; +fp t211; +fp t239; +fp t201; +fp t238; +fp t223; +fp t213; +fp t237; +fp t200; +fp t236; +fp t203; +fp t235; +fp t202; +fp t234; +fp t224; +fp t216; +fp t221; +fp t233; +fp t225; +fp t217; +fp t232; +fp t218; +fp t231; +fp t230; +fp t220; +fp t229; +fp t228; +fp t215; +fp t210; +fp t207; +fp t206; +fp t205; +fp t204; + t240 = x*y; + t227 = y*y; + t219 = RATIONAL(-3.0,8.0)*t227; + t226 = x*x; + t222 = RATIONAL(1.0,40.0); + t211 = t222*t226; + t239 = t219+t211; + t201 = RATIONAL(1.0,20.0)*t240; + t238 = t201+RATIONAL(9.0,40.0)*y; + t223 = RATIONAL(-1.0,40.0); + t213 = t223*t226; + t237 = t219+t213; + t200 = RATIONAL(3.0,20.0)*t240; + t236 = t200+RATIONAL(7.0,40.0)*y; + t203 = RATIONAL(-3.0,20.0)*t240; + t235 = t203+RATIONAL(13.0,40.0)*y; + t202 = RATIONAL(-1.0,20.0)*t240; + t234 = t202+RATIONAL(11.0,40.0)*y; + t224 = RATIONAL(3.0,40.0); + t216 = t224*t226; + t221 = RATIONAL(-1.0,8.0)*t227; + t233 = t216+t221; + t225 = RATIONAL(-3.0,40.0); + t217 = t225*t226; + t232 = t217+t221; + t218 = RATIONAL(1.0,8.0)*t227; + t231 = t218+t217; + t230 = t216+t218; + t220 = RATIONAL(3.0,8.0)*t227; + t229 = t220+t213; + t228 = t220+t211; + t215 = RATIONAL(-1.0,50.0)*x; + t210 = RATIONAL(2.0,25.0)*x; + t207 = RATIONAL(-9.0,100.0)*x; + t206 = RATIONAL(-1.0,100.0)*x; + t205 = RATIONAL(7.0,100.0)*x; + t204 = RATIONAL(-13.0,100.0)*x; + coeffs_dy->coeff_m1_m1 = RATIONAL(6.0,25.0)*x+RATIONAL(-109.0,1200.0)+ +t232+t235; + coeffs_dy->coeff_0_m1 = t215+RATIONAL(-223.0,1200.0)+t233+t234; + coeffs_dy->coeff_p1_m1 = RATIONAL(-157.0,1200.0)+t204+t233+t238; + coeffs_dy->coeff_p2_m1 = t207+RATIONAL(89.0,1200.0)+t232+t236; + coeffs_dy->coeff_m1_0 = RATIONAL(-23.0,40.0)*y+RATIONAL(-31.0,400.0)+t200 ++t215+t229; + coeffs_dy->coeff_0_0 = RATIONAL(-21.0,40.0)*y+RATIONAL(-57.0,400.0)+ +RATIONAL(-1.0,25.0)*x+t201+t228; + coeffs_dy->coeff_p1_0 = t202+RATIONAL(-19.0,40.0)*y+RATIONAL(-63.0,400.0) ++t206+t228; + coeffs_dy->coeff_p2_0 = t203+RATIONAL(-49.0,400.0)+RATIONAL(-17.0,40.0)*y ++t205+t229; + coeffs_dy->coeff_m1_p1 = t204+RATIONAL(111.0,400.0)+t236+t239; + coeffs_dy->coeff_0_p1 = t206+RATIONAL(117.0,400.0)+t237+t238; + coeffs_dy->coeff_p1_p1 = RATIONAL(3.0,50.0)*x+RATIONAL(103.0,400.0)+t234+ +t237; + coeffs_dy->coeff_p2_p1 = t210+RATIONAL(69.0,400.0)+t235+t239; + coeffs_dy->coeff_m1_p2 = t207+t203+RATIONAL(-131.0,1200.0)+t224*y+t230; + coeffs_dy->coeff_0_p2 = t205+t202+t222*y+RATIONAL(43.0,1200.0)+t231; + coeffs_dy->coeff_p1_p2 = t201+RATIONAL(37.0,1200.0)+t223*y+t210+t231; + coeffs_dy->coeff_p2_p2 = t225*y+t200+RATIONAL(-3.0,50.0)*x+RATIONAL( +-149.0,1200.0)+t230; diff --git a/src/GeneralizedPolynomial-Uniform/Lagrange/2d.coeffs/2d.cube.order3.smooth0/coeffs-dyy.compute.c b/src/GeneralizedPolynomial-Uniform/Lagrange/2d.coeffs/2d.cube.order3.smooth0/coeffs-dyy.compute.c new file mode 100644 index 0000000..6bbdde7 --- /dev/null +++ b/src/GeneralizedPolynomial-Uniform/Lagrange/2d.coeffs/2d.cube.order3.smooth0/coeffs-dyy.compute.c @@ -0,0 +1,40 @@ +fp t276; +fp t290; +fp t277; +fp t289; +fp t278; +fp t288; +fp t275; +fp t287; +fp t282; +fp t281; +fp t280; +fp t279; + t276 = RATIONAL(-1.0,20.0)*x; + t290 = t276+RATIONAL(11.0,40.0); + t277 = RATIONAL(3.0,20.0)*x; + t289 = RATIONAL(7.0,40.0)+t277; + t278 = RATIONAL(-3.0,20.0)*x; + t288 = t278+RATIONAL(13.0,40.0); + t275 = RATIONAL(1.0,20.0)*x; + t287 = t275+RATIONAL(9.0,40.0); + t282 = RATIONAL(1.0,4.0)*y; + t281 = RATIONAL(-1.0,4.0)*y; + t280 = RATIONAL(3.0,4.0)*y; + t279 = RATIONAL(-3.0,4.0)*y; + coeffs_dyy->coeff_m1_m1 = t281+t288; + coeffs_dyy->coeff_0_m1 = t281+t290; + coeffs_dyy->coeff_p1_m1 = t281+t287; + coeffs_dyy->coeff_p2_m1 = t281+t289; + coeffs_dyy->coeff_m1_0 = t277+t280+RATIONAL(-23.0,40.0); + coeffs_dyy->coeff_0_0 = t280+t275+RATIONAL(-21.0,40.0); + coeffs_dyy->coeff_p1_0 = t280+RATIONAL(-19.0,40.0)+t276; + coeffs_dyy->coeff_p2_0 = RATIONAL(-17.0,40.0)+t280+t278; + coeffs_dyy->coeff_m1_p1 = t279+t289; + coeffs_dyy->coeff_0_p1 = t279+t287; + coeffs_dyy->coeff_p1_p1 = t279+t290; + coeffs_dyy->coeff_p2_p1 = t279+t288; + coeffs_dyy->coeff_m1_p2 = t278+t282+RATIONAL(3.0,40.0); + coeffs_dyy->coeff_0_p2 = t282+RATIONAL(1.0,40.0)+t276; + coeffs_dyy->coeff_p1_p2 = t275+t282+RATIONAL(-1.0,40.0); + coeffs_dyy->coeff_p2_p2 = t277+RATIONAL(-3.0,40.0)+t282; 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 deleted file mode 100644 index cfe391a..0000000 --- a/src/GeneralizedPolynomial-Uniform/Lagrange/2d.coeffs/2d.cube.order4.smooth0/coeff-I.compute.c +++ /dev/null @@ -1,243 +0,0 @@ -fp t417, - t418, - t468, - t414, - t411, - t361, - t416, - t373, - t435, - t446, - t467, - t410, - t367, - t368, - t437, - t466, - t465, - t441, - t464, - t463, - t395, - t445, - t400, - t403, - t366, - t462, - t397, - t404, - t362, - t401, - t461, - t389, - t409, - t316, - t460, - t402, - t312, - t384, - t459, - t412, - t443, - t407, - t309, - t458, - t295, - t399, - t442, - t457, - t406, - t405, - t319, - t456, - t444, - t370, - t455, - t315, - t454, - t325, - t453, - t358, - t398, - t294, - t452, - t451, - t413, - t364, - t450, - t357, - t380, - t381, - t449, - t332, - t415, - t356, - t448, - t447, - t440, - t360, - t439, - t438, - t436, - t390, - t408, - t363, - t392, - t434, - t386, - t374, - t433, - t372, - t432, - t431, - t317, - t430, - t429, - t326, - t428, - t427, - t385, - t426, - t425, - t388, - t424, - t423, - t394, - t393, - t391, - t387; - t417 = x*x; - t418 = y*y; - t468 = t418+t417; - t414 = t417*x; - t411 = RATIONAL(1.0,60.0); - t361 = t411*t414; - t416 = y*t418; - t373 = t411*t416; - t435 = x*t373+y*t361; - t446 = x*y; - t467 = RATIONAL(-11.0,150.0)*t446+t435; - t410 = RATIONAL(-1.0,60.0); - t367 = t410*t416; - t368 = t410*t414; - t437 = x*t367+y*t368; - t466 = RATIONAL(37.0,300.0)*t446+t437; - t465 = t437+RATIONAL(11.0,150.0)*t446; - t441 = t417*t418; - t464 = RATIONAL(-1.0,98.0)*t441+RATIONAL(-99.0,1225.0); - t463 = RATIONAL(-37.0,300.0)*t446+t435; - t395 = RATIONAL(-31.0,210.0); - t445 = x*t418; - t400 = RATIONAL(1.0,140.0); - t403 = RATIONAL(1.0,30.0); - t366 = t403*t414; - t462 = t400*t445+t395*x+t366; - t397 = RATIONAL(31.0,210.0); - t404 = RATIONAL(-1.0,30.0); - t362 = t404*t414; - t401 = RATIONAL(-1.0,140.0); - t461 = t362+t401*t445+t397*x; - t389 = RATIONAL(31.0,420.0); - t409 = RATIONAL(-1.0,35.0); - t316 = t409*t445; - t460 = t389*x+t316+t368; - t402 = RATIONAL(1.0,35.0); - t312 = t402*t445; - t384 = RATIONAL(-31.0,420.0); - t459 = t312+t384*x+t361; - t412 = RATIONAL(1.0,84.0); - t443 = y*t417; - t407 = RATIONAL(-1.0,70.0); - t309 = t407*t443; - t458 = t309+t412*y+t373; - t295 = RATIONAL(21.0,200.0)*t446; - t399 = RATIONAL(1.0,120.0); - t442 = x*t416; - t457 = t295+t399*t442+y*t362; - t406 = RATIONAL(-1.0,84.0); - t405 = RATIONAL(1.0,70.0); - t319 = t405*t443; - t456 = t406*y+t319+t367; - t444 = t414*y; - t370 = t404*t416; - t455 = t295+t399*t444+x*t370; - t315 = t405*t445; - t454 = t406*x+t315+t368; - t325 = t407*t445; - t453 = t361+t412*x+t325; - t358 = t403*t416; - t398 = RATIONAL(-1.0,120.0); - t294 = RATIONAL(-21.0,200.0)*t446; - t452 = t398*t444+t294+x*t358; - t451 = y*t366+t398*t442+t294; - t413 = t417*t417; - t364 = t404*t413; - t450 = RATIONAL(-24.0,1225.0)+t364+RATIONAL(1.0,196.0)*t441+t468*RATIONAL -(181.0,1470.0); - t357 = t399*t413; - t380 = RATIONAL(71.0,5880.0); - t381 = RATIONAL(113.0,735.0); - t449 = t357+t380*t417+t381*t418+t464; - t332 = RATIONAL(1.0,49.0)*t441; - t415 = t418*t418; - t356 = t399*t415; - t448 = RATIONAL(51.0,1225.0)+t357+t332+t356+t468*RATIONAL(-289.0,5880.0); - t447 = t380*t418+t364+t381*t417+t356+t464; - t440 = RATIONAL(1.0,98.0)*t441+RATIONAL(246.0,1225.0); - t360 = t404*t415; - t439 = t360+t370; - t438 = t358+t360; - t436 = RATIONAL(-1.0,49.0)*t441+RATIONAL(96.0,1225.0); - t390 = RATIONAL(83.0,735.0); - t408 = RATIONAL(1.0,20.0); - t363 = t408*t413; - t392 = RATIONAL(-53.0,196.0); - t434 = t360+t390*t418+t363+t392*t417+t440; - t386 = RATIONAL(-41.0,196.0); - t374 = RATIONAL(191.0,5880.0); - t433 = t386*t417+t374*t418+t363+t356+t436; - t372 = t408*t415; - t432 = t357+t374*t417+t386*t418+t372+t436; - t431 = t364+t392*t418+t390*t417+t372+t440; - t317 = t402*t443; - t430 = t384*y+t317+t373+t448; - t429 = t395*y+t400*t443+t438+t450; - t326 = t409*t443; - t428 = t389*y+t367+t326+t448; - t427 = t397*y+t401*t443+t439+t450; - t385 = RATIONAL(11.0,105.0); - t426 = t385*x+t315+t362+t447; - t425 = t385*y+t319+t439+t449; - t388 = RATIONAL(-11.0,105.0); - t424 = t388*y+t309+t438+t449; - t423 = t388*x+t366+t325+t447; - t394 = RATIONAL(17.0,420.0); - t393 = RATIONAL(17.0,105.0); - t391 = RATIONAL(-17.0,105.0); - t387 = RATIONAL(-17.0,420.0); - coeff_I_m2_m2 = t428+t460+t467; - coeff_I_m1_m2 = t423+t456+t457; - coeff_I_0_m2 = t317+t367+t387*y+t433; - coeff_I_p1_m2 = t426+t451+t456; - coeff_I_p2_m2 = t428+t459+t465; - coeff_I_m2_m1 = t424+t454+t455; - coeff_I_m1_m1 = t429+t462+t466; - coeff_I_0_m1 = t319+t391*y+t358+t434; - coeff_I_p1_m1 = t429+t461+t463; - coeff_I_p2_m1 = t424+t452+t453; - coeff_I_m2_0 = t368+t312+t387*x+t432; - coeff_I_m1_0 = t366+t315+t391*x+t431; - coeff_I_0_0 = RATIONAL(541.0,1225.0)+t363+t372+t332+t468*RATIONAL(-57.0, -196.0); - coeff_I_p1_0 = t325+t362+t393*x+t431; - coeff_I_p2_0 = t394*x+t316+t361+t432; - coeff_I_m2_p1 = t425+t452+t454; - coeff_I_m1_p1 = t427+t462+t463; - coeff_I_0_p1 = t309+t370+t393*y+t434; - coeff_I_p1_p1 = t427+t461+t466; - coeff_I_p2_p1 = t425+t453+t455; - coeff_I_m2_p2 = t430+t460+t465; - coeff_I_m1_p2 = t423+t451+t458; - coeff_I_0_p2 = t394*y+t326+t373+t433; - coeff_I_p1_p2 = t426+t457+t458; - coeff_I_p2_p2 = t430+t459+t467; 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 deleted file mode 100644 index e56ee74..0000000 --- a/src/GeneralizedPolynomial-Uniform/Lagrange/2d.coeffs/2d.cube.order4.smooth0/coeff-dx.compute.c +++ /dev/null @@ -1,155 +0,0 @@ -fp t538, - t507, - t539, - t537, - t559, - t573, - t560, - t572, - t514, - t554, - t571, - t526, - t536, - t516, - t570, - t531, - t517, - t569, - t481, - t568, - t561, - t494, - t567, - t495, - t566, - t525, - t506, - t565, - t510, - t564, - t511, - t563, - t562, - t475, - t527, - t509, - t558, - t513, - t557, - t529, - t512, - t556, - t555, - t518, - t478, - t553, - t479, - t552, - t551, - t550, - t490, - t549, - t548, - t547, - t546, - t545, - t544, - t543, - t542, - t486, - t482, - t477, - t474, - t471, - t469; - t538 = x*x; - t507 = RATIONAL(1.0,20.0)*t538; - t539 = y*y; - t537 = y*t539; - t559 = y*t507+RATIONAL(1.0,60.0)*t537; - t573 = RATIONAL(-37.0,300.0)*y+t559; - t560 = x*t539; - t572 = RATIONAL(1.0,98.0)*t560+RATIONAL(181.0,735.0)*x; - t514 = RATIONAL(-1.0,20.0)*t538; - t554 = y*t514+RATIONAL(-1.0,60.0)*t537; - t571 = RATIONAL(37.0,300.0)*y+t554; - t526 = RATIONAL(1.0,30.0); - t536 = x*t538; - t516 = t526*t536; - t570 = RATIONAL(71.0,2940.0)*x+t516; - t531 = RATIONAL(-1.0,35.0); - t517 = t531*t539; - t569 = RATIONAL(31.0,420.0)+t514+t517; - t481 = RATIONAL(2.0,49.0)*t560; - t568 = RATIONAL(-289.0,2940.0)*x+t481+t516; - t561 = y*t538; - t494 = RATIONAL(-21.0,200.0)*y; - t567 = RATIONAL(-1.0,40.0)*t561+t494+t526*t537; - t495 = RATIONAL(21.0,200.0)*y; - t566 = RATIONAL(1.0,40.0)*t561+t495+RATIONAL(-1.0,30.0)*t537; - t525 = RATIONAL(1.0,35.0); - t506 = t525*t539; - t565 = RATIONAL(-31.0,420.0)+t507+t506; - t510 = RATIONAL(-1.0,10.0)*t538; - t564 = y*t510+t495+RATIONAL(1.0,120.0)*t537; - t511 = RATIONAL(1.0,10.0)*t538; - t563 = y*t511+t494+RATIONAL(-1.0,120.0)*t537; - t562 = x*y; - t475 = RATIONAL(-1.0,49.0)*t560; - t527 = RATIONAL(1.0,70.0); - t509 = t527*t539; - t558 = t475+t509; - t513 = RATIONAL(-2.0,15.0)*t536; - t557 = t510+t513; - t529 = RATIONAL(-1.0,70.0); - t512 = t529*t539; - t556 = t475+t512; - t555 = t511+t513; - t518 = RATIONAL(1.0,5.0)*t536; - t478 = RATIONAL(-2.0,49.0)*t560; - t553 = t518+RATIONAL(-41.0,98.0)*x+t478; - t479 = RATIONAL(1.0,49.0)*t560; - t552 = t518+RATIONAL(-53.0,98.0)*x+t479; - t551 = RATIONAL(191.0,2940.0)*x+t478+t516; - t550 = RATIONAL(166.0,735.0)*x+t479+t513; - t490 = RATIONAL(226.0,735.0)*x; - t549 = RATIONAL(11.0,105.0)+t490+t557+t558; - t548 = t490+RATIONAL(-11.0,105.0)+t555+t556; - t547 = RATIONAL(-31.0,210.0)+RATIONAL(1.0,140.0)*t539+t555+t572; - t546 = RATIONAL(31.0,210.0)+RATIONAL(-1.0,140.0)*t539+t557+t572; - t545 = RATIONAL(1.0,84.0)+t507+t556+t570; - t544 = RATIONAL(-11.0,150.0)*y+t559+t568; - t543 = RATIONAL(-1.0,84.0)+t514+t558+t570; - t542 = RATIONAL(11.0,150.0)*y+t554+t568; - t486 = t525*t562; - t482 = RATIONAL(-2.0,35.0)*t562; - t477 = RATIONAL(2.0,35.0)*t562; - t474 = t531*t562; - t471 = t527*t562; - t469 = t529*t562; - coeff_dx_m2_m2 = t482+t544+t569; - coeff_dx_m1_m2 = t486+t548+t564; - coeff_dx_0_m2 = t477+t553; - coeff_dx_p1_m2 = t486+t549+t563; - coeff_dx_p2_m2 = t482+t542+t565; - coeff_dx_m2_m1 = t474+t543+t566; - coeff_dx_m1_m1 = t471+t547+t571; - coeff_dx_0_m1 = t486+t552; - coeff_dx_p1_m1 = t471+t546+t573; - coeff_dx_p2_m1 = t474+t545+t567; - coeff_dx_m2_0 = t514+RATIONAL(-17.0,420.0)+t506+t551; - coeff_dx_m1_0 = t511+RATIONAL(-17.0,105.0)+t509+t550; - coeff_dx_0_0 = t481+t518+RATIONAL(-57.0,98.0)*x; - coeff_dx_p1_0 = RATIONAL(17.0,105.0)+t510+t512+t550; - coeff_dx_p2_0 = t507+RATIONAL(17.0,420.0)+t517+t551; - coeff_dx_m2_p1 = t486+t543+t567; - coeff_dx_m1_p1 = t469+t547+t573; - coeff_dx_0_p1 = t474+t552; - coeff_dx_p1_p1 = t469+t546+t571; - coeff_dx_p2_p1 = t486+t545+t566; - coeff_dx_m2_p2 = t477+t542+t569; - coeff_dx_m1_p2 = t474+t548+t563; - coeff_dx_0_p2 = t482+t553; - coeff_dx_p1_p2 = t474+t549+t564; - coeff_dx_p2_p2 = t477+t544+t565; 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 deleted file mode 100644 index 3ac1e8c..0000000 --- a/src/GeneralizedPolynomial-Uniform/Lagrange/2d.coeffs/2d.cube.order4.smooth0/coeff-dxx.compute.c +++ /dev/null @@ -1,107 +0,0 @@ -fp t713, - t683, - t708, - t712, - t687, - t730, - t696, - t729, - t728, - t690, - t727, - t726, - t697, - t686, - t698, - t691, - t693, - t725, - t699, - t724, - t682, - t723, - t722, - t695, - t721, - t720, - t689, - t719, - t718, - t717, - t716, - t678, - t715, - t676, - t714, - t694, - t685, - t681, - t680, - t679, - t677; - t713 = y*y; - t683 = RATIONAL(-1.0,49.0)*t713; - t708 = RATIONAL(1.0,10.0); - t712 = x*x; - t687 = t708*t712; - t730 = RATIONAL(71.0,2940.0)+t683+t687; - t696 = RATIONAL(-2.0,5.0)*t712; - t729 = t696+RATIONAL(226.0,735.0)+t683; - t728 = t696+RATIONAL(181.0,735.0)+RATIONAL(1.0,98.0)*t713; - t690 = RATIONAL(2.0,49.0)*t713; - t727 = t690+RATIONAL(-289.0,2940.0)+t687; - t726 = x*y; - t697 = RATIONAL(1.0,5.0)*x; - t686 = RATIONAL(-1.0,10.0)*x; - t698 = RATIONAL(-1.0,5.0)*x; - t691 = t708*x; - t693 = RATIONAL(-2.0,49.0)*t713; - t725 = t693+RATIONAL(191.0,2940.0)+t687; - t699 = RATIONAL(3.0,5.0)*t712; - t724 = t693+t699+RATIONAL(-41.0,98.0); - t682 = RATIONAL(1.0,49.0)*t713; - t723 = t699+RATIONAL(-53.0,98.0)+t682; - t722 = t696+RATIONAL(166.0,735.0)+t682; - t695 = RATIONAL(1.0,35.0)*y; - t721 = t695+t730; - t720 = RATIONAL(-1.0,70.0)*y+t728; - t689 = RATIONAL(-1.0,35.0)*y; - t719 = t689+t730; - t718 = t698+t729; - t717 = RATIONAL(1.0,70.0)*y+t728; - t716 = t697+t729; - t678 = y*t686; - t715 = t678+t727; - t676 = y*t691; - t714 = t676+t727; - t694 = RATIONAL(2.0,35.0)*y; - t685 = RATIONAL(-2.0,35.0)*y; - t681 = y*t697; - t680 = y*t698; - t679 = RATIONAL(1.0,20.0)*t726; - t677 = RATIONAL(-1.0,20.0)*t726; - coeff_dxx_m2_m2 = t686+t685+t714; - coeff_dxx_m1_m2 = t680+t695+t716; - coeff_dxx_0_m2 = t694+t724; - coeff_dxx_p1_m2 = t695+t681+t718; - coeff_dxx_p2_m2 = t691+t685+t715; - coeff_dxx_m2_m1 = t679+t686+t719; - coeff_dxx_m1_m1 = t697+t678+t717; - coeff_dxx_0_m1 = t695+t723; - coeff_dxx_p1_m1 = t698+t676+t717; - coeff_dxx_p2_m1 = t691+t677+t719; - coeff_dxx_m2_0 = t686+t725; - coeff_dxx_m1_0 = t697+t722; - coeff_dxx_0_0 = RATIONAL(-57.0,98.0)+t690+t699; - coeff_dxx_p1_0 = t698+t722; - coeff_dxx_p2_0 = t691+t725; - coeff_dxx_m2_p1 = t677+t686+t721; - coeff_dxx_m1_p1 = t697+t676+t720; - coeff_dxx_0_p1 = t689+t723; - coeff_dxx_p1_p1 = t678+t698+t720; - coeff_dxx_p2_p1 = t691+t679+t721; - coeff_dxx_m2_p2 = t694+t686+t715; - coeff_dxx_m1_p2 = t689+t681+t716; - coeff_dxx_0_p2 = t685+t724; - coeff_dxx_p1_p2 = t680+t689+t718; - coeff_dxx_p2_p2 = t694+t691+t714; 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 deleted file mode 100644 index 0505d86..0000000 --- a/src/GeneralizedPolynomial-Uniform/Lagrange/2d.coeffs/2d.cube.order4.smooth0/coeff-dxy.compute.c +++ /dev/null @@ -1,119 +0,0 @@ -fp t779, - t778, - t797, - t772, - t774, - t796, - t795, - t770, - t739, - t766, - t752, - t794, - t742, - t759, - t793, - t792, - t733, - t791, - t790, - t789, - t788, - t787, - t735, - t786, - t732, - t785, - t773, - t777, - t784, - t783, - t782, - t781, - t780, - t776, - t771, - t768, - t767, - t758, - t757, - t755, - t750, - t743, - t738, - t737, - t736, - t734, - t731; - t779 = y*y; - t778 = x*x; - t797 = t778+t779; - t772 = RATIONAL(1.0,10.0); - t774 = RATIONAL(-1.0,40.0); - t796 = t772*t778+t774*t779; - t795 = t774*t778+t772*t779; - t770 = RATIONAL(-1.0,35.0); - t739 = t770*y; - t766 = RATIONAL(1.0,35.0); - t752 = t766*x; - t794 = t739+t752; - t742 = t770*x; - t759 = t766*y; - t793 = t742+t759; - t792 = x*y; - t733 = RATIONAL(-2.0,49.0)*t792; - t791 = t733+RATIONAL(-21.0,200.0); - t790 = t797*RATIONAL(-1.0,20.0); - t789 = t797*RATIONAL(1.0,20.0); - t788 = t733+RATIONAL(21.0,200.0); - t787 = t752+t759+t791; - t735 = RATIONAL(4.0,49.0)*t792; - t786 = t735+RATIONAL(11.0,150.0)+t790; - t732 = RATIONAL(1.0,49.0)*t792; - t785 = t732+RATIONAL(37.0,300.0)+t790; - t773 = RATIONAL(-1.0,10.0); - t777 = RATIONAL(1.0,40.0); - t784 = t773*t779+t777*t778+t788; - t783 = t735+RATIONAL(-11.0,150.0)+t789; - t782 = t732+RATIONAL(-37.0,300.0)+t789; - t781 = t773*t778+t777*t779+t788; - t780 = t739+t742+t791; - t776 = RATIONAL(-2.0,35.0); - t771 = RATIONAL(2.0,35.0); - t768 = RATIONAL(-1.0,70.0); - t767 = RATIONAL(1.0,70.0); - t758 = t771*y; - t757 = t768*y; - t755 = t776*x; - t750 = t771*x; - t743 = t767*x; - t738 = t767*y; - t737 = t776*y; - t736 = t768*x; - t734 = RATIONAL(-4.0,49.0)*t792; - t731 = RATIONAL(2.0,49.0)*t792; - coeff_dxy_m2_m2 = t755+t737+t783; - coeff_dxy_m1_m2 = t781+t794; - coeff_dxy_0_m2 = t750+t734; - coeff_dxy_p1_m2 = t787+t796; - coeff_dxy_p2_m2 = t758+t755+t786; - coeff_dxy_m2_m1 = t784+t793; - coeff_dxy_m1_m1 = t738+t743+t785; - coeff_dxy_0_m1 = t752+t731; - coeff_dxy_p1_m1 = t757+t743+t782; - coeff_dxy_p2_m1 = t780+t795; - coeff_dxy_m2_0 = t734+t758; - coeff_dxy_m1_0 = t731+t759; - coeff_dxy_0_0 = t735; - coeff_dxy_p1_0 = t731+t739; - coeff_dxy_p2_0 = t734+t737; - coeff_dxy_m2_p1 = t787+t795; - coeff_dxy_m1_p1 = t738+t736+t782; - coeff_dxy_0_p1 = t731+t742; - coeff_dxy_p1_p1 = t757+t736+t785; - coeff_dxy_p2_p1 = t784+t794; - coeff_dxy_m2_p2 = t750+t737+t786; - coeff_dxy_m1_p2 = t780+t796; - coeff_dxy_0_p2 = t755+t734; - coeff_dxy_p1_p2 = t781+t793; - coeff_dxy_p2_p2 = t758+t750+t783; 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 deleted file mode 100644 index f4baa3a..0000000 --- a/src/GeneralizedPolynomial-Uniform/Lagrange/2d.coeffs/2d.cube.order4.smooth0/coeff-dy.compute.c +++ /dev/null @@ -1,147 +0,0 @@ -fp t644, - t609, - t630, - t643, - t616, - t675, - t636, - t621, - t622, - t674, - t661, - t601, - t641, - t673, - t631, - t600, - t672, - t613, - t671, - t662, - t642, - t611, - t670, - t587, - t669, - t619, - t590, - t668, - t634, - t612, - t667, - t666, - t610, - t665, - t632, - t618, - t664, - t663, - t660, - t659, - t574, - t658, - t586, - t657, - t623, - t656, - t655, - t654, - t653, - t652, - t651, - t650, - t649, - t648, - t647, - t592, - t584, - t583, - t582, - t578, - t576; - t644 = y*y; - t609 = RATIONAL(1.0,20.0)*t644; - t630 = RATIONAL(1.0,35.0); - t643 = x*x; - t616 = t630*t643; - t675 = t609+t616+RATIONAL(-31.0,420.0); - t636 = RATIONAL(-1.0,35.0); - t621 = t636*t643; - t622 = RATIONAL(-1.0,20.0)*t644; - t674 = t621+t622+RATIONAL(31.0,420.0); - t661 = x*t644; - t601 = RATIONAL(21.0,200.0)*x; - t641 = x*t643; - t673 = RATIONAL(1.0,40.0)*t661+t601+RATIONAL(-1.0,30.0)*t641; - t631 = RATIONAL(1.0,30.0); - t600 = RATIONAL(-21.0,200.0)*x; - t672 = t600+RATIONAL(-1.0,40.0)*t661+t631*t641; - t613 = RATIONAL(1.0,10.0)*t644; - t671 = RATIONAL(1.0,140.0)*t643+t613+RATIONAL(-31.0,210.0); - t662 = y*t643; - t642 = y*t644; - t611 = RATIONAL(-2.0,15.0)*t642; - t670 = RATIONAL(181.0,735.0)*y+t611+RATIONAL(1.0,98.0)*t662; - t587 = RATIONAL(-1.0,49.0)*t662; - t669 = RATIONAL(226.0,735.0)*y+t587+t611; - t619 = t631*t642; - t590 = RATIONAL(2.0,49.0)*t662; - t668 = RATIONAL(-289.0,2940.0)*y+t619+t590; - t634 = RATIONAL(-1.0,70.0); - t612 = t634*t643; - t667 = t613+t612+RATIONAL(-11.0,105.0); - t666 = RATIONAL(71.0,2940.0)*y+t587+t619; - t610 = RATIONAL(-1.0,10.0)*t644; - t665 = RATIONAL(31.0,210.0)+RATIONAL(-1.0,140.0)*t643+t610; - t632 = RATIONAL(1.0,70.0); - t618 = t632*t643; - t664 = t610+t618+RATIONAL(11.0,105.0); - t663 = x*y; - t660 = x*t622+RATIONAL(-1.0,60.0)*t641; - t659 = x*t609+RATIONAL(1.0,60.0)*t641; - t574 = RATIONAL(1.0,49.0)*t662; - t658 = RATIONAL(166.0,735.0)*y+t611+t574; - t586 = RATIONAL(-2.0,49.0)*t662; - t657 = RATIONAL(191.0,2940.0)*y+t619+t586; - t623 = RATIONAL(1.0,5.0)*t642; - t656 = RATIONAL(-41.0,98.0)*y+t623+t586; - t655 = RATIONAL(-53.0,98.0)*y+t623+t574; - t654 = RATIONAL(37.0,300.0)*x+t660+t670; - t653 = RATIONAL(11.0,150.0)*x+t660+t668; - t652 = RATIONAL(-37.0,300.0)*x+t659+t670; - t651 = RATIONAL(-1.0,84.0)+t618+t622+t666; - t650 = x*t610+t601+RATIONAL(1.0,120.0)*t641+t669; - t649 = RATIONAL(1.0,84.0)+t609+t612+t666; - t648 = RATIONAL(-11.0,150.0)*x+t659+t668; - t647 = t600+RATIONAL(-1.0,120.0)*t641+x*t613+t669; - t592 = RATIONAL(-2.0,35.0)*t663; - t584 = t630*t663; - t583 = t636*t663; - t582 = RATIONAL(2.0,35.0)*t663; - t578 = t634*t663; - t576 = t632*t663; - coeff_dy_m2_m2 = t592+t648+t674; - coeff_dy_m1_m2 = t583+t651+t673; - coeff_dy_0_m2 = t616+RATIONAL(-17.0,420.0)+t622+t657; - coeff_dy_p1_m2 = t584+t651+t672; - coeff_dy_p2_m2 = t582+t653+t674; - coeff_dy_m2_m1 = t584+t650+t667; - coeff_dy_m1_m1 = t576+t654+t671; - coeff_dy_0_m1 = t618+t613+RATIONAL(-17.0,105.0)+t658; - coeff_dy_p1_m1 = t578+t652+t671; - coeff_dy_p2_m1 = t583+t647+t667; - coeff_dy_m2_0 = t582+t656; - coeff_dy_m1_0 = t584+t655; - coeff_dy_0_0 = RATIONAL(-57.0,98.0)*y+t590+t623; - coeff_dy_p1_0 = t583+t655; - coeff_dy_p2_0 = t592+t656; - coeff_dy_m2_p1 = t584+t647+t664; - coeff_dy_m1_p1 = t576+t652+t665; - coeff_dy_0_p1 = t612+t610+RATIONAL(17.0,105.0)+t658; - coeff_dy_p1_p1 = t578+t654+t665; - coeff_dy_p2_p1 = t583+t650+t664; - coeff_dy_m2_p2 = t592+t653+t675; - coeff_dy_m1_p2 = t583+t649+t672; - coeff_dy_0_p2 = t621+RATIONAL(17.0,420.0)+t609+t657; - coeff_dy_p1_p2 = t584+t649+t673; - coeff_dy_p2_p2 = t582+t648+t675; 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 deleted file mode 100644 index 9e72985..0000000 --- a/src/GeneralizedPolynomial-Uniform/Lagrange/2d.coeffs/2d.cube.order4.smooth0/coeff-dyy.compute.c +++ /dev/null @@ -1,115 +0,0 @@ -fp t834, - t813, - t854, - t853, - t830, - t835, - t814, - t852, - t816, - t851, - t850, - t815, - t808, - t818, - t820, - t821, - t849, - t848, - t810, - t847, - t804, - t846, - t819, - t845, - t844, - t843, - t842, - t841, - t801, - t840, - t798, - t839, - t838, - t837, - t836, - t817, - t812, - t811, - t809, - t806, - t805, - t803, - t802, - t800, - t799; - t834 = x*x; - t813 = RATIONAL(-1.0,49.0)*t834; - t854 = t813+RATIONAL(226.0,735.0); - t853 = RATIONAL(1.0,98.0)*t834+RATIONAL(181.0,735.0); - t830 = RATIONAL(1.0,10.0); - t835 = y*y; - t814 = t830*t835; - t852 = t813+t814+RATIONAL(71.0,2940.0); - t816 = RATIONAL(2.0,49.0)*t834; - t851 = t814+t816+RATIONAL(-289.0,2940.0); - t850 = x*y; - t815 = RATIONAL(-1.0,10.0)*y; - t808 = t830*y; - t818 = RATIONAL(-1.0,5.0)*y; - t820 = RATIONAL(1.0,5.0)*y; - t821 = RATIONAL(-2.0,5.0)*t835; - t849 = t818+t821; - t848 = t820+t821; - t810 = RATIONAL(-2.0,49.0)*t834; - t847 = t810+t814+RATIONAL(191.0,2940.0); - t804 = RATIONAL(1.0,49.0)*t834; - t846 = t804+t821+RATIONAL(166.0,735.0); - t819 = RATIONAL(3.0,5.0)*t835; - t845 = t810+t819+RATIONAL(-41.0,98.0); - t844 = t804+t819+RATIONAL(-53.0,98.0); - t843 = t849+t853; - t842 = t848+t853; - t841 = t849+t854; - t801 = x*t808; - t840 = t801+t851; - t798 = x*t815; - t839 = t798+t851; - t838 = t808+t852; - t837 = t848+t854; - t836 = t815+t852; - t817 = RATIONAL(-2.0,35.0)*x; - t812 = RATIONAL(-1.0,35.0)*x; - t811 = RATIONAL(1.0,35.0)*x; - t809 = RATIONAL(2.0,35.0)*x; - t806 = RATIONAL(1.0,70.0)*x; - t805 = RATIONAL(-1.0,70.0)*x; - t803 = x*t818; - t802 = x*t820; - t800 = RATIONAL(1.0,20.0)*t850; - t799 = RATIONAL(-1.0,20.0)*t850; - coeff_dyy_m2_m2 = t815+t817+t840; - coeff_dyy_m1_m2 = t800+t812+t836; - coeff_dyy_0_m2 = t815+t847; - coeff_dyy_p1_m2 = t799+t811+t836; - coeff_dyy_p2_m2 = t815+t809+t839; - coeff_dyy_m2_m1 = t803+t811+t837; - coeff_dyy_m1_m1 = t806+t798+t842; - coeff_dyy_0_m1 = t820+t846; - coeff_dyy_p1_m1 = t805+t801+t842; - coeff_dyy_p2_m1 = t802+t812+t837; - coeff_dyy_m2_0 = t809+t845; - coeff_dyy_m1_0 = t811+t844; - coeff_dyy_0_0 = RATIONAL(-57.0,98.0)+t819+t816; - coeff_dyy_p1_0 = t812+t844; - coeff_dyy_p2_0 = t817+t845; - coeff_dyy_m2_p1 = t802+t811+t841; - coeff_dyy_m1_p1 = t806+t801+t843; - coeff_dyy_0_p1 = t818+t846; - coeff_dyy_p1_p1 = t805+t798+t843; - coeff_dyy_p2_p1 = t812+t803+t841; - coeff_dyy_m2_p2 = t817+t808+t839; - coeff_dyy_m1_p2 = t812+t799+t838; - coeff_dyy_0_p2 = t808+t847; - coeff_dyy_p1_p2 = t811+t800+t838; - coeff_dyy_p2_p2 = t808+t809+t840; diff --git a/src/GeneralizedPolynomial-Uniform/Lagrange/2d.coeffs/2d.cube.order4.smooth0/coeffs-I.compute.c b/src/GeneralizedPolynomial-Uniform/Lagrange/2d.coeffs/2d.cube.order4.smooth0/coeffs-I.compute.c new file mode 100644 index 0000000..51bec8d --- /dev/null +++ b/src/GeneralizedPolynomial-Uniform/Lagrange/2d.coeffs/2d.cube.order4.smooth0/coeffs-I.compute.c @@ -0,0 +1,239 @@ +fp t411; +fp t412; +fp t464; +fp t414; +fp t404; +fp t356; +fp t413; +fp t360; +fp t433; +fp t442; +fp t463; +fp t407; +fp t365; +fp t364; +fp t436; +fp t462; +fp t435; +fp t461; +fp t438; +fp t389; +fp t398; +fp t406; +fp t362; +fp t460; +fp t401; +fp t368; +fp t395; +fp t397; +fp t459; +fp t439; +fp t359; +fp t458; +fp t370; +fp t291; +fp t396; +fp t440; +fp t457; +fp t405; +fp t322; +fp t409; +fp t456; +fp t402; +fp t400; +fp t327; +fp t455; +fp t383; +fp t454; +fp t311; +fp t453; +fp t390; +fp t307; +fp t452; +fp t451; +fp t384; +fp t408; +fp t318; +fp t450; +fp t399; +fp t296; +fp t449; +fp t448; +fp t387; +fp t403; +fp t330; +fp t447; +fp t416; +fp t354; +fp t415; +fp t355; +fp t314; +fp t446; +fp t369; +fp t358; +fp t445; +fp t376; +fp t379; +fp t444; +fp t443; +fp t441; +fp t437; +fp t434; +fp t410; +fp t366; +fp t394; +fp t372; +fp t432; +fp t371; +fp t431; +fp t391; +fp t393; +fp t430; +fp t429; +fp t428; +fp t306; +fp t427; +fp t426; +fp t316; +fp t425; +fp t424; +fp t423; +fp t422; +fp t421; +fp t392; +fp t386; +fp t385; +fp t382; + t411 = x*x; + t412 = y*y; + t464 = t411+t412; + t414 = t411*x; + t404 = RATIONAL(-1.0,60.0); + t356 = t404*t414; + t413 = y*t412; + t360 = t404*t413; + t433 = y*t356+x*t360; + t442 = x*y; + t463 = RATIONAL(11.0,150.0)*t442+t433; + t407 = RATIONAL(1.0,60.0); + t365 = t407*t414; + t364 = t407*t413; + t436 = x*t364+y*t365; + t462 = RATIONAL(-11.0,150.0)*t442+t436; + t435 = t411*t412; + t461 = RATIONAL(-99.0,1225.0)+RATIONAL(-1.0,98.0)*t435; + t438 = y*t411; + t389 = RATIONAL(31.0,210.0); + t398 = RATIONAL(-1.0,140.0); + t406 = RATIONAL(-1.0,30.0); + t362 = t406*t413; + t460 = t389*y+t398*t438+t362; + t401 = RATIONAL(1.0,30.0); + t368 = t401*t413; + t395 = RATIONAL(-31.0,210.0); + t397 = RATIONAL(1.0,140.0); + t459 = t368+t395*y+t397*t438; + t439 = x*t412; + t359 = t401*t414; + t458 = t395*x+t397*t439+t359; + t370 = t406*t414; + t291 = RATIONAL(21.0,200.0)*t442; + t396 = RATIONAL(1.0,120.0); + t440 = x*t413; + t457 = t291+t396*t440+y*t370; + t405 = RATIONAL(-1.0,70.0); + t322 = t405*t438; + t409 = RATIONAL(1.0,84.0); + t456 = t322+t409*y+t364; + t402 = RATIONAL(-1.0,84.0); + t400 = RATIONAL(1.0,70.0); + t327 = t400*t439; + t455 = t327+t356+t402*x; + t383 = RATIONAL(-11.0,105.0); + t454 = t368+t322+t383*y; + t311 = t405*t439; + t453 = t409*x+t311+t365; + t390 = RATIONAL(11.0,105.0); + t307 = t400*t438; + t452 = t307+t390*y+t362; + t451 = t307+t402*y+t360; + t384 = RATIONAL(31.0,420.0); + t408 = RATIONAL(-1.0,35.0); + t318 = t408*t439; + t450 = t384*x+t356+t318; + t399 = RATIONAL(-1.0,120.0); + t296 = RATIONAL(-21.0,200.0)*t442; + t449 = y*t359+t399*t440+t296; + t448 = t389*x+t370+t398*t439; + t387 = RATIONAL(-31.0,420.0); + t403 = RATIONAL(1.0,35.0); + t330 = t403*t439; + t447 = t330+t387*x+t365; + t416 = t412*t412; + t354 = t396*t416; + t415 = t411*t411; + t355 = t396*t415; + t314 = RATIONAL(1.0,49.0)*t435; + t446 = RATIONAL(51.0,1225.0)+t314+t354+t355+t464*RATIONAL(-289.0,5880.0); + t369 = t406*t416; + t358 = t406*t415; + t445 = RATIONAL(1.0,196.0)*t435+t369+RATIONAL(-24.0,1225.0)+t358+t464* +RATIONAL(181.0,1470.0); + t376 = RATIONAL(71.0,5880.0); + t379 = RATIONAL(113.0,735.0); + t444 = t376*t412+t379*t411+t354+t358+t461; + t443 = t369+t376*t411+t379*t412+t355+t461; + t441 = t414*y; + t437 = RATIONAL(-1.0,49.0)*t435+RATIONAL(96.0,1225.0); + t434 = RATIONAL(1.0,98.0)*t435+RATIONAL(246.0,1225.0); + t410 = RATIONAL(1.0,20.0); + t366 = t410*t416; + t394 = RATIONAL(-41.0,196.0); + t372 = RATIONAL(191.0,5880.0); + t432 = t366+t372*t411+t355+t394*t412+t437; + t371 = t410*t415; + t431 = t371+t394*t411+t372*t412+t354+t437; + t391 = RATIONAL(83.0,735.0); + t393 = RATIONAL(-53.0,196.0); + t430 = t369+t391*t412+t371+t393*t411+t434; + t429 = t366+t393*t412+t391*t411+t358+t434; + t428 = RATIONAL(-37.0,300.0)*t442+t436+t445; + t306 = t403*t438; + t427 = t306+t387*y+t364+t446; + t426 = RATIONAL(37.0,300.0)*t442+t433+t445; + t316 = t408*t438; + t425 = t384*y+t360+t316+t446; + t424 = t390*x+t370+t327+t444; + t423 = t383*x+t311+t359+t444; + t422 = x*t368+t296+t399*t441+t443; + t421 = t396*t441+t291+x*t362+t443; + t392 = RATIONAL(-17.0,105.0); + t386 = RATIONAL(17.0,105.0); + t385 = RATIONAL(17.0,420.0); + t382 = RATIONAL(-17.0,420.0); + coeffs_I->coeff_m2_m2 = t425+t450+t462; + coeffs_I->coeff_m1_m2 = t423+t451+t457; + coeffs_I->coeff_0_m2 = t306+t382*y+t360+t431; + coeffs_I->coeff_p1_m2 = t424+t449+t451; + coeffs_I->coeff_p2_m2 = t425+t447+t463; + coeffs_I->coeff_m2_m1 = t421+t454+t455; + coeffs_I->coeff_m1_m1 = t426+t458+t459; + coeffs_I->coeff_0_m1 = t392*y+t368+t307+t430; + coeffs_I->coeff_p1_m1 = t428+t448+t459; + coeffs_I->coeff_p2_m1 = t422+t453+t454; + coeffs_I->coeff_m2_0 = t356+t330+t382*x+t432; + coeffs_I->coeff_m1_0 = t359+t392*x+t327+t429; + coeffs_I->coeff_0_0 = RATIONAL(541.0,1225.0)+t366+t314+t371+t464*RATIONAL +(-57.0,196.0); + coeffs_I->coeff_p1_0 = t311+t370+t386*x+t429; + coeffs_I->coeff_p2_0 = t318+t365+t385*x+t432; + coeffs_I->coeff_m2_p1 = t422+t452+t455; + coeffs_I->coeff_m1_p1 = t428+t458+t460; + coeffs_I->coeff_0_p1 = t386*y+t362+t322+t430; + coeffs_I->coeff_p1_p1 = t426+t448+t460; + coeffs_I->coeff_p2_p1 = t421+t452+t453; + coeffs_I->coeff_m2_p2 = t427+t450+t463; + coeffs_I->coeff_m1_p2 = t423+t449+t456; + coeffs_I->coeff_0_p2 = t316+t364+t385*y+t431; + coeffs_I->coeff_p1_p2 = t424+t456+t457; + coeffs_I->coeff_p2_p2 = t427+t447+t462; diff --git a/src/GeneralizedPolynomial-Uniform/Lagrange/2d.coeffs/2d.cube.order4.smooth0/coeffs-dx.compute.c b/src/GeneralizedPolynomial-Uniform/Lagrange/2d.coeffs/2d.cube.order4.smooth0/coeffs-dx.compute.c new file mode 100644 index 0000000..295e9c5 --- /dev/null +++ b/src/GeneralizedPolynomial-Uniform/Lagrange/2d.coeffs/2d.cube.order4.smooth0/coeffs-dx.compute.c @@ -0,0 +1,155 @@ +fp t532; +fp t503; +fp t533; +fp t534; +fp t551; +fp t569; +fp t511; +fp t550; +fp t568; +fp t557; +fp t567; +fp t522; +fp t535; +fp t507; +fp t566; +fp t481; +fp t565; +fp t526; +fp t512; +fp t564; +fp t505; +fp t491; +fp t563; +fp t506; +fp t490; +fp t562; +fp t556; +fp t561; +fp t524; +fp t501; +fp t560; +fp t559; +fp t558; +fp t508; +fp t555; +fp t554; +fp t482; +fp t525; +fp t509; +fp t553; +fp t521; +fp t504; +fp t552; +fp t475; +fp t549; +fp t467; +fp t548; +fp t514; +fp t547; +fp t546; +fp t545; +fp t544; +fp t486; +fp t543; +fp t542; +fp t541; +fp t540; +fp t539; +fp t538; +fp t483; +fp t479; +fp t473; +fp t469; +fp t466; +fp t465; + t532 = x*x; + t503 = RATIONAL(1.0,20.0)*t532; + t533 = y*y; + t534 = y*t533; + t551 = y*t503+RATIONAL(1.0,60.0)*t534; + t569 = t551+RATIONAL(-37.0,300.0)*y; + t511 = RATIONAL(-1.0,20.0)*t532; + t550 = RATIONAL(-1.0,60.0)*t534+y*t511; + t568 = t550+RATIONAL(37.0,300.0)*y; + t557 = x*t533; + t567 = RATIONAL(1.0,98.0)*t557+RATIONAL(181.0,735.0)*x; + t522 = RATIONAL(1.0,30.0); + t535 = x*t532; + t507 = t522*t535; + t566 = RATIONAL(71.0,2940.0)*x+t507; + t481 = RATIONAL(2.0,49.0)*t557; + t565 = RATIONAL(-289.0,2940.0)*x+t481+t507; + t526 = RATIONAL(-1.0,35.0); + t512 = t526*t533; + t564 = t511+t512+RATIONAL(31.0,420.0); + t505 = RATIONAL(-1.0,10.0)*t532; + t491 = RATIONAL(21.0,200.0)*y; + t563 = y*t505+t491+RATIONAL(1.0,120.0)*t534; + t506 = RATIONAL(1.0,10.0)*t532; + t490 = RATIONAL(-21.0,200.0)*y; + t562 = t490+RATIONAL(-1.0,120.0)*t534+y*t506; + t556 = y*t532; + t561 = RATIONAL(1.0,40.0)*t556+t491+RATIONAL(-1.0,30.0)*t534; + t524 = RATIONAL(1.0,35.0); + t501 = t524*t533; + t560 = t503+RATIONAL(-31.0,420.0)+t501; + t559 = RATIONAL(-1.0,40.0)*t556+t490+t522*t534; + t558 = x*y; + t508 = RATIONAL(-2.0,15.0)*t535; + t555 = t506+t508; + t554 = t505+t508; + t482 = RATIONAL(-1.0,49.0)*t557; + t525 = RATIONAL(-1.0,70.0); + t509 = t525*t533; + t553 = t482+t509; + t521 = RATIONAL(1.0,70.0); + t504 = t521*t533; + t552 = t504+t482; + t475 = RATIONAL(-2.0,49.0)*t557; + t549 = t475+RATIONAL(191.0,2940.0)*x+t507; + t467 = RATIONAL(1.0,49.0)*t557; + t548 = t508+RATIONAL(166.0,735.0)*x+t467; + t514 = RATIONAL(1.0,5.0)*t535; + t547 = t475+t514+RATIONAL(-41.0,98.0)*x; + t546 = t467+RATIONAL(-53.0,98.0)*x+t514; + t545 = t511+RATIONAL(-1.0,84.0)+t552+t566; + t544 = t503+RATIONAL(1.0,84.0)+t553+t566; + t486 = RATIONAL(226.0,735.0)*x; + t543 = t486+RATIONAL(-11.0,105.0)+t553+t555; + t542 = RATIONAL(1.0,140.0)*t533+RATIONAL(-31.0,210.0)+t555+t567; + t541 = t486+RATIONAL(11.0,105.0)+t552+t554; + t540 = RATIONAL(-1.0,140.0)*t533+RATIONAL(31.0,210.0)+t554+t567; + t539 = RATIONAL(-11.0,150.0)*y+t551+t565; + t538 = RATIONAL(11.0,150.0)*y+t550+t565; + t483 = t524*t558; + t479 = RATIONAL(-2.0,35.0)*t558; + t473 = RATIONAL(2.0,35.0)*t558; + t469 = t526*t558; + t466 = t521*t558; + t465 = t525*t558; + coeffs_dx->coeff_m2_m2 = t479+t539+t564; + coeffs_dx->coeff_m1_m2 = t483+t543+t563; + coeffs_dx->coeff_0_m2 = t473+t547; + coeffs_dx->coeff_p1_m2 = t483+t541+t562; + coeffs_dx->coeff_p2_m2 = t479+t538+t560; + coeffs_dx->coeff_m2_m1 = t469+t545+t561; + coeffs_dx->coeff_m1_m1 = t466+t542+t568; + coeffs_dx->coeff_0_m1 = t483+t546; + coeffs_dx->coeff_p1_m1 = t466+t540+t569; + coeffs_dx->coeff_p2_m1 = t469+t544+t559; + coeffs_dx->coeff_m2_0 = t501+t511+RATIONAL(-17.0,420.0)+t549; + coeffs_dx->coeff_m1_0 = t504+RATIONAL(-17.0,105.0)+t506+t548; + coeffs_dx->coeff_0_0 = t514+RATIONAL(-57.0,98.0)*x+t481; + coeffs_dx->coeff_p1_0 = t509+RATIONAL(17.0,105.0)+t505+t548; + coeffs_dx->coeff_p2_0 = RATIONAL(17.0,420.0)+t512+t503+t549; + coeffs_dx->coeff_m2_p1 = t483+t545+t559; + coeffs_dx->coeff_m1_p1 = t465+t542+t569; + coeffs_dx->coeff_0_p1 = t469+t546; + coeffs_dx->coeff_p1_p1 = t465+t540+t568; + coeffs_dx->coeff_p2_p1 = t483+t544+t561; + coeffs_dx->coeff_m2_p2 = t473+t538+t564; + coeffs_dx->coeff_m1_p2 = t469+t543+t562; + coeffs_dx->coeff_0_p2 = t479+t547; + coeffs_dx->coeff_p1_p2 = t469+t541+t563; + coeffs_dx->coeff_p2_p2 = t473+t539+t560; diff --git a/src/GeneralizedPolynomial-Uniform/Lagrange/2d.coeffs/2d.cube.order4.smooth0/coeffs-dxx.compute.c b/src/GeneralizedPolynomial-Uniform/Lagrange/2d.coeffs/2d.cube.order4.smooth0/coeffs-dxx.compute.c new file mode 100644 index 0000000..c1b3b65 --- /dev/null +++ b/src/GeneralizedPolynomial-Uniform/Lagrange/2d.coeffs/2d.cube.order4.smooth0/coeffs-dxx.compute.c @@ -0,0 +1,107 @@ +fp t704; +fp t708; +fp t685; +fp t709; +fp t688; +fp t726; +fp t683; +fp t695; +fp t725; +fp t724; +fp t723; +fp t722; +fp t679; +fp t692; +fp t694; +fp t687; +fp t681; +fp t721; +fp t693; +fp t720; +fp t686; +fp t719; +fp t718; +fp t682; +fp t717; +fp t678; +fp t716; +fp t684; +fp t715; +fp t714; +fp t713; +fp t712; +fp t689; +fp t711; +fp t710; +fp t677; +fp t676; +fp t675; +fp t674; +fp t673; +fp t672; + t704 = RATIONAL(1.0,10.0); + t708 = x*x; + t685 = t704*t708; + t709 = y*y; + t688 = RATIONAL(2.0,49.0)*t709; + t726 = t685+t688+RATIONAL(-289.0,2940.0); + t683 = RATIONAL(-1.0,49.0)*t709; + t695 = RATIONAL(-2.0,5.0)*t708; + t725 = t683+t695+RATIONAL(226.0,735.0); + t724 = RATIONAL(1.0,98.0)*t709+t695+RATIONAL(181.0,735.0); + t723 = t683+t685+RATIONAL(71.0,2940.0); + t722 = x*y; + t679 = t704*x; + t692 = RATIONAL(-1.0,5.0)*x; + t694 = RATIONAL(1.0,5.0)*x; + t687 = RATIONAL(-1.0,10.0)*x; + t681 = RATIONAL(-2.0,49.0)*t709; + t721 = t681+t685+RATIONAL(191.0,2940.0); + t693 = RATIONAL(3.0,5.0)*t708; + t720 = t681+t693+RATIONAL(-41.0,98.0); + t686 = RATIONAL(1.0,49.0)*t709; + t719 = t686+t695+RATIONAL(166.0,735.0); + t718 = t686+t693+RATIONAL(-53.0,98.0); + t682 = RATIONAL(2.0,35.0)*y; + t717 = t682+t726; + t678 = RATIONAL(-1.0,35.0)*y; + t716 = t678+t725; + t684 = RATIONAL(1.0,35.0)*y; + t715 = t684+t725; + t714 = t679+t723; + t713 = t687+t723; + t712 = RATIONAL(-1.0,70.0)*y+t724; + t689 = RATIONAL(-2.0,35.0)*y; + t711 = t689+t726; + t710 = RATIONAL(1.0,70.0)*y+t724; + t677 = y*t692; + t676 = y*t694; + t675 = RATIONAL(1.0,20.0)*t722; + t674 = y*t687; + t673 = RATIONAL(-1.0,20.0)*t722; + t672 = y*t679; + coeffs_dxx->coeff_m2_m2 = t672+t687+t711; + coeffs_dxx->coeff_m1_m2 = t677+t694+t715; + coeffs_dxx->coeff_0_m2 = t682+t720; + coeffs_dxx->coeff_p1_m2 = t676+t692+t715; + coeffs_dxx->coeff_p2_m2 = t674+t679+t711; + coeffs_dxx->coeff_m2_m1 = t675+t678+t713; + coeffs_dxx->coeff_m1_m1 = t694+t674+t710; + coeffs_dxx->coeff_0_m1 = t684+t718; + coeffs_dxx->coeff_p1_m1 = t692+t672+t710; + coeffs_dxx->coeff_p2_m1 = t678+t673+t714; + coeffs_dxx->coeff_m2_0 = t687+t721; + coeffs_dxx->coeff_m1_0 = t694+t719; + coeffs_dxx->coeff_0_0 = t693+RATIONAL(-57.0,98.0)+t688; + coeffs_dxx->coeff_p1_0 = t692+t719; + coeffs_dxx->coeff_p2_0 = t679+t721; + coeffs_dxx->coeff_m2_p1 = t684+t673+t713; + coeffs_dxx->coeff_m1_p1 = t672+t694+t712; + coeffs_dxx->coeff_0_p1 = t678+t718; + coeffs_dxx->coeff_p1_p1 = t692+t674+t712; + coeffs_dxx->coeff_p2_p1 = t684+t675+t714; + coeffs_dxx->coeff_m2_p2 = t687+t674+t717; + coeffs_dxx->coeff_m1_p2 = t694+t676+t716; + coeffs_dxx->coeff_0_p2 = t689+t720; + coeffs_dxx->coeff_p1_p2 = t677+t692+t716; + coeffs_dxx->coeff_p2_p2 = t679+t672+t717; diff --git a/src/GeneralizedPolynomial-Uniform/Lagrange/2d.coeffs/2d.cube.order4.smooth0/coeffs-dxy.compute.c b/src/GeneralizedPolynomial-Uniform/Lagrange/2d.coeffs/2d.cube.order4.smooth0/coeffs-dxy.compute.c new file mode 100644 index 0000000..dd9a747 --- /dev/null +++ b/src/GeneralizedPolynomial-Uniform/Lagrange/2d.coeffs/2d.cube.order4.smooth0/coeffs-dxy.compute.c @@ -0,0 +1,119 @@ +fp t774; +fp t775; +fp t793; +fp t763; +fp t747; +fp t765; +fp t732; +fp t792; +fp t740; +fp t753; +fp t791; +fp t767; +fp t771; +fp t790; +fp t789; +fp t788; +fp t787; +fp t786; +fp t731; +fp t785; +fp t784; +fp t728; +fp t783; +fp t730; +fp t782; +fp t768; +fp t773; +fp t781; +fp t780; +fp t779; +fp t778; +fp t777; +fp t776; +fp t770; +fp t766; +fp t764; +fp t762; +fp t755; +fp t754; +fp t751; +fp t746; +fp t745; +fp t742; +fp t736; +fp t733; +fp t729; +fp t727; + t774 = x*x; + t775 = y*y; + t793 = t774+t775; + t763 = RATIONAL(1.0,35.0); + t747 = t763*x; + t765 = RATIONAL(-1.0,35.0); + t732 = t765*y; + t792 = t747+t732; + t740 = t765*x; + t753 = t763*y; + t791 = t740+t753; + t767 = RATIONAL(1.0,10.0); + t771 = RATIONAL(-1.0,40.0); + t790 = t767*t775+t771*t774; + t789 = t767*t774+t771*t775; + t788 = x*y; + t787 = t793*RATIONAL(-1.0,20.0); + t786 = t793*RATIONAL(1.0,20.0); + t731 = RATIONAL(-2.0,49.0)*t788; + t785 = t731+RATIONAL(21.0,200.0); + t784 = t731+RATIONAL(-21.0,200.0); + t728 = RATIONAL(1.0,49.0)*t788; + t783 = t728+RATIONAL(-37.0,300.0)+t786; + t730 = RATIONAL(4.0,49.0)*t788; + t782 = t730+RATIONAL(-11.0,150.0)+t786; + t768 = RATIONAL(-1.0,10.0); + t773 = RATIONAL(1.0,40.0); + t781 = t773*t775+t768*t774+t785; + t780 = t768*t775+t773*t774+t785; + t779 = t740+t732+t784; + t778 = t730+RATIONAL(11.0,150.0)+t787; + t777 = t747+t753+t784; + t776 = t728+RATIONAL(37.0,300.0)+t787; + t770 = RATIONAL(-2.0,35.0); + t766 = RATIONAL(2.0,35.0); + t764 = RATIONAL(-1.0,70.0); + t762 = RATIONAL(1.0,70.0); + t755 = t762*y; + t754 = t770*y; + t751 = t770*x; + t746 = t766*y; + t745 = t766*x; + t742 = t764*y; + t736 = t762*x; + t733 = t764*x; + t729 = RATIONAL(-4.0,49.0)*t788; + t727 = RATIONAL(2.0,49.0)*t788; + coeffs_dxy->coeff_m2_m2 = t754+t751+t782; + coeffs_dxy->coeff_m1_m2 = t781+t792; + coeffs_dxy->coeff_0_m2 = t729+t745; + coeffs_dxy->coeff_p1_m2 = t777+t789; + coeffs_dxy->coeff_p2_m2 = t746+t751+t778; + coeffs_dxy->coeff_m2_m1 = t780+t791; + coeffs_dxy->coeff_m1_m1 = t755+t736+t776; + coeffs_dxy->coeff_0_m1 = t727+t747; + coeffs_dxy->coeff_p1_m1 = t742+t736+t783; + coeffs_dxy->coeff_p2_m1 = t779+t790; + coeffs_dxy->coeff_m2_0 = t746+t729; + coeffs_dxy->coeff_m1_0 = t727+t753; + coeffs_dxy->coeff_0_0 = t730; + coeffs_dxy->coeff_p1_0 = t732+t727; + coeffs_dxy->coeff_p2_0 = t754+t729; + coeffs_dxy->coeff_m2_p1 = t777+t790; + coeffs_dxy->coeff_m1_p1 = t733+t755+t783; + coeffs_dxy->coeff_0_p1 = t727+t740; + coeffs_dxy->coeff_p1_p1 = t733+t742+t776; + coeffs_dxy->coeff_p2_p1 = t780+t792; + coeffs_dxy->coeff_m2_p2 = t745+t754+t778; + coeffs_dxy->coeff_m1_p2 = t779+t789; + coeffs_dxy->coeff_0_p2 = t729+t751; + coeffs_dxy->coeff_p1_p2 = t781+t791; + coeffs_dxy->coeff_p2_p2 = t745+t746+t782; diff --git a/src/GeneralizedPolynomial-Uniform/Lagrange/2d.coeffs/2d.cube.order4.smooth0/coeffs-dy.compute.c b/src/GeneralizedPolynomial-Uniform/Lagrange/2d.coeffs/2d.cube.order4.smooth0/coeffs-dy.compute.c new file mode 100644 index 0000000..a1b62b5 --- /dev/null +++ b/src/GeneralizedPolynomial-Uniform/Lagrange/2d.coeffs/2d.cube.order4.smooth0/coeffs-dy.compute.c @@ -0,0 +1,147 @@ +fp t638; +fp t617; +fp t637; +fp t640; +fp t655; +fp t671; +fp t613; +fp t656; +fp t670; +fp t658; +fp t586; +fp t627; +fp t639; +fp t614; +fp t669; +fp t587; +fp t608; +fp t668; +fp t667; +fp t666; +fp t596; +fp t657; +fp t665; +fp t631; +fp t618; +fp t664; +fp t629; +fp t611; +fp t663; +fp t597; +fp t662; +fp t607; +fp t630; +fp t606; +fp t661; +fp t626; +fp t612; +fp t605; +fp t660; +fp t659; +fp t572; +fp t654; +fp t581; +fp t653; +fp t619; +fp t652; +fp t651; +fp t650; +fp t649; +fp t648; +fp t647; +fp t646; +fp t645; +fp t644; +fp t643; +fp t584; +fp t582; +fp t577; +fp t576; +fp t573; +fp t571; + t638 = y*y; + t617 = RATIONAL(-1.0,20.0)*t638; + t637 = x*x; + t640 = x*t637; + t655 = x*t617+RATIONAL(-1.0,60.0)*t640; + t671 = t655+RATIONAL(37.0,300.0)*x; + t613 = RATIONAL(1.0,20.0)*t638; + t656 = x*t613+RATIONAL(1.0,60.0)*t640; + t670 = RATIONAL(-37.0,300.0)*x+t656; + t658 = y*t637; + t586 = RATIONAL(2.0,49.0)*t658; + t627 = RATIONAL(1.0,30.0); + t639 = y*t638; + t614 = t627*t639; + t669 = RATIONAL(-289.0,2940.0)*y+t586+t614; + t587 = RATIONAL(-1.0,49.0)*t658; + t608 = RATIONAL(-2.0,15.0)*t639; + t668 = t587+RATIONAL(226.0,735.0)*y+t608; + t667 = t587+t614+RATIONAL(71.0,2940.0)*y; + t666 = RATIONAL(1.0,98.0)*t658+RATIONAL(181.0,735.0)*y+t608; + t596 = RATIONAL(-21.0,200.0)*x; + t657 = x*t638; + t665 = t627*t640+RATIONAL(-1.0,40.0)*t657+t596; + t631 = RATIONAL(-1.0,35.0); + t618 = t631*t637; + t664 = t617+t618+RATIONAL(31.0,420.0); + t629 = RATIONAL(1.0,35.0); + t611 = t629*t637; + t663 = t611+t613+RATIONAL(-31.0,420.0); + t597 = RATIONAL(21.0,200.0)*x; + t662 = RATIONAL(1.0,40.0)*t657+RATIONAL(-1.0,30.0)*t640+t597; + t607 = RATIONAL(1.0,10.0)*t638; + t630 = RATIONAL(-1.0,70.0); + t606 = t630*t637; + t661 = RATIONAL(-11.0,105.0)+t607+t606; + t626 = RATIONAL(1.0,70.0); + t612 = t626*t637; + t605 = RATIONAL(-1.0,10.0)*t638; + t660 = t612+t605+RATIONAL(11.0,105.0); + t659 = x*y; + t572 = RATIONAL(1.0,49.0)*t658; + t654 = t572+RATIONAL(166.0,735.0)*y+t608; + t581 = RATIONAL(-2.0,49.0)*t658; + t653 = RATIONAL(191.0,2940.0)*y+t581+t614; + t619 = RATIONAL(1.0,5.0)*t639; + t652 = t581+t619+RATIONAL(-41.0,98.0)*y; + t651 = t572+t619+RATIONAL(-53.0,98.0)*y; + t650 = RATIONAL(-11.0,150.0)*x+t656+t669; + t649 = RATIONAL(1.0,140.0)*t637+t607+RATIONAL(-31.0,210.0)+t666; + t648 = x*t607+t596+RATIONAL(-1.0,120.0)*t640+t668; + t647 = x*t605+t597+RATIONAL(1.0,120.0)*t640+t668; + t646 = RATIONAL(-1.0,140.0)*t637+t605+RATIONAL(31.0,210.0)+t666; + t645 = t612+t617+RATIONAL(-1.0,84.0)+t667; + t644 = t613+t606+RATIONAL(1.0,84.0)+t667; + t643 = RATIONAL(11.0,150.0)*x+t655+t669; + t584 = RATIONAL(-2.0,35.0)*t659; + t582 = t631*t659; + t577 = t629*t659; + t576 = RATIONAL(2.0,35.0)*t659; + t573 = t626*t659; + t571 = t630*t659; + coeffs_dy->coeff_m2_m2 = t584+t650+t664; + coeffs_dy->coeff_m1_m2 = t582+t645+t662; + coeffs_dy->coeff_0_m2 = t611+t617+RATIONAL(-17.0,420.0)+t653; + coeffs_dy->coeff_p1_m2 = t577+t645+t665; + coeffs_dy->coeff_p2_m2 = t576+t643+t664; + coeffs_dy->coeff_m2_m1 = t577+t647+t661; + coeffs_dy->coeff_m1_m1 = t573+t649+t671; + coeffs_dy->coeff_0_m1 = RATIONAL(-17.0,105.0)+t607+t612+t654; + coeffs_dy->coeff_p1_m1 = t571+t649+t670; + coeffs_dy->coeff_p2_m1 = t582+t648+t661; + coeffs_dy->coeff_m2_0 = t576+t652; + coeffs_dy->coeff_m1_0 = t577+t651; + coeffs_dy->coeff_0_0 = RATIONAL(-57.0,98.0)*y+t619+t586; + coeffs_dy->coeff_p1_0 = t582+t651; + coeffs_dy->coeff_p2_0 = t584+t652; + coeffs_dy->coeff_m2_p1 = t577+t648+t660; + coeffs_dy->coeff_m1_p1 = t573+t646+t670; + coeffs_dy->coeff_0_p1 = RATIONAL(17.0,105.0)+t606+t605+t654; + coeffs_dy->coeff_p1_p1 = t571+t646+t671; + coeffs_dy->coeff_p2_p1 = t582+t647+t660; + coeffs_dy->coeff_m2_p2 = t584+t643+t663; + coeffs_dy->coeff_m1_p2 = t582+t644+t665; + coeffs_dy->coeff_0_p2 = t618+RATIONAL(17.0,420.0)+t613+t653; + coeffs_dy->coeff_p1_p2 = t577+t644+t662; + coeffs_dy->coeff_p2_p2 = t576+t650+t663; diff --git a/src/GeneralizedPolynomial-Uniform/Lagrange/2d.coeffs/2d.cube.order4.smooth0/coeffs-dyy.compute.c b/src/GeneralizedPolynomial-Uniform/Lagrange/2d.coeffs/2d.cube.order4.smooth0/coeffs-dyy.compute.c new file mode 100644 index 0000000..29eef2a --- /dev/null +++ b/src/GeneralizedPolynomial-Uniform/Lagrange/2d.coeffs/2d.cube.order4.smooth0/coeffs-dyy.compute.c @@ -0,0 +1,107 @@ +fp t831; +fp t817; +fp t830; +fp t848; +fp t808; +fp t826; +fp t811; +fp t847; +fp t846; +fp t813; +fp t845; +fp t844; +fp t814; +fp t816; +fp t800; +fp t812; +fp t807; +fp t843; +fp t815; +fp t842; +fp t801; +fp t841; +fp t840; +fp t839; +fp t838; +fp t806; +fp t837; +fp t809; +fp t836; +fp t835; +fp t834; +fp t833; +fp t832; +fp t810; +fp t805; +fp t799; +fp t798; +fp t797; +fp t796; +fp t795; +fp t794; + t831 = y*y; + t817 = RATIONAL(-2.0,5.0)*t831; + t830 = x*x; + t848 = RATIONAL(1.0,98.0)*t830+t817+RATIONAL(181.0,735.0); + t808 = RATIONAL(-1.0,49.0)*t830; + t826 = RATIONAL(1.0,10.0); + t811 = t826*t831; + t847 = t808+t811+RATIONAL(71.0,2940.0); + t846 = t808+t817+RATIONAL(226.0,735.0); + t813 = RATIONAL(2.0,49.0)*t830; + t845 = t811+t813+RATIONAL(-289.0,2940.0); + t844 = x*y; + t814 = RATIONAL(-1.0,5.0)*y; + t816 = RATIONAL(1.0,5.0)*y; + t800 = t826*y; + t812 = RATIONAL(-1.0,10.0)*y; + t807 = RATIONAL(-2.0,49.0)*t830; + t843 = t807+t811+RATIONAL(191.0,2940.0); + t815 = RATIONAL(3.0,5.0)*t831; + t842 = t807+t815+RATIONAL(-41.0,98.0); + t801 = RATIONAL(1.0,49.0)*t830; + t841 = t801+t815+RATIONAL(-53.0,98.0); + t840 = t801+RATIONAL(166.0,735.0)+t817; + t839 = t800+t845; + t838 = t812+t845; + t806 = RATIONAL(1.0,35.0)*x; + t837 = t806+t847; + t809 = RATIONAL(-1.0,35.0)*x; + t836 = t809+t847; + t835 = t814+t846; + t834 = t816+t846; + t833 = RATIONAL(1.0,70.0)*x+t848; + t832 = RATIONAL(-1.0,70.0)*x+t848; + t810 = RATIONAL(-2.0,35.0)*x; + t805 = RATIONAL(2.0,35.0)*x; + t799 = x*t814; + t798 = x*t816; + t797 = RATIONAL(1.0,20.0)*t844; + t796 = x*t812; + t795 = RATIONAL(-1.0,20.0)*t844; + t794 = x*t800; + coeffs_dyy->coeff_m2_m2 = t810+t794+t838; + coeffs_dyy->coeff_m1_m2 = t797+t812+t836; + coeffs_dyy->coeff_0_m2 = t812+t843; + coeffs_dyy->coeff_p1_m2 = t812+t795+t837; + coeffs_dyy->coeff_p2_m2 = t796+t805+t838; + coeffs_dyy->coeff_m2_m1 = t799+t806+t834; + coeffs_dyy->coeff_m1_m1 = t796+t816+t833; + coeffs_dyy->coeff_0_m1 = t816+t840; + coeffs_dyy->coeff_p1_m1 = t794+t816+t832; + coeffs_dyy->coeff_p2_m1 = t809+t798+t834; + coeffs_dyy->coeff_m2_0 = t805+t842; + coeffs_dyy->coeff_m1_0 = t806+t841; + coeffs_dyy->coeff_0_0 = t813+t815+RATIONAL(-57.0,98.0); + coeffs_dyy->coeff_p1_0 = t809+t841; + coeffs_dyy->coeff_p2_0 = t810+t842; + coeffs_dyy->coeff_m2_p1 = t806+t798+t835; + coeffs_dyy->coeff_m1_p1 = t794+t814+t833; + coeffs_dyy->coeff_0_p1 = t814+t840; + coeffs_dyy->coeff_p1_p1 = t814+t796+t832; + coeffs_dyy->coeff_p2_p1 = t799+t809+t835; + coeffs_dyy->coeff_m2_p2 = t796+t810+t839; + coeffs_dyy->coeff_m1_p2 = t795+t800+t836; + coeffs_dyy->coeff_0_p2 = t800+t843; + coeffs_dyy->coeff_p1_p2 = t800+t797+t837; + coeffs_dyy->coeff_p2_p2 = t794+t805+t839; -- cgit v1.2.3