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.

  Matt


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