Sorry about our late reply. Isn't VecGetValues similar? Note that VecGetValuesSection currently leaks a mutable pointer and really shouldn't exist like this without a matching VecRestoreValuesSection.
PetscCall(VecGetArray(v, &baseArray)); *values = &baseArray[s->atlasOff[p]]; PetscCall(VecRestoreArray(v, &baseArray)); Can you make a draft MR with what you've tried thus far? Blaise Bourdin <bour...@mcmaster.ca> writes: > Hi, > > I’ve been trying to add a fortran90 binding for VecGetValuesSection, without > success, and could not find any function with a similar prototype to learn > from. Can anybody help? > > Regards, > Blaise > > -- > Professor, Department of Mathematics & Statistics > Hamilton Hall room 409A, McMaster University > 1280 Main Street West, Hamilton, Ontario L8S 4K1, Canada > https://www.math.mcmaster.ca/bourdin | +1 (905) 525 9140 ext. 27243