Block sparse matrix with repeated values

Barry Smith bsmith at mcs.anl.gov
Mon Oct 8 12:23:25 CDT 2007


  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
> >
> >
> >
> 
> 
> 




More information about the petsc-users mailing list