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
<<attachment: for_gsl.zip>>
