add support for CarpetInterp2 to InterpToArray

git-svn-id: http://svn.cactuscode.org/arrangements/CactusNumerical/InterpToArray/trunk@27 6ca9aeac-0e4f-0410-a746-fe4df63e9d0c

19 files changed, 166 insertions, 0 deletions

diff --git a/test/test.ccl b/test/test.ccl
new file mode 100644
index 0000000..7202173
--- /dev/null
+++ b/test/test.ccl
@@ -0,0 +1,5 @@
+# Carpet output contains empty lines whose number depends the number of processes
+TEST wavetoy2
+{
+  NPROCS 2
+}
diff --git a/test/wavetoy2.par b/test/wavetoy2.par
new file mode 100644
index 0000000..a32dddf
--- /dev/null
+++ b/test/wavetoy2.par
@@ -0,0 +1,88 @@
+!DESC "Test InterpToArray with WaveToy initial data"
+
+ActiveThorns = "
+        Boundary
+        CartGrid3D
+        CoordBase
+        IDScalarWaveC
+        CarpetIOASCII
+        CarpetIOBasic
+        CarpetIOScalar
+        IOUtil
+        InterpToArray
+        LocalReduce
+        Carpet
+        CarpetLib
+        GSl
+        LoopControl
+        InitBase
+        CarpetInterp2
+        CarpetReduce
+        CarpetSlab
+        SymBase
+        WaveToyC
+"
+
+Cactus::cctk_itlast = 0
+
+driver::global_nx  = 15
+driver::global_ny  = 15
+driver::global_nz  = 15
+driver::ghost_size =  2
+
+grid::type = "byspacing"
+grid::dxyz = 0.6
+
+IDScalarWave::initial_data = "Gaussian"
+IDScalarWave::sigma        = 2.8
+IDScalarWave::radius       = 5.0
+
+
+
+InterpToArray::use_carpetinterp2  = "yes"
+InterpToArray::carpetinterp2_interpolator_order = 2
+
+InterpToArray::narrays2d          = 1
+InterpToArray::array2d_vars[0]    = "wavetoy::phi"
+InterpToArray::array2d_x0         = -1.0
+InterpToArray::array2d_y0         = -1.0
+InterpToArray::array2d_z0         = -1.0
+InterpToArray::array2d_dx_i       =  0.5
+InterpToArray::array2d_dy_j       =  0.4
+InterpToArray::array2d_dz_j       =  0.1
+InterpToArray::array2d_npoints_i  =  3
+InterpToArray::array2d_npoints_j  =  4
+
+InterpToArray::nparrays3d         = 1
+InterpToArray::parray3d_vars[0]   = "wavetoy::phi"
+InterpToArray::parray3d_x0        = -1.0
+InterpToArray::parray3d_y0        = -1.0
+InterpToArray::parray3d_z0        = -1.0
+InterpToArray::parray3d_dx_i      =  0.5
+InterpToArray::parray3d_dy_j      =  0.5
+InterpToArray::parray3d_dz_k      =  0.5
+InterpToArray::parray3d_npoints_i =  3
+InterpToArray::parray3d_npoints_j =  4
+InterpToArray::parray3d_npoints_k =  5
+
+
+
+IO::out_dir       = $parfile
+IO::parfile_write = "no"
+IO::out_fileinfo  = "none"
+
+IOScalar::outScalar_reductions = "norm1 norm2 minimum maximum norm_inf"
+IOScalar::outScalar_every      = 1
+IOScalar::outScalar_vars       = "
+        InterpToArray::arrays2d[0]
+        InterpToArray::parrays3d[0]
+"
+
+IOASCII::out1D_every = 1
+IOASCII::compact_format = "yes"
+IOASCII::output_ghost_points = "no"
+IOASCII::out1D_vars  = "
+        InterpToArray::scalars
+        InterpToArray::arrays2d
+        InterpToArray::parrays3d
+"
diff --git a/test/wavetoy2/arrays2d[0].maximum.asc b/test/wavetoy2/arrays2d[0].maximum.asc
new file mode 100644
index 0000000..a0f9244
--- /dev/null
+++ b/test/wavetoy2/arrays2d[0].maximum.asc
@@ -0,0 +1,3 @@
+# Scalar ASCII output created by CarpetIOScalar
+#
+0 0 0.255929043803271
diff --git a/test/wavetoy2/arrays2d[0].minimum.asc b/test/wavetoy2/arrays2d[0].minimum.asc
new file mode 100644
index 0000000..eb12060
--- /dev/null
+++ b/test/wavetoy2/arrays2d[0].minimum.asc
@@ -0,0 +1,3 @@
+# Scalar ASCII output created by CarpetIOScalar
+#
+0 0 0.0981997271447794
diff --git a/test/wavetoy2/arrays2d[0].norm1.asc b/test/wavetoy2/arrays2d[0].norm1.asc
new file mode 100644
index 0000000..8845a35
--- /dev/null
+++ b/test/wavetoy2/arrays2d[0].norm1.asc
@@ -0,0 +1,3 @@
+# Scalar ASCII output created by CarpetIOScalar
+#
+0 0 0.162223918469515
diff --git a/test/wavetoy2/arrays2d[0].norm2.asc b/test/wavetoy2/arrays2d[0].norm2.asc
new file mode 100644
index 0000000..69182e8
--- /dev/null
+++ b/test/wavetoy2/arrays2d[0].norm2.asc
@@ -0,0 +1,3 @@
+# Scalar ASCII output created by CarpetIOScalar
+#
+0 0 0.168525955789019
diff --git a/test/wavetoy2/arrays2d[0].norm_inf.asc b/test/wavetoy2/arrays2d[0].norm_inf.asc
new file mode 100644
index 0000000..a0f9244
--- /dev/null
+++ b/test/wavetoy2/arrays2d[0].norm_inf.asc
@@ -0,0 +1,3 @@
+# Scalar ASCII output created by CarpetIOScalar
+#
+0 0 0.255929043803271
diff --git a/test/wavetoy2/arrays2d[0].x.asc b/test/wavetoy2/arrays2d[0].x.asc
new file mode 100644
index 0000000..d821c81
--- /dev/null
+++ b/test/wavetoy2/arrays2d[0].x.asc
@@ -0,0 +1,6 @@
+# 1D ASCII output created by CarpetIOASCII
+#
+0	0 0	0	0.255929043803271
+0	1 0	0	0.20945834089817
+0	2 0	0	0.193918918322488
+
diff --git a/test/wavetoy2/arrays2d[0].y.asc b/test/wavetoy2/arrays2d[0].y.asc
new file mode 100644
index 0000000..a7fd7ef
--- /dev/null
+++ b/test/wavetoy2/arrays2d[0].y.asc
@@ -0,0 +1,7 @@
+# 1D ASCII output created by CarpetIOASCII
+#
+0	0 0	0	0.255929043803271
+0	0 1	0	0.204490368041785
+0	0 2	0	0.173853957866224
+0	0 3	0	0.164387756102711
+
diff --git a/test/wavetoy2/arrays2d[0].z.asc b/test/wavetoy2/arrays2d[0].z.asc
new file mode 100644
index 0000000..98bbcc6
--- /dev/null
+++ b/test/wavetoy2/arrays2d[0].z.asc
@@ -0,0 +1,3 @@
+# 1D ASCII output created by CarpetIOASCII
+#
+0	0 0	0	0.255929043803271
diff --git a/test/wavetoy2/parrays3d[0].d.asc b/test/wavetoy2/parrays3d[0].d.asc
new file mode 100644
index 0000000..d537add
--- /dev/null
+++ b/test/wavetoy2/parrays3d[0].d.asc
@@ -0,0 +1,5 @@
+# 1D ASCII output created by CarpetIOASCII
+#
+0	0 0 0	0	0 0 0	0.255929043803271
+0	1 1 1	0	1 1 1	0.111155027559738
+0	2 2 2	0	2 2 2	0.0575070381476306
diff --git a/test/wavetoy2/parrays3d[0].maximum.asc b/test/wavetoy2/parrays3d[0].maximum.asc
new file mode 100644
index 0000000..a0f9244
--- /dev/null
+++ b/test/wavetoy2/parrays3d[0].maximum.asc
@@ -0,0 +1,3 @@
+# Scalar ASCII output created by CarpetIOScalar
+#
+0 0 0.255929043803271
diff --git a/test/wavetoy2/parrays3d[0].minimum.asc b/test/wavetoy2/parrays3d[0].minimum.asc
new file mode 100644
index 0000000..b779824
--- /dev/null
+++ b/test/wavetoy2/parrays3d[0].minimum.asc
@@ -0,0 +1,3 @@
+# Scalar ASCII output created by CarpetIOScalar
+#
+0 0 0.0575070381476306
diff --git a/test/wavetoy2/parrays3d[0].norm1.asc b/test/wavetoy2/parrays3d[0].norm1.asc
new file mode 100644
index 0000000..d3cbda7
--- /dev/null
+++ b/test/wavetoy2/parrays3d[0].norm1.asc
@@ -0,0 +1,3 @@
+# Scalar ASCII output created by CarpetIOScalar
+#
+0 0 0.147238120771799
diff --git a/test/wavetoy2/parrays3d[0].norm2.asc b/test/wavetoy2/parrays3d[0].norm2.asc
new file mode 100644
index 0000000..bb2d29f
--- /dev/null
+++ b/test/wavetoy2/parrays3d[0].norm2.asc
@@ -0,0 +1,3 @@
+# Scalar ASCII output created by CarpetIOScalar
+#
+0 0 0.154432972668436
diff --git a/test/wavetoy2/parrays3d[0].norm_inf.asc b/test/wavetoy2/parrays3d[0].norm_inf.asc
new file mode 100644
index 0000000..a0f9244
--- /dev/null
+++ b/test/wavetoy2/parrays3d[0].norm_inf.asc
@@ -0,0 +1,3 @@
+# Scalar ASCII output created by CarpetIOScalar
+#
+0 0 0.255929043803271
diff --git a/test/wavetoy2/parrays3d[0].x.asc b/test/wavetoy2/parrays3d[0].x.asc
new file mode 100644
index 0000000..729d4a3
--- /dev/null
+++ b/test/wavetoy2/parrays3d[0].x.asc
@@ -0,0 +1,6 @@
+# 1D ASCII output created by CarpetIOASCII
+#
+0	0 0 0	0	0.255929043803271
+0	1 0 0	0	0.20945834089817
+0	2 0 0	0	0.193918918322488
+
diff --git a/test/wavetoy2/parrays3d[0].y.asc b/test/wavetoy2/parrays3d[0].y.asc
new file mode 100644
index 0000000..282ade5
--- /dev/null
+++ b/test/wavetoy2/parrays3d[0].y.asc
@@ -0,0 +1,7 @@
+# 1D ASCII output created by CarpetIOASCII
+#
+0	0 0 0	0	0.255929043803271
+0	0 1 0	0	0.20945834089817
+0	0 2 0	0	0.193918918322488
+0	0 3 0	0	0.20945834089817
+
diff --git a/test/wavetoy2/parrays3d[0].z.asc b/test/wavetoy2/parrays3d[0].z.asc
new file mode 100644
index 0000000..e9c5ced
--- /dev/null
+++ b/test/wavetoy2/parrays3d[0].z.asc
@@ -0,0 +1,9 @@
+# 1D ASCII output created by CarpetIOASCII
+#
+0	0 0 0	0	0.255929043803271
+0	0 0 1	0	0.20945834089817
+0	0 0 2	0	0.193918918322488
+
+0	0 0 3	0	0.20945834089817
+0	0 0 4	0	0.255929043803271
+