<div dir="ltr"><div class="gmail_extra"><div class="gmail_quote">On Sat, Aug 2, 2014 at 11:17 AM, Sun, Hui <span dir="ltr"><<a href="mailto:hus003@ucsd.edu" target="_blank">hus003@ucsd.edu</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
I do have my code set up based on ex70. The reason I do this is that I need to use fieldsplit and schur complement. I can't find other examples doing that, except ex55 and ex70, however I have a hard time get ex55 running.<br>

<br>
I don't understand what you said about "This example uses a simple decomposition, not what you use above". What do you mean by "simple decomposition"?<br></blockquote><div><br></div><div>I use FieldSplit and schur on SNES ex62.</div>
<div><br></div><div>   Matt</div><div> </div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
________________________________________<br>
From: Jed Brown [<a href="mailto:jed@jedbrown.org">jed@jedbrown.org</a>]<br>
Sent: Friday, August 01, 2014 11:15 PM<br>
To: Sun, Hui; <a href="mailto:petsc-users@mcs.anl.gov">petsc-users@mcs.anl.gov</a><br>
Subject: Re: [petsc-users] Some possible bugs with PETSC_VIEWER_BINARY_?<br>
<br>
"Sun, Hui" <<a href="mailto:hus003@ucsd.edu">hus003@ucsd.edu</a>> writes:<br>
> Here you can see I have 4 processors running at the same time. I read in the output from MATLAB using the<br>
><br>
> command [A,B,C,D,x,b,y] = PetscBinaryRead('binaryoutput').<br>
><br>
><br>
> And I visualize the variable y by the commands:<br>
><br>
> mesh(reshape(y(20001:30000),100,100))<br>
><br>
> mesh(reshape(y(10001:20000),100,100))<br>
><br>
> mesh(reshape(y(1:10000),100,100))<br>
><br>
><br>
> The output plots are different from what I get if I use 1 processor.<br>
<br>
This example uses a simple decomposition, not what you use above.<br>
<br>
Note that this example was contributed and we do not endorse it as<br>
recommended practice.  If your code is already set up very much like<br>
this, it may be worth looking at, but don't write new code in this way<br>
because it is algorithmically constraining.<br>
</blockquote></div><br><br clear="all"><div><br></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
</div></div>