[petsc-dev] Fortran binding for VecGetValuesSection

Jed Brown jed at jedbrown.org
Sun Jun 26 00:20:38 CDT 2022

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 <bourdin at 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

More information about the petsc-dev mailing list