BLAS is safe (numpy.dot releases GIL during call). The LAPACK guys also try to maintain threadsafety, but as far as I know this is not 100% reliable at least in the older versions, as some of their routines use fortran-allocated temp arrays, which gfortran makes static if too big (-> fun times). So maybe better leave `nogil` off here to synchronize it here...