[petsc-users] Writing a domain decomposition code with PETSc

[Åsmund Ervik]

Hi again,

Regarding the routine DMDAVecGetArrayF90(): I think this warrants it's
own manual page, in addition to the mention on the DMDAVecGetArray()
manpage.

A specific point that confused me is that one should *not* allocate
the pointer array going into this routine. This was far from obvious
from the existing manpage, I only realised it after staring at
dm/ex11f90.F.

Åsmund


On 07. okt. 2013 17:55, Jed Brown wrote:
> Åsmund Ervik <asmund.ervik at ntnu.no> writes:
>> Thanks for this comment; I have been suspecting this is a bad
>> idea, nice to get confirmation. I will follow your advice.
> 
> Cool.
> 
>> Other than that, does the mapping from local to global using
>> Fortran assumed-shape arrays make sense, or is there a better way
>> of doing it?
> 
> It looks okay to me.
>