efficient matrix block assignment

Matthew Knepley knepley at gmail.com
Thu Aug 6 12:45:35 CDT 2009

On Thu, Aug 6, 2009 at 12:29 PM, Umut Tabak <u.tabak at tudelft.nl> wrote:

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

No, both "cans" are correct. You can do both things.


> Thanks for the clarifications,
> Umut

What most experimenters take for granted before they begin their experiments
is infinitely more interesting than any results to which their experiments
-- Norbert Wiener
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20090806/f7650a37/attachment.htm>

More information about the petsc-users mailing list