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.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 ++++++ 12 files changed, 136 insertions(+), 136 deletions(-) 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 (limited to 'src/GeneralizedPolynomial-Uniform/Lagrange/2d.coeffs/2d.cube.order2.smooth0') 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; -- cgit v1.2.3