[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