From 15c98b3abee2525b24ca89932256ff2cda63c9bc Mon Sep 17 00:00:00 2001
From: jthorn <jthorn@df1f8a13-aa1d-4dd4-9681-27ded5b42416>
Date: Sun, 18 Aug 2002 15:12:49 +0000
Subject: This commit reorganizes the Maple scripts which generate the
 interpolation coefficients, and the coefficient files themselves, to properly
 support multiple interpolation operators, and adds some skeleton support for
 Hermite interpolation.  I will add full support for Hermite in a future
 checkin.

git-svn-id: http://svn.cactuscode.org/arrangements/CactusBase/LocalInterp/trunk@81 df1f8a13-aa1d-4dd4-9681-27ded5b42416
---
 .../2d.cube.order3.smooth0/coeff-dx.compute.c      | 95 ++++++++++++++++++++++
 1 file changed, 95 insertions(+)
 create mode 100644 src/GeneralizedPolynomial-Uniform/Lagrange/2d.coeffs/2d.cube.order3.smooth0/coeff-dx.compute.c

(limited to 'src/GeneralizedPolynomial-Uniform/Lagrange/2d.coeffs/2d.cube.order3.smooth0/coeff-dx.compute.c')

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
new file mode 100644
index 0000000..4650764
--- /dev/null
+++ b/src/GeneralizedPolynomial-Uniform/Lagrange/2d.coeffs/2d.cube.order3.smooth0/coeff-dx.compute.c
@@ -0,0 +1,95 @@
+fp t201,
+   t164,
+   t200,
+   t163,
+   t199,
+   t183,
+   t188,
+   t173,
+   t187,
+   t179,
+   t198,
+   t161,
+   t197,
+   t162,
+   t196,
+   t185,
+   t176,
+   t195,
+   t181,
+   t194,
+   t193,
+   t186,
+   t177,
+   t182,
+   t192,
+   t180,
+   t191,
+   t184,
+   t172,
+   t190,
+   t189,
+   t178,
+   t175,
+   t170,
+   t168,
+   t166,
+   t165;
+      t201 = x*y;
+      t164 = RATIONAL(-3.0,20.0)*t201;
+      t200 = RATIONAL(13.0,40.0)*x+t164;
+      t163 = RATIONAL(3.0,20.0)*t201;
+      t199 = RATIONAL(7.0,40.0)*x+t163;
+      t183 = RATIONAL(1.0,40.0);
+      t188 = y*y;
+      t173 = t183*t188;
+      t187 = x*x;
+      t179 = RATIONAL(-3.0,8.0)*t187;
+      t198 = t173+t179;
+      t161 = RATIONAL(1.0,20.0)*t201;
+      t197 = RATIONAL(9.0,40.0)*x+t161;
+      t162 = RATIONAL(-1.0,20.0)*t201;
+      t196 = RATIONAL(11.0,40.0)*x+t162;
+      t185 = RATIONAL(-1.0,40.0);
+      t176 = t185*t188;
+      t195 = t176+t179;
+      t181 = RATIONAL(3.0,8.0)*t187;
+      t194 = t173+t181;
+      t193 = t176+t181;
+      t186 = RATIONAL(-3.0,40.0);
+      t177 = t186*t188;
+      t182 = RATIONAL(-1.0,8.0)*t187;
+      t192 = t177+t182;
+      t180 = RATIONAL(1.0,8.0)*t187;
+      t191 = t177+t180;
+      t184 = RATIONAL(3.0,40.0);
+      t172 = t184*t188;
+      t190 = t172+t180;
+      t189 = t172+t182;
+      t178 = RATIONAL(2.0,25.0)*y;
+      t175 = RATIONAL(-1.0,50.0)*y;
+      t170 = RATIONAL(7.0,100.0)*y;
+      t168 = RATIONAL(-9.0,100.0)*y;
+      t166 = RATIONAL(-1.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+t200;
+      coeff_dx_0_m1 = RATIONAL(-31.0,400.0)+RATIONAL(-23.0,40.0)*x+t175+t163+
+t193;
+      coeff_dx_p1_m1 = RATIONAL(111.0,400.0)+t165+t198+t199;
+      coeff_dx_p2_m1 = t168+t164+RATIONAL(-131.0,1200.0)+t184*x+t190;
+      coeff_dx_m1_0 = RATIONAL(-223.0,1200.0)+t175+t189+t196;
+      coeff_dx_0_0 = RATIONAL(-57.0,400.0)+RATIONAL(-21.0,40.0)*x+RATIONAL(-1.0
+,25.0)*y+t161+t194;
+      coeff_dx_p1_0 = RATIONAL(117.0,400.0)+t166+t195+t197;
+      coeff_dx_p2_0 = RATIONAL(43.0,1200.0)+t183*x+t162+t170+t191;
+      coeff_dx_m1_p1 = RATIONAL(-157.0,1200.0)+t165+t189+t197;
+      coeff_dx_0_p1 = RATIONAL(-63.0,400.0)+t166+t162+RATIONAL(-19.0,40.0)*x+
+t194;
+      coeff_dx_p1_p1 = RATIONAL(3.0,50.0)*y+RATIONAL(103.0,400.0)+t195+t196;
+      coeff_dx_p2_p1 = t185*x+t161+RATIONAL(37.0,1200.0)+t178+t191;
+      coeff_dx_m1_p2 = t168+RATIONAL(89.0,1200.0)+t192+t199;
+      coeff_dx_0_p2 = t164+t170+RATIONAL(-17.0,40.0)*x+RATIONAL(-49.0,400.0)+
+t193;
+      coeff_dx_p1_p2 = RATIONAL(69.0,400.0)+t178+t198+t200;
+      coeff_dx_p2_p2 = t163+RATIONAL(-149.0,1200.0)+t186*x+RATIONAL(-3.0,50.0)*
+y+t190;
-- 
cgit v1.2.3