Re: [petsc-users] PETSC ERROR in DMGetLocalBoundingBox?

2023-04-30 Thread Matthew Knepley
On Sun, Apr 30, 2023 at 2:57 AM Danyang Su  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: *
> *Date: *Friday, March 17, 2023 at 11:02 AM
> *To: *'Matthew Knepley' 
> *Cc: *
> *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 
> *Sent:* March 4, 2023 4:38 PM
> *To:* Matthew Knepley 
> *Cc:* petsc-users@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.2000E+01
>
> Total volume of simulation domain   0.2000E+01
>
>
>
>
>
> -> The output from PETSc 3.17 and earlier, no error
>
> input file: stedvs.dat
>
>
>
> ------------
>
> global control parameters
>
> 
>
>
>
> Total volume of simulation domain   0.2000E+01
>
> Total volume of simulation domain   0.2000E+01
>
>
>
>
>
> Thanks,
>
>
>
> Danyang
>
> *From: *Matthew Knepley 
> *Date: *Friday, March 3, 2023 at 8:58 PM
> *To: *
> *Cc: *
> *Subject: *Re: [petsc-users] PETSC ERROR in DMGetLocalBoundingBox?
>
>
>
> On Sat, Mar 4, 2023 at 1:35 AM  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)
>
> CHKER

Re: [petsc-users] PETSC ERROR in DMGetLocalBoundingBox?

2023-04-30 Thread Danyang Su
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.

 

Regards,

 

Danyang

 

From: 
Date: Friday, March 17, 2023 at 11:02 AM
To: 'Matthew Knepley' 
Cc: 
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  
Sent: March 4, 2023 4:38 PM
To: Matthew Knepley 
Cc: petsc-users@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.2000E+01

Total volume of simulation domain   0.2000E+01

 

 

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

input file: stedvs.dat

 



global control parameters



 

Total volume of simulation domain   0.2000E+01

Total volume of simulation domain   0.2000E+01

 

 

Thanks,

 

Danyang

From: Matthew Knepley 
Date: Friday, March 3, 2023 at 8:58 PM
To: 
Cc: 
Subject: Re: [petsc-users] PETSC ERROR in DMGetLocalBoundingBox?

 

On Sat, Mar 4, 2023 at 1:35 AM  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 --down

Re: [petsc-users] PETSC ERROR in DMGetLocalBoundingBox?

2023-03-17 Thread danyang.su
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  
Sent: March 4, 2023 4:38 PM
To: Matthew Knepley 
Cc: petsc-users@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.2000E+01

Total volume of simulation domain   0.2000E+01

 

 

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

input file: stedvs.dat

 



global control parameters



 

Total volume of simulation domain   0.2000E+01

Total volume of simulation domain   0.2000E+01

 

 

Thanks,

 

Danyang

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

 

On Sat, Mar 4, 2023 at 1:35 AM mailto: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]P

Re: [petsc-users] PETSC ERROR in DMGetLocalBoundingBox?

2023-03-04 Thread Danyang Su
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.2000E+01

Total volume of simulation domain   0.2000E+01

 

 

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

input file: stedvs.dat

 



global control parameters



 

Total volume of simulation domain   0.2000E+01

Total volume of simulation domain   0.2000E+01

 

 

Thanks,

 

Danyang

From: Matthew Knepley 
Date: Friday, March 3, 2023 at 8:58 PM
To: 
Cc: 
Subject: Re: [petsc-users] PETSC ERROR in DMGetLocalBoundingBox?

 

On Sat, Mar 4, 2023 at 1:35 AM  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

Re: [petsc-users] PETSC ERROR in DMGetLocalBoundingBox?

2023-03-03 Thread Matthew Knepley
On Sat, Mar 4, 2023 at 1:35 AM  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-users] PETSC ERROR in DMGetLocalBoundingBox?

2023-03-03 Thread danyang.su
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. 

 

!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_d
dmethod.F90:2140

 

Any suggestion on this?

 

Thanks,

 

Danyang