<html dir="ltr">
<head>
<meta http-equiv="Content-Type" content="text/html; charset=iso-8859-1">
<style type="text/css" id="owaParaStyle"></style>
</head>
<body fpstyle="1" ocsi="0">
<div style="direction: ltr;font-family: Tahoma;color: #000000;font-size: 10pt;">That's true. But I'm using finite difference.
<div><br>
</div>
<div>Best,</div>
<div>Hui</div>
<div><br>
<div style="font-family: Times New Roman; color: #000000; font-size: 16px">
<hr tabindex="-1">
<div id="divRpF750021" style="direction: ltr;"><font face="Tahoma" size="2" color="#000000"><b>From:</b> Matthew Knepley [knepley@gmail.com]<br>
<b>Sent:</b> Saturday, August 02, 2014 9:51 AM<br>
<b>To:</b> Sun, Hui<br>
<b>Cc:</b> Jed Brown; petsc-users@mcs.anl.gov<br>
<b>Subject:</b> Re: [petsc-users] Some possible bugs with PETSC_VIEWER_BINARY_?<br>
</font><br>
</div>
<div></div>
<div>
<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" target="_blank">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" target="_blank">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" target="_blank">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>
</div>
</div>
</div>
</div>
</body>
</html>