Multiplication of numerical matrices



Are there optimized matrix multiplication routines
for specific datatypes (like in Maple)?

For example, multiplication of two 1000 by 1000 zero
matrices takes few minutes in both Maxima and Maple.
If we specify the option datatype=float[8] in Maple,
multiplication will take only a second.