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.