MatGetArrayF90 returns 2d array

[Matt Knepley]

I cannot see the logic in getting the array when we already allow  
getting the aij structures. I cannot condone a function which changes  
it's meaning for every matrix type.

   Matt

 From the phone

On Jun 30, 2009, at 6:37 PM, Jed Brown <jed at 59A2.org> wrote:

> 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
> interface.
>
> Jed
>