<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 Rolnd and Zach, </p>
<p>Thanks for your explanation. My domain size is <span>(-540,-540,-540)x(540,540,540), so the <span>interpolation point is outside of the domain as you said.</span></span></p>
<p>For using <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 <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;">out_of_bounds_xyz
in the parameter file</span> ?</span></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 <zachetie@gmail.com><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> </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="m_5683372223795461865m_1823707245693416716gmail_signature" data-smartmail="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">* * *</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 <<a href="mailto:rhaas@illinois.edu" target="_blank">rhaas@illinois.edu</a>> 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>
"DoOneRK4StepForParticleTracerET") 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>
> Hello,<br>
> <br>
> 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> .<br>
> <br>
> The simulation runs well until t~5300 and an error occurs. The error message is like this:<br>
> <br>
> <br>
> WARNING level 1 from host nid02210 process 0<br>
> while executing schedule bin particle_tracerET, routine particle_tracerET::DoOneRK4StepForParticleTracerET<br>
> in thorn AEILocalInterp, file /global/u1/c/chlin/et/Cactus/arrangements/Numerical/AEILocalInterp/src/Lagrange-tensor-product/../template.c:1109:<br>
> -> <br>
> CCTK_InterpLocalUniform():<br>
> interpolation point is either outside the grid,<br>
> or inside but too close to the grid boundary!<br>
> (this may be caused by a global interpolation with<br>
> driver::ghost_size too small)<br>
> 0-origin interpolation point number pt=0 of N_interp_points=1<br>
> interpolation point (x,y,z)=(-157.725,-576.349,-99.2027)<br>
> grid x_min(delta_x)x_max = -576(18)576<br>
> grid y_min(delta_y)y_max = -576(18)54<br>
> grid z_min(delta_z)z_max = -576(18)54<br>
> <br>
> WARNING level 1 from host nid02210 process 0<br>
> while executing schedule bin particle_tracerET, routine particle_tracerET::DoOneRK4StepForParticleTracerET<br>
> 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>
> -> interpolation screwed up <br>
> WARNING level 1 from host nid02210 process 1<br>
> while executing schedule bin particle_tracerET, routine particle_tracerET::DoOneRK4StepForParticleTracerET<br>
> 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>
> -> interpolation screwed up <br>
> srun: error: nid02211: tasks 2-3: Exited with exit code 1<br>
> srun: Terminating job step 21753671.0<br>
> srun: error: nid02210: tasks 0-1: Exited with exit code 1<br>
> <br>
> <br>
> I think maybe t can be avoided by decreasing the timestep so that the tracer particle would not move too fast to make the interpolation screwed up but I am not sure. So what I want to ask is :<br>
> <br>
> 1. It this error can be avoided by decreasing the timestep ?<br>
> 2. If so, how can I decrease the timestep in parameter file (or in other places)?<br>
> <br>
> <br>
> Thank you.<br>
> Best regards,<br>
> 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>
</body>
</html>