Stavros Macrakis wrote:
>
> I assume you plan to improve the error returns to give a
> human-understandable error rather than a number (and to suppress the
> Fortran-style cruft -- it would be nice to suppress the compilation
It comes directly from the translation of the Fortran code. I have
thought about doing something about it, but wasn't sure exactly what.
It wouldn't be hard to have the function that prints this do nothing at
all. (I think the original slatec code wanted to cause a core dump!)
>
> I don't know much about numeric integration routines; should nintegrate be
> using other, more robust ones in some of the above cases?
nintegrate uses quad_qags or quad_qagi. These are both adaptive
Gauss-Kronrod integration routines. I think they're reasonably robust,
but obviously not perfect because they don't really analyze the
integrand in any meaningful way.
Ray