aboutsummaryrefslogtreecommitdiff
diff options
context:
space:
mode:
authorErik Schnetter <schnetter@aei.mpg.de>2005-01-27 10:47:00 +0000
committerErik Schnetter <schnetter@aei.mpg.de>2005-01-27 10:47:00 +0000
commitff45e76f759e8e6d666e5867654fb148ab8ad7c2 (patch)
tree3f7e171ff0a4c95a567077f69a0c8c749b7a0d85
parent2d14b56f07f70752181e24039b57c673b5a2b0d8 (diff)
Carpet: set cGH identity
darcs-hash:20050127104759-891bb-5a7b05cb10a33dbdaa5afe9bf854ada9838b9161.gz
-rw-r--r--Carpet/Carpet/src/SetupGH.cc24
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);