Sajid,
This seems like a good fix.
Barry
> On Jul 27, 2020, at 11:34 PM, Sajid Ali <[email protected]>
> wrote:
>
> 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 <http://s-sajid-ali.github.io/>
