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

[Smith, Barry F.]

> On Nov 1, 2019, at 4:50 PM, Zhang, Junchao via petsc-users <petsc-users at mcs.anl.gov> wrote:
> 
> I know nothing about Vec FFTW,

  You are lucky :-)

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