Hi list,
I need to multiply a triangular matrix by a normal (dense) matrix.
As far as I can tell ublas does not take advantage of the knowledge that (roughly) half the elements of the triangular matrix are zero. I've tried prod() and axpy_prod().
Is there a way to do this efficiently with Boot.Ublas or do I have to roll my own product?
Also, is there a way to detect that a matrix expression will produce a triangular 'shape'?
Thanks in advance!