Barry Smith <[email protected]> writes: > On Sep 29, 2013, at 2:52 PM, Jed Brown <[email protected]> wrote: > >> Barry Smith <[email protected]> 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.
pgpKzOpOtVb2q.pgp
Description: PGP signature
