[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