diff options
author | Erik Schnetter <schnetter@aei.mpg.de> | 2005-01-27 10:47:00 +0000 |
---|---|---|
committer | Erik Schnetter <schnetter@aei.mpg.de> | 2005-01-27 10:47:00 +0000 |
commit | ff45e76f759e8e6d666e5867654fb148ab8ad7c2 (patch) | |
tree | 3f7e171ff0a4c95a567077f69a0c8c749b7a0d85 | |
parent | 2d14b56f07f70752181e24039b57c673b5a2b0d8 (diff) |
Carpet: set cGH identity
darcs-hash:20050127104759-891bb-5a7b05cb10a33dbdaa5afe9bf854ada9838b9161.gz
-rw-r--r-- | Carpet/Carpet/src/SetupGH.cc | 24 |
1 files changed, 22 insertions, 2 deletions
diff --git a/Carpet/Carpet/src/SetupGH.cc b/Carpet/Carpet/src/SetupGH.cc index 12f036208..22548fb54 100644 --- a/Carpet/Carpet/src/SetupGH.cc +++ b/Carpet/Carpet/src/SetupGH.cc @@ -361,6 +361,8 @@ namespace Carpet { Waypoint ("Setting up the grid hierarchy"); + cgh->identity = strdup (""); + // Processor information Output ("Carpet is running on %d processors", CCTK_nProcs(cgh)); @@ -529,7 +531,7 @@ namespace Carpet { // Adapt for convergence level rvect const spacing - = base_spacing * ipow (CCTK_REAL (convergence_factor), basemglevel); + = base_spacing * ipow ((CCTK_REAL) convergence_factor, basemglevel); // Calculate global number of grid points // SW note this and other examples break the encapsulation of vect class @@ -791,6 +793,24 @@ namespace Carpet { int ierr = CCTK_GroupData (group, &gp); assert (!ierr); + if (gp.compact) { + char * const groupname = CCTK_GroupName (group); + CCTK_VWarn (0, __LINE__, __FILE__, CCTK_THORNSTRING, + "The group \"%s\" has COMPACT=1. Compact groups are not yet supported", + groupname); + free (groupname); + } + +#if 0 + if (gp.contiguous) { + char * const groupname = CCTK_GroupName (group); + CCTK_VWarn (0, __LINE__, __FILE__, CCTK_THORNSTRING, + "The group \"%s\" has CONTIGUOUS=1. Contiguous groups are not yet supported", + groupname); + free (groupname); + } +#endif + switch (gp.grouptype) { case CCTK_GF: { @@ -893,7 +913,7 @@ namespace Carpet { + 2*ghostsizes[d]; assert (sizes[d] >= 0); } - + const ivect alb(0); const ivect aub(sizes-1); const ivect astr(1); |