Re: [Bug-gsl] bug with gsl_integration

[Rhys Ulerich]

> I have found a possible issue with the gsl_integration_qng
> and gsl_integration_qag functions.
>
> ...
>
> mex intTwoK_mex.c -I/usr/local/include/ -lmwlapack -lmwblas -lgsl
>

Can you reproduce the issue without using MEX or Matlab?

Looking at your code, the wrap_intfunc and intfunc are relatively simple.
The mex Function has quite a bit of behavior that might segfault.

I am not saying your code is the culprit, but removing that much code from
the debug would be helpful. Just understanding the ultimate arguments of
the GSL function call would help.

- Rhys

>