efficient matrix block assignment

Umut Tabak u.tabak at tudelft.nl
Thu Aug 6 12:29:24 CDT 2009

Barry Smith wrote:
> On Aug 6, 2009, at 8:55 AM, Umut Tabak wrote:
>    The "blocks" in the matrix are always square. When you call 
> MatSetValues() you can put in a rectangular block of values. When you 
> call MatSetValuesBlock() you can put in a rectangular set of blocks.
>    Barry
Hi Barry,

The blocks issue is clear from my tries. But you use two successive 
"can" for rectangular matrices, I guess the second can is "can not" for 
MatSetValuesBlock(), right?

Thanks for the clarifications,

More information about the petsc-users mailing list