SAMRAI
SAMRAI copied to clipboard
Multiblock in restart
If it is not a restart case, TimeRefinementIntegrator::initializeHierarchy will call PatchHierarchy::makeCoarsestLevel and then call adjustMultiblockPatchLevelBoundaries at the coarsest level, and Patch will have correct PatchGeometry in all following simulations. In restart case, adjustMultiblockPatchLevelBoundaries is never called during initialization. Patch will have incorrect PatchGeometry in multiblock case (all face, edge, node boundaries will be regarded as physical boundaries), unless I called adjustMultiblockPatchLevelBoundaries in my own code from coarsest level to finest level.