From: Gunter Winkler (guwi17_at_[hidden])
Date: 2008-03-05 17:44:28
Nico Galoppo wrote:
> 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.
IMO the best solution would be to provide a custom operator like
matrix_matrix_prod which knows how to compute element (i,j) of the result.