The LinBox matrix-vector multiplication is vectorMul (Vector1 &w, const Matrix &A, const Vector2 &v) // w <- Av vectorAxpyin (Vector1 &y, const Matrix &A, const Vector2 &x) // y <- y+Ax
I don't see a rowvector-matrix mul routine, but then that must be avoided in C storage order (row-major). -- To post to this group, send an email to sage-devel@googlegroups.com To unsubscribe from this group, send an email to sage-devel+unsubscr...@googlegroups.com For more options, visit this group at http://groups.google.com/group/sage-devel URL: http://www.sagemath.org