more stuff

git-svn-id: http://svn.cactuscode.org/arrangements/CactusPUGH/PUGHSlab/trunk@11 10716dce-81a3-4424-a2c8-48026a0d3035

1 files changed, 65 insertions, 0 deletions

diff --git a/src/TestSlab2D.c b/src/TestSlab2D.c
new file mode 100644
index 0000000..0dbe3de
--- /dev/null
+++ b/src/TestSlab2D.c
@@ -0,0 +1,65 @@
+
+#include <stdio.h>
+#include <stdlib.h>
+#include <math.h>  
+
+#include "cctk.h"
+#include "cctk_Parameters.h"
+#include "cctk_Arguments.h"
+
+#include "Hyperslab.h"
+#include "CactusPUGH/PUGH/src/include/pugh.h"
+#include "CactusBase/IOUtil/src/ioGH.h"
+
+#define MAX_DIM 3
+#define MAX_FACE 6
+#define EPS 1e-10
+#define BAD -42
+
+#define ABS(a) ((a)<0 ? -(a) : (a))
+
+
+void TestSlab2D(CCTK_ARGUMENTS) {  
+  DECLARE_CCTK_ARGUMENTS
+
+    int gs[3], di1[3], di2[3]; /* global start/direction (2times) */
+  int ds[2];                   /* downsampling */
+  int ls[3], le[3];        /* local start/end */
+  int dsize[1];          /* number of points */
+  int nprocs,iproc;        /* number of procs */
+
+  int l;
+  int *index1, *index2;
+
+  int vindex, vdim, vtimelvl;
+  CCTK_REAL *dptr;
+
+  /*$vindex   = CCTK_VarIndex("hyperslab::tz");$*/
+  vindex   = CCTK_VarIndex("grid::z");
+  vtimelvl = CCTK_NumTimeLevelsFromVarI (vindex) - 1;
+  vdim     = CCTK_GroupDimFromVarI(vindex);
+
+  gs[0]=1;
+  gs[1]=1;
+  gs[2]=1;
+ 
+  di1[0]= 0;
+  di1[1]= 0;
+  di1[2]= 1;
+
+  di2[0]= 0;
+  di2[1]= 1;
+  di2[2]= 0;
+
+
+  ds[0]=1;
+  ds[1]=2;
+
+  printf("Calling CollectData2D \n");
+  CollectLocalData2D(cctkGH, vindex, vtimelvl, vdim, 
+		gs, di1, di2, ds, 145, -1, 
+		&dptr, index1, index2, dsize);
+
+  printf("TestSlab: %d %d\n",dsize[0], dsize[2]);
+
+}