<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">
<p>Dear Roland,</p>
<p>Yes, I did use the parameter file from gallery page and use unequal mass initial data from Lorene.</p>
<p>I just check the initial location of the binary star is (-34.59 km , 45.41 km).However you said it<span> makes the hard-coded assumption that if one of the stars is at (x,y,z) the other one is at (-x,-y,z) .How should I set <span>CarpetRegrid2::position_x_1
and <span>CarpetRegrid2::position_x_2 ?</span></span></span></p>
<p><span><span><span>Thank you.</span></span></span></p>
<p><span><span><span>Best regards,</span></span></span></p>
<p><span><span><span>Chia-Hui</span></span></span></p>
</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月15日 上午 02:13:06<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,<br>
<br>
thank you for the attached file.<br>
<br>
Yes, the error is related to the grid structure:<br>
<br>
WARNING level 0 from host nid05920 process 0<br>
while executing schedule bin (none), routine (no thorn)::(no routine)<br>
in thorn CarpetLib, file /global/u1/c/chlin/Cactus/arrangements/Carpet/CarpetLib/src/gh.cc:178:<br>
-> The refinement hierarchy is not properly nested.<br>
<br>
This is for the run where you used the parameter file from the gallery<br>
page but used unequal mass initial data from LORENE?<br>
<br>
I am asking because, given that the parameter file specifies<br>
Rotating180 symmetry, the evolution will be as if equal mass NS were<br>
given in the initial data file (it will duplicate the right hand NS and<br>
ignore the left hand one). So the error that you saw should not be due<br>
to the unequal mass per-se. It could be if the location of the stars<br>
changed initially (eg for unequal mass the more massive one will be<br>
closer to the origin than the other one). This still should not give<br>
you this particular error (it would be wrong in other ways though)<br>
since the parameter file uses CarpetRegrid2::ensure_proper_nesting =<br>
"yes" which should avoid this particular error.<br>
<br>
Looking at the plot from your initial email it may be that there was an<br>
error in the initial setups as it looks to me as if the neutron stars<br>
have actually drifted apart. At least the gallery examples has the NS<br>
initially at ~+/-15 units from the origin but in you plot it looks like<br>
much more.<br>
<br>
Wen you added the new ID you must make sure to adjust the initial<br>
location in the parfile in the lines reading:<br>
<br>
CarpetRegrid2::position_x_1 = 15.1875<br>
CarpetRegrid2::position_x_2 = -15.1875<br>
<br>
The binary black hole example does not have the restriction of<br>
assuming 180 degree symmetry, no.<br>
<br>
Yours,<br>
Roland<br>
<br>
> Dear Roland,<br>
> 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.<br>
> For not able to track multiple stars, does it occur in binary blackhole case(<a href="http://einsteintoolkit.org/gallery/bbh/index.html">http://einsteintoolkit.org/gallery/bbh/index.html</a>) ?<br>
> Best regards,<br>
> Chia-Hui<br>
> ________________________________<br>
> 寄件者: Roland Haas <rhaas@illinois.edu><br>
> 寄件日期: 2018年5月13日 下午 09:28:51<br>
> 收件者: 林家暉<br>
> 副本: Einstein Toolkit Users<br>
> 主旨: Re: [Users] problems of reproducing the simulation of binary neutron star with different mass<br>
> <br>
> 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>
<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>