[petsc-users] PETSC ERROR in DMGetLocalBoundingBox?

Fri Mar 17 13:02:33 CDT 2023

Hi Matt,

I am following up to check if you can reproduce the problem on your side. 

Thanks and have a great weekend,

Danyang

From: Danyang Su <danyang.su at gmail.com> 
Sent: March 4, 2023 4:38 PM
To: Matthew Knepley <knepley at gmail.com>
Cc: petsc-users at mcs.anl.gov
Subject: Re: [petsc-users] PETSC ERROR in DMGetLocalBoundingBox?

Hi Matt,

Attached is the source code and example. I have deleted most of the unused source code but it is still a bit length. Sorry about that. The errors come after DMGetLocalBoundingBox and DMGetBoundingBox.

-> To compile the code

Please type 'make exe' and the executable file petsc_bounding will be created under the same folder.

-> To test the code

Please go to fold 'test' and type 'mpiexec -n 1 ../petsc_bounding'.

-> The output from PETSc 3.18, error information

input file: stedvs.dat

------------------------------------------------------------------------

global control parameters

------------------------------------------------------------------------

[0]PETSC ERROR: --------------------- Error Message --------------------------------------------------------------

[0]PETSC ERROR: Corrupt argument: https://petsc.org/release/faq/#valgrind

[0]PETSC ERROR: Object already free: Parameter # 1

[0]PETSC ERROR: See https://petsc.org/release/faq/ for trouble shooting.

[0]PETSC ERROR: Petsc Release Version 3.18.3, Dec 28, 2022 

[0]PETSC ERROR: ../petsc_bounding on a linux-gnu-dbg named starblazer by dsu Sat Mar  4 16:20:51 2023

[0]PETSC ERROR: Configure options --with-cc=gcc --with-cxx=g++ --with-fc=gfortran --download-mpich --download-scalapack --download-parmetis --download-metis --download-mumps --download-ptscotch --download-chaco --download-fblaslapack --download-hypre --download-superlu_dist --download-hdf5=yes --download-ctetgen --download-zlib --download-pnetcdf --download-cmake --with-hdf5-fortran-bindings --with-debugging=1

[0]PETSC ERROR: #1 VecGetArrayRead() at /home/dsu/Soft/petsc/petsc-3.18.3/src/vec/vec/interface/rvector.c:1928

[0]PETSC ERROR: #2 DMGetLocalBoundingBox() at /home/dsu/Soft/petsc/petsc-3.18.3/src/dm/interface/dmcoordinates.c:897

[0]PETSC ERROR: #3 /home/dsu/Work/bug-check/petsc_bounding/src/solver_ddmethod.F90:1920

Total volume of simulation domain   0.20000000E+01

Total volume of simulation domain   0.20000000E+01

-> The output from PETSc 3.17 and earlier, no error

input file: stedvs.dat

------------------------------------------------------------------------

global control parameters

------------------------------------------------------------------------

Total volume of simulation domain   0.20000000E+01

Total volume of simulation domain   0.20000000E+01

Thanks,

Danyang

From: Matthew Knepley <knepley at gmail.com <mailto:knepley at gmail.com> >
Date: Friday, March 3, 2023 at 8:58 PM
To: <danyang.su at gmail.com <mailto:danyang.su at gmail.com> >
Cc: <petsc-users at mcs.anl.gov <mailto:petsc-users at mcs.anl.gov> >
Subject: Re: [petsc-users] PETSC ERROR in DMGetLocalBoundingBox?

On Sat, Mar 4, 2023 at 1:35 AM <danyang.su at gmail.com <mailto:danyang.su at gmail.com> > wrote:

Hi All,

I get a very strange error after upgrading PETSc version to 3.18.3, indicating some object is already free. The error is begin and does not crash the code. There is no error before PETSc 3.17.5 versions.

We have changed the way coordinates are handled in order to support higher order coordinate fields. Is it possible

to send something that we can run that has this error? It could be on our end, but it could also be that you are

destroying a coordinate vector accidentally.

  Thanks,

     Matt

        !Check coordinates

        call DMGetCoordinateDM(dmda_flow%da,cda,ierr)

        CHKERRQ(ierr)

        call DMGetCoordinates(dmda_flow%da,gc,ierr)

        CHKERRQ(ierr)

        call DMGetLocalBoundingBox(dmda_flow%da,lmin,lmax,ierr)

        CHKERRQ(ierr)

        call DMGetBoundingBox(dmda_flow%da,gmin,gmax,ierr)

        CHKERRQ(ierr)

[0]PETSC ERROR: --------------------- Error Message --------------------------------------------------------------

[0]PETSC ERROR: Corrupt argument: https://petsc.org/release/faq/#valgrind

[0]PETSC ERROR: Object already free: Parameter # 1

[0]PETSC ERROR: See https://petsc.org/release/faq/ for trouble shooting.

[0]PETSC ERROR: Petsc Release Version 3.18.3, Dec 28, 2022

[0]PETSC ERROR: ../min3p-hpc-mpi-petsc-3.18.3 on a linux-gnu-dbg named starblazer by dsu Fri Mar  3 16:26:03 2023

[0]PETSC ERROR: Configure options --with-cc=gcc --with-cxx=g++ --with-fc=gfortran --download-mpich --download-scalapack --download-parmetis --download-metis --download-mumps --download-ptscotch --download-chaco --download-fblaslapack --download-hypre --download-superlu_dist --download-hdf5=yes --download-ctetgen --download-zlib --download-pnetcdf --download-cmake --with-hdf5-fortran-bindings --with-debugging=1

[0]PETSC ERROR: #1 VecGetArrayRead() at /home/dsu/Soft/petsc/petsc-3.18.3/src/vec/vec/interface/rvector.c:1928

[0]PETSC ERROR: #2 DMGetLocalBoundingBox() at /home/dsu/Soft/petsc/petsc-3.18.3/src/dm/interface/dmcoordinates.c:897

[0]PETSC ERROR: #3 /home/dsu/Work/min3p-dbs-backup/src/project/makefile_p/../../solver/solver_ddmethod.F90:2140

Any suggestion on this?

Thanks,

Danyang

-- 

What most experimenters take for granted before they begin their experiments is infinitely more interesting than any results to which their experiments lead.
-- Norbert Wiener

https://www.cse.buffalo.edu/~knepley/ <http://www.cse.buffalo.edu/~knepley/> 

-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20230317/72c13909/attachment.html>