<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 found the tracking problem is very obvious that the canter of mass moves toward right. </p>
<p>Therefore,I have tried the thorn you mentioned (<a href="https://bitbucket.org/zach_etienne/wvuthorns_diagnostics" class="x_OWAAutoLink" id="LPlnk114222">https://bitbucket.org/zach_etienne/wvuthorns_diagnostics</a>).</p>
<p>However I met some problem. I added following lines in <span style="font-size:12pt">nsnstohmns.th</span></p>
<p><span style="font-size:12pt"><br>
</span></p>
<div># wvuthorns_diagnostics</div>
<div>!TARGET = $ARR</div>
<div>!TYPE = git</div>
<div>!URL = https://bitbucket.org/zach_etienne/wvuthorns_diagnostics.git</div>
<div>!REPO_PATH= $2</div>
<div>!CHECKOUT =</div>
<div>wvuthorns_diagnostics/VolumeIntegrals_GRMHD</div>
<br>
<div>and following lines in <span style="font-size:12pt">nsnstohmns.par</span></div>
<div><span style="font-size:12pt"><br>
</span></div>
<div>
<div>ActiveThorns = "VolumeIntegrals_GRMHD"</div>
<div>VolumeIntegrals_GRMHD::volintegral_sphere__center_x_initial[101] ="0.0"</div>
<div>VolumeIntegrals_GRMHD::volintegral_sphere__center_y_initial[101] ="0.0"</div>
<div>VolumeIntegrals_GRMHD::volintegral_sphere__center_z_initial[101] ="0.0"</div>
<br>
</div>
<div>then I ran the command in order to obtain the thorn <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">VolumeIntegrals_GRMHD</span> :</div>
<div></div>
<span>./GetComponents </span><span>thornlists/nsnstohmns.th</span>
<div><br>
<div>but the following error occurs</div>
<div><br>
<div>Error: Thorn VolumeIntegrals_GRMHD not found</div>
<div>Activation failed - 1 errors in activation sequence</div>
<div>WARNING level 0 from host nid02749 process 0</div>
<div> while executing schedule bin (none), routine (no thorn)::(no routine)</div>
<div> in thorn Cactus, file /global/u1/c/chlin/Cactus/src/main/SetParams.c:93:</div>
<div> -> CCTKi_SetParameter: Error at line 585 in parameter file /scratch1/scratchdirs/chlin/simulations/nsns_09_063_thorn_test5/output-0000/nsnstohmns.par while activating thorns</div>
</div>
</div>
<div><br>
</div>
<div>I wonder that I missed some steps since it is my first time to add a new thorn. Is there any possible solution?</div>
<div>Thank you.</div>
<div>Best regards,</div>
<div>Chia-Hui</div>
<div><br>
</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月15日 下午 08:16:12<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>
> I do not quite understand what causes ,in your example,<br>
> <br>
> CarpetRegrid2::position_x_1 = 42.<br>
> CarpetRegrid2::position_x_2 = -17<br>
> <br>
> to become fail.<br>
This has to do with the way that the thorn GRHydro_Analysis that looks<br>
for the location of the star works. It first determines the location on<br>
the grid with the largest density, then it computes a center of mass<br>
integral in a region around that point. This computed center of mass if<br>
used as the location of the star and to move boxes. There is code in<br>
thorn NSTracker that uses this single found location and moves two<br>
refined regions to (x,y,z) and (-x,-y,z) where (x,y,z) is the found<br>
location of the star.<br>
<br>
There are two reasons why this does not work well with unequal masses:<br>
<br>
1. it only looks for a single star and bases that star's location on<br>
the highest density<br>
2. NSTracker makes the hard-coded assumption that the two regions that<br>
cover the star are mirror immages of each other<br>
<br>
In your case since the star at -17 (always assuming you actually have<br>
set up a grid that has no symmetries) is the heavier one (since it is<br>
closer to the center of rotation) and thus has the higher central<br>
density, therefore GRHydro_Analysis will use it as the star to track<br>
and thus will find the star location to be (-17,0,0) when it looks for<br>
the star.<br>
<br>
> In my case, if the initial location of the binary star is (-34.59<br>
> km , 45.41 km) , I should set<br>
> <br>
> CarpetRegrid2::position_x_1 = -34.59/1.48=-23.37<br>
> CarpetRegrid2::position_x_2 = 45.41/1.48=30.68<br>
> <br>
> right? <br>
Yes that is correct.<br>
<br>
> Would this soon be set to (-23.37,23.37)? Because the mass difference<br>
> is not so much, I am wondering whether there exists some maximum value<br>
> of mass ratio which can tolerate small mass difference.<br>
> <br>
> By the way ,the link <a href="https://trac.einsteintoolkit.org/ticket/1974">https://trac.einsteintoolkit.org/ticket/1974</a><br>
> seems not work. I would like to check it up.<br>
Thank you for pointing this out. It seems that there is an endless<br>
redirection loop in the trac website already for<br>
<a href="https://trac.einsteintoolkit.org">https://trac.einsteintoolkit.org</a><br>
<br>
Yours,<br>
Roland<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>