Ingoing ReggeWheelerRadial at large radius

[LeanneDurkan]

The ingoing solution for the ReggeWheelerRadial function struggles at very large radii.

For Example:

[omega]=N[Sqrt[1/100^3], 32]
ReggeWheelerRadial[2, 2,[omega]][[1]][100000]

and

[omega]=N[Sqrt[1/10^3], 32]
ReggeWheelerRadial[2, 2,[omega] ][[1]][10000]

both fail, producing an error similar to:

-0.*10^254 + 0.*10^254 I

However with greatly increased precision, ceases to produce an error.

For Example:

[omega]=N[Sqrt[1/100^3], 100]
ReggeWheelerRadial[2, 2,[omega]][[1]][100000]

will produce a numerical answer with only a few significant figures.

[barrywardell]

This is a fundamental problem as the ingoing solution is hard to compute accurately at large radius. It may be possible to improve the accuracy by using a different MST representation.