source code (no blas)