<div style="line-height:1.7;color:#000000;font-size:14px;font-family:arial"><DIV>Hello,Matt<BR>    Thanks again.<BR>    What you said("Use -ksp_type preonly -pc_type lu") means that I should use a direct solver.And I am puzzled.<BR>    (1)I think this problem (single processor is ok and more than one are problematic)is typical,but I cannot <BR>find useful infomation both in the manual and the petsc-users.So can you give me some explanation for this <BR>phenomenon.<BR>    (2)When I use -pc_type lu,the error message is "Matrix format mpiaij does not have a built-in PETSc LU!".<BR>    (3)The manual says that "This approach prevents the user from the unexpected surprise of having a <BR>corrupted matrix after a linear solve"(P72).So do you mean the problem of mine is here?<BR>    (4)Generally I want to use the iterative approach other than the direct method.<BR>    As a beginne
 r,I have to learn more about the manual.And also I think an experienced senior will be a more<BR> useful help for me because I have gotten stuck on this problem for more than a week.<BR>     Thanks.                                Jim<BR><BR></DIV>
<DIV></DIV>
<DIV id="divNeteaseMailCard"></DIV>
<DIV><BR></DIV>>At 2012-05-26 20:17:42,"Matthew Knepley" <knepley@gmail.com> wrote:<BR>
<BLOCKQUOTE style="BORDER-LEFT: #ccc 1px solid; MARGIN: 0px 0px 0px 0.8ex; PADDING-LEFT: 1ex" id="isReplyContent">>Use -ksp_type preonly -pc_type lu. 
<DIV><BR></DIV>
<DIV> >  Matt</DIV></BLOCKQUOTE></div><br><br><span title="neteasefooter"><span id="netease_mail_footer"></span></span>