Many people doing "serious" matrix stuff, rely on BLAS, adding LaPack /ATLAS(normal matrices) orUMFPACK(sparse matrices) for more advanced math. The reason is that this code is well-tested, stable, reliable, and quite fast. Furthermore, you can buy them directly from a vendor (e.g. ...