<html>
<head>
<meta http-equiv="Content-Type" content="text/html; charset=big5">
<meta name="Generator" content="Microsoft Exchange Server">
<!-- converted from text --><style><!-- .EmailQuote { margin-left: 1pt; padding-left: 4pt; border-left: #800000 2px solid; } --></style>
</head>
<body>
<meta content="text/html; charset=UTF-8">
<style type="text/css" style="">
<!--
p
{margin-top:0;
margin-bottom:0}
-->
</style>
<div dir="ltr">
<div id="x_divtagdefaultwrapper" dir="ltr" style="font-size:12pt; color:#000000; font-family:Calibri,Helvetica,sans-serif">
Dear Roland,
<div>Thank you for your reply.<br>
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.</div>
<div>For not able to track multiple stars, does it occur in binary blackhole case(<a href="http://einsteintoolkit.org/gallery/bbh/index.html" class="x_OWAAutoLink" id="LPlnk912279">http://einsteintoolkit.org/gallery/bbh/index.html</a>) ?</div>
<div>Best regards,</div>
<div>Chia-Hui</div>
</div>
<hr tabindex="-1" style="display:inline-block; width:98%">
<div id="x_divRplyFwdMsg" dir="ltr"><font face="Calibri, sans-serif" color="#000000" style="font-size:11pt"><b>寄件者:</b> Roland Haas <rhaas@illinois.edu><br>
<b>寄件日期:</b> 2018年5月13日 下午 09:28:51<br>
<b>收件者:</b> 林家暉<br>
<b>副本:</b> Einstein Toolkit Users<br>
<b>主旨:</b> Re: [Users] problems of reproducing the simulation of binary neutron star with different mass</font>
<div> </div>
</div>
</div>
<font size="2"><span style="font-size:10pt;">
<div class="PlainText">Hello Chia-Hui Lin,<br>
<br>
changing the central enthalpy seems like the right thing to do to<br>
change masses in LORENE since it prescribes the central enthalpy of<br>
each star.<br>
<br>
The backtrace that you showed and the error message is unfortunately<br>
not detailed enough to make a guess as to what could have gone wrong.<br>
This output is produced by a call to CCTK_Error() however the<br>
CCTK_Error function takes an error string as an argument which is not<br>
present in the lines that you have copied. It is likely in the *.err<br>
file that the simulation produced and you should try and include that<br>
file as an attachment in your email.<br>
<br>
If you have non-equal mass stars then you have to remove the 180 degree<br>
rotation symmetry from the parameter files. Ie you have to change the<br>
lines (these are scattered throughout throughout the par file)<br>
<br>
ActiveThorns = "RotatingSymmetry180 ReflectionSymmetry"<br>
CoordBase::xmin = 0.00<br>
CoordBase::boundary_shiftout_x_lower = 1<br>
CarpetRegrid2::symmetry_rotating180 = "yes"<br>
Hydro_Analysis::Hydro_Analysis_rho_max_loc_only_positive_x = "true"<br>
Hydro_Analysis::Hydro_Analysis_average_multiple_maxima_locations = "yes"<br>
<br>
to<br>
<br>
ActiveThorns = "ReflectionSymmetry"<br>
CoordBase::xmin = 396.00<br>
CoordBase::boundary_shiftout_x_lower = 0<br>
CarpetRegrid2::symmetry_rotating180 = "no"<br>
Hydro_Analysis::Hydro_Analysis_rho_max_loc_only_positive_x = "false"<br>
Hydro_Analysis::Hydro_Analysis_average_multiple_maxima_locations = "no"<br>
<br>
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.<br>
<br>
Yours,<br>
Roland<br>
<br>
> Dear whom it may concern,<br>
> <br>
> I am trying to reproduce the simulation of binary neutron star (<a href="http://einsteintoolkit.org/gallery/bns/index.html">http://einsteintoolkit.org/gallery/bns/index.html</a>) 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:<br>
> <br>
> 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:<br>
> <br>
> <br>
> <br>
> cactus_sim: /global/u1/c/chlin/Cactus/arrangements/Carpet/Carpet/src/helpers.cc:275: int Carpet::Abort(const _cGH *, int): Assertion `0' failed.<br>
> <br>
> cactus_sim: /global/u1/c/chlin/Cactus/arrangements/Carpet/Carpet/src/helpers.cc:275: int Carpet::Abort(const _cGH *, int): Assertion `0' failed.<br>
> <br>
> Rank 0 with PID 47930 received signal 6<br>
> <br>
> Writing backtrace to nsnstohmns/backtrace.0.txt<br>
> <br>
> Rank 1 with PID 47931 received signal 6<br>
> <br>
> Writing backtrace to nsnstohmns/backtrace.1.txt<br>
> <br>
> cactus_sim: /global/u1/c/chlin/Cactus/arrangements/Carpet/Carpet/src/helpers.cc:275: int Carpet::Abort(const _cGH *, int): Assertion `0' failed.<br>
> <br>
> cactus_sim: /global/u1/c/chlin/Cactus/arrangements/Carpet/Carpet/src/helpers.cc:275: int Carpet::Abort(const _cGH *, int): Assertion `0' failed.<br>
> <br>
> Rank 2 with PID 13239 received signal 6<br>
> <br>
> Rank 3 with PID 13240 received signal 6<br>
> <br>
> Writing backtrace to nsnstohmns/backtrace.2.txt<br>
> <br>
> Writing backtrace to nsnstohmns/backtrace.3.txt<br>
> <br>
> srun: error: nid05930: tasks 2-3: Aborted<br>
> <br>
> srun: Terminating job step 9192104.0<br>
> <br>
> srun: error: nid05920: tasks 0-1: Aborted<br>
> <br>
> <br>
> The code is aborted without completing the whole process until the merger. How can I perform the complete process?<br>
> <br>
> 2. When I use visit to do the visualization ,I found the file rho.xy.h5 only has half of the whole simulation as :<br>
> <br>
> [cid:e66a53d5-d406-4599-997f-00e90a864dc0]<br>
> <br>
> <br>
> <br>
> <br>
> <br>
> 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 ?<br>
> <br>
> <br>
> Thank you.<br>
> <br>
> Best regards,<br>
> <br>
> Chia-Hui Lin<br>
<br>
<br>
<br>
-- <br>
My email is as private as my paper mail. I therefore support encrypting<br>
and signing email messages. Get my PGP key from <a href="http://pgp.mit.edu">http://pgp.mit.edu</a> .<br>
</div>
</span></font>
</body>
</html>