Kernel crash when evaluating a TeukolskyRadialFunction outside its domain

[barrywardell]

When evaluating a TeukolskyRadialFunction (constructed using the NumericalIntegration method) outside of the domain on which it has been computed the Mathematica kernel crashes. This appears to be triggering a bug in Mathematica's InterpolatingFunction, which causes the kernel to crash when an InterpolatingFunction is evaluated outside of its domain.

Here is a simple example:

R = TeukolskyRadial[2, 2, -2, 0, 0.08,   Method -> {"NumericalIntegration", "Domain" -> <|"In" -> {3, 5}, "Up" -> {5, 11}|>}];

We get a kernel crash with R["In"][2.9] and messages plus a kernel crash with R["Up"][11.1]. On the other hand, R["In"][5.1] and R["Up"][4.9] both work fine.

[barrywardell]

This has been resolved by fixing issue #40.