[petsc-users] Size of off-diagonal matrix

John Fettig john.fettig at gmail.com
Thu Jun 23 09:47:11 CDT 2011


On Thu, Jun 23, 2011 at 10:25 AM, Jed Brown <jed at 59a2.org> wrote:
> On Thu, Jun 23, 2011 at 14:53, John Fettig <john.fettig at gmail.com> wrote:
>>  Is there
>> anyplace in the documentation that would help me understand this?
>> Alternatively, where would I look in mpiaij.c?
>
> The MPIAIJ format works with blocks that are in SeqAIJ format. In most
> cases, it's not necessary to be familiar with the internal storage format.
> Are you just curious how it's implemented, or is there something you are
> trying to do that depends on the storage format?

I'm trying to use MatCreateMPIAIJWithSplitArrays, but it isn't clear
to me how the off-diagonal portion is supposed to be constructed.

John


More information about the petsc-users mailing list