From 575b95f668a005466d8eeb1ecedc565df3adbb72 Mon Sep 17 00:00:00 2001 From: lanfer Date: Fri, 5 May 2000 08:45:59 +0000 Subject: more stuff git-svn-id: http://svn.cactuscode.org/arrangements/CactusPUGH/PUGHSlab/trunk@11 10716dce-81a3-4424-a2c8-48026a0d3035 --- src/TestSlab2D.c | 65 ++++++++++++++++++++++++++++++++++++++++++++++++++++++++ 1 file changed, 65 insertions(+) create mode 100644 src/TestSlab2D.c (limited to 'src/TestSlab2D.c') 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 +#include +#include + +#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]); + +} -- cgit v1.2.3