[petsc-users] VecDuplicate for FFTW-Vec causes VecDestroy to fail conditionally on VecLoad

[Zhang, Junchao]

I know nothing about Vec FFTW, but if you can provide hdf5 files in your test, I will see if I can reproduce it.
--Junchao Zhang


On Fri, Nov 1, 2019 at 2:08 PM Sajid Ali via petsc-users <petsc-users at mcs.anl.gov<mailto:petsc-users at mcs.anl.gov>> wrote:
Hi PETSc-developers,

I'm unable to debug a crash with VecDestroy that seems to depend only on whether or not a VecLoad was performed on a vector that was generated by duplicating one generated by MatCreateVecsFFTW.

I'm attaching two examples ex1.c and ex2.c. The first one just creates vectors aligned as per FFTW layout, duplicates one of them and destroys all at the end. A bug related to this was fixed sometime between the 3.11 release and 3.12 release. I've tested this code with the versions 3.11.1 and 3.12.1 and as expected it runs with no issues for 3.12.1 and fails with 3.11.1.

Now, the second one just adds a few lines which load a vector from memory to the duplicated vector before destroying all. For some reason, this code fails for both 3.11.1 and 3.12.1 versions. I'm lost as to what may cause this error and would appreciate any help in how to debug this. Thanks in advance for the help!

PS: I've attached the two codes, ex1.c/ex2.c, the log files for both make and run and finally a bash script that was run to compile/log and control the version of petsc used.


--
Sajid Ali
Applied Physics
Northwestern University
s-sajid-ali.github.io<http://s-sajid-ali.github.io>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20191101/79910655/attachment.html>