Add example and test case

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

26 files changed, 225 insertions, 0 deletions

diff --git a/test/wavetoy.par b/test/wavetoy.par
new file mode 100644
index 0000000..f1425bd
--- /dev/null
+++ b/test/wavetoy.par
@@ -0,0 +1,83 @@
+!DESC "Test InterpToArray with WaveToy initial data"
+
+ActiveThorns = "
+        Boundary
+        CartGrid3D
+        CoordBase
+        IDScalarWaveC
+        IOASCII
+        IOBasic
+        IOUtil
+        InterpToArray
+        LocalInterp
+        LocalReduce
+        PUGH
+        PUGHInterp
+        PUGHReduce
+        PUGHSlab
+        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::interpolator_name = "uniform cartesian"
+
+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"
+
+IOBasic::outScalar_reductions = "norm1 norm2 minimum maximum norm_inf"
+IOBasic::outScalar_every      = 1
+IOBasic::outScalar_vars       = "
+        wavetoy::phi
+        InterpToArray::arrays2d[0]
+        InterpToArray::parrays3d[0]
+"
+
+IOASCII::out1D_every = 1
+IOASCII::out1D_vars  = "
+        wavetoy::phi
+        InterpToArray::scalars
+        InterpToArray::arrays2d
+        InterpToArray::parrays3d
+"
diff --git a/test/wavetoy/arrays2d[0]_maximum.xg b/test/wavetoy/arrays2d[0]_maximum.xg
new file mode 100644
index 0000000..5267316
--- /dev/null
+++ b/test/wavetoy/arrays2d[0]_maximum.xg
@@ -0,0 +1,2 @@
+"arrays2d[0] v time
+0.0000000000000	0.2559290438033
diff --git a/test/wavetoy/arrays2d[0]_minimum.xg b/test/wavetoy/arrays2d[0]_minimum.xg
new file mode 100644
index 0000000..a319a06
--- /dev/null
+++ b/test/wavetoy/arrays2d[0]_minimum.xg
@@ -0,0 +1,2 @@
+"arrays2d[0] v time
+0.0000000000000	0.0981997271448
diff --git a/test/wavetoy/arrays2d[0]_norm1.xg b/test/wavetoy/arrays2d[0]_norm1.xg
new file mode 100644
index 0000000..d727340
--- /dev/null
+++ b/test/wavetoy/arrays2d[0]_norm1.xg
@@ -0,0 +1,2 @@
+"arrays2d[0] v time
+0.0000000000000	0.1622239184695
diff --git a/test/wavetoy/arrays2d[0]_norm2.xg b/test/wavetoy/arrays2d[0]_norm2.xg
new file mode 100644
index 0000000..ccfb470
--- /dev/null
+++ b/test/wavetoy/arrays2d[0]_norm2.xg
@@ -0,0 +1,2 @@
+"arrays2d[0] v time
+0.0000000000000	0.1685259557890
diff --git a/test/wavetoy/arrays2d[0]_norm_inf.xg b/test/wavetoy/arrays2d[0]_norm_inf.xg
new file mode 100644
index 0000000..5267316
--- /dev/null
+++ b/test/wavetoy/arrays2d[0]_norm_inf.xg
@@ -0,0 +1,2 @@
+"arrays2d[0] v time
+0.0000000000000	0.2559290438033
diff --git a/test/wavetoy/arrays2d[0]_x_[0].xg b/test/wavetoy/arrays2d[0]_x_[0].xg
new file mode 100644
index 0000000..273e4dc
--- /dev/null
+++ b/test/wavetoy/arrays2d[0]_x_[0].xg
@@ -0,0 +1,6 @@
+
+
+"Time = 0.0000000000000
+0		0.2559290438033
+1		0.2094583408982
+2		0.1939189183225
diff --git a/test/wavetoy/arrays2d[0]_y_[0].xg b/test/wavetoy/arrays2d[0]_y_[0].xg
new file mode 100644
index 0000000..06c162c
--- /dev/null
+++ b/test/wavetoy/arrays2d[0]_y_[0].xg
@@ -0,0 +1,7 @@
+
+
+"Time = 0.0000000000000
+0		0.2559290438033
+1		0.2044903680418
+2		0.1738539578662
+3		0.1643877561027
diff --git a/test/wavetoy/parrays3d[0]_3D_diagonal.xg b/test/wavetoy/parrays3d[0]_3D_diagonal.xg
new file mode 100644
index 0000000..0c61eb0
--- /dev/null
+++ b/test/wavetoy/parrays3d[0]_3D_diagonal.xg
@@ -0,0 +1,6 @@
+
+
+"Time = 0.0000000000000
+0		0.2559290438033
+1		0.1111550275597
+2		0.0575070381476
diff --git a/test/wavetoy/parrays3d[0]_maximum.xg b/test/wavetoy/parrays3d[0]_maximum.xg
new file mode 100644
index 0000000..c29e9fc
--- /dev/null
+++ b/test/wavetoy/parrays3d[0]_maximum.xg
@@ -0,0 +1,2 @@
+"parrays3d[0] v time
+0.0000000000000	0.2559290438033
diff --git a/test/wavetoy/parrays3d[0]_minimum.xg b/test/wavetoy/parrays3d[0]_minimum.xg
new file mode 100644
index 0000000..c814a33
--- /dev/null
+++ b/test/wavetoy/parrays3d[0]_minimum.xg
@@ -0,0 +1,2 @@
+"parrays3d[0] v time
+0.0000000000000	0.0575070381476
diff --git a/test/wavetoy/parrays3d[0]_norm1.xg b/test/wavetoy/parrays3d[0]_norm1.xg
new file mode 100644
index 0000000..5a96903
--- /dev/null
+++ b/test/wavetoy/parrays3d[0]_norm1.xg
@@ -0,0 +1,2 @@
+"parrays3d[0] v time
+0.0000000000000	0.1472381207718
diff --git a/test/wavetoy/parrays3d[0]_norm2.xg b/test/wavetoy/parrays3d[0]_norm2.xg
new file mode 100644
index 0000000..491c016
--- /dev/null
+++ b/test/wavetoy/parrays3d[0]_norm2.xg
@@ -0,0 +1,2 @@
+"parrays3d[0] v time
+0.0000000000000	0.1544329726684
diff --git a/test/wavetoy/parrays3d[0]_norm_inf.xg b/test/wavetoy/parrays3d[0]_norm_inf.xg
new file mode 100644
index 0000000..c29e9fc
--- /dev/null
+++ b/test/wavetoy/parrays3d[0]_norm_inf.xg
@@ -0,0 +1,2 @@
+"parrays3d[0] v time
+0.0000000000000	0.2559290438033
diff --git a/test/wavetoy/parrays3d[0]_x_[0][0].xg b/test/wavetoy/parrays3d[0]_x_[0][0].xg
new file mode 100644
index 0000000..273e4dc
--- /dev/null
+++ b/test/wavetoy/parrays3d[0]_x_[0][0].xg
@@ -0,0 +1,6 @@
+
+
+"Time = 0.0000000000000
+0		0.2559290438033
+1		0.2094583408982
+2		0.1939189183225
diff --git a/test/wavetoy/parrays3d[0]_y_[0][0].xg b/test/wavetoy/parrays3d[0]_y_[0][0].xg
new file mode 100644
index 0000000..3ec649e
--- /dev/null
+++ b/test/wavetoy/parrays3d[0]_y_[0][0].xg
@@ -0,0 +1,7 @@
+
+
+"Time = 0.0000000000000
+0		0.2559290438033
+1		0.2094583408982
+2		0.1939189183225
+3		0.2094583408982
diff --git a/test/wavetoy/parrays3d[0]_z_[0][0].xg b/test/wavetoy/parrays3d[0]_z_[0][0].xg
new file mode 100644
index 0000000..9454b68
--- /dev/null
+++ b/test/wavetoy/parrays3d[0]_z_[0][0].xg
@@ -0,0 +1,8 @@
+
+
+"Time = 0.0000000000000
+0		0.2559290438033
+1		0.2094583408982
+2		0.1939189183225
+3		0.2094583408982
+4		0.2559290438033
diff --git a/test/wavetoy/phi_3D_diagonal.xg b/test/wavetoy/phi_3D_diagonal.xg
new file mode 100644
index 0000000..7261bd4
--- /dev/null
+++ b/test/wavetoy/phi_3D_diagonal.xg
@@ -0,0 +1,18 @@
+
+
+"Time = 0.0000000000000
+-6.7549981495186		0.6751235936513
+-5.7157676649773		0.9367420144351
+-4.6765371804360		0.9867432223551
+-3.6373066958946		0.7891067225602
+-2.5980762113533		0.4790874265193
+-1.5588457268120		0.2208214420284
+-0.5196152422707		0.0772707537752
+0.5196152422707		0.0772707537752
+1.5588457268120		0.2208214420284
+2.5980762113533		0.4790874265193
+3.6373066958946		0.7891067225602
+4.6765371804360		0.9867432223551
+5.7157676649773		0.9367420144351
+6.7549981495186		0.6751235936513
+7.7942286340599		0.3693975497431
diff --git a/test/wavetoy/phi_maximum.xg b/test/wavetoy/phi_maximum.xg
new file mode 100644
index 0000000..09ef046
--- /dev/null
+++ b/test/wavetoy/phi_maximum.xg
@@ -0,0 +1,2 @@
+"phi v time
+0.0000000000000	0.9999198827161
diff --git a/test/wavetoy/phi_minimum.xg b/test/wavetoy/phi_minimum.xg
new file mode 100644
index 0000000..dfff009
--- /dev/null
+++ b/test/wavetoy/phi_minimum.xg
@@ -0,0 +1,2 @@
+"phi v time
+0.0000000000000	0.0772707537752
diff --git a/test/wavetoy/phi_norm1.xg b/test/wavetoy/phi_norm1.xg
new file mode 100644
index 0000000..8c182d5
--- /dev/null
+++ b/test/wavetoy/phi_norm1.xg
@@ -0,0 +1,2 @@
+"phi v time
+0.0000000000000	0.8148785107206
diff --git a/test/wavetoy/phi_norm2.xg b/test/wavetoy/phi_norm2.xg
new file mode 100644
index 0000000..b618b94
--- /dev/null
+++ b/test/wavetoy/phi_norm2.xg
@@ -0,0 +1,2 @@
+"phi v time
+0.0000000000000	0.8424968585241
diff --git a/test/wavetoy/phi_norm_inf.xg b/test/wavetoy/phi_norm_inf.xg
new file mode 100644
index 0000000..09ef046
--- /dev/null
+++ b/test/wavetoy/phi_norm_inf.xg
@@ -0,0 +1,2 @@
+"phi v time
+0.0000000000000	0.9999198827161
diff --git a/test/wavetoy/phi_x_[7][7].xg b/test/wavetoy/phi_x_[7][7].xg
new file mode 100644
index 0000000..d993f96
--- /dev/null
+++ b/test/wavetoy/phi_x_[7][7].xg
@@ -0,0 +1,18 @@
+
+
+"Time = 0.0000000000000
+-3.9000000000000		0.8624762718586
+-3.3000000000000		0.6998161919677
+-2.7000000000000		0.5192102694005
+-2.1000000000000		0.3529089299372
+-1.5000000000000		0.2208214420284
+-0.9000000000000		0.1292593546971
+-0.3000000000000		0.0772707537752
+0.3000000000000		0.0772707537752
+0.9000000000000		0.1292593546971
+1.5000000000000		0.2208214420284
+2.1000000000000		0.3529089299372
+2.7000000000000		0.5192102694005
+3.3000000000000		0.6998161919677
+3.9000000000000		0.8624762718586
+4.5000000000000		0.9710346043534
diff --git a/test/wavetoy/phi_y_[7][7].xg b/test/wavetoy/phi_y_[7][7].xg
new file mode 100644
index 0000000..d993f96
--- /dev/null
+++ b/test/wavetoy/phi_y_[7][7].xg
@@ -0,0 +1,18 @@
+
+
+"Time = 0.0000000000000
+-3.9000000000000		0.8624762718586
+-3.3000000000000		0.6998161919677
+-2.7000000000000		0.5192102694005
+-2.1000000000000		0.3529089299372
+-1.5000000000000		0.2208214420284
+-0.9000000000000		0.1292593546971
+-0.3000000000000		0.0772707537752
+0.3000000000000		0.0772707537752
+0.9000000000000		0.1292593546971
+1.5000000000000		0.2208214420284
+2.1000000000000		0.3529089299372
+2.7000000000000		0.5192102694005
+3.3000000000000		0.6998161919677
+3.9000000000000		0.8624762718586
+4.5000000000000		0.9710346043534
diff --git a/test/wavetoy/phi_z_[7][7].xg b/test/wavetoy/phi_z_[7][7].xg
new file mode 100644
index 0000000..d993f96
--- /dev/null
+++ b/test/wavetoy/phi_z_[7][7].xg
@@ -0,0 +1,18 @@
+
+
+"Time = 0.0000000000000
+-3.9000000000000		0.8624762718586
+-3.3000000000000		0.6998161919677
+-2.7000000000000		0.5192102694005
+-2.1000000000000		0.3529089299372
+-1.5000000000000		0.2208214420284
+-0.9000000000000		0.1292593546971
+-0.3000000000000		0.0772707537752
+0.3000000000000		0.0772707537752
+0.9000000000000		0.1292593546971
+1.5000000000000		0.2208214420284
+2.1000000000000		0.3529089299372
+2.7000000000000		0.5192102694005
+3.3000000000000		0.6998161919677
+3.9000000000000		0.8624762718586
+4.5000000000000		0.9710346043534