initial checkin of new LocalInterp thorn

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

14 files changed, 132 insertions, 0 deletions

diff --git a/src/GeneralizedPolynomial-Uniform/2d.coeffs/2d.cube.size3/coeff-I.dcl.c b/src/GeneralizedPolynomial-Uniform/2d.coeffs/2d.cube.size3/coeff-I.dcl.c
new file mode 100644
index 0000000..3f4b0cb
--- /dev/null
+++ b/src/GeneralizedPolynomial-Uniform/2d.coeffs/2d.cube.size3/coeff-I.dcl.c
@@ -0,0 +1,9 @@
+fp coeff_I_m1_m1,
+   coeff_I_0_m1,
+   coeff_I_p1_m1,
+   coeff_I_m1_0,
+   coeff_I_0_0,
+   coeff_I_p1_0,
+   coeff_I_m1_p1,
+   coeff_I_0_p1,
+   coeff_I_p1_p1;
diff --git a/src/GeneralizedPolynomial-Uniform/2d.coeffs/2d.cube.size3/coeff-dx.dcl.c b/src/GeneralizedPolynomial-Uniform/2d.coeffs/2d.cube.size3/coeff-dx.dcl.c
new file mode 100644
index 0000000..3efed8a
--- /dev/null
+++ b/src/GeneralizedPolynomial-Uniform/2d.coeffs/2d.cube.size3/coeff-dx.dcl.c
@@ -0,0 +1,9 @@
+fp coeff_dx_m1_m1,
+   coeff_dx_0_m1,
+   coeff_dx_p1_m1,
+   coeff_dx_m1_0,
+   coeff_dx_0_0,
+   coeff_dx_p1_0,
+   coeff_dx_m1_p1,
+   coeff_dx_0_p1,
+   coeff_dx_p1_p1;
diff --git a/src/GeneralizedPolynomial-Uniform/2d.coeffs/2d.cube.size3/coeff-dxx.dcl.c b/src/GeneralizedPolynomial-Uniform/2d.coeffs/2d.cube.size3/coeff-dxx.dcl.c
new file mode 100644
index 0000000..2c62a96
--- /dev/null
+++ b/src/GeneralizedPolynomial-Uniform/2d.coeffs/2d.cube.size3/coeff-dxx.dcl.c
@@ -0,0 +1,9 @@
+fp coeff_dxx_m1_m1,
+   coeff_dxx_0_m1,
+   coeff_dxx_p1_m1,
+   coeff_dxx_m1_0,
+   coeff_dxx_0_0,
+   coeff_dxx_p1_0,
+   coeff_dxx_m1_p1,
+   coeff_dxx_0_p1,
+   coeff_dxx_p1_p1;
diff --git a/src/GeneralizedPolynomial-Uniform/2d.coeffs/2d.cube.size3/coeff-dxy.dcl.c b/src/GeneralizedPolynomial-Uniform/2d.coeffs/2d.cube.size3/coeff-dxy.dcl.c
new file mode 100644
index 0000000..0b6a45b
--- /dev/null
+++ b/src/GeneralizedPolynomial-Uniform/2d.coeffs/2d.cube.size3/coeff-dxy.dcl.c
@@ -0,0 +1,9 @@
+fp coeff_dxy_m1_m1,
+   coeff_dxy_0_m1,
+   coeff_dxy_p1_m1,
+   coeff_dxy_m1_0,
+   coeff_dxy_0_0,
+   coeff_dxy_p1_0,
+   coeff_dxy_m1_p1,
+   coeff_dxy_0_p1,
+   coeff_dxy_p1_p1;
diff --git a/src/GeneralizedPolynomial-Uniform/2d.coeffs/2d.cube.size3/coeff-dy.dcl.c b/src/GeneralizedPolynomial-Uniform/2d.coeffs/2d.cube.size3/coeff-dy.dcl.c
new file mode 100644
index 0000000..1aebea8
--- /dev/null
+++ b/src/GeneralizedPolynomial-Uniform/2d.coeffs/2d.cube.size3/coeff-dy.dcl.c
@@ -0,0 +1,9 @@
+fp coeff_dy_m1_m1,
+   coeff_dy_0_m1,
+   coeff_dy_p1_m1,
+   coeff_dy_m1_0,
+   coeff_dy_0_0,
+   coeff_dy_p1_0,
+   coeff_dy_m1_p1,
+   coeff_dy_0_p1,
+   coeff_dy_p1_p1;
diff --git a/src/GeneralizedPolynomial-Uniform/2d.coeffs/2d.cube.size3/coeff-dyy.dcl.c b/src/GeneralizedPolynomial-Uniform/2d.coeffs/2d.cube.size3/coeff-dyy.dcl.c
new file mode 100644
index 0000000..52b76c9
--- /dev/null
+++ b/src/GeneralizedPolynomial-Uniform/2d.coeffs/2d.cube.size3/coeff-dyy.dcl.c
@@ -0,0 +1,9 @@
+fp coeff_dyy_m1_m1,
+   coeff_dyy_0_m1,
+   coeff_dyy_p1_m1,
+   coeff_dyy_m1_0,
+   coeff_dyy_0_0,
+   coeff_dyy_p1_0,
+   coeff_dyy_m1_p1,
+   coeff_dyy_0_p1,
+   coeff_dyy_p1_p1;
diff --git a/src/GeneralizedPolynomial-Uniform/2d.coeffs/2d.cube.size3/data-var.assign.c b/src/GeneralizedPolynomial-Uniform/2d.coeffs/2d.cube.size3/data-var.assign.c
new file mode 100644
index 0000000..bb68189
--- /dev/null
+++ b/src/GeneralizedPolynomial-Uniform/2d.coeffs/2d.cube.size3/data-var.assign.c
@@ -0,0 +1,9 @@
+data_m1_m1 = DATA(-1,-1);
+data_0_m1 = DATA(0,-1);
+data_p1_m1 = DATA(1,-1);
+data_m1_0 = DATA(-1,0);
+data_0_0 = DATA(0,0);
+data_p1_0 = DATA(1,0);
+data_m1_p1 = DATA(-1,1);
+data_0_p1 = DATA(0,1);
+data_p1_p1 = DATA(1,1);
diff --git a/src/GeneralizedPolynomial-Uniform/2d.coeffs/2d.cube.size3/data-var.dcl.c b/src/GeneralizedPolynomial-Uniform/2d.coeffs/2d.cube.size3/data-var.dcl.c
new file mode 100644
index 0000000..77a064e
--- /dev/null
+++ b/src/GeneralizedPolynomial-Uniform/2d.coeffs/2d.cube.size3/data-var.dcl.c
@@ -0,0 +1,9 @@
+fp data_m1_m1,
+   data_0_m1,
+   data_p1_m1,
+   data_m1_0,
+   data_0_0,
+   data_p1_0,
+   data_m1_p1,
+   data_0_p1,
+   data_p1_p1;
diff --git a/src/GeneralizedPolynomial-Uniform/2d.coeffs/2d.cube.size3/interp-I.compute.c b/src/GeneralizedPolynomial-Uniform/2d.coeffs/2d.cube.size3/interp-I.compute.c
new file mode 100644
index 0000000..b9249fa
--- /dev/null
+++ b/src/GeneralizedPolynomial-Uniform/2d.coeffs/2d.cube.size3/interp-I.compute.c
@@ -0,0 +1,10 @@
+result =
+	coeff_I_m1_m1*data_m1_m1
+	+ coeff_I_0_m1*data_0_m1
+	+ coeff_I_p1_m1*data_p1_m1
+	+ coeff_I_m1_0*data_m1_0
+	+ coeff_I_0_0*data_0_0
+	+ coeff_I_p1_0*data_p1_0
+	+ coeff_I_m1_p1*data_m1_p1
+	+ coeff_I_0_p1*data_0_p1
+	+ coeff_I_p1_p1*data_p1_p1;
diff --git a/src/GeneralizedPolynomial-Uniform/2d.coeffs/2d.cube.size3/interp-dx.compute.c b/src/GeneralizedPolynomial-Uniform/2d.coeffs/2d.cube.size3/interp-dx.compute.c
new file mode 100644
index 0000000..f7d5f3a
--- /dev/null
+++ b/src/GeneralizedPolynomial-Uniform/2d.coeffs/2d.cube.size3/interp-dx.compute.c
@@ -0,0 +1,10 @@
+result =
+	coeff_dx_m1_m1*data_m1_m1
+	+ coeff_dx_0_m1*data_0_m1
+	+ coeff_dx_p1_m1*data_p1_m1
+	+ coeff_dx_m1_0*data_m1_0
+	+ coeff_dx_0_0*data_0_0
+	+ coeff_dx_p1_0*data_p1_0
+	+ coeff_dx_m1_p1*data_m1_p1
+	+ coeff_dx_0_p1*data_0_p1
+	+ coeff_dx_p1_p1*data_p1_p1;
diff --git a/src/GeneralizedPolynomial-Uniform/2d.coeffs/2d.cube.size3/interp-dxx.compute.c b/src/GeneralizedPolynomial-Uniform/2d.coeffs/2d.cube.size3/interp-dxx.compute.c
new file mode 100644
index 0000000..9764dbc
--- /dev/null
+++ b/src/GeneralizedPolynomial-Uniform/2d.coeffs/2d.cube.size3/interp-dxx.compute.c
@@ -0,0 +1,10 @@
+result =
+	coeff_dxx_m1_m1*data_m1_m1
+	+ coeff_dxx_0_m1*data_0_m1
+	+ coeff_dxx_p1_m1*data_p1_m1
+	+ coeff_dxx_m1_0*data_m1_0
+	+ coeff_dxx_0_0*data_0_0
+	+ coeff_dxx_p1_0*data_p1_0
+	+ coeff_dxx_m1_p1*data_m1_p1
+	+ coeff_dxx_0_p1*data_0_p1
+	+ coeff_dxx_p1_p1*data_p1_p1;
diff --git a/src/GeneralizedPolynomial-Uniform/2d.coeffs/2d.cube.size3/interp-dxy.compute.c b/src/GeneralizedPolynomial-Uniform/2d.coeffs/2d.cube.size3/interp-dxy.compute.c
new file mode 100644
index 0000000..5c9faa7
--- /dev/null
+++ b/src/GeneralizedPolynomial-Uniform/2d.coeffs/2d.cube.size3/interp-dxy.compute.c
@@ -0,0 +1,10 @@
+result =
+	coeff_dxy_m1_m1*data_m1_m1
+	+ coeff_dxy_0_m1*data_0_m1
+	+ coeff_dxy_p1_m1*data_p1_m1
+	+ coeff_dxy_m1_0*data_m1_0
+	+ coeff_dxy_0_0*data_0_0
+	+ coeff_dxy_p1_0*data_p1_0
+	+ coeff_dxy_m1_p1*data_m1_p1
+	+ coeff_dxy_0_p1*data_0_p1
+	+ coeff_dxy_p1_p1*data_p1_p1;
diff --git a/src/GeneralizedPolynomial-Uniform/2d.coeffs/2d.cube.size3/interp-dy.compute.c b/src/GeneralizedPolynomial-Uniform/2d.coeffs/2d.cube.size3/interp-dy.compute.c
new file mode 100644
index 0000000..d4eb9ca
--- /dev/null
+++ b/src/GeneralizedPolynomial-Uniform/2d.coeffs/2d.cube.size3/interp-dy.compute.c
@@ -0,0 +1,10 @@
+result =
+	coeff_dy_m1_m1*data_m1_m1
+	+ coeff_dy_0_m1*data_0_m1
+	+ coeff_dy_p1_m1*data_p1_m1
+	+ coeff_dy_m1_0*data_m1_0
+	+ coeff_dy_0_0*data_0_0
+	+ coeff_dy_p1_0*data_p1_0
+	+ coeff_dy_m1_p1*data_m1_p1
+	+ coeff_dy_0_p1*data_0_p1
+	+ coeff_dy_p1_p1*data_p1_p1;
diff --git a/src/GeneralizedPolynomial-Uniform/2d.coeffs/2d.cube.size3/interp-dyy.compute.c b/src/GeneralizedPolynomial-Uniform/2d.coeffs/2d.cube.size3/interp-dyy.compute.c
new file mode 100644
index 0000000..6ee2445
--- /dev/null
+++ b/src/GeneralizedPolynomial-Uniform/2d.coeffs/2d.cube.size3/interp-dyy.compute.c
@@ -0,0 +1,10 @@
+result =
+	coeff_dyy_m1_m1*data_m1_m1
+	+ coeff_dyy_0_m1*data_0_m1
+	+ coeff_dyy_p1_m1*data_p1_m1
+	+ coeff_dyy_m1_0*data_m1_0
+	+ coeff_dyy_0_0*data_0_0
+	+ coeff_dyy_p1_0*data_p1_0
+	+ coeff_dyy_m1_p1*data_m1_p1
+	+ coeff_dyy_0_p1*data_0_p1
+	+ coeff_dyy_p1_p1*data_p1_p1;