Boost logo

Ublas :

From: Nico Galoppo (ngaloppo_at_[hidden])
Date: 2008-03-04 22:12:47


Paul,

thanks for the tip. Unfortunately, the matrices A and B that I'm
dealing with are not monomial. As a matter of fact, A itself is
constructed as a low-rank multiplication: A = U1 * U2' with v1 and v2
both r-by-3 matrices.

--nico

On Tue, Mar 4, 2008 at 6:40 PM, Paul C. Leopardi
<paul.leopardi_at_[hidden]> wrote:
> On Wed, 5 Mar 2008, Nico Galoppo wrote:
> > I'm currently computing the Kronecker tensor product of two matrices A
> > and B, both r-by-r to compute the Kronecker tensor product, i.e. in
> > matrix parlance:
> >
> > C = A(:) * B(:)';
> >
> > I'm currently doing this with BLAS and ublas in a way that I think is
> > pretty dangerous/hacky, but it was the only way I could come up with.
> > If anyone out there has better options, please let me know.
>
> Nico,
> GluCat just computes Kronecker product one element at a time.
> http://glucat.cvs.sourceforge.net/glucat/glucat/glucat/matrix_imp.h?revision=1.14&view=markup
> This works well for GluCat because quite often one of the matrices is monomial
> (only one non-zero per row and column).
> Best, Paul
> _______________________________________________
> ublas mailing list
> ublas_at_[hidden]
> http://lists.boost.org/mailman/listinfo.cgi/ublas
>

-- 
Nico Galoppo :: http://www.ngaloppo.org