<div dir="ltr">On Thu, Sep 26, 2013 at 1:15 PM, Mark F. Adams <span dir="ltr"><<a href="mailto:mfadams@lbl.gov" target="_blank">mfadams@lbl.gov</a>></span> wrote:<br><div class="gmail_extra"><div class="gmail_quote">
<blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left-width:1px;border-left-color:rgb(204,204,204);border-left-style:solid;padding-left:1ex">><br>
> We have ghost cells in TS ex11, and I think we do them the correct way :) We have a function you register that fills them up, and we<br>
> loop over them explicitly.<br>
><br>
<br>
Well damn Matt, this looks like the right way to do FV in PETSc, regardless of the nice BC abstraction, but am I going to get be seduced and crash into this island and turn into stone?  i.e., can stupid people deal with your Sections and Plexes and crap (Barry?).<br>
</blockquote><div><br></div><div>So the important thing here is the organization of the boundary functions. You call</div><div><br></div><div>  ModelBoundaryRegister(bdFunc, bdMarker)</div><div><br></div><div>where the function is</div>
<div><br></div><div>  PhysicsBoundary_Advect_Inflow(Model mod, PetscReal time, const PetscReal *centroid, const PetscReal *normal, const PetscScalar *xInterior, PetscScalar *xGhost, void *ctx)</div><div><br></div><div>Then we loop over all the mesh pieces marked with each boundary marker, and call the function on those faces. You get the interior</div>
<div>value, and make the ghost accordingly. I think you would do exactly the same thing in your simple example below. The only optimization</div><div>would be in handling the topology, which you do not see anyway (unless you are a masochist).</div>
<div><br></div><div>   Matt</div><div><br></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left-width:1px;border-left-color:rgb(204,204,204);border-left-style:solid;padding-left:1ex">
I'm thinking of a much simpler problem (2D, 9-point stencil, constant coefficient Laplacian) and I want to add a convergence study wrapper around a code like this.  So have an outer loop that refines grids, gets errors, and prints convergence rates.<br>

<br>
Perhaps I could morph this example into a simpler example that shows how to do FV in the canonically blessed way in PETSc, abstracted BC functions, and a convergence test for FMG.  This might be a worthy addition.<br>
<br>
I suppose I could keep TS and SNES even though I just plan on looking at a linear elliptic solve for now.<br>
<br>
Any suggestions?<br>
<span class=""><font color="#888888"><br>
Mark</font></span></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>