Invalidate whole field for CHECK=4

boutproject · Dec 4, 2023 · e27e1bb · e27e1bb
1 parent 37e20c0
commit e27e1bb
Show file tree

Hide file tree

Showing 3 changed files with 9 additions and 3 deletions.
diff --git a/src/field/field2d.cxx b/src/field/field2d.cxx
@@ -102,7 +102,9 @@ Field2D& Field2D::allocate() {
       ny = fieldmesh->LocalNy;
     }
     data.reallocate(nx * ny);
-#if CHECK > 2
+#if CHECK > 3
+    *this = BoutNaN;
+#elif CHECK > 2
     invalidateGuards(*this);
 #endif
   } else {

diff --git a/src/field/field3d.cxx b/src/field/field3d.cxx
@@ -116,7 +116,9 @@ Field3D& Field3D::allocate() {
       nz = fieldmesh->LocalNz;
     }
     data.reallocate(nx * ny * nz);
-#if CHECK > 2
+#if CHECK > 3
+    *this = BoutNaN;
+#elif CHECK > 2
     invalidateGuards(*this);
 #endif
   } else {

diff --git a/src/field/fieldperp.cxx b/src/field/fieldperp.cxx
@@ -66,7 +66,9 @@ FieldPerp& FieldPerp::allocate() {
       nz = fieldmesh->LocalNz;
     }
     data.reallocate(nx * nz);
-#if CHECK > 2
+#if CHECK > 3
+    *this = BoutNaN;
+#elif CHECK > 2
     invalidateGuards(*this);
 #endif
   } else {