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 <rh...@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 .

Attachment: nsns_09_063.err
Description: nsns_09_063.err

_______________________________________________
Users mailing list
Users@einsteintoolkit.org
http://lists.einsteintoolkit.org/mailman/listinfo/users

Reply via email to