You probably just need to add "matrix->close()" somewhere in your code.
Might need to put print statements in the code to figure out where the error is... I can never get this one to show up in the debugger properly. David On 08/21/2012 10:38 AM, DYMITROWSKA Magdalena wrote: > Hello all libmesh community, > > I am a new user of our home made libmesh-based code for solving diphase flow > in porous media. > At the moment I have troubles with switching from libmesh-0.7.1 with > petsc-3.0.0 to libmesh-0.7.2 with petsc-3.2. > I have seen a similar problem discussed in January: > http://sourceforge.net/mailarchive/forum.php?thread_name=4F1D9C36.9050603%40purdue.edu&forum_name=libmesh-users > however, I must admit, that I do not understand the solution: > http://sourceforge.net/mailarchive/forum.php?thread_name=4F1DA2D3.2040802%40purdue.edu&forum_name=libmesh-users > It seems that I have to modify the initialisations of jacobian and of > residual, but I am not sure where to look for this, in my application or in > libmesh ? > I give the run time error I get at the bottom of my message. > I would greatly appreciate Your help, > best regards > > Magdalena Dymitrowska > Research Engineer > IRSN > France > > > > > > > Time step 1 ; T = 0.31536 > 0 SNES Function norm 1.390060969904e+86 > [0]PETSC ERROR: --------------------- Error Message > ------------------------------------ > [0]PETSC ERROR: Object is in wrong state! > [0]PETSC ERROR: Matrix must be set first! > [0]PETSC ERROR: > ------------------------------------------------------------------------ > [0]PETSC ERROR: Petsc Release Version 3.2.0, Patch 6, Wed Jan 11 09:28:45 CST > 2012 > [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: ./forge_new on a gnu-linux named BERIS-7 by dymitrow Tue Aug > 21 13:46:08 2012 > [0]PETSC ERROR: Libraries linked from > /produits/freeware/petsc-3.2-p6_x86_64/lib > [0]PETSC ERROR: Configure run at Wed Feb 8 16:42:20 2012 > [0]PETSC ERROR: Configure options --download-f-blas-lapack=1 --with-mpi=1 > --download-mpich=/tmp/mpich2-1.4.1p1.tar.gz --with-shared-libraries=1 > --with-scalapack=1 --download-scalapack=/tmp/SCALAPACK-1.7.tar.gz > --with-blacs=1 --download-blacs=/tmp/blacs-dev.tar.gz --with-mumps=1 > --download-mumps=/tmp/MUMPS_4.10.0.tar.gz --with-parmetis=1 > --download-parmetis=1 --with-superlu_dist=1 --with-ptscotch=1 > --download-ptscotch=/tmp/scotch_5.1.11_esmumps.tar.gz > --download-superlu_dist=/tmp/superlu_dist_2.5.tar.gz > --download-f-blas-lapack=1 --docleanup=1 --with-debugging=0 --with-pthread=1 > --with-pthreadclasses=1 --prefix=/produits/freeware/petsc-3.2-p6_x86_64 > [0]PETSC ERROR: > ------------------------------------------------------------------------ > [0]PETSC ERROR: PCSetUp() line 799 in > /produits/sources/petsc/petsc-3.2-p6/src/ksp/pc/interface/precon.c > [0]PETSC ERROR: KSPSetUp() line 260 in > /produits/sources/petsc/petsc-3.2-p6/src/ksp/ksp/interface/itfunc.c > [0]PETSC ERROR: KSPSolve() line 379 in > /produits/sources/petsc/petsc-3.2-p6/src/ksp/ksp/interface/itfunc.c > [0]PETSC ERROR: SNES_KSPSolve() line 3396 in > /produits/sources/petsc/petsc-3.2-p6/src/snes/interface/snes.c > [0]PETSC ERROR: SNESSolve_LS() line 190 in > /produits/sources/petsc/petsc-3.2-p6/src/snes/impls/ls/ls.c > [0]PETSC ERROR: SNESSolve() line 2676 in > /produits/sources/petsc/petsc-3.2-p6/src/snes/interface/snes.c > [0]PETSC ERROR: User provided function() line 381 in > "unknowndirectory/"src/solvers/petsc_nonlinear_solver.C > application called MPI_Abort(comm=0x84000000, 73) - process 0 > [unset]: aborting job: > application called MPI_Abort(comm=0x84000000, 73) - process 0 > ------------------------------------------------------------------------------ > Live Security Virtual Conference > Exclusive live event will cover all the ways today's security and > threat landscape has changed and how IT managers can respond. Discussions > will include endpoint security, mobile security and the latest in malware > threats. http://www.accelacomm.com/jaw/sfrnl04242012/114/50122263/ > _______________________________________________ > Libmesh-users mailing list > [email protected] > https://lists.sourceforge.net/lists/listinfo/libmesh-users ------------------------------------------------------------------------------ Live Security Virtual Conference Exclusive live event will cover all the ways today's security and threat landscape has changed and how IT managers can respond. Discussions will include endpoint security, mobile security and the latest in malware threats. http://www.accelacomm.com/jaw/sfrnl04242012/114/50122263/ _______________________________________________ Libmesh-users mailing list [email protected] https://lists.sourceforge.net/lists/listinfo/libmesh-users
