Seg. fault using Data Fitting Akima spline

Hi,

I am encountering a problem with the Akima spline that causes my code to seg. fault.

Debugging the code with optimisations and debug symbols, I can repeat the problem. It occurs when dfdInterpolate1D is called with a value that is ~7e-18 within the upper boundary of the spline.

I realise that this value at double precision is near as dammit on the upper boundary, and it is a relatively rare occurrence in my code but it sadly it occasionally occurs.

Is there some underlying floating-point tolerance within the interpolation routine that might be causing this issue?

At the moment I can put in some boundary checking to catch values within tolerance to work around the problem, but ideally I would like to avoid this as it will slow down the code!

Thanks,

Ewan