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
