Re: [Bug-gsl] bug with gsl_integration

[Johannes Hendriks]

Dear sir/madam,

I would like to close this thread as it has been resolved. The problem
ended up being the combination of the default error handler along with the
mex compiler. The default error handler would do use 'printf' to output
error message to screen and abort the function. Unfortunately this would
cause matlab to crash as output to display requires mexPrintf when using
mex compiler.

So I recognise this is not a bug in your functions. However, maybe it would
be good to highlight somewhere that if mex is being used then the default
error handler needs to be disabled / replaced. I am sure there are quite a
few academic and research users who might come across this problem.

Regards,
Johannes

On Wed, Oct 17, 2018 at 1:18 PM Johannes Hendriks <
address@hidden> wrote:

> Dear sir/madam
>
> I have found a possible issue with the gsl_integration_qng
> and gsl_integration_qag functions. Under particular conditions when these
> functions fail to meet the specified relative or absolute tolerance,
> instead of failing gracefully they are causing a seg fault.
>
> I am running on osx 10.12.6 and compiling using the matlab mex compiling.
> I have matlab 2017b, and am compiling with the command
>
> mex intTwoK_mex.c -I/usr/local/include/ -lmwlapack -lmwblas -lgsl
>
> I have attached a zip file containing a matlab script that will cause the
> integration routine to fault.
>
> Let me know if you require any further infomation.
>
> Regards,
> Johannes Hendriks
>
>
>