Hello Jimmy,
Sorry for the delay in responding.
The assert error you see (stride == other.stride) is triggered if Carpet tries to copy (but not interpolate) data between different refinement levels.
This should never happen, and the assert is there to catch bugs.
Unfortunately without more context it is not easy to say what may be the root cause.
The first step in trying to deduce the cause would be to take a look at the backtrace (in the files backtrace.*.txt) to see the call tree that lead to the assert failure.
Do you those files and could attach them to an email to this list?
The only thing in your parfile that is slightly unusual (but should work and I do not see how it would directly cause issues) is use of the
CarpetRegrid2::radius_x_3[5] = 18.0 CarpetRegrid2::radius_y_3[5] = 18.0 CarpetRegrid2::radius_z_3[5] = 15.0
rather than a single radius.
When does the error actually happen? Right at the beginning (and in that case, would you be fine to share the full parfile) or only after a long simulation time and to reproducible?
Yours, Roland
Hi, An error reported by the Thorn Carpet occurred when I run a bns merger simulation in the unequal mass case when applying a certain grid structure, which reads:
Cactus/arrangements/Carpet/CarpetLib/src/bboxset2.hh:262: CarpetLib::bboxset2::bboxset<T, D> CarpetLib::bboxset2::bboxset<T, D>::binary_operator(const F&, const CarpetLib::bboxset2::bboxset<T, D>D>&) const [with F = CarpetLib::bboxset2::bboxset<T, D>D>D>::operator&(const CarpetLib::bboxset2::bboxset<T, D>&) const D>D>D>[with T = int; int D = 3]::<lambda(const bboxset1&, const D>D>D>bboxset1&)>; T = int; int D = 3]: Assertion `all(stride == D>D>D>other.stride)' failed. D>D>D>~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ However, I D>D>D>don't think the grid structure I applied in my simulation is D>D>D>peculiar:CoordBase::xmin = -400 CoordBase::xmax = 400 CoordBase::ymin = -400 CoordBase::ymax = 400 CoordBase::zmin = 0 CoordBase::zmax = 400
CoordBase::spacing = "numcells" CoordBase::ncells_x = 40 CoordBase::ncells_y = 40 CoordBase::ncells_z = 20
...
CarpetRegrid2::num_centres = 3 CarpetRegrid2::regrid_every = 128 CarpetRegrid2::snap_to_coarse = "yes" CarpetRegrid2::freeze_unaligned_levels = "yes" CarpetRegrid2::freeze_unaligned_parent_levels = "yes"
# ----------- Region 1 -------------------- CarpetRegrid2::active_1 = "yes" CarpetRegrid2::num_levels_1 = 7 CarpetRegrid2::position_x_1 = -17.4316 CarpetRegrid2::position_y_1 = 0.0 CarpetRegrid2::radius_1[1] = 240.0 CarpetRegrid2::radius_1[2] = 120.0 CarpetRegrid2::radius_1[3] = 60.0 CarpetRegrid2::radius_1[4] = 30.0 CarpetRegrid2::radius_1[5] = 20.0 CarpetRegrid2::radius_1[6] = 15.0
# ----------- Region 2 -------------------- CarpetRegrid2::active_2 = "yes" CarpetRegrid2::num_levels_2 = 7 CarpetRegrid2::position_x_2 = 23.1912 CarpetRegrid2::position_y_2 = -0.0 CarpetRegrid2::radius_2[1] = 240.0 CarpetRegrid2::radius_2[2] = 120.0 CarpetRegrid2::radius_2[3] = 60.0 CarpetRegrid2::radius_2[4] = 30.0 CarpetRegrid2::radius_2[5] = 18.0 CarpetRegrid2::radius_2[6] = 14.0
# ----------- Region 3 -------------------- CarpetRegrid2::active_3 = "yes" CarpetRegrid2::num_levels_3 = 3 CarpetRegrid2::position_x_3 = 0
CarpetRegrid2::radius_3[1] = 240.0 CarpetRegrid2::radius_3[2] = 120.0 CarpetRegrid2::radius_3[3] = 66.0 CarpetRegrid2::radius_3[4] = 33.0 CarpetRegrid2::radius_x_3[5] = 18.0 CarpetRegrid2::radius_y_3[5] = 18.0 CarpetRegrid2::radius_z_3[5] = 15.0 CarpetRegrid2::radius_3[6] = 12.0 CarpetRegrid2::radius_3[7] = 8.0 # -----------------------------------------
Though the resulotion is relatively low for test. I am wondering what the error indicates and I'm looking forward to addressing it with your help. Thank you in advance, Jimmy _______________________________________________ Users mailing list Users@einsteintoolkit.org https://urldefense.com/v3/__http://lists.einsteintoolkit.org/mailman/listinfo/users__;!!DZ3fjg!7J7ZLqS8Exs7ZbQTYF0L9KyJpqqXT7bswbvhND87I74UxEOlXSH9N9G4josTAHijv2Q584yrgN0eIFmqKYx3QQ$