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 at anu.edu.au> 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 >
