<div dir="ltr"><div dir="ltr">On Wed, Dec 21, 2022 at 1:34 PM Matteo Semplice <<a href="mailto:matteo.semplice@uninsubria.it">matteo.semplice@uninsubria.it</a>> wrote:<br></div><div class="gmail_quote"><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex">Dear Matt,<br>
<br>
     I am trying to translate the local nCell numbers returned by swarm <br>
particle location routines into the global (i,j,k) indices in the dmda <br>
and I have run into some oddities in the 3d case. I have dug into the <br>
private_DMDALocatePointsIS_2D_Regular and <br>
private_DMDALocatePointsIS_3D_Regular routines and, using 14 points per <br>
direction and 1 ghost for the DMDA, on 2 cpus:<br>
<br>
- the 2d routine thinks to have 13x6 Q1 elements on cpu0 and 13x7 on cpu <br>
1 (which looks correct to me as 6+7 makes 13 elements in the y direction)<br>
<br>
- the 2d routine thinks to have 13x13x6 elements on cpu0 and 13x13x8 on <br>
cpu1 (which looks one too many)<br>
<br>
Removing lines 155-157 from src/dm/impls/da/dageometry.c (patch <br>
attached) seems to fix some things in my code. I have not studied in <br>
depth the code of private_DMDALocatePointsIS_3D_Regular, but I have seen <br>
that the left shift of xs,ys,zs are already performed at lines 137-139.<br>
<br>
Could you have a look, please?<br>
<br>
(This looked more appropriate for petsc-dev than petsc-users, but I am <br>
not a subscriber, so could you please reply to me directly as well?)<br></blockquote><div><br></div><div>Dang! You can see it was fixed in 2D but not in 3D. Thanks for finding this. I put it in this branch</div><div><br></div><div>  knepley/fix-da-point-location</div><div><br></div><div>which is the last one I mailed you.</div><div><br></div><div>   Matt</div><div> </div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex">
Thanks<br>
<br>
     Matteo<br>
</blockquote></div><br clear="all"><div><br></div>-- <br><div dir="ltr" class="gmail_signature"><div dir="ltr"><div><div dir="ltr"><div><div dir="ltr"><div>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><br></div><div><a href="http://www.cse.buffalo.edu/~knepley/" target="_blank">https://www.cse.buffalo.edu/~knepley/</a><br></div></div></div></div></div></div></div></div>