Are the matrices poorly conditioned for multiplication? Finally, are BLAS/LINPACK/... or other matrix libraries available to you? Very few people need to optimize machine-level matrix multiplication.

I've run into a problem with my project with both the speed and accuracy with a home brew matrix class using doubles. I'm inverting a matrix and multiplying it by a vector and transposed vector ...