<html>
<head>
<meta http-equiv="Content-Type" content="text/html; charset=big5">
</head>
<body>
<style type="text/css" style="display:none;"><!-- P {margin-top:0;margin-bottom:0;} --></style>
<div id="divtagdefaultwrapper" style="font-size:12pt;color:#000000;font-family:Calibri,Helvetica,sans-serif;" dir="ltr">
<p>Dear Zach,</p>
<p>I see. Thank you !</p>
<p>Best regards,</p>
<p>Chia-Hui</p>
</div>
<hr style="display:inline-block;width:98%" tabindex="-1">
<div id="divRplyFwdMsg" dir="ltr"><font face="Calibri, sans-serif" style="font-size:11pt" color="#000000"><b>寄件者:</b> Zach Etienne &lt;zachetie@gmail.com&gt;<br>
<b>寄件日期:</b> 2019年5月31日 上午 12:47:00<br>
<b>收件者:</b> 林家暉<br>
<b>副本:</b> Haas, Roland; Einstein Toolkit Users<br>
<b>主旨:</b> Re: [Users] interpolation error</font>
<div>&nbsp;</div>
</div>
<div>
<div dir="ltr">
<div dir="ltr">Dear Chia-Hui,
<div><br>
</div>
<div><span style="color:rgb(0,0,0);font-family:Calibri,Helvetica,sans-serif,Helvetica,EmojiFont,&quot;Apple Color Emoji&quot;,&quot;Segoe UI Emoji&quot;,NotoColorEmoji,&quot;Segoe UI Symbol&quot;,&quot;Android Emoji&quot;,EmojiSymbols;font-size:16px">out_of_bounds_xyz is denoted as&nbsp;</span>STEERABLE=ALWAYS.
 Therefore you can change the value at any time in your parameter file, even after a checkpoint.<br>
</div>
<div><br clear="all">
<div>
<div dir="ltr" class="gmail_signature">
<div dir="ltr">
<div>
<div dir="ltr">
<div>
<div dir="ltr">
<div>
<div dir="ltr">
<div dir="ltr">
<div dir="ltr">
<div style="font-size:12.8px">-Zach</div>
<div style="font-size:12.8px"><br>
</div>
<span style="font-size:12.8px">*&nbsp;&nbsp;&nbsp;&nbsp; *&nbsp;&nbsp;&nbsp;&nbsp; *</span><br style="font-size:12.8px">
<span style="font-size:12.8px">Prof. Zachariah Etienne</span><br style="font-size:12.8px">
<div style="font-size:12.8px">West Virginia University</div>
<div><font color="#0000ee"><u><a href="https://math.wvu.edu/~zetienne/" target="_blank">https://math.wvu.edu/~zetienne/</a></u></font></div>
<div><a href="https://blackholesathome.net" style="font-size:12.8px" target="_blank">https://blackholesathome.net</a><br>
</div>
</div>
</div>
</div>
</div>
</div>
</div>
</div>
</div>
</div>
</div>
</div>
<br>
</div>
</div>
</div>
<br>
<div class="gmail_quote">
<div dir="ltr" class="gmail_attr">On Thu, May 30, 2019 at 11:43 AM 林家暉 &lt;<a href="mailto:r06222015@ntu.edu.tw">r06222015@ntu.edu.tw</a>&gt; wrote:<br>
</div>
<blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex">
<div>
<div id="gmail-m_-8578809389745933503divtagdefaultwrapper" style="font-size:12pt;color:rgb(0,0,0);font-family:Calibri,Helvetica,sans-serif" dir="ltr">
<p>Dear Rolnd and Zach,&nbsp;</p>
<p>Thanks for your explanation. My domain size is <span>(-540,-540,-540)x(540,540,540), so the&nbsp;<span>interpolation point is outside of the domain as you said.</span></span></p>
<p>For using&nbsp;<span>out_of_bounds_xyz, do I need to rerun the simulation from the beginning or I can run it from the last checkpoint after adding&nbsp;<span style="font-family:Calibri,Helvetica,sans-serif,Helvetica,EmojiFont,&quot;Apple Color Emoji&quot;,&quot;Segoe UI Emoji&quot;,NotoColorEmoji,&quot;Segoe UI Symbol&quot;,&quot;Android Emoji&quot;,EmojiSymbols;font-size:16px">out_of_bounds_xyz
 in the parameter file</span>&nbsp;?</span></p>
<p>Best regards,</p>
<p>Chia-Hui</p>
</div>
<hr style="display:inline-block;width:98%">
<div id="gmail-m_-8578809389745933503divRplyFwdMsg" dir="ltr"><font face="Calibri, sans-serif" style="font-size:11pt" color="#000000"><b>寄件者:</b> Zach Etienne &lt;<a href="mailto:zachetie@gmail.com" target="_blank">zachetie@gmail.com</a>&gt;<br>
<b>寄件日期:</b> 2019年5月30日 下午 11:03:03<br>
<b>收件者:</b> Haas, Roland<br>
<b>副本:</b> 林家暉; Einstein Toolkit Users<br>
<b>主旨:</b> Re: [Users] interpolation error</font>
<div>&nbsp;</div>
</div>
<div>
<div dir="ltr">Hi Chia-Hui,
<div><br>
</div>
<div>particle_tracerET already has a parameter for this: out_of_bounds_xyz</div>
<div><br>
</div>
<div>I follow the general rule of setting it to about 80% of my (coarsest) grid size.</div>
<div><br clear="all">
<div>
<div dir="ltr" class="gmail-m_-8578809389745933503m_5683372223795461865m_1823707245693416716gmail_signature">
<div dir="ltr">
<div>
<div dir="ltr">
<div>
<div dir="ltr">
<div>
<div dir="ltr">
<div dir="ltr">
<div dir="ltr">
<div style="font-size:12.8px">-Zach</div>
<div style="font-size:12.8px"><br>
</div>
<span style="font-size:12.8px">*&nbsp;&nbsp;&nbsp;&nbsp; *&nbsp;&nbsp;&nbsp;&nbsp; *</span><br style="font-size:12.8px">
<span style="font-size:12.8px">Prof. Zachariah Etienne</span><br style="font-size:12.8px">
<div style="font-size:12.8px">West Virginia University</div>
<div><font color="#0000ee"><u><a href="https://math.wvu.edu/~zetienne/" target="_blank">https://math.wvu.edu/~zetienne/</a></u></font></div>
<div><a href="https://blackholesathome.net" style="font-size:12.8px" target="_blank">https://blackholesathome.net</a><br>
</div>
</div>
</div>
</div>
</div>
</div>
</div>
</div>
</div>
</div>
</div>
</div>
<br>
</div>
</div>
<br>
<div class="gmail_quote">
<div dir="ltr" class="gmail_attr">On Thu, May 30, 2019 at 10:58 AM Haas, Roland &lt;<a href="mailto:rhaas@illinois.edu" target="_blank">rhaas@illinois.edu</a>&gt; wrote:<br>
</div>
<blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex">
Hello Chia-Hui,<br>
<br>
the error message indicates that the code tried to interpolation a<br>
quantity to location (x,y,z)=(-157.725,-576.349,-99.2027) but that the<br>
grid only has a size of (xmin,ymin,zmin)x(xmax,ymax,zmax) =<br>
(-576,-576,-576)x(576,576,576) (where I made a guess for the ymax and<br>
zmax values given that the output only shows the shape of the grid on<br>
one MPI rank).<br>
<br>
So the interpolation point is outside of the domain.<br>
<br>
Depending on what caused the interpolation this could mean eg that a<br>
tracer particle (given the function name<br>
&quot;DoOneRK4StepForParticleTracerET&quot;) tried to leave the domain. This is<br>
something that needsto be fixed in the particle_tracerET thorn,<br>
probably by making it deactivate particles that get too close to the<br>
edge of the domain.<br>
<br>
Yours,<br>
Roland<br>
<br>
&gt; Hello,<br>
&gt; <br>
&gt; I am using IllinoisGRMHD with the test with nsns_test.par in <a href="https://bitbucket.org/zach_etienne/wvuthorns_diagnostics/src/master/NSNS_parameter_files/" rel="noreferrer" target="_blank">
https://bitbucket.org/zach_etienne/wvuthorns_diagnostics/src/master/NSNS_parameter_files/</a>&nbsp; .<br>
&gt; <br>
&gt; The simulation runs well until t~5300 and an error occurs. The error message is like this:<br>
&gt; <br>
&gt; <br>
&gt; WARNING level 1 from host nid02210 process 0<br>
&gt;&nbsp; &nbsp;while executing schedule bin particle_tracerET, routine particle_tracerET::DoOneRK4StepForParticleTracerET<br>
&gt;&nbsp; &nbsp;in thorn AEILocalInterp, file /global/u1/c/chlin/et/Cactus/arrangements/Numerical/AEILocalInterp/src/Lagrange-tensor-product/../template.c:1109:<br>
&gt;&nbsp; &nbsp;-&gt;&nbsp; <br>
&gt;&nbsp; &nbsp; CCTK_InterpLocalUniform():<br>
&gt;&nbsp; &nbsp; &nbsp; &nbsp; &nbsp;interpolation point is either outside the grid,<br>
&gt;&nbsp; &nbsp; &nbsp; &nbsp; &nbsp;or inside but too close to the grid boundary!<br>
&gt;&nbsp; &nbsp; &nbsp; &nbsp; &nbsp;(this may be caused by a global interpolation with<br>
&gt;&nbsp; &nbsp; &nbsp; &nbsp; &nbsp; driver::ghost_size too small)<br>
&gt;&nbsp; &nbsp; &nbsp; &nbsp; &nbsp;0-origin interpolation point number pt=0 of N_interp_points=1<br>
&gt;&nbsp; &nbsp; &nbsp; &nbsp; &nbsp;interpolation point (x,y,z)=(-157.725,-576.349,-99.2027)<br>
&gt;&nbsp; &nbsp; &nbsp; &nbsp; &nbsp;grid x_min(delta_x)x_max = -576(18)576<br>
&gt;&nbsp; &nbsp; &nbsp; &nbsp; &nbsp;grid y_min(delta_y)y_max = -576(18)54<br>
&gt;&nbsp; &nbsp; &nbsp; &nbsp; &nbsp;grid z_min(delta_z)z_max = -576(18)54<br>
&gt; <br>
&gt; WARNING level 1 from host nid02210 process 0<br>
&gt;&nbsp; &nbsp;while executing schedule bin particle_tracerET, routine particle_tracerET::DoOneRK4StepForParticleTracerET<br>
&gt;&nbsp; &nbsp;in thorn particle_tracerET, file /global/u1/c/chlin/et/Cactus/arrangements/WVUThorns_Diagnostics/particle_tracerET/src/Interpolate_velocities_at_particle_positions.C:131:<br>
&gt;&nbsp; &nbsp;-&gt; interpolation screwed up&nbsp; <br>
&gt; WARNING level 1 from host nid02210 process 1<br>
&gt;&nbsp; &nbsp;while executing schedule bin particle_tracerET, routine particle_tracerET::DoOneRK4StepForParticleTracerET<br>
&gt;&nbsp; &nbsp;in thorn particle_tracerET, file /global/u1/c/chlin/et/Cactus/arrangements/WVUThorns_Diagnostics/particle_tracerET/src/Interpolate_velocities_at_particle_positions.C:131:<br>
&gt;&nbsp; &nbsp;-&gt; interpolation screwed up&nbsp; <br>
&gt; srun: error: nid02211: tasks 2-3: Exited with exit code 1<br>
&gt; srun: Terminating job step 21753671.0<br>
&gt; srun: error: nid02210: tasks 0-1: Exited with exit code 1<br>
&gt; <br>
&gt; <br>
&gt; I think maybe t can be avoided by decreasing the timestep so that the tracer particle would not move too fast to make the&nbsp; interpolation screwed up but I am not sure. So what I want to ask is :<br>
&gt; <br>
&gt;&nbsp; &nbsp;1.&nbsp; It this error can be avoided by decreasing the timestep ?<br>
&gt;&nbsp; &nbsp;2.&nbsp; If so, how can I&nbsp; decrease the timestep in parameter file (or in other places)?<br>
&gt; <br>
&gt; <br>
&gt; Thank you.<br>
&gt; Best regards,<br>
&gt; Chia-Hui<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" rel="noreferrer" target="_blank">
http://pgp.mit.edu</a> .<br>
_______________________________________________<br>
Users mailing list<br>
<a href="mailto:Users@einsteintoolkit.org" target="_blank">Users@einsteintoolkit.org</a><br>
<a href="http://lists.einsteintoolkit.org/mailman/listinfo/users" rel="noreferrer" target="_blank">http://lists.einsteintoolkit.org/mailman/listinfo/users</a><br>
</blockquote>
</div>
</div>
</div>
</blockquote>
</div>
</div>
</body>
</html>