[petsc-dev] Calls to DMCreateFieldIS missing from fieldsplit.c

[Jungho Lee]

DMCreateFieldIS, etc., in fieldsplit.c, are back in. Had anyone
rigorously test that part of the code before it went away, though? It
fails in the PetscFree call (line 237):


(serial version of runex28_3)
$ ./ex28 -da_grid_x 20 -snes_converged_reason -snes_monitor_short
-ksp_monitor_short -problem_type 2 -snes_mf_operator -pack_dm_mat_type
aij -pc_type fieldsplit -pc_fieldsplit_type additive
-fieldsplit_u_ksp_type gmres -fieldsplit_k_pc_type jacobi
  0 SNES Function norm 20.0067
[0]PETSC ERROR: PetscTrFreeDefault() called from
PCFieldSplitSetDefaults() line 237 in
/Users/jungholee/projects/petsc-dev-same/src/ksp/pc/impls/fieldsplit/fieldsplit.c
[0]PETSC ERROR: Block at address 0x7fff5fbfa5a0 is corrupted; cannot free;
may be block not allocated with PetscMalloc()
[0]PETSC ERROR: --------------------- Error Message
------------------------------------
[0]PETSC ERROR: Memory corruption!
[0]PETSC ERROR: Bad location or corrupted memory!
[0]PETSC ERROR:
------------------------------------------------------------------------
[0]PETSC ERROR: Petsc Development HG revision:
25d401b609b5244f31ee51f342bdb0f3d0fdb309  HG Date: Thu Mar 15 23:26:29
2012 -0500
[0]PETSC ERROR: See docs/changes/index.html for recent updates.
[0]PETSC ERROR: See docs/faq.html for hints about trouble shooting.
[0]PETSC ERROR: See docs/index.html for manual pages.
[0]PETSC ERROR:
------------------------------------------------------------------------
[0]PETSC ERROR: ./ex28 on a arch-c named mcswl179.mcs.anl.gov by
jungholee Fri Mar 16 11:55:11 2012
[0]PETSC ERROR: Libraries linked from
/Users/jungholee/projects/petsc-dev-same/arch-c/lib
[0]PETSC ERROR: Configure run at Fri Mar 16 00:04:21 2012
[0]PETSC ERROR: Configure options --with-superlu_dist-dir=/opt/local
--with-superlu-dir=/opt/local --with-hypre-dir=/opt/local
--with-sowing-dir=/opt/local --with-c2html-dir=/opt/local
--with-parmetis-dir=/opt/local --with-metis-dir=/opt/local
--with-shared-libraries
[0]PETSC ERROR:
------------------------------------------------------------------------
[0]PETSC ERROR: PetscTrFreeDefault() line 287 in
/Users/jungholee/projects/petsc-dev-same/src/sys/memory/mtr.c
[0]PETSC ERROR: PCFieldSplitSetDefaults() line 237 in
/Users/jungholee/projects/petsc-dev-same/src/ksp/pc/impls/fieldsplit/fieldsplit.c
[0]PETSC ERROR: PCSetUp_FieldSplit() line 366 in
/Users/jungholee/projects/petsc-dev-same/src/ksp/pc/impls/fieldsplit/fieldsplit.c
[0]PETSC ERROR: PCSetUp() line 832 in
/Users/jungholee/projects/petsc-dev-same/src/ksp/pc/interface/precon.c
[0]PETSC ERROR: KSPSetUp() line 261 in
/Users/jungholee/projects/petsc-dev-same/src/ksp/ksp/interface/itfunc.c
[0]PETSC ERROR: KSPSolve() line 385 in
/Users/jungholee/projects/petsc-dev-same/src/ksp/ksp/interface/itfunc.c
[0]PETSC ERROR: SNES_KSPSolve() line 4057 in
/Users/jungholee/projects/petsc-dev-same/src/snes/interface/snes.c
[0]PETSC ERROR: SNESSolve_LS() line 193 in
/Users/jungholee/projects/petsc-dev-same/src/snes/impls/ls/ls.c
[0]PETSC ERROR: SNESSolve() line 3277 in
/Users/jungholee/projects/petsc-dev-same/src/snes/interface/snes.c
[0]PETSC ERROR: main() line 432 in src/snes/examples/tutorials/ex28.c
application called MPI_Abort(MPI_COMM_WORLD, 1) - process 0

On Thu, Mar 15, 2012 at 1:13 PM, Dmitry Karpeev <karpeev at mcs.anl.gov> wrote:
> Okay, thanks.
> Dmitry.