<div dir="ltr"><div class="gmail_extra"><div class="gmail_quote">On Fri, Jun 26, 2015 at 12:27 PM, Gianluca Meneghello <span dir="ltr"><<a href="mailto:gianmail@gmail.com" target="_blank">gianmail@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 dir="ltr">Dear all,<div><br></div><div>I would like to solve a PDE discretized on a nonuniform --- but rectangular --- grid and I wanted to use the DM coordinates vector to compute the metric terms by finite difference. The only alternative I see is to recompute the coordinates (and then the metric terms) at every function and jacobian evaluation call.</div><div><br></div><div>My question is what is the best (or even correct) way to provide the coordinates to the newly created da objects. Is there anything like a DMDASetNonUniformCoordinates to which to pass a function computing the coordinates? As far as I can tell the fine grid coordinates are currently linearly interpolated from the coarse grid ones.</div></div></blockquote><div><br></div><div>We do not have that function, however you can operate on the coordinate vector the same as you would any other vector, namely</div><div>you call DMDAVecGetArrayDOF() on the coordinate vector, write your nested loops over the vertices, and set coords[k][j][i][d] to</div><div>the d-th coordinate of the point.</div><div><br></div><div>You can ask for the interpolator from the coordinate DM to project these to coarse levels, or just set them directly.</div><div><br></div><div>Does this make sense?</div><div><br></div><div>  Thanks,</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"><div dir="ltr"><div>Please also let me thank you for your great work: it has been and it currently is of enormous help.</div><div><br></div><div>Best</div><span class="HOEnZb"><font color="#888888"><div><br></div><div>Gianluca</div></font></span></div>
</blockquote></div><br><br clear="all"><div><br></div>-- <br><div class="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>