[petsc-users] DMDAVecRestoreArrayF90 error

[Karl Rupp]

Hi,

what you describe looks a lot like memory corruption. Does your code run 
cleanly through valgrind?

Best regards,
Karli


On 07/01/2016 08:33 AM, TAY wee-beng wrote:
> Hi,
>
> I had problems with DMDAVecRestoreArrayF90 last time when I used an old
> version of Intel Fortran compiler. It works fine in gfortran and new
> version of Intel compiler. It was determined as a bug by the PETSc team.
>
> To use in old version of Intel, I had to use -O1 instead of -O3 -ipo in
> subroutines when DMDAVecRestoreArrayF90 is called.
>
> Recently, my cluster was reset and all files were deleted. I upload my
> files and compile my code again. However, this time, with a new version
> of Intel, I got segmentation error with DMDAVecRestoreArrayF90. Changing
> to -O1 works. But I thought it was working fine. So maybe I need to
> check with my admin if the ver before and after reset are the same.
>
> Another thing was using PETSc ver 3.6.4 and 3.7.2. Using v3.6.4
> (compiled with -O3 -ipo), it encountered segmentation err right from the
> code start, when DMDAVecRestoreArrayF90  is called.
>
> However, for v3.7.2 (compiled with -O3 -ipo), it only happened during
> the 2nd time step, when I need to use DMDAVecRestoreArrayF90 in order to
> use KSP to solve the linear equation later on.
>
> So I wonder why v3.7.2 can get pass the 1st time step w/o problem and
> getting the right answer and only give errors at the 2nd time step.
>
> Any explanation for this?
>
> Btw, same result whether MPI is used or not.
>
>
>