
<html dir="ltr">

<head>

<meta http-equiv="Content-Type" content="text/html; charset=iso-8859-1">

<style id="owaParaStyle" type="text/css">P {margin-top:0;margin-bottom:0;}</style>

</head>

<body ocsi="0" fpstyle="1">

<div style="direction: ltr;font-family: Tahoma;color: #000000;font-size: 10pt;">Hi All,<br>

<br>

<br>

I'm having a problem with the spooles direct solver. I have it working, but unfortunately it seems to be lacking robustness. I'm dealing with the stiffness matrices of a FE program, where some iterations converge, and for others I get an error. It might be

 important to mention I'm dealing with ill-conditioned matrices (condition number ~10^10). Any suggestions to solve this problem?<br>

<br>

<br>

thanks.<br>

<br>

[0]PETSC ERROR: ------------------------------------------------------------------------<br>

[0]PETSC ERROR: Caught signal number 11 SEGV: Segmentation Violation, probably memory access out of range<br>

[0]PETSC ERROR: Try option -start_in_debugger or -on_error_attach_debugger<br>

[0]PETSC ERROR: or see http://www.mcs.anl.gov/petsc/documentation/faq.html#valgrind[0]PETSC ERROR: or try http://valgrind.org on GNU/linux and Apple Mac OS X to find memory corruption errors<br>

[0]PETSC ERROR: likely location of problem given in stack below<br>

[0]PETSC ERROR: ---------------------  Stack Frames ------------------------------------<br>

[0]PETSC ERROR: Note: The EXACT line numbers in the stack are not available,<br>

[0]PETSC ERROR:       INSTEAD the line number of the start of the function<br>

[0]PETSC ERROR:       is given.<br>

[0]PETSC ERROR: [0] MatFactorNumeric_SeqSpooles line 132 src/mat/impls/aij/seq/spooles/spooles.c<br>

[0]PETSC ERROR: [0] MatCholeskyFactorNumeric line 2950 src/mat/interface/matrix.c<br>

[0]PETSC ERROR: [0] PCSetUp_Cholesky line 92 src/ksp/pc/impls/factor/cholesky/cholesky.c<br>

[0]PETSC ERROR: [0] PCSetUp line 810 src/ksp/pc/interface/precon.c<br>

[0]PETSC ERROR: [0] KSPSetUp line 182 src/ksp/ksp/interface/itfunc.c<br>

[0]PETSC ERROR: [0] KSPSolve line 351 src/ksp/ksp/interface/itfunc.c<br>

<br>

</div>

</body>

</html>