<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>