[Users] A Carpet error occurred in an unequal mass bns simulation

[Roland Haas]

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:
> ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ 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>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 at einsteintoolkit.org
> https://urldefense.com/v3/__http://lists.einsteintoolkit.org/mailman/listinfo/users__;!!DZ3fjg!7J7ZLqS8Exs7ZbQTYF0L9KyJpqqXT7bswbvhND87I74UxEOlXSH9N9G4josTAHijv2Q584yrgN0eIFmqKYx3QQ$ 


-- 
My email is as private as my paper mail. I therefore support encrypting
and signing email messages. Get my PGP key from http://pgp.mit.edu .
-------------- next part --------------
A non-text attachment was scrubbed...
Name: not available
Type: application/pgp-signature
Size: 833 bytes
Desc: OpenPGP digital signature
Url : http://lists.einsteintoolkit.org/pipermail/users/attachments/20220818/493049de/attachment.bin