On 1/23/2012 12:32 AM, Jed Brown wrote:
> On Sun, Jan 22, 2012 at 23:13, Michael Povolotskyi
> <[email protected] <mailto:[email protected]>> wrote:
>
> Dear libmesh developers and users,
> I'd like to upgrade my code to use libmesh-0.7.2 with petsc3.2.
>
> My application uses nonlinear newton solver.
> Unfortunately with the new version of libmesh and petsc I get the
> following error:
> [0]PETSC ERROR: --------------------- Error Message
> ------------------------------------
> [0]PETSC ERROR: Object is in wrong state!
> [0]PETSC ERROR: Matrix must be set first!
> [0]PETSC ERROR:
> ------------------------------------------------------------------------
>
>
> Please send the entire error message (preferably using a debug build,
> otherwise run in a debugger so you can get a stack trace, you can run
> with -on_error_attach_debugger, then type "bt" when the debugging
> xterm comes up).
>
>
> I searched over the web and found a discussion about this issue:
>
> http://www.mail-archive.com/[email protected]/msg03530.html
>
> I tried to follow this thread, but could not find a solution.
> Could you, please, tell me how should I change my application in order
> to have the matrix in a correct state?
> Thank you,
> Michael
>
>
>
> ------------------------------------------------------------------------------
> Try before you buy = See our experts in action!
> The most comprehensive online learning library for Microsoft
> developers
> is just $99.99! Visual Studio, SharePoint, SQL - plus HTML5, CSS3,
> MVC3,
> Metro Style Apps, more. Free future releases when you subscribe now!
> http://p.sf.net/sfu/learndevnow-dev2
> _______________________________________________
> Libmesh-users mailing list
> [email protected]
> <mailto:[email protected]>
> https://lists.sourceforge.net/lists/listinfo/libmesh-users
>
>
Here it is ( I put a break point at PetscError):
Breakpoint 1, 0x00002aaab247b8c0 in PetscError(int, int, char const*,
char const *, char const*, int, PetscErrorType, char const*, ...) ()
from
/apps/rhel5/petsc-3.2-p5/64/nemesis-mkl-intel-12.0.084/linux/lib/libpets
c.so
(gdb) where
#0 0x00002aaab247b8c0 in PetscError(int, int, char const*, char const*,
char co nst*, int, PetscErrorType, char const*, ...) ()
from
/apps/rhel5/petsc-3.2-p5/64/nemesis-mkl-intel-12.0.084/linux/lib/libpets
c.so
#1 0x00002aaab26eb80d in PCSetUp(_p_PC*) ()
from
/apps/rhel5/petsc-3.2-p5/64/nemesis-mkl-intel-12.0.084/linux/lib/libpets
c.so
#2 0x00002aaab2558492 in KSPSetUp(_p_KSP*) ()
from
/apps/rhel5/petsc-3.2-p5/64/nemesis-mkl-intel-12.0.084/linux/lib/libpets
c.so
#3 0x00002aaab2554880 in KSPSolve(_p_KSP*, _p_Vec*, _p_Vec*) ()
from
/apps/rhel5/petsc-3.2-p5/64/nemesis-mkl-intel-12.0.084/linux/lib/libpets
c.so
#4 0x00002aaab27b2024 in SNES_KSPSolve(_p_SNES*, _p_KSP*, _p_Vec*,
_p_Vec*) ()
from
/apps/rhel5/petsc-3.2-p5/64/nemesis-mkl-intel-12.0.084/linux/lib/libpets
c.so
#5 0x00002aaab25694fb in SNESSolve_LS(_p_SNES*) ()
from
/apps/rhel5/petsc-3.2-p5/64/nemesis-mkl-intel-12.0.084/linux/lib/libpets
c.so
#6 0x00002aaab27ae1f6 in SNESSolve(_p_SNES*, _p_Vec*, _p_Vec*) ()
from
/apps/rhel5/petsc-3.2-p5/64/nemesis-mkl-intel-12.0.084/linux/lib/libpets
c.so
#7 0x00002aaab07b816f in
libMesh::PetscNonlinearSolver<double>::solve(libMesh:: ---Type
<return> to continue, or q <return> to quit---
SparseMatrix<double>&, libMesh::NumericVector<double>&,
libMesh::NumericVector<d ouble>&, double, unsigned int) ()
from
/home/ba01/u120/mpovolot/app-nemo/NEMO/libs/libmesh/libmesh/lib/x86_64-u
nknown-linux-gnu_opt/libmesh.so
#8 0x00002aaab07f47dd in libMesh::NonlinearImplicitSystem::solve() ()
from
/home/ba01/u120/mpovolot/app-nemo/NEMO/libs/libmesh/libmesh/lib/x86_64-u
nknown-linux-gnu_opt/libmesh.so
------------------------------------------------------------------------------
Try before you buy = See our experts in action!
The most comprehensive online learning library for Microsoft developers
is just $99.99! Visual Studio, SharePoint, SQL - plus HTML5, CSS3, MVC3,
Metro Style Apps, more. Free future releases when you subscribe now!
http://p.sf.net/sfu/learndevnow-dev2
_______________________________________________
Libmesh-users mailing list
[email protected]
https://lists.sourceforge.net/lists/listinfo/libmesh-users
