diff options
| author | tradke <tradke@94b1c47f-dcfd-45ef-a468-0854c0e9e350> | 2002-10-23 09:08:31 +0000 |
|---|---|---|
| committer | tradke <tradke@94b1c47f-dcfd-45ef-a468-0854c0e9e350> | 2002-10-23 09:08:31 +0000 |
| commit | 2c294a2c35524f87c480668a04469cb7355adc2f (patch) | |
| tree | f5597d04d498e8962957034091332825f7230ff2 /src | |
| parent | bdc9e64dc21109384094c8281022037471080d9b (diff) | |
Fix to output complex variables in xgraph format.
This closes PR CactusBase/1264.
git-svn-id: http://svn.cactuscode.org/arrangements/CactusBase/IOASCII/trunk@134 94b1c47f-dcfd-45ef-a468-0854c0e9e350
Diffstat (limited to 'src')
| -rw-r--r-- | src/Write1D.c | 24 |
1 files changed, 15 insertions, 9 deletions
diff --git a/src/Write1D.c b/src/Write1D.c index 46cb6b9..f4d73d5 100644 --- a/src/Write1D.c +++ b/src/Write1D.c @@ -23,6 +23,7 @@ #include "cctk.h" #include "cctk_Parameters.h" #include "Hyperslab.h" +#include "util_Table.h" #include "CactusBase/IOUtil/src/ioGH.h" #include "CactusBase/IOUtil/src/ioutil_AdvertisedFiles.h" #include "ioASCIIGH.h" @@ -46,20 +47,21 @@ CCTK_FILEVERSION(CactusBase_IOASCII_Write1D_c) const CCTK_REAL *_coord = (const CCTK_REAL *) hdata[1]; \ \ \ - for (_h = hstart; _h < hstride * hsize; _h += hstride) \ + for (_h = hstart; _h < hsize; _h++) \ { \ if (_coord) \ { \ - fprintf (file, fmt, (double) _coord[_h], (c_type) _hdata[_h]); \ + fprintf (file, fmt, (double) _coord[_h], \ + (c_type) _hdata[_h * hstride]); \ } \ else \ { \ - fprintf (file, fmt, _h, (c_type) _hdata[_h]); \ + fprintf (file, fmt, _h, (c_type) _hdata[_h * hstride]); \ } \ \ if (is_complex && hstride == 1) \ { \ - fprintf (file, format[2], (c_type) _hdata[_h+1]); \ + fprintf (file, format[2], (c_type) _hdata[_h + 1]); \ } \ fputc ('\n', file); \ } \ @@ -128,7 +130,7 @@ int IOASCII_Write1D (const cGH *GH, int vindex, const char *alias) { asciiioGH *myGH; int do_dir[4], coord_index[3]; - int i, dir, myproc, groupindex, have_coords; + int i, dir, myproc, groupindex, have_coords, mapping, table; #if 0 int lower, upper; #endif @@ -142,7 +144,7 @@ int IOASCII_Write1D (const cGH *GH, int vindex, const char *alias) CCTK_REAL offset; CCTK_REAL coord_lower[3]; CCTK_INT *origin, *direction; - CCTK_INT mapping, hsize, extent; + CCTK_INT hsize, extent; CCTK_INT vindices[2]; void *hdata[2]; DECLARE_CCTK_PARAMETERS @@ -330,16 +332,20 @@ fprintf (stderr, "%d/%d\n", origin[dir], extent); memset (origin, 0, gdata.dim * sizeof (CCTK_INT)); } - - + table = out_physical_boundaries ? + Util_TableCreateFromString ("with_ghostzones = 1") : -1; mapping = Hyperslab_DefineGlobalMappingByIndex (GH, vindex, 1, direction, origin, &extent, NULL, /* downsample */ - -1, /* table handle */ + table, NULL /* conversion fn */, &hsize); + if (table >= 0) + { + Util_TableDestroy (table); + } if (mapping < 0) { CCTK_VWarn (1, __LINE__, __FILE__, CCTK_THORNSTRING, |