<div class="gmail_quote">On Mon, Feb 21, 2011 at 21:39, Lisandro Dalcin <span dir="ltr"><<a href="mailto:dalcinl@gmail.com">dalcinl@gmail.com</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex;">
<div id=":2ol">Perhaps I'm missing something obvious, but DA refine/coarsen with<br>
stencil_width=0 for nproc>1 is failing, see the example below</div></blockquote></div><br><div>Lisandro and I chatted about this earlier. The problem is that getting interpolation will fail because a coarse basis function always crosses the process boundary. You can have the locality if you use Q0 interpolation (maybe that should be set as default for stencil_width=0?), but that will fail even in serial unless you set periodic (which can be used to emulate cell-centered).</div>
<div><br></div><div>There is a degenerate case that I think should work in parallel, but does not: Q0 interpolation and periodic refinement. The constraint equations I wrote do not have a solution in this case. I'm not sure if it requires handling as a special case or not. </div>
<div><br></div><div>I put the consistency check in because I didn't like needing to wait for a confusing error in DMGetInterpolation when trying to coarsen in an incompatible way. But perhaps it is useful to refine or coarsen with no intention of having an interpolation operator, in which case the check can be removed.</div>