[ublas] triangular matrix product

Subject: [ublas] triangular matrix product
From: dariomt_at_[hidden]
Date: 2012-01-16 05:44:34

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!