> 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 >