Boost logo

Ublas :

From: Paul C. Leopardi (paul.leopardi_at_[hidden])
Date: 2008-03-04 18:40:25


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