[petsc-users] PETSC ERROR in DMGetLocalBoundingBox?

Sat Mar 4 18:38:15 CST 2023

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>
Date: Friday, March 3, 2023 at 8:58 PM
To: <danyang.su at gmail.com>
Cc: <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> 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/

-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20230304/e908b7e4/attachment-0001.html>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: petsc_bounding_check.tar.gz
Type: application/x-gzip
Size: 77551 bytes
Desc: not available
URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20230304/e908b7e4/attachment-0001.gz>