<html>
<head>
<meta content="text/html; charset=utf-8" http-equiv="Content-Type">
</head>
<body bgcolor="#FFFFFF" text="#000000">
<br>
<br>
<div class="moz-cite-prefix">On 20/07/15 22:51, Matthew Knepley
wrote:<br>
</div>
<blockquote
cite="mid:CAMYG4GmJ2COQJqx_1vdBEfP6iwn=mrA=oxnBSkH3AeG6Q2dzzA@mail.gmail.com"
type="cite">
<meta http-equiv="Content-Type" content="text/html; charset=utf-8">
<div dir="ltr">
<div class="gmail_extra">
<div class="gmail_quote"><br>
<div>Thats a bug. CreateGhostCells() makes an entirely new
DM, and the SetNatural is a piece of state which did not
get carried over. I will</div>
<div>fix it.</div>
</div>
</div>
</div>
</blockquote>
<br>
OK, that's what I suspected. Thanks.<br>
<blockquote
cite="mid:CAMYG4GmJ2COQJqx_1vdBEfP6iwn=mrA=oxnBSkH3AeG6Q2dzzA@mail.gmail.com"
type="cite">
<div dir="ltr">
<div class="gmail_extra">
<div class="gmail_quote">
<div><br>
</div>
<div>I have a question about this strategy. The design
philosophy for the PETSc discretization stuff is to make
as much</div>
<div>mesh-independent as possible. I usually think of
boundary and initial conditions as analytic functions that
are</div>
<div>applied to a given geometric region on the mesh, and
then projected into a function space once I choose my</div>
<div>discretization. Why doesn't that work in this case?</div>
</div>
</div>
</div>
</blockquote>
<br>
Yeah, I noticed that's the way all the example problems work, and
I'd do it that way too if I could.<br>
<br>
However my initial conditions aren't defined by analytic functions.
Often they are taken from the end result of a previous model run, so
they're defined as values on particular cells. In most such cases I
should be able to load them in from HDF5 files, but even then I
gather I'll still need this natural ordering stuff to make them end
up on the right cells in the mesh, especially as the number of
processors won't necessarily be the same between runs.<br>
<br>
- Adrian<br>
<pre class="moz-signature" cols="72">--
Dr Adrian Croucher
Senior Research Fellow
Department of Engineering Science
University of Auckland, New Zealand
email: <a class="moz-txt-link-abbreviated" href="mailto:a.croucher@auckland.ac.nz">a.croucher@auckland.ac.nz</a>
tel: +64 (0)9 923 84611
</pre>
</body>
</html>