lanl/phoebus

Interpolated cached geometry not working in 1D with particles

Closed this issue · 1 comments

When interpolating cached geometry to particle positions for geodesic updates, in 1D I see segfaults because even for a particle in the domain in X2, the geometry caching code will try to access indices at j=-1. Some additional checking is needed for lower dimensions, I think

Will be resolved by #112