reorganize common to Lagrange and Hermite interpolators

(i.e. code depending only on molecule size/shape)
--> now there are separate *functions* rather than just code fragments
--> ../template.c should hopefully compile in finite time.....


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

1 files changed, 16 insertions, 0 deletions

diff --git a/src/GeneralizedPolynomial-Uniform/common/2d.cube.size4/store-coeffs.c b/src/GeneralizedPolynomial-Uniform/common/2d.cube.size4/store-coeffs.c
new file mode 100644
index 0000000..804500c
--- /dev/null
+++ b/src/GeneralizedPolynomial-Uniform/common/2d.cube.size4/store-coeffs.c
@@ -0,0 +1,16 @@
+COEFF(-1,-1) = factor * coeffs->coeff_m1_m1;
+COEFF(0,-1) = factor * coeffs->coeff_0_m1;
+COEFF(1,-1) = factor * coeffs->coeff_p1_m1;
+COEFF(2,-1) = factor * coeffs->coeff_p2_m1;
+COEFF(-1,0) = factor * coeffs->coeff_m1_0;
+COEFF(0,0) = factor * coeffs->coeff_0_0;
+COEFF(1,0) = factor * coeffs->coeff_p1_0;
+COEFF(2,0) = factor * coeffs->coeff_p2_0;
+COEFF(-1,1) = factor * coeffs->coeff_m1_p1;
+COEFF(0,1) = factor * coeffs->coeff_0_p1;
+COEFF(1,1) = factor * coeffs->coeff_p1_p1;
+COEFF(2,1) = factor * coeffs->coeff_p2_p1;
+COEFF(-1,2) = factor * coeffs->coeff_m1_p2;
+COEFF(0,2) = factor * coeffs->coeff_0_p2;
+COEFF(1,2) = factor * coeffs->coeff_p1_p2;
+COEFF(2,2) = factor * coeffs->coeff_p2_p2;