Comment by khawkins

I was actually unfamiliar with ATLAS, but I am familiar with Eigen, which I use whenever implementing matrix operations in C++: http://eigen.tuxfamily.org

With so many available options, you should never really be implementing your own matrix multiplication.