Fix numerical stability problems with RaySolver

org-arl / AcousticRayTracers.jl

Differentiable acoustic ray tracers

MIT License

6 stars 0 forks source link

Fix numerical stability problems with RaySolver #5

Open mchitre opened 3 years ago

mchitre commented 3 years ago

Observed occasional numerical stability problems in some cases:

long range and complicated bathymetry
complicated sound speed profile with :smooth interpolation

In most cases these are resolved by using a lower tolerance or a stiffer solver, but not always. For now :smooth interpolation on SSP should be considered "broken".

mchitre commented 3 years ago

The problem with :smooth SSP seems to do with q becoming negative. Perhaps the curvature in SSP is leading to caustics that are not dealt with correctly. Related to org-arl/AcousticRayTracers.jl#6 ?