[petsc-users] PETSC ERROR in DMGetLocalBoundingBox?

Matthew Knepley knepley at gmail.com
Sun Apr 30 09:40:59 CDT 2023 

On Sun, Apr 30, 2023 at 2:57 AM Danyang Su <danyang.su at gmail.com> wrote:

> Hi Matt,
>
>
>
> Just let you know that the error problem in DMGetLocalBoundingBox seems
> fixed in the latest dev version. I didn’t catch the error information any
> more.
>

Sorry, I forgot to mail back. Thanks for reporting that.

   Matt


>
>
> Regards,
>
>
>
> Danyang
>
>
>
> *From: *<danyang.su at gmail.com>
> *Date: *Friday, March 17, 2023 at 11:02 AM
> *To: *'Matthew Knepley' <knepley at gmail.com>
> *Cc: *<petsc-users at mcs.anl.gov>
> *Subject: *RE: [petsc-users] PETSC ERROR in DMGetLocalBoundingBox?
>
>
>
> 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>
> *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/
> <http://www.cse.buffalo.edu/~knepley/>
>


-- 
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/20230430/7c087399/attachment.html>

More information about the petsc-users mailing list