[petsc-users] VecLoad into a SubVector ?

[Sajid Ali]

Hi Barry/Matt,

The fix to this bug would be to disable replacearray op on a subvector. I
modified the source code for vecio.c forcing VecLoad_HDF5 to always perform
an array copy and the above test passes for both binary and hdf5 viewers in
serial and parallel.

I can open a PR that adds a line Z->ops->replacearray = NULL; at line 1286
in the rvector.c file if one of you can confirm that the above logic is
correct. The example attached in the last email could be used as a test for
the same if necessary.

Thank You,
Sajid Ali | PhD Candidate
Applied Physics
Northwestern University
s-sajid-ali.github.io
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20200727/89875cdc/attachment.html>