modify/remove the assertions about no hidden dimension since this cas…

…e is working in ERF
AMReX-Codes · Oct 3, 2024 · 7fbe6e0 · 7fbe6e0
1 parent 8349789
commit 7fbe6e0
Showing 1 changed file with 5 additions and 3 deletions.
diff --git a/Src/LinearSolvers/MLMG/AMReX_MLCellLinOp.H b/Src/LinearSolvers/MLMG/AMReX_MLCellLinOp.H
@@ -536,7 +536,11 @@ MLCellLinOpT<MF>::setLevelBC (int amrlev, const MF* a_levelbcdata, const MF* rob
     {
         if (this->needsCoarseDataForBC())
         {
-            AMREX_ALWAYS_ASSERT(!this->hasHiddenDimension());
+            // AMREX_ALWAYS_ASSERT(!this->hasHiddenDimension());
+            if (this->hasHiddenDimension()) {
+                int hidden_dir = this->hiddenDirection();
+                AMREX_ALWAYS_ASSERT(this->m_coarse_data_crse_ratio[hidden_dir] == 1);
+            }
             br_ref_ratio = this->m_coarse_data_crse_ratio.allGT(0) ? this->m_coarse_data_crse_ratio : IntVect(2);
             if (m_crse_sol_br[amrlev] == nullptr && br_ref_ratio.allGT(0))
             {
@@ -1946,8 +1950,6 @@ MLCellLinOpT<MF>::computeVolInv () const
         m_volinv[amrlev].resize(this->NMGLevels(amrlev));
     }
 
-    AMREX_ASSERT(this->m_coarse_fine_bc_type == LinOpBCType::Dirichlet || ! this->hasHiddenDimension());
-
     // We don't need to compute for every level
 
     auto f = [&] (int amrlev, int mglev) {