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

5 files changed, 80 insertions, 0 deletions

diff --git a/src/GeneralizedPolynomial-Uniform/common/2d.cube.size4/coeffs-dcl.h b/src/GeneralizedPolynomial-Uniform/common/2d.cube.size4/coeffs-dcl.h
new file mode 100644
index 0000000..beb0825
--- /dev/null
+++ b/src/GeneralizedPolynomial-Uniform/common/2d.cube.size4/coeffs-dcl.h
@@ -0,0 +1,16 @@
+fp coeff_m1_m1;
+fp coeff_0_m1;
+fp coeff_p1_m1;
+fp coeff_p2_m1;
+fp coeff_m1_0;
+fp coeff_0_0;
+fp coeff_p1_0;
+fp coeff_p2_0;
+fp coeff_m1_p1;
+fp coeff_0_p1;
+fp coeff_p1_p1;
+fp coeff_p2_p1;
+fp coeff_m1_p2;
+fp coeff_0_p2;
+fp coeff_p1_p2;
+fp coeff_p2_p2;
diff --git a/src/GeneralizedPolynomial-Uniform/common/2d.cube.size4/data-dcl.h b/src/GeneralizedPolynomial-Uniform/common/2d.cube.size4/data-dcl.h
new file mode 100644
index 0000000..a55365f
--- /dev/null
+++ b/src/GeneralizedPolynomial-Uniform/common/2d.cube.size4/data-dcl.h
@@ -0,0 +1,16 @@
+fp data_m1_m1;
+fp data_0_m1;
+fp data_p1_m1;
+fp data_p2_m1;
+fp data_m1_0;
+fp data_0_0;
+fp data_p1_0;
+fp data_p2_0;
+fp data_m1_p1;
+fp data_0_p1;
+fp data_p1_p1;
+fp data_p2_p1;
+fp data_m1_p2;
+fp data_0_p2;
+fp data_p1_p2;
+fp data_p2_p2;
diff --git a/src/GeneralizedPolynomial-Uniform/common/2d.cube.size4/evaluate-molecule.c b/src/GeneralizedPolynomial-Uniform/common/2d.cube.size4/evaluate-molecule.c
new file mode 100644
index 0000000..03e306b
--- /dev/null
+++ b/src/GeneralizedPolynomial-Uniform/common/2d.cube.size4/evaluate-molecule.c
@@ -0,0 +1,16 @@
+    coeffs->coeff_m1_m1*data->data_m1_m1
+  + coeffs->coeff_0_m1*data->data_0_m1
+  + coeffs->coeff_p1_m1*data->data_p1_m1
+  + coeffs->coeff_p2_m1*data->data_p2_m1
+  + coeffs->coeff_m1_0*data->data_m1_0
+  + coeffs->coeff_0_0*data->data_0_0
+  + coeffs->coeff_p1_0*data->data_p1_0
+  + coeffs->coeff_p2_0*data->data_p2_0
+  + coeffs->coeff_m1_p1*data->data_m1_p1
+  + coeffs->coeff_0_p1*data->data_0_p1
+  + coeffs->coeff_p1_p1*data->data_p1_p1
+  + coeffs->coeff_p2_p1*data->data_p2_p1
+  + coeffs->coeff_m1_p2*data->data_m1_p2
+  + coeffs->coeff_0_p2*data->data_0_p2
+  + coeffs->coeff_p1_p2*data->data_p1_p2
+  + coeffs->coeff_p2_p2*data->data_p2_p2;
diff --git a/src/GeneralizedPolynomial-Uniform/common/2d.cube.size4/fetch-data.c b/src/GeneralizedPolynomial-Uniform/common/2d.cube.size4/fetch-data.c
new file mode 100644
index 0000000..ee17aca
--- /dev/null
+++ b/src/GeneralizedPolynomial-Uniform/common/2d.cube.size4/fetch-data.c
@@ -0,0 +1,16 @@
+data->data_m1_m1 = DATA(-1,-1);
+data->data_0_m1 = DATA(0,-1);
+data->data_p1_m1 = DATA(1,-1);
+data->data_p2_m1 = DATA(2,-1);
+data->data_m1_0 = DATA(-1,0);
+data->data_0_0 = DATA(0,0);
+data->data_p1_0 = DATA(1,0);
+data->data_p2_0 = DATA(2,0);
+data->data_m1_p1 = DATA(-1,1);
+data->data_0_p1 = DATA(0,1);
+data->data_p1_p1 = DATA(1,1);
+data->data_p2_p1 = DATA(2,1);
+data->data_m1_p2 = DATA(-1,2);
+data->data_0_p2 = DATA(0,2);
+data->data_p1_p2 = DATA(1,2);
+data->data_p2_p2 = DATA(2,2);
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;