MatGetArrayF90 returns 2d array

Jed Brown jed at
Tue Jun 30 18:37:25 CDT 2009

Matthew Knepley wrote:
> I thought the idea was that MatGetArray() never applies to a sparse
> matrix.  No other sparse format supports this, does it?

That's not true at all, but the result is implementation-dependent.  For
example, the array for AIJ is different from the array for BAIJ.  For
this reason, you shouldn't be calling MatGetArray unless you know the
matrix type, but of course the F90 interface should agree with the C


-------------- next part --------------
A non-text attachment was scrubbed...
Name: signature.asc
Type: application/pgp-signature
Size: 260 bytes
Desc: OpenPGP digital signature
URL: <>

More information about the petsc-users mailing list