<div dir="ltr">Hi,<div><br></div><div>I am trying to change my problem to using periodic boundary conditions.</div><div><br></div><div>However, when I use DMDASetUniformCoordinates on the DA, the spacing changes.</div><div><br></div><div>This is due to an additional point e.g. in dm/impls/da/gr1.c</div><div><br></div><div><div>else if (dim == 2) {</div><div>    if (bx == DM_BOUNDARY_PERIODIC) hx = (xmax-xmin)/(M);</div><div>    else hx = (xmax-xmin)/(M-1);</div><div>    if (by == DM_BOUNDARY_PERIODIC) hy = (ymax-ymin)/(N);</div><div>    else hy = (ymax-ymin)/(N-1);</div></div><div><br></div><div>I don't understand the logic here, since xmin an xmax refer to the physical domain, how does changing to a periodic BC change the discretization ?</div><div><br></div><div><br></div><div>Could someone clarify or point to a reference ?</div><div><br></div><div>Thanks</div></div>