<div style="line-height:1.7;color:#000000;font-size:14px;font-family:arial"><DIV>Hello Matt<BR>    What you said that "I suspect that you do not understand that only part of a Vec is stored on each MPI process" may be the reason.And I have spent lots of time testing this,especially get results from Vec x (the result of Ax=b)to the main program variable.  And I will have a try again.<BR>    Besides, I have a puzzle about the process (mpiexec -n 2).As you know,I only use PETSc in the subroutine(a function named PETSCSOLVE,may be taken as a tool function) and the main program call it for many times. So,where should I call "PetscInitialize" and "PetscFinalize",in the main program or in the subroutine PETSCSOLVE?<BR>    Thanks.                                                  Jim<BR><BR></DIV>
<DIV></DIV>
<DIV id="divNeteaseMailCard"></DIV>
<DIV><BR></DIV>ÔÚ 2012-05-22 19:33:10£¬"Matthew Knepley" <knepley@gmail.com> Ð´µÀ£º<BR>
<BLOCKQUOTE style="BORDER-LEFT: #ccc 1px solid; MARGIN: 0px 0px 0px 0.8ex; PADDING-LEFT: 1ex" id="isReplyContent">On Tue, May 22, 2012 at 5:17 AM, w_ang_temp <SPAN dir="ltr"><<A href="mailto:w_ang_temp@163.com" target="_blank">w_ang_temp@163.com</A>></SPAN> wrote:<BR>
<DIV class="gmail_quote">
<BLOCKQUOTE style="BORDER-LEFT: #ccc 1px solid; MARGIN: 0px 0px 0px 0.8ex; PADDING-LEFT: 1ex" class="gmail_quote">
<DIV style="LINE-HEIGHT: 1.7; FONT-FAMILY: arial; FONT-SIZE: 14px">
<DIV>Hello</DIV>
<DIV>     I am a beginner of PETSc.I find the petsc-users lists is very useful, however I do not know how to ask a question in it. So can you tell me? And below is what I want to ask.<BR>    I am solving Ax=b with PETSc. And I only use it in a fortran subroutine to get x ,then the main function of my program use x to do other things. The question is that: when the process is 1(mpiexec -n 1),the result of my program is ok, but when the process is more than 1 the result is wrong.<BR></DIV></DIV></BLOCKQUOTE>
<DIV><BR></DIV>
<DIV>"Wrong" is not an acceptable description. I suspect that you do not understand that only part of a Vec</DIV>
<DIV>is stored on each MPI process. First, get the book Using MPI, and then read the PETSc manual</DIV>
<DIV>section on vectors.</DIV>
<DIV><BR></DIV>
<DIV>  Thanks,</DIV>
<DIV><BR></DIV>
<DIV>     Matt</DIV>
<DIV> </DIV>
<BLOCKQUOTE style="BORDER-LEFT: #ccc 1px solid; MARGIN: 0px 0px 0px 0.8ex; PADDING-LEFT: 1ex" class="gmail_quote">
<DIV style="LINE-HEIGHT: 1.7; FONT-FAMILY: arial; FONT-SIZE: 14px">
<DIV>    So where could be wrong? Can you give me some hints on how I can debug? Thanks.</DIV>
<DIV>                                  Jim</DIV></DIV></BLOCKQUOTE></DIV>-- <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></BLOCKQUOTE></div><br><br><span title="neteasefooter"><span id="netease_mail_footer"></span></span>