<div dir="ltr"><div class="gmail_extra"><div class="gmail_quote">On Sat, Jun 25, 2016 at 8:03 PM, Barry Smith <span dir="ltr"><<a href="mailto:bsmith@mcs.anl.gov" target="_blank">bsmith@mcs.anl.gov</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><br>
 On Jun 25, 2016, at 7:56 PM, Matthew Knepley <<a href="mailto:knepley@gmail.com">knepley@gmail.com</a>> wrote:<br>
><br>
><br>
> On Sat, Jun 25, 2016 at 7:18 PM, Xiangdong <<a href="mailto:epscodes@gmail.com">epscodes@gmail.com</a>> wrote:<br>
> Maybe my previous email is not clear. What is the relationship between the parameter x in DMDASetOverlap and the parameter s in DMDACreate2d? Does x need to be smaller than s? Are they completely independent?<br>
><br>
> Yes. The Overlaps are only used in decomposing the DA into small blocks for things like block GS and NASM.<br>
<br>
The manual page for DMDASetOverlap() is pitiful, it could at least indicate what this overlap thing is. Even the sentence above has far more information than the manual page but is still completely unclear.<br></blockquote><div><br></div><div>I agree with you. However, I think most of what Peter did here should be canned. I am trying to rewrite</div><div>SNES completely to take these things into account (he ignores BCs most of the time, does not interface</div><div>well with unstructured or SAMR, etc.), but it got bogged down just like me.</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>
>  Thanks,<br>
><br>
>     Matt<br>
><br>
> PetscErrorCode  DMDASetOverlap(DM da,PetscInt x,PetscInt y,PetscInt z)<br>
><br>
> PetscErrorCode  DMDACreate2d(MPI_Comm comm,DMBoundaryType bx,DMBoundaryType by,DMDAStencilType stencil_type, PetscInt M,PetscInt N,PetscInt m,PetscInt n,PetscInt dof,PetscInt s,const PetscInt lx[],const PetscInt ly[],DM *da)<br>
><br>
> Thanks.<br>
><br>
> Xiangdong<br>
><br>
><br>
> On Fri, Jun 24, 2016 at 5:23 PM, Barry Smith <<a href="mailto:bsmith@mcs.anl.gov">bsmith@mcs.anl.gov</a>> wrote:<br>
><br>
>    What arguments are you referring to?<br>
><br>
> PetscErrorCode  DMDACreate2d(MPI_Comm comm,DMBoundaryType bx,DMBoundaryType by,DMDAStencilType stencil_type,<br>
>                           PetscInt M,PetscInt N,PetscInt m,PetscInt n,PetscInt dof,PetscInt s,const PetscInt lx[],const PetscInt ly[],DM *da)<br>
><br>
> the argument s is the stencil width and depends on the finite difference stencil you want to use. There is no da_overlap argument.<br>
><br>
>   Barry<br>
><br>
> > On Jun 24, 2016, at 4:17 PM, Xiangdong <<a href="mailto:epscodes@gmail.com">epscodes@gmail.com</a>> wrote:<br>
> ><br>
> > Hello everyone,<br>
> ><br>
> > I have a quick question. What is the relationship between da_overlap and the width parameter in dmdacreate2d?<br>
> ><br>
> > From the examples, it seems that width=1 works with da_overlap=4. Could you please explain the relationship or point me to a reference page?<br>
> ><br>
> > Thanks.<br>
> ><br>
> > Xiangdong<br>
><br>
><br>
><br>
><br>
<span class="HOEnZb"><font color="#888888">><br>
> --<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>
<br>
</font></span></blockquote></div><br><br clear="all"><div><br></div>-- <br><div class="gmail_signature" data-smartmail="gmail_signature">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>