diff --git a/src/Utilities/iterateIdealH.cc b/src/Utilities/iterateIdealH.cc index 9b4a1a66c..3746da220 100644 --- a/src/Utilities/iterateIdealH.cc +++ b/src/Utilities/iterateIdealH.cc @@ -102,18 +102,18 @@ iterateIdealH(DataBase& dataBase, } } - // Check if we're using ASPH and radialOnly. If so we'll switch to fixShape for the iteration. - auto radialOnly = false; - ASPHSmoothingScale* asphPkg = nullptr; - for (auto* pkg: packages) { - asphPkg = dynamic_cast*>(pkg); - if (asphPkg != nullptr and asphPkg->radialOnly()) { - radialOnly = true; - asphPkg->radialOnly(false); - asphPkg->fixShape(true); - break; - } - } + // // Check if we're using ASPH and radialOnly. If so we'll switch to fixShape for the iteration. + // auto radialOnly = false; + // ASPHSmoothingScale* asphPkg = nullptr; + // for (auto* pkg: packages) { + // asphPkg = dynamic_cast*>(pkg); + // if (asphPkg != nullptr and asphPkg->radialOnly()) { + // radialOnly = true; + // asphPkg->radialOnly(false); + // asphPkg->fixShape(true); + // break; + // } + // } // Build a list of flags to indicate which nodes have been completed. auto flagNodeDone = dataBase.newFluidFieldList(0, "node completed"); @@ -265,12 +265,12 @@ iterateIdealH(DataBase& dataBase, for (auto* boundaryPtr: range(boundaries.begin(), boundaries.end())) boundaryPtr->applyFieldListGhostBoundary(m); for (auto* boundaryPtr: range(boundaries.begin(), boundaries.end())) boundaryPtr->finalizeGhostBoundary(); - // Restore ASPH radialOnly choice if necessary - if (radialOnly) { - CHECK(asphPkg != nullptr); - asphPkg->radialOnly(true); - asphPkg->fixShape(false); - } + // // Restore ASPH radialOnly choice if necessary + // if (radialOnly) { + // CHECK(asphPkg != nullptr); + // asphPkg->radialOnly(true); + // asphPkg->fixShape(false); + // } // Report the final timing. const auto t1 = clock();