Floating Point Matrix Multiplication

I want to multiply a 56x56 matrix with a 56x1 matrix in floating point. There is a altfp_matrix_mult megafunction which does not compute this quickly enough. I am looking for ideas to implement ...