aboutsummaryrefslogtreecommitdiff
path: root/src/GeneralizedPolynomial-Uniform/common/2d.cube.size4
diff options
context:
space:
mode:
Diffstat (limited to 'src/GeneralizedPolynomial-Uniform/common/2d.cube.size4')
-rw-r--r--src/GeneralizedPolynomial-Uniform/common/2d.cube.size4/coeff-I.dcl.c16
-rw-r--r--src/GeneralizedPolynomial-Uniform/common/2d.cube.size4/coeff-I.store.c16
-rw-r--r--src/GeneralizedPolynomial-Uniform/common/2d.cube.size4/coeff-dx.dcl.c16
-rw-r--r--src/GeneralizedPolynomial-Uniform/common/2d.cube.size4/coeff-dx.store.c16
-rw-r--r--src/GeneralizedPolynomial-Uniform/common/2d.cube.size4/coeff-dxx.dcl.c16
-rw-r--r--src/GeneralizedPolynomial-Uniform/common/2d.cube.size4/coeff-dxx.store.c16
-rw-r--r--src/GeneralizedPolynomial-Uniform/common/2d.cube.size4/coeff-dxy.dcl.c16
-rw-r--r--src/GeneralizedPolynomial-Uniform/common/2d.cube.size4/coeff-dxy.store.c16
-rw-r--r--src/GeneralizedPolynomial-Uniform/common/2d.cube.size4/coeff-dy.dcl.c16
-rw-r--r--src/GeneralizedPolynomial-Uniform/common/2d.cube.size4/coeff-dy.store.c16
-rw-r--r--src/GeneralizedPolynomial-Uniform/common/2d.cube.size4/coeff-dyy.dcl.c16
-rw-r--r--src/GeneralizedPolynomial-Uniform/common/2d.cube.size4/coeff-dyy.store.c16
-rw-r--r--src/GeneralizedPolynomial-Uniform/common/2d.cube.size4/data-var.assign.c16
-rw-r--r--src/GeneralizedPolynomial-Uniform/common/2d.cube.size4/data-var.dcl.c16
-rw-r--r--src/GeneralizedPolynomial-Uniform/common/2d.cube.size4/interp-I.compute.c17
-rw-r--r--src/GeneralizedPolynomial-Uniform/common/2d.cube.size4/interp-dx.compute.c17
-rw-r--r--src/GeneralizedPolynomial-Uniform/common/2d.cube.size4/interp-dxx.compute.c17
-rw-r--r--src/GeneralizedPolynomial-Uniform/common/2d.cube.size4/interp-dxy.compute.c17
-rw-r--r--src/GeneralizedPolynomial-Uniform/common/2d.cube.size4/interp-dy.compute.c17
-rw-r--r--src/GeneralizedPolynomial-Uniform/common/2d.cube.size4/interp-dyy.compute.c17
20 files changed, 326 insertions, 0 deletions
diff --git a/src/GeneralizedPolynomial-Uniform/common/2d.cube.size4/coeff-I.dcl.c b/src/GeneralizedPolynomial-Uniform/common/2d.cube.size4/coeff-I.dcl.c
new file mode 100644
index 0000000..d35d81b
--- /dev/null
+++ b/src/GeneralizedPolynomial-Uniform/common/2d.cube.size4/coeff-I.dcl.c
@@ -0,0 +1,16 @@
+fp coeff_I_m1_m1,
+ coeff_I_0_m1,
+ coeff_I_p1_m1,
+ coeff_I_p2_m1,
+ coeff_I_m1_0,
+ coeff_I_0_0,
+ coeff_I_p1_0,
+ coeff_I_p2_0,
+ coeff_I_m1_p1,
+ coeff_I_0_p1,
+ coeff_I_p1_p1,
+ coeff_I_p2_p1,
+ coeff_I_m1_p2,
+ coeff_I_0_p2,
+ coeff_I_p1_p2,
+ coeff_I_p2_p2;
diff --git a/src/GeneralizedPolynomial-Uniform/common/2d.cube.size4/coeff-I.store.c b/src/GeneralizedPolynomial-Uniform/common/2d.cube.size4/coeff-I.store.c
new file mode 100644
index 0000000..729b8ea
--- /dev/null
+++ b/src/GeneralizedPolynomial-Uniform/common/2d.cube.size4/coeff-I.store.c
@@ -0,0 +1,16 @@
+COEFF(-1,-1) = coeff_I_m1_m1;
+COEFF(0,-1) = coeff_I_0_m1;
+COEFF(1,-1) = coeff_I_p1_m1;
+COEFF(2,-1) = coeff_I_p2_m1;
+COEFF(-1,0) = coeff_I_m1_0;
+COEFF(0,0) = coeff_I_0_0;
+COEFF(1,0) = coeff_I_p1_0;
+COEFF(2,0) = coeff_I_p2_0;
+COEFF(-1,1) = coeff_I_m1_p1;
+COEFF(0,1) = coeff_I_0_p1;
+COEFF(1,1) = coeff_I_p1_p1;
+COEFF(2,1) = coeff_I_p2_p1;
+COEFF(-1,2) = coeff_I_m1_p2;
+COEFF(0,2) = coeff_I_0_p2;
+COEFF(1,2) = coeff_I_p1_p2;
+COEFF(2,2) = coeff_I_p2_p2;
diff --git a/src/GeneralizedPolynomial-Uniform/common/2d.cube.size4/coeff-dx.dcl.c b/src/GeneralizedPolynomial-Uniform/common/2d.cube.size4/coeff-dx.dcl.c
new file mode 100644
index 0000000..edfc9bb
--- /dev/null
+++ b/src/GeneralizedPolynomial-Uniform/common/2d.cube.size4/coeff-dx.dcl.c
@@ -0,0 +1,16 @@
+fp coeff_dx_m1_m1,
+ coeff_dx_0_m1,
+ coeff_dx_p1_m1,
+ coeff_dx_p2_m1,
+ coeff_dx_m1_0,
+ coeff_dx_0_0,
+ coeff_dx_p1_0,
+ coeff_dx_p2_0,
+ coeff_dx_m1_p1,
+ coeff_dx_0_p1,
+ coeff_dx_p1_p1,
+ coeff_dx_p2_p1,
+ coeff_dx_m1_p2,
+ coeff_dx_0_p2,
+ coeff_dx_p1_p2,
+ coeff_dx_p2_p2;
diff --git a/src/GeneralizedPolynomial-Uniform/common/2d.cube.size4/coeff-dx.store.c b/src/GeneralizedPolynomial-Uniform/common/2d.cube.size4/coeff-dx.store.c
new file mode 100644
index 0000000..c1bdc39
--- /dev/null
+++ b/src/GeneralizedPolynomial-Uniform/common/2d.cube.size4/coeff-dx.store.c
@@ -0,0 +1,16 @@
+COEFF(-1,-1) = factor * coeff_dx_m1_m1;
+COEFF(0,-1) = factor * coeff_dx_0_m1;
+COEFF(1,-1) = factor * coeff_dx_p1_m1;
+COEFF(2,-1) = factor * coeff_dx_p2_m1;
+COEFF(-1,0) = factor * coeff_dx_m1_0;
+COEFF(0,0) = factor * coeff_dx_0_0;
+COEFF(1,0) = factor * coeff_dx_p1_0;
+COEFF(2,0) = factor * coeff_dx_p2_0;
+COEFF(-1,1) = factor * coeff_dx_m1_p1;
+COEFF(0,1) = factor * coeff_dx_0_p1;
+COEFF(1,1) = factor * coeff_dx_p1_p1;
+COEFF(2,1) = factor * coeff_dx_p2_p1;
+COEFF(-1,2) = factor * coeff_dx_m1_p2;
+COEFF(0,2) = factor * coeff_dx_0_p2;
+COEFF(1,2) = factor * coeff_dx_p1_p2;
+COEFF(2,2) = factor * coeff_dx_p2_p2;
diff --git a/src/GeneralizedPolynomial-Uniform/common/2d.cube.size4/coeff-dxx.dcl.c b/src/GeneralizedPolynomial-Uniform/common/2d.cube.size4/coeff-dxx.dcl.c
new file mode 100644
index 0000000..b5add16
--- /dev/null
+++ b/src/GeneralizedPolynomial-Uniform/common/2d.cube.size4/coeff-dxx.dcl.c
@@ -0,0 +1,16 @@
+fp coeff_dxx_m1_m1,
+ coeff_dxx_0_m1,
+ coeff_dxx_p1_m1,
+ coeff_dxx_p2_m1,
+ coeff_dxx_m1_0,
+ coeff_dxx_0_0,
+ coeff_dxx_p1_0,
+ coeff_dxx_p2_0,
+ coeff_dxx_m1_p1,
+ coeff_dxx_0_p1,
+ coeff_dxx_p1_p1,
+ coeff_dxx_p2_p1,
+ coeff_dxx_m1_p2,
+ coeff_dxx_0_p2,
+ coeff_dxx_p1_p2,
+ coeff_dxx_p2_p2;
diff --git a/src/GeneralizedPolynomial-Uniform/common/2d.cube.size4/coeff-dxx.store.c b/src/GeneralizedPolynomial-Uniform/common/2d.cube.size4/coeff-dxx.store.c
new file mode 100644
index 0000000..e65fc6f
--- /dev/null
+++ b/src/GeneralizedPolynomial-Uniform/common/2d.cube.size4/coeff-dxx.store.c
@@ -0,0 +1,16 @@
+COEFF(-1,-1) = factor * coeff_dxx_m1_m1;
+COEFF(0,-1) = factor * coeff_dxx_0_m1;
+COEFF(1,-1) = factor * coeff_dxx_p1_m1;
+COEFF(2,-1) = factor * coeff_dxx_p2_m1;
+COEFF(-1,0) = factor * coeff_dxx_m1_0;
+COEFF(0,0) = factor * coeff_dxx_0_0;
+COEFF(1,0) = factor * coeff_dxx_p1_0;
+COEFF(2,0) = factor * coeff_dxx_p2_0;
+COEFF(-1,1) = factor * coeff_dxx_m1_p1;
+COEFF(0,1) = factor * coeff_dxx_0_p1;
+COEFF(1,1) = factor * coeff_dxx_p1_p1;
+COEFF(2,1) = factor * coeff_dxx_p2_p1;
+COEFF(-1,2) = factor * coeff_dxx_m1_p2;
+COEFF(0,2) = factor * coeff_dxx_0_p2;
+COEFF(1,2) = factor * coeff_dxx_p1_p2;
+COEFF(2,2) = factor * coeff_dxx_p2_p2;
diff --git a/src/GeneralizedPolynomial-Uniform/common/2d.cube.size4/coeff-dxy.dcl.c b/src/GeneralizedPolynomial-Uniform/common/2d.cube.size4/coeff-dxy.dcl.c
new file mode 100644
index 0000000..853b503
--- /dev/null
+++ b/src/GeneralizedPolynomial-Uniform/common/2d.cube.size4/coeff-dxy.dcl.c
@@ -0,0 +1,16 @@
+fp coeff_dxy_m1_m1,
+ coeff_dxy_0_m1,
+ coeff_dxy_p1_m1,
+ coeff_dxy_p2_m1,
+ coeff_dxy_m1_0,
+ coeff_dxy_0_0,
+ coeff_dxy_p1_0,
+ coeff_dxy_p2_0,
+ coeff_dxy_m1_p1,
+ coeff_dxy_0_p1,
+ coeff_dxy_p1_p1,
+ coeff_dxy_p2_p1,
+ coeff_dxy_m1_p2,
+ coeff_dxy_0_p2,
+ coeff_dxy_p1_p2,
+ coeff_dxy_p2_p2;
diff --git a/src/GeneralizedPolynomial-Uniform/common/2d.cube.size4/coeff-dxy.store.c b/src/GeneralizedPolynomial-Uniform/common/2d.cube.size4/coeff-dxy.store.c
new file mode 100644
index 0000000..442cfe8
--- /dev/null
+++ b/src/GeneralizedPolynomial-Uniform/common/2d.cube.size4/coeff-dxy.store.c
@@ -0,0 +1,16 @@
+COEFF(-1,-1) = factor * coeff_dxy_m1_m1;
+COEFF(0,-1) = factor * coeff_dxy_0_m1;
+COEFF(1,-1) = factor * coeff_dxy_p1_m1;
+COEFF(2,-1) = factor * coeff_dxy_p2_m1;
+COEFF(-1,0) = factor * coeff_dxy_m1_0;
+COEFF(0,0) = factor * coeff_dxy_0_0;
+COEFF(1,0) = factor * coeff_dxy_p1_0;
+COEFF(2,0) = factor * coeff_dxy_p2_0;
+COEFF(-1,1) = factor * coeff_dxy_m1_p1;
+COEFF(0,1) = factor * coeff_dxy_0_p1;
+COEFF(1,1) = factor * coeff_dxy_p1_p1;
+COEFF(2,1) = factor * coeff_dxy_p2_p1;
+COEFF(-1,2) = factor * coeff_dxy_m1_p2;
+COEFF(0,2) = factor * coeff_dxy_0_p2;
+COEFF(1,2) = factor * coeff_dxy_p1_p2;
+COEFF(2,2) = factor * coeff_dxy_p2_p2;
diff --git a/src/GeneralizedPolynomial-Uniform/common/2d.cube.size4/coeff-dy.dcl.c b/src/GeneralizedPolynomial-Uniform/common/2d.cube.size4/coeff-dy.dcl.c
new file mode 100644
index 0000000..e27b5ff
--- /dev/null
+++ b/src/GeneralizedPolynomial-Uniform/common/2d.cube.size4/coeff-dy.dcl.c
@@ -0,0 +1,16 @@
+fp coeff_dy_m1_m1,
+ coeff_dy_0_m1,
+ coeff_dy_p1_m1,
+ coeff_dy_p2_m1,
+ coeff_dy_m1_0,
+ coeff_dy_0_0,
+ coeff_dy_p1_0,
+ coeff_dy_p2_0,
+ coeff_dy_m1_p1,
+ coeff_dy_0_p1,
+ coeff_dy_p1_p1,
+ coeff_dy_p2_p1,
+ coeff_dy_m1_p2,
+ coeff_dy_0_p2,
+ coeff_dy_p1_p2,
+ coeff_dy_p2_p2;
diff --git a/src/GeneralizedPolynomial-Uniform/common/2d.cube.size4/coeff-dy.store.c b/src/GeneralizedPolynomial-Uniform/common/2d.cube.size4/coeff-dy.store.c
new file mode 100644
index 0000000..fa3aefd
--- /dev/null
+++ b/src/GeneralizedPolynomial-Uniform/common/2d.cube.size4/coeff-dy.store.c
@@ -0,0 +1,16 @@
+COEFF(-1,-1) = factor * coeff_dy_m1_m1;
+COEFF(0,-1) = factor * coeff_dy_0_m1;
+COEFF(1,-1) = factor * coeff_dy_p1_m1;
+COEFF(2,-1) = factor * coeff_dy_p2_m1;
+COEFF(-1,0) = factor * coeff_dy_m1_0;
+COEFF(0,0) = factor * coeff_dy_0_0;
+COEFF(1,0) = factor * coeff_dy_p1_0;
+COEFF(2,0) = factor * coeff_dy_p2_0;
+COEFF(-1,1) = factor * coeff_dy_m1_p1;
+COEFF(0,1) = factor * coeff_dy_0_p1;
+COEFF(1,1) = factor * coeff_dy_p1_p1;
+COEFF(2,1) = factor * coeff_dy_p2_p1;
+COEFF(-1,2) = factor * coeff_dy_m1_p2;
+COEFF(0,2) = factor * coeff_dy_0_p2;
+COEFF(1,2) = factor * coeff_dy_p1_p2;
+COEFF(2,2) = factor * coeff_dy_p2_p2;
diff --git a/src/GeneralizedPolynomial-Uniform/common/2d.cube.size4/coeff-dyy.dcl.c b/src/GeneralizedPolynomial-Uniform/common/2d.cube.size4/coeff-dyy.dcl.c
new file mode 100644
index 0000000..ee37cfe
--- /dev/null
+++ b/src/GeneralizedPolynomial-Uniform/common/2d.cube.size4/coeff-dyy.dcl.c
@@ -0,0 +1,16 @@
+fp coeff_dyy_m1_m1,
+ coeff_dyy_0_m1,
+ coeff_dyy_p1_m1,
+ coeff_dyy_p2_m1,
+ coeff_dyy_m1_0,
+ coeff_dyy_0_0,
+ coeff_dyy_p1_0,
+ coeff_dyy_p2_0,
+ coeff_dyy_m1_p1,
+ coeff_dyy_0_p1,
+ coeff_dyy_p1_p1,
+ coeff_dyy_p2_p1,
+ coeff_dyy_m1_p2,
+ coeff_dyy_0_p2,
+ coeff_dyy_p1_p2,
+ coeff_dyy_p2_p2;
diff --git a/src/GeneralizedPolynomial-Uniform/common/2d.cube.size4/coeff-dyy.store.c b/src/GeneralizedPolynomial-Uniform/common/2d.cube.size4/coeff-dyy.store.c
new file mode 100644
index 0000000..3f29fbf
--- /dev/null
+++ b/src/GeneralizedPolynomial-Uniform/common/2d.cube.size4/coeff-dyy.store.c
@@ -0,0 +1,16 @@
+COEFF(-1,-1) = factor * coeff_dyy_m1_m1;
+COEFF(0,-1) = factor * coeff_dyy_0_m1;
+COEFF(1,-1) = factor * coeff_dyy_p1_m1;
+COEFF(2,-1) = factor * coeff_dyy_p2_m1;
+COEFF(-1,0) = factor * coeff_dyy_m1_0;
+COEFF(0,0) = factor * coeff_dyy_0_0;
+COEFF(1,0) = factor * coeff_dyy_p1_0;
+COEFF(2,0) = factor * coeff_dyy_p2_0;
+COEFF(-1,1) = factor * coeff_dyy_m1_p1;
+COEFF(0,1) = factor * coeff_dyy_0_p1;
+COEFF(1,1) = factor * coeff_dyy_p1_p1;
+COEFF(2,1) = factor * coeff_dyy_p2_p1;
+COEFF(-1,2) = factor * coeff_dyy_m1_p2;
+COEFF(0,2) = factor * coeff_dyy_0_p2;
+COEFF(1,2) = factor * coeff_dyy_p1_p2;
+COEFF(2,2) = factor * coeff_dyy_p2_p2;
diff --git a/src/GeneralizedPolynomial-Uniform/common/2d.cube.size4/data-var.assign.c b/src/GeneralizedPolynomial-Uniform/common/2d.cube.size4/data-var.assign.c
new file mode 100644
index 0000000..d0316f1
--- /dev/null
+++ b/src/GeneralizedPolynomial-Uniform/common/2d.cube.size4/data-var.assign.c
@@ -0,0 +1,16 @@
+data_m1_m1 = DATA(-1,-1);
+data_0_m1 = DATA(0,-1);
+data_p1_m1 = DATA(1,-1);
+data_p2_m1 = DATA(2,-1);
+data_m1_0 = DATA(-1,0);
+data_0_0 = DATA(0,0);
+data_p1_0 = DATA(1,0);
+data_p2_0 = DATA(2,0);
+data_m1_p1 = DATA(-1,1);
+data_0_p1 = DATA(0,1);
+data_p1_p1 = DATA(1,1);
+data_p2_p1 = DATA(2,1);
+data_m1_p2 = DATA(-1,2);
+data_0_p2 = DATA(0,2);
+data_p1_p2 = DATA(1,2);
+data_p2_p2 = DATA(2,2);
diff --git a/src/GeneralizedPolynomial-Uniform/common/2d.cube.size4/data-var.dcl.c b/src/GeneralizedPolynomial-Uniform/common/2d.cube.size4/data-var.dcl.c
new file mode 100644
index 0000000..abbfb5d
--- /dev/null
+++ b/src/GeneralizedPolynomial-Uniform/common/2d.cube.size4/data-var.dcl.c
@@ -0,0 +1,16 @@
+fp data_m1_m1,
+ data_0_m1,
+ data_p1_m1,
+ data_p2_m1,
+ data_m1_0,
+ data_0_0,
+ data_p1_0,
+ data_p2_0,
+ data_m1_p1,
+ data_0_p1,
+ data_p1_p1,
+ data_p2_p1,
+ data_m1_p2,
+ data_0_p2,
+ data_p1_p2,
+ data_p2_p2;
diff --git a/src/GeneralizedPolynomial-Uniform/common/2d.cube.size4/interp-I.compute.c b/src/GeneralizedPolynomial-Uniform/common/2d.cube.size4/interp-I.compute.c
new file mode 100644
index 0000000..cf13e49
--- /dev/null
+++ b/src/GeneralizedPolynomial-Uniform/common/2d.cube.size4/interp-I.compute.c
@@ -0,0 +1,17 @@
+result =
+ coeff_I_m1_m1*data_m1_m1
+ + coeff_I_0_m1*data_0_m1
+ + coeff_I_p1_m1*data_p1_m1
+ + coeff_I_p2_m1*data_p2_m1
+ + coeff_I_m1_0*data_m1_0
+ + coeff_I_0_0*data_0_0
+ + coeff_I_p1_0*data_p1_0
+ + coeff_I_p2_0*data_p2_0
+ + coeff_I_m1_p1*data_m1_p1
+ + coeff_I_0_p1*data_0_p1
+ + coeff_I_p1_p1*data_p1_p1
+ + coeff_I_p2_p1*data_p2_p1
+ + coeff_I_m1_p2*data_m1_p2
+ + coeff_I_0_p2*data_0_p2
+ + coeff_I_p1_p2*data_p1_p2
+ + coeff_I_p2_p2*data_p2_p2;
diff --git a/src/GeneralizedPolynomial-Uniform/common/2d.cube.size4/interp-dx.compute.c b/src/GeneralizedPolynomial-Uniform/common/2d.cube.size4/interp-dx.compute.c
new file mode 100644
index 0000000..69f9f56
--- /dev/null
+++ b/src/GeneralizedPolynomial-Uniform/common/2d.cube.size4/interp-dx.compute.c
@@ -0,0 +1,17 @@
+result =
+ coeff_dx_m1_m1*data_m1_m1
+ + coeff_dx_0_m1*data_0_m1
+ + coeff_dx_p1_m1*data_p1_m1
+ + coeff_dx_p2_m1*data_p2_m1
+ + coeff_dx_m1_0*data_m1_0
+ + coeff_dx_0_0*data_0_0
+ + coeff_dx_p1_0*data_p1_0
+ + coeff_dx_p2_0*data_p2_0
+ + coeff_dx_m1_p1*data_m1_p1
+ + coeff_dx_0_p1*data_0_p1
+ + coeff_dx_p1_p1*data_p1_p1
+ + coeff_dx_p2_p1*data_p2_p1
+ + coeff_dx_m1_p2*data_m1_p2
+ + coeff_dx_0_p2*data_0_p2
+ + coeff_dx_p1_p2*data_p1_p2
+ + coeff_dx_p2_p2*data_p2_p2;
diff --git a/src/GeneralizedPolynomial-Uniform/common/2d.cube.size4/interp-dxx.compute.c b/src/GeneralizedPolynomial-Uniform/common/2d.cube.size4/interp-dxx.compute.c
new file mode 100644
index 0000000..9ce20a1
--- /dev/null
+++ b/src/GeneralizedPolynomial-Uniform/common/2d.cube.size4/interp-dxx.compute.c
@@ -0,0 +1,17 @@
+result =
+ coeff_dxx_m1_m1*data_m1_m1
+ + coeff_dxx_0_m1*data_0_m1
+ + coeff_dxx_p1_m1*data_p1_m1
+ + coeff_dxx_p2_m1*data_p2_m1
+ + coeff_dxx_m1_0*data_m1_0
+ + coeff_dxx_0_0*data_0_0
+ + coeff_dxx_p1_0*data_p1_0
+ + coeff_dxx_p2_0*data_p2_0
+ + coeff_dxx_m1_p1*data_m1_p1
+ + coeff_dxx_0_p1*data_0_p1
+ + coeff_dxx_p1_p1*data_p1_p1
+ + coeff_dxx_p2_p1*data_p2_p1
+ + coeff_dxx_m1_p2*data_m1_p2
+ + coeff_dxx_0_p2*data_0_p2
+ + coeff_dxx_p1_p2*data_p1_p2
+ + coeff_dxx_p2_p2*data_p2_p2;
diff --git a/src/GeneralizedPolynomial-Uniform/common/2d.cube.size4/interp-dxy.compute.c b/src/GeneralizedPolynomial-Uniform/common/2d.cube.size4/interp-dxy.compute.c
new file mode 100644
index 0000000..5130c7b
--- /dev/null
+++ b/src/GeneralizedPolynomial-Uniform/common/2d.cube.size4/interp-dxy.compute.c
@@ -0,0 +1,17 @@
+result =
+ coeff_dxy_m1_m1*data_m1_m1
+ + coeff_dxy_0_m1*data_0_m1
+ + coeff_dxy_p1_m1*data_p1_m1
+ + coeff_dxy_p2_m1*data_p2_m1
+ + coeff_dxy_m1_0*data_m1_0
+ + coeff_dxy_0_0*data_0_0
+ + coeff_dxy_p1_0*data_p1_0
+ + coeff_dxy_p2_0*data_p2_0
+ + coeff_dxy_m1_p1*data_m1_p1
+ + coeff_dxy_0_p1*data_0_p1
+ + coeff_dxy_p1_p1*data_p1_p1
+ + coeff_dxy_p2_p1*data_p2_p1
+ + coeff_dxy_m1_p2*data_m1_p2
+ + coeff_dxy_0_p2*data_0_p2
+ + coeff_dxy_p1_p2*data_p1_p2
+ + coeff_dxy_p2_p2*data_p2_p2;
diff --git a/src/GeneralizedPolynomial-Uniform/common/2d.cube.size4/interp-dy.compute.c b/src/GeneralizedPolynomial-Uniform/common/2d.cube.size4/interp-dy.compute.c
new file mode 100644
index 0000000..f671033
--- /dev/null
+++ b/src/GeneralizedPolynomial-Uniform/common/2d.cube.size4/interp-dy.compute.c
@@ -0,0 +1,17 @@
+result =
+ coeff_dy_m1_m1*data_m1_m1
+ + coeff_dy_0_m1*data_0_m1
+ + coeff_dy_p1_m1*data_p1_m1
+ + coeff_dy_p2_m1*data_p2_m1
+ + coeff_dy_m1_0*data_m1_0
+ + coeff_dy_0_0*data_0_0
+ + coeff_dy_p1_0*data_p1_0
+ + coeff_dy_p2_0*data_p2_0
+ + coeff_dy_m1_p1*data_m1_p1
+ + coeff_dy_0_p1*data_0_p1
+ + coeff_dy_p1_p1*data_p1_p1
+ + coeff_dy_p2_p1*data_p2_p1
+ + coeff_dy_m1_p2*data_m1_p2
+ + coeff_dy_0_p2*data_0_p2
+ + coeff_dy_p1_p2*data_p1_p2
+ + coeff_dy_p2_p2*data_p2_p2;
diff --git a/src/GeneralizedPolynomial-Uniform/common/2d.cube.size4/interp-dyy.compute.c b/src/GeneralizedPolynomial-Uniform/common/2d.cube.size4/interp-dyy.compute.c
new file mode 100644
index 0000000..112e958
--- /dev/null
+++ b/src/GeneralizedPolynomial-Uniform/common/2d.cube.size4/interp-dyy.compute.c
@@ -0,0 +1,17 @@
+result =
+ coeff_dyy_m1_m1*data_m1_m1
+ + coeff_dyy_0_m1*data_0_m1
+ + coeff_dyy_p1_m1*data_p1_m1
+ + coeff_dyy_p2_m1*data_p2_m1
+ + coeff_dyy_m1_0*data_m1_0
+ + coeff_dyy_0_0*data_0_0
+ + coeff_dyy_p1_0*data_p1_0
+ + coeff_dyy_p2_0*data_p2_0
+ + coeff_dyy_m1_p1*data_m1_p1
+ + coeff_dyy_0_p1*data_0_p1
+ + coeff_dyy_p1_p1*data_p1_p1
+ + coeff_dyy_p2_p1*data_p2_p1
+ + coeff_dyy_m1_p2*data_m1_p2
+ + coeff_dyy_0_p2*data_0_p2
+ + coeff_dyy_p1_p2*data_p1_p2
+ + coeff_dyy_p2_p2*data_p2_p2;