On Mon, 14 Jun 2010 07:23:20 -0500, Matthew Knepley <knepley at gmail.com> wrote: > Is this a new function? Or do you mean just VecGetArray?
This introduces VecRestoreArrayRead, but the array is still not marked const. http://petsc.cs.iit.edu/petsc/petsc-dev/rev/bf752c718a45 It actually has different behavior than the usual VecRestoreArray, but I would prefer a symmetric interface with a VecGetArrayRead that is equivalent to VecGetArray except for a const array. This is as much an issue about symmetry and grepability as not littering the code with (PetscScalar**) casts. Jed
