<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>I do not quite understand what causes ,in your example,</p>
<p><span style="color:rgb(33,33,33); font-family:"Microsoft JHengHei UI","Microsoft JHengHei",微軟正黑體,MingLiU,細明體,sans-serif,serif,EmojiFont; font-size:13.3333px">CarpetRegrid2::position_x_1 = 42.</span><br style="color:rgb(33,33,33); font-family:"Microsoft JHengHei UI","Microsoft JHengHei",微軟正黑體,MingLiU,細明體,sans-serif,serif,EmojiFont; font-size:13.3333px">
<span style="color:rgb(33,33,33); font-family:"Microsoft JHengHei UI","Microsoft JHengHei",微軟正黑體,MingLiU,細明體,sans-serif,serif,EmojiFont; font-size:13.3333px">CarpetRegrid2::position_x_2 = -17</span><br>
</p>
<p><span style="color:rgb(33,33,33); font-family:"Microsoft JHengHei UI","Microsoft JHengHei",微軟正黑體,MingLiU,細明體,sans-serif,serif,EmojiFont; font-size:13.3333px"><span style="font-family:Calibri,Helvetica,sans-serif,Helvetica,EmojiFont,"Apple Color Emoji","Segoe UI Emoji",NotoColorEmoji,"Segoe UI Symbol","Android Emoji",EmojiSymbols; font-size:16px">to
become fail. In my case, if the <span style="font-family:Calibri,Helvetica,sans-serif,serif,EmojiFont; font-size:16px">initial location of the binary star is (-34.59 km , 45.41 km) ,I should set </span></span><br>
</span></p>
<p><span style="color:rgb(33,33,33); font-family:"Microsoft JHengHei UI","Microsoft JHengHei",微軟正黑體,MingLiU,細明體,sans-serif,serif,EmojiFont; font-size:13.3333px"><span style="font-family:Calibri,Helvetica,sans-serif,Helvetica,EmojiFont,"Apple Color Emoji","Segoe UI Emoji",NotoColorEmoji,"Segoe UI Symbol","Android Emoji",EmojiSymbols; font-size:16px"><span style="font-family:Calibri,Helvetica,sans-serif,serif,EmojiFont; font-size:16px"><span style="color:rgb(33,33,33); font-family:"Microsoft JHengHei UI","Microsoft JHengHei",微軟正黑體,MingLiU,細明體,sans-serif,serif,EmojiFont; font-size:13.3333px">CarpetRegrid2::position_x_1
= -34.59/1.48=-23.37</span><br style="color:rgb(33,33,33); font-family:"Microsoft JHengHei UI","Microsoft JHengHei",微軟正黑體,MingLiU,細明體,sans-serif,serif,EmojiFont; font-size:13.3333px">
<span style="color:rgb(33,33,33); font-family:"Microsoft JHengHei UI","Microsoft JHengHei",微軟正黑體,MingLiU,細明體,sans-serif,serif,EmojiFont; font-size:13.3333px">CarpetRegrid2::position_x_2 = 45.41/1.48=30.68</span><br>
</span></span></span></p>
<p><span style="color:rgb(33,33,33); font-family:Calibri,Helvetica,sans-serif,serif,EmojiFont; font-size:16px">right? Would this soon be set to (-23.37,23.37)? Because the mass difference is not so much, I am wondering whether there exists some maximum value
of mass ratio which can tolerate small mass difference.</span><br>
</p>
<p><span style="color:rgb(33,33,33); font-family:Calibri,Helvetica,sans-serif,serif,EmojiFont; font-size:16px">By the way ,the link <a href="https://trac.einsteintoolkit.org/ticket/1974" class="x_OWAAutoLink" id="LPlnk155431">https://trac.einsteintoolkit.org/ticket/1974</a> seems not
work. I would like to check it up.</span></p>
<p><span style="color:rgb(33,33,33); font-family:Calibri,Helvetica,sans-serif,serif,EmojiFont; font-size:16px">Thanks,</span></p>
<p><font color="#212121" face="Calibri, Helvetica, sans-serif, serif, EmojiFont">Best regards,</font></p>
<p><font color="#212121" face="Calibri, Helvetica, sans-serif, serif, EmojiFont">Chia-Hui</font></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:46:15<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>
ok that would make the distance from the origin ~38 units (Msun)<br>
which is certainly not what the parameter file is designed to handle.<br>
<br>
The code only flips the x and y coordinates, not z as well so (x,y,z)<br>
and (-x,-y,z) thought this will not make a difference for you as the<br>
center of the stars is at z=0 anyway.<br>
<br>
You have to edit the .par file and there are two lines in there that<br>
read:<br>
<br>
CarpetRegrid2::position_x_1 = 15.1875<br>
<br>
and<br>
<br>
CarpetRegrid2::position_x_2 = -15.1875<br>
<br>
You have to change 15.1875 to the initial x position of one of the stars (measured in solar masses ie about 1.48 km) and the -15.1875 to that of the other one.<br>
<br>
You then have to restart your simulation from the beginning. Note that I expect this to still fail since even if you initially give the code:<br>
<br>
CarpetRegrid2::position_x_1 = 42.<br>
CarpetRegrid2::position_x_2 = -17<br>
<br>
it will very soon set the centers to be<br>
<br>
x_1 = 17, x_2 = -17 <br>
<br>
since the star at x_2 = -17 is the heavier one (and thus has larger<br>
central density).<br>
<br>
There is a thorn that Zach Etienne is proposing for inclusion in the<br>
Einstein Toolkit: <a href="https://trac.einsteintoolkit.org/ticket/1974">https://trac.einsteintoolkit.org/ticket/1974</a> that he<br>
says can track two unequal mass stars. The thorn in question would seem<br>
to be VolumeIntegrals_GRMHD .<br>
<br>
Yours,<br>
Roland<br>
<br>
> Dear Roland,<br>
> <br>
> Yes, I did use the parameter file from gallery page and use unequal mass initial data from Lorene.<br>
> <br>
> I just check the initial location of the binary star is (-34.59 km , 45.41 km).However you said 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) .How should I set CarpetRegrid2::position_x_1 and CarpetRegrid2::position_x_2
?<br>
> <br>
> Thank you.<br>
> <br>
> Best regards,<br>
> <br>
> Chia-Hui<br>
> <br>
> ________________________________<br>
> 寄件者: Roland Haas <rhaas@illinois.edu><br>
> 寄件日期: 2018年5月15日 上午 02:13:06<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,<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>
<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>