Re: [ublas] Kronecker tensor product

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