[petsc-users] DMDAVecGetArrayReadF90

[Adrian Croucher]

hi, I'm using a SNES to solve a nonlinear problem on a DMDA, using Fortran.

In my function to be minimized I need an array on the solution vector. 
If I use DMDAVecGetArrayF90(), I get an 'object is in wrong state' error 
because the solution vector has been locked read-only by the SNES.

If I'm using DMPLex I use VecGetArrayReadF90() to get a read-only array, 
but for DMDA there doesn't seem to be a corresponding 
DMDAVecGetArrayReadF90(). Is there something else I should use instead?

- Adrian

-- 
Dr Adrian Croucher
Senior Research Fellow
Department of Engineering Science
University of Auckland, New Zealand
email: a.croucher at auckland.ac.nz
tel: +64 (0)9 923 84611