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 <knep...@gmail.com> Date: Friday, March 3, 2023 at 8:58 PM To: <danyang...@gmail.com> Cc: <petsc-users@mcs.anl.gov> Subject: Re: [petsc-users] PETSC ERROR in DMGetLocalBoundingBox? On Sat, Mar 4, 2023 at 1:35 AM <danyang...@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/
petsc_bounding_check.tar.gz
Description: GNU Zip compressed data
