On 11/30/2010 6:03 AM, Anton Voropaev wrote:
> 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.
>
> _______________________________________________
> Maxima mailing list
> Maxima at math.utexas.edu
> http://www.math.utexas.edu/mailman/listinfo/maxima
This sounds like a job for LAPACK, which Raymond Toy has imported into
Maxima.
Are those zero matrices double precision? or complex zeros? Is it
diagonal? :)
I don't see that the variety of matrix multiplication possibilities
is exposed to the user in Maxima -- you may have to call the routines from
Lisp, which means fussing with the data representation, and prefixing
the program names with "?".
So I think the answer is "yes" but it is not completely integrated into
the system at the moment.
RJF