Hi Chia-Hui,
particle_tracerET already has a parameter for this: out_of_bounds_xyz
I follow the general rule of setting it to about 80% of my (coarsest) grid size.
-Zach
* * * Prof. Zachariah Etienne West Virginia University *https://math.wvu.edu/~zetienne/ https://math.wvu.edu/~zetienne/* https://blackholesathome.net
On Thu, May 30, 2019 at 10:58 AM Haas, Roland rhaas@illinois.edu wrote:
Hello Chia-Hui,
the error message indicates that the code tried to interpolation a quantity to location (x,y,z)=(-157.725,-576.349,-99.2027) but that the grid only has a size of (xmin,ymin,zmin)x(xmax,ymax,zmax) = (-576,-576,-576)x(576,576,576) (where I made a guess for the ymax and zmax values given that the output only shows the shape of the grid on one MPI rank).
So the interpolation point is outside of the domain.
Depending on what caused the interpolation this could mean eg that a tracer particle (given the function name "DoOneRK4StepForParticleTracerET") tried to leave the domain. This is something that needsto be fixed in the particle_tracerET thorn, probably by making it deactivate particles that get too close to the edge of the domain.
Yours, Roland
Hello,
I am using IllinoisGRMHD with the test with nsns_test.par in
https://bitbucket.org/zach_etienne/wvuthorns_diagnostics/src/master/NSNS_par... .
The simulation runs well until t~5300 and an error occurs. The error
message is like this:
WARNING level 1 from host nid02210 process 0 while executing schedule bin particle_tracerET, routine
particle_tracerET::DoOneRK4StepForParticleTracerET
in thorn AEILocalInterp, file
/global/u1/c/chlin/et/Cactus/arrangements/Numerical/AEILocalInterp/src/Lagrange-tensor-product/../template.c:1109:
-> CCTK_InterpLocalUniform(): interpolation point is either outside the grid, or inside but too close to the grid boundary! (this may be caused by a global interpolation with driver::ghost_size too small) 0-origin interpolation point number pt=0 of N_interp_points=1 interpolation point (x,y,z)=(-157.725,-576.349,-99.2027) grid x_min(delta_x)x_max = -576(18)576 grid y_min(delta_y)y_max = -576(18)54 grid z_min(delta_z)z_max = -576(18)54
WARNING level 1 from host nid02210 process 0 while executing schedule bin particle_tracerET, routine
particle_tracerET::DoOneRK4StepForParticleTracerET
in thorn particle_tracerET, file
/global/u1/c/chlin/et/Cactus/arrangements/WVUThorns_Diagnostics/particle_tracerET/src/Interpolate_velocities_at_particle_positions.C:131:
-> interpolation screwed up WARNING level 1 from host nid02210 process 1 while executing schedule bin particle_tracerET, routine
particle_tracerET::DoOneRK4StepForParticleTracerET
in thorn particle_tracerET, file
/global/u1/c/chlin/et/Cactus/arrangements/WVUThorns_Diagnostics/particle_tracerET/src/Interpolate_velocities_at_particle_positions.C:131:
-> interpolation screwed up srun: error: nid02211: tasks 2-3: Exited with exit code 1 srun: Terminating job step 21753671.0 srun: error: nid02210: tasks 0-1: Exited with exit code 1
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 :
- It this error can be avoided by decreasing the timestep ?
- If so, how can I decrease the timestep in parameter file (or in
other places)?
Thank you. Best regards, Chia-Hui
-- My email is as private as my paper mail. I therefore support encrypting and signing email messages. Get my PGP key from http://pgp.mit.edu . _______________________________________________ Users mailing list Users@einsteintoolkit.org http://lists.einsteintoolkit.org/mailman/listinfo/users