Block sparse matrix with repeated values

Tim Sears tim.sears at anu.edu.au
Mon Oct 8 04:49:20 CDT 2007


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
>




More information about the petsc-users mailing list