[Users] problems of reproducing the simulation of binary neutron star with different mass

[林家暉]

Dear Roland,
Thank you for your reply.
The attached file is .err file corresponding to the simulation of binary star with different central enthalpy.  Can you have a look at it ?Thank you.
For not able to track multiple stars, does it occur in binary blackhole case(http://einsteintoolkit.org/gallery/bbh/index.html) ?
Best regards,
Chia-Hui
________________________________
寄件者: Roland Haas <rhaas at illinois.edu>
寄件日期: 2018年5月13日 下午 09:28:51
收件者: 林家暉
副本: Einstein Toolkit Users
主旨: Re: [Users] problems of reproducing the simulation of binary neutron star with different mass

Hello Chia-Hui Lin,

changing the central enthalpy seems like the right thing to do to
change masses in LORENE since it prescribes the central enthalpy of
each star.

The backtrace that you showed and the error message is unfortunately
not detailed enough to make a guess as to what could have gone wrong.
This output is produced by a call to CCTK_Error() however the
CCTK_Error function takes an error string as an argument which is not
present in the lines that you have copied. It is likely in the *.err
file that the simulation produced and you should try and include that
file as an attachment in your email.

If you have non-equal mass stars then you have to remove the 180 degree
rotation symmetry from the parameter files. Ie you have to change the
lines (these are scattered throughout throughout the par file)

ActiveThorns = "RotatingSymmetry180 ReflectionSymmetry"
CoordBase::xmin =    0.00
CoordBase::boundary_shiftout_x_lower = 1
CarpetRegrid2::symmetry_rotating180    = "yes"
Hydro_Analysis::Hydro_Analysis_rho_max_loc_only_positive_x = "true"
Hydro_Analysis::Hydro_Analysis_average_multiple_maxima_locations = "yes"

to

ActiveThorns = "ReflectionSymmetry"
CoordBase::xmin = 396.00
CoordBase::boundary_shiftout_x_lower = 0
CarpetRegrid2::symmetry_rotating180    = "no"
Hydro_Analysis::Hydro_Analysis_rho_max_loc_only_positive_x = "false"
Hydro_Analysis::Hydro_Analysis_average_multiple_maxima_locations = "no"

Please note though that currently the thorn NSTracker and HydroAnalysis *cannot* properly track multiple stars. It makes the hard-coded assumption that if one of the stars is at (x,y,z) the other one is at (-x,-y,z) which is only true for equal mass neutron stars. This will mean that one of the stars will eventually leave the box that is supposed to track it.

Yours,
Roland

> Dear whom it may concern,
>
> I am trying to reproduce the simulation of binary neutron star (http://einsteintoolkit.org/gallery/bns/index.html) with different mass. I use Lorene code to produce the initial data . The parameter files I used is locate in /Lorene/Codes/Bin_star/Parameters/Polytrope_irrot/Test_GR . And I changed  central enthalpy of each star in par_init.d to get different mass (Is there other way I can get different mass?). And then I simulate with the ET. However I met some problems:
>
> 1.I requested 48 cores with 10 hours to run the code , however the code terminated after running ~4 hours. And the end of .err file reads:
>
>
>
> cactus_sim: /global/u1/c/chlin/Cactus/arrangements/Carpet/Carpet/src/helpers.cc:275: int Carpet::Abort(const _cGH *, int): Assertion `0' failed.
>
> cactus_sim: /global/u1/c/chlin/Cactus/arrangements/Carpet/Carpet/src/helpers.cc:275: int Carpet::Abort(const _cGH *, int): Assertion `0' failed.
>
> Rank 0 with PID 47930 received signal 6
>
> Writing backtrace to nsnstohmns/backtrace.0.txt
>
> Rank 1 with PID 47931 received signal 6
>
> Writing backtrace to nsnstohmns/backtrace.1.txt
>
> cactus_sim: /global/u1/c/chlin/Cactus/arrangements/Carpet/Carpet/src/helpers.cc:275: int Carpet::Abort(const _cGH *, int): Assertion `0' failed.
>
> cactus_sim: /global/u1/c/chlin/Cactus/arrangements/Carpet/Carpet/src/helpers.cc:275: int Carpet::Abort(const _cGH *, int): Assertion `0' failed.
>
> Rank 2 with PID 13239 received signal 6
>
> Rank 3 with PID 13240 received signal 6
>
> Writing backtrace to nsnstohmns/backtrace.2.txt
>
> Writing backtrace to nsnstohmns/backtrace.3.txt
>
> srun: error: nid05930: tasks 2-3: Aborted
>
> srun: Terminating job step 9192104.0
>
> srun: error: nid05920: tasks 0-1: Aborted
>
>
> The code is aborted without completing the whole process until the merger. How can I perform the complete process?
>
> 2. When I use visit to do the visualization ,I found the file rho.xy.h5 only has half of the whole simulation as :
>
> [cid:e66a53d5-d406-4599-997f-00e90a864dc0]
>
>
>
>
>
> The tutorial(scripts.tar.gz<http://einsteintoolkit.org/gallery/bns/scripts.tar.gz>) tells that I should use symmetric reflection to perform the whole visualization. However I cannot use this method since I simulated binary star with different mass . So there is no symmetry. How can I get the whole picture instead of only half ?
>
>
> Thank you.
>
> Best regards,
>
> Chia-Hui Lin



--
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 --------------
An HTML attachment was scrubbed...
URL: http://lists.einsteintoolkit.org/pipermail/users/attachments/20180514/5fb95be3/attachment-0001.html 
-------------- next part --------------
A non-text attachment was scrubbed...
Name: nsns_09_063.err
Type: application/octet-stream
Size: 47955 bytes
Desc: nsns_09_063.err
Url : http://lists.einsteintoolkit.org/pipermail/users/attachments/20180514/5fb95be3/attachment-0001.obj