Formating

CMakeLists.txt

@@ -86,4 +86,5 @@ add_subdirectory(src)
 if (AthenaPK_ENABLE_TESTING)
   include(CTest)
   add_subdirectory(tst/regression)
-endif()
+endif()
+

src/eos/adiabatic_glmmhd.cpp

@@ -51,12 +51,6 @@ void AdiabaticGLMMHDEOS::ConservedToPrimitive(MeshData<Real> *md) const {
         auto &prim = prim_pack(b);
         // auto &nu = entropy_pack(b);
 
-        // Getting the global indexing
-
-        auto pmb = md->GetBlockData(b)->GetBlockPointer();
-        auto pm = pmb->pmy_mesh;
-        auto hydro_pkg = pmb->packages.Get("Hydro");
-
         return this_on_device.ConsToPrim(cons, prim, nhydro, nscalars, k, j, i);
       });
 }

src/eos/adiabatic_glmmhd.hpp

@@ -93,9 +93,9 @@ class AdiabaticGLMMHDEOS : public EquationOfState {
 
     // Let's apply floors explicitly, i.e., by default floor will be disabled (<=0)
     // and the code will fail if a negative density is encountered.
-    // PARTHENON_REQUIRE(u_d > 0.0 || density_floor_ > 0.0,
-    //                  "Got negative density. Consider enabling first-order flux "
-    //                  "correction or setting a reasonable density floor.");
+    PARTHENON_REQUIRE(u_d > 0.0 || density_floor_ > 0.0,
+                      "Got negative density. Consider enabling first-order flux "
+                      "correction or setting a reasonable density floor.");
 
     // apply density floor, without changing momentum or energy
     u_d = (u_d > density_floor_) ? u_d : density_floor_;
@@ -134,10 +134,10 @@ class AdiabaticGLMMHDEOS : public EquationOfState {
 
     // Let's apply floors explicitly, i.e., by default floor will be disabled (<=0)
     // and the code will fail if a negative pressure is encountered.
-    // PARTHENON_REQUIRE(w_p > 0.0 || pressure_floor_ > 0.0 || e_floor_ > 0.0,
-    //                  "Got negative pressure. Consider enabling first-order flux "
-    //                  "correction or setting a reasonable pressure or temperature
-    //                  floor.");
+    PARTHENON_REQUIRE(
+        w_p > 0.0 || pressure_floor_ > 0.0 || e_floor_ > 0.0,
+        "Got negative pressure. Consider enabling first-order flux "
+        "correction or setting a reasonable pressure or temperature floor.");
 
     // Pressure floor (if present) takes precedence over temperature floor
     if ((pressure_floor_ > 0.0) && (w_p < pressure_floor_)) {