Unfortunately we do not have code specifically for this purpose. One would need to write a MATSHELL matrix that then uses repeated multiplications against another Mat (representing the B presumably).
Barry On Mon, 8 Oct 2007, 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 > > > > > > > > >
