Proposed interface to dgemm



As mentioned in another recent message, maxima doesn't have an interface
to dgemm (lapack's general matrix multiplication routine).  In that
message, I supplied a simple (undocumented) interface.  Here is
"official" proposed interface.

dgemm(a, b, options)

where options is a list of maxima keyword arguments:

c    - matrix c
transpose_a - true if a is to be transposed
transpose_b - true if b is to be transposed
alpha - defaults to 1
beta - defaults to 0

dgemm computes alpha*(op(a) . op(b)) + beta*c where op is either the
identity or transpose depending on the value of transpose_a/transpose_b.

Ray