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.

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