[petsc-users] problems with MatCreateSeqAIJWithArrays

Jed Brown jedbrown at mcs.anl.gov
Tue Nov 29 06:47:19 CST 2011


On Tue, Nov 29, 2011 at 06:12, Nisoli Isaia <isaia.nisoli at gmail.com> wrote:

> Sorry, I mixed up things.
> It is probably simply
>
> rowindices[size]=rowindices[size-1];
>

That gives you an empty row.


Some algorithms require diagonal entries, so it's sometimes better to
preallocate them and put an explicit zero there than to skip them.
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20111129/9af6a4d5/attachment.htm>


More information about the petsc-users mailing list