From: Nico Galoppo (ngaloppo_at_[hidden])
Date: 2008-03-04 22:12:47
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.
On Tue, Mar 4, 2008 at 6:40 PM, Paul C. Leopardi
> 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.
> GluCat just computes Kronecker product one element at a time.
> 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
-- Nico Galoppo :: http://www.ngaloppo.org