I've updated the comment in the code
The code below "cheats" from PETSc style
1) VecRestoreArray() is called immediately after VecGetArray() and the
array values are still accessed; the reason for the immediate
restore is that Vec operations are done on some of the vectors dur
What is the reason for the fakery with VecGet/RestoreArray()? If there is no
chance of
overwriting, why restore immediately? If there is, why doesn't it happen?
Matt
--
What most experimenters take for granted before they begin their experiments
is infinitely more interesting than any results