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: ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ exec: 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>&) const [with F = CarpetLib::bboxset2::bboxset<T, D>::operator&(const CarpetLib::bboxset2::bboxset<T, D>&) const [with T = int; int D = 3]::<lambda(const bboxset1&, const bboxset1&)>; T = int; int D = 3]: Assertion `all(stride == other.stride)' failed. ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
However, I don't think the grid structure I applied in my simulation is 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