Thanks Matt,
Instead of constructing the nested vector by VecCreateNest(), I use
VecRestoreSubVector() to solve this issue.
But I got problems for field splitting method,
I use the following options:
#PETSc Option Table entries:
-fieldsplit_0_ksp_converged_reason true
-fieldsplit_0_ksp_rtol 1e-12
-fieldsplit_0_ksp_type cg
-fieldsplit_0_pc_type ilu
-fieldsplit_1_ksp_converged_reason true
-fieldsplit_1_ksp_rtol 1e-12
-fieldsplit_1_ksp_type cg
-fieldsplit_1_pc_type hypre
-ksp_converged_reason
-ksp_rtol 1e-12
-ksp_type fgmres
-ksp_view
-pc_fieldsplit_schur_fact_type upper
-pc_fieldsplit_schur_precondition seflp
-pc_type fieldsplit
#End of PETSc Option Table entries
There are some options that were not used.
WARNING! There are options you set that were not used!
WARNING! could be spelling mistake, etc!
Option left: name:-fieldsplit_0_ksp_converged_reason value: true
Option left: name:-fieldsplit_1_ksp_converged_reason value: true
Option left: name:-pc_fieldsplit_schur_fact_type value: upper
Option left: name:-pc_fieldsplit_schur_precondition value: seflp
Kind Regards,
Shidi
On 2020-01-14 11:53, Matthew Knepley wrote:
It says that memory is being overwritten somewhere. You can track this
down using valgrind,
as it suggests in the error message.
Thanks,
Matt
On Tue, Jan 14, 2020 at 5:44 AM Y. Shidi <[email protected]> wrote:
Dear developers,
I have a 2x2 nested matrix and the corresponding nested vector.
When I running the code with field splitting, it gets the following
errors:
[0]PETSC ERROR: PetscTrFreeDefault() called from
VecRestoreArray_Nest()
line 678 in
/home/ys453/Sources/petsc/src/vec/vec/impls/nest/vecnest.c
[0]PETSC ERROR: Block at address 0x3f95f60 is corrupted; cannot
free;
may be block not allocated with PetscMalloc()
[0]PETSC ERROR: --------------------- Error Message
--------------------------------------------------------------
[0]PETSC ERROR: Memory corruption:
http://www.mcs.anl.gov/petsc/documentation/installation.html#valgrind
[0]PETSC ERROR: Bad location or corrupted memory
[0]PETSC ERROR: See
http://www.mcs.anl.gov/petsc/documentation/faq.html
for trouble shooting.
[0]PETSC ERROR: Petsc Release Version 3.9.3, unknown
[0]PETSC ERROR: 2DPetscSpuriousTest on a arch-linux2-c-debug named
merlin by ys453 Tue Jan 14 10:36:53 2020
[0]PETSC ERROR: Configure options --download-scalapack
--download-mumps
--download-parmetis --download-metis --download-ptscotch
--download-superlu_dist --download-hypre
[0]PETSC ERROR: #1 PetscTrFreeDefault() line 269 in
/home/ys453/Sources/petsc/src/sys/memory/mtr.c
[0]PETSC ERROR: #2 VecRestoreArray_Nest() line 678 in
/home/ys453/Sources/petsc/src/vec/vec/impls/nest/vecnest.c
[0]PETSC ERROR: #3 VecRestoreArrayRead() line 1835 in
/home/ys453/Sources/petsc/src/vec/vec/interface/rvector.c
[0]PETSC ERROR: #4 VecRestoreArrayPair() line 511 in
/home/ys453/Sources/petsc/include/petscvec.h
[0]PETSC ERROR: #5 VecScatterBegin_SSToSS() line 671 in
/home/ys453/Sources/petsc/src/vec/vscat/impls/vscat.c
[0]PETSC ERROR: #6 VecScatterBegin() line 1779 in
/home/ys453/Sources/petsc/src/vec/vscat/impls/vscat.c
[0]PETSC ERROR: #7 PCApply_FieldSplit() line 1010 in
/home/ys453/Sources/petsc/src/ksp/pc/impls/fieldsplit/fieldsplit.c
[0]PETSC ERROR: #8 PCApply() line 457 in
/home/ys453/Sources/petsc/src/ksp/pc/interface/precon.c
[0]PETSC ERROR: #9 KSP_PCApply() line 276 in
/home/ys453/Sources/petsc/include/petsc/private/kspimpl.h
[0]PETSC ERROR: #10 KSPFGMRESCycle() line 166 in
/home/ys453/Sources/petsc/src/ksp/ksp/impls/gmres/fgmres/fgmres.c
[0]PETSC ERROR: #11 KSPSolve_FGMRES() line 291 in
/home/ys453/Sources/petsc/src/ksp/ksp/impls/gmres/fgmres/fgmres.c
[0]PETSC ERROR: #12 KSPSolve() line 669 in
/home/ys453/Sources/petsc/src/ksp/ksp/interface/itfunc.c
I am not sure why it happens.
Thank you for your time.
Kind Regards,
Shidi
--
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/ [1]
Links:
------
[1] http://www.cse.buffalo.edu/~knepley/
