[Date Prev][Date Next][Thread Prev][Thread Next][Date Index][Thread Index]

Re: Block sparse matrix with repeated values



The kron here merely involves the identity matrix, so any solution that promotes reusing the memory holding B would be useful.
Hopefully ending with a call to MatMult.
Thanks,
Tim
On Oct 8, 2007, at 7:32 PM, Matthew Knepley wrote:


We have no format that handles Kronecker products explicitly. However,
these generally reduce to one-dimensional versions, and are maintained
unassembled I think. However, maybe Barry has a better suggestion.

  Matt

On 10/8/07, Tim Sears <tim.sears@xxxxxxxxxx> wrote:
I wish to perform the matrix-vector multiplication Ax, where A = (I
'kron' B), and B is a rectangular matrix.
'kron' is the Kronecker product of two matrices.
Clearly A contains many repeated values. Is there a PETSC format that
can store this matrix efficiently?
If so how does one populate it?
Cheers,
--Tim





--
What most experimenters take for granted before they begin their
experiments is infinitely more interesting than any results to which
their experiments lead.
-- Norbert Wiener