[Users] A Carpet error occurred in an unequal mass bns simulation
白济民
beki-cat at sjtu.edu.cn
Mon Aug 15 11:21:37 CDT 2022
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
More information about the Users
mailing list