[petsc-users] MatMPIBAIJSetPreallocationCSR i and j indices

[Jed Brown]

Barry Smith <bsmith at mcs.anl.gov> writes:

> On Sep 29, 2013, at 2:52 PM, Jed Brown <jedbrown at mcs.anl.gov> wrote:
>
>> Barry Smith <bsmith at mcs.anl.gov> writes:
>>> 
>>>    Also there is no way to get this flag into MatCreateMPIBAIJWithArrays()
>> 
>> The user calls MatSetOption before making that call and changes it
>> afterward if they want.
>
>    Huh?  I can't do
>
>    MatSetOption(?, MAT_ROW_ORIENTED,something)
>    MatCreateMPIBAIJWithArrays(……, &mat);
>
>    because mat doesn't exist before the MatCreate… () call. Now
>    obviously they could call the several routines instead of
>    MatCreateMPIBAIJWithArrays()

Exactly.

>    but then we still need to decide what format
>    MatCreateMPIBAIJWithArrays() should use since the user cannot
>    choose from several. I conclude the best thing to do is to just use
>    the default column oriented as the only possibility.  

Agreed, but I would do that by calling
MatSetOption(A,MAT_ROW_ORIENTED,PETSC_FALSE) inside
MatCreateMPIBAIJWithArrays and restoring it afterward, so that
MatMPIBAIJSetPreallocationCSR does not change that option.
-------------- next part --------------
A non-text attachment was scrubbed...
Name: not available
Type: application/pgp-signature
Size: 835 bytes
Desc: not available
URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20130929/94f526e5/attachment.pgp>