1) The solution norm you report is NaN, which indicates memory corruption in your code<br><br>2) The number of allocated nonzeros is really suspicious. There are only 2 nonzeros per row?<br><br> Matt<br><br><div class="gmail_quote">
On Tue, Jan 5, 2010 at 12:01 PM, <span dir="ltr"><<a href="mailto:hxie@umn.edu">hxie@umn.edu</a>></span> wrote:<br><blockquote class="gmail_quote" style="border-left: 1px solid rgb(204, 204, 204); margin: 0pt 0pt 0pt 0.8ex; padding-left: 1ex;">
<br>
The information for superlu_dist <br>
Linear solve converged due to CONVERGED_ITS iterations 1<br>
KSP Object:<br>
type: preonly<br>
maximum iterations=1000, initial guess is zero<br>
tolerances: relative=1e-06, absolute=1e-10, divergence=100<br>
left preconditioning<br>
PC Object:<br>
type: lu<br>
LU: out-of-place factorization<br>
matrix ordering: natural<br>
LU: tolerance for zero pivot 1e-12<br>
LU: factor fill ratio needed 0<br>
Factored matrix follows<br>
Matrix Object:<br>
type=mpiaij, rows=291896, cols=291896<br>
package used to perform factorization: superlu_dist<br>
total: nonzeros=0, allocated nonzeros=583792<br>
SuperLU_DIST run parameters:<br>
Process grid nprow 2 x npcol 2 Equilibrate matrix TRUE Matrix input mode 1 Replace tiny pivots TRUE Use iterative refinement FALSE Processors in row 2 col partition 2 Row permutation LargeDiag Column permutation MMD_AT_PLUS_A<br>
Parallel symbolic factorization FALSE Repeated factorization SamePattern<br>
linear system matrix followed by preconditioner matrix:<br>
Matrix Object:<br>
type=mpiaij, rows=291896, cols=291896<br>
total: nonzeros=33670280, allocated nonzeros=36641300<br>
using I-node (on process 0) routines: found 21607 nodes, limit used is 5<br>
Matrix Object:<br>
type=mpiaij, rows=291896, cols=291896<br>
total: nonzeros=33670280, allocated nonzeros=36641300<br>
using I-node (on process 0) routines: found 21607 nodes, limit used is 5<br>
<br>
WARNING: Tolerance not met by convergence norms<br>
-----------------------------------------------------------------------------<br>
Solution norm = NaN<br>
Total CPU time = 39.0 min<br>
<br>
For mumps, it takes only 2min to get some result "Solution norm = 0.255E+03".<br>
<br>
<br>
</blockquote></div><br><br clear="all"><br>-- <br>What most experimenters take for granted before they begin their experiments is infinitely more interesting than any results to which their experiments lead.<br>-- Norbert Wiener<br>