Certain input birth rate landscapes result in SingularException. Changing the order for the interpolation fixes the issue, for example via: interp = PolyharmonicInterpolation.PolyharmonicInterpolator(X,Y,3) in ploidyMovement.jl. May want to add this as user parameter with a default setting.