matrix.test + 3 friends crash when linked against libmkl

[Peter]

matrix.test, kernel_pca.test, regression.test and svd.test fail when linking against libmkl.

For failure in matrix.test I've traced down the abortion to this function call in operator*=(const Matrix&)

gsl_blas_dgemm(CblasNoTrans, CblasNoTrans, 1.0, m_, other.m_,
               0.0, blas_result_);

Is this a problem in yat, in libmkl's implementation of cblas or in GSL's wrapper around cblas, i.e., gsl_cblas* functions?

[Peter]

For regression.test I've traced the error to gsl_multifit_linear which is called in MultiDimensional::fit, which in turn is called from Polynomial::fit.

[Peter]

Replying to peter:

For regression.test I've traced the error to gsl_multifit_linear which is called in MultiDimensional::fit, which in turn is called from Polynomial::fit.

And that function calls gsl_linalg_SV_decomp_mod (A, QSI, Q, S, xt); so this error is possibly the same as in SVD.

[Peter]

I no longer have access to the system on which the failure was detected.