On 02/13/2012 06:19 PM, Mark Wiebe wrote: > It might be nice to turn the matrix class into a short class hierarchy, > something like this: > > class MatrixBase > class DenseMatrix(MatrixBase) > class TriangularMatrix(MatrixBase) # Maybe a few variations of > upper/lower triangular and whether the diagonal is stored > class SymmetricMatrix(MatrixBase) > > These other matrix classes could use packed storage, and could call the > specific optimized BLAS/LAPACK functions to get higher performance when > it is known the matrix is triangular or symmetric. I'm not sure whether > this affects the discussion of the matrix * and \ operators, but it's a > possibility to consider.
I've been working on exactly this (+ some more) in January, and will be continuing to in the months to come. (Can write more tomorrow if anybody's interested -- or email me directly as I don't have a 0.1 release to show yet -- got to go now) Dag _______________________________________________ NumPy-Discussion mailing list NumPy-Discussion@scipy.org http://mail.scipy.org/mailman/listinfo/numpy-discussion
