#797: N13CarpetRegrid211snap_coarseE does not hold after being enforced in cell centering ---------------------+------------------------------------------------------ Reporter: rhaas | Owner: rhaas Type: defect | Status: review Priority: major | Milestone: Component: Carpet | Version: Resolution: | Keywords: ---------------------+------------------------------------------------------
Comment (by rhaas):
With the patch applied as is, it aborts with an assert(): {{{ arrangements/Carpet/CarpetRegrid2/src/property.cc:286: ibset CarpetRegrid2::snap_coarse::snapped_regions(const gh&, const dh&, const CarpetRegrid2::level_boundary&, const std::vector<bboxset<int, 3> >&, int): Assertion `original_region <= base' failed }}} if I then comment out the assert(), it no longer fails but the grid structure created is incorrect (only a single box on refinement level 1 when there should be two).