<HTML><BODY style="word-wrap: break-word; -khtml-nbsp-mode: space; -khtml-line-break: after-white-space; ">Hi All,<DIV><BR class="khtml-block-placeholder"></DIV><DIV>Is anyone aware of PETSc particle tracking code? I'm looking for something that will do Lagrangian particle tracking in parallel for 2D fluid flow on a regular cartesian mesh.</DIV><DIV><BR class="khtml-block-placeholder"></DIV><DIV>If this doesn't exist, I may try to create it---quick and dirty style. My thinking was to a master-slaves parallelism: keep a complete directory of the particle locations on node 0 then send the entire directory to each node after each timestep to do an update. Each node updates any points that fall into it's subdomain, flags those particle entries, then returns the entire directory to node 0. Node 0 merges the updates into the master directory. </DIV><DIV><BR class="khtml-block-placeholder"></DIV><DIV>This isn't scalable but it is simpler to program. Does anyone have any suggestions about strategy, implementation, PETSc tools that would facilitate this, etc?</DIV><DIV><BR class="khtml-block-placeholder"></DIV><DIV>Cheers</DIV><DIV>Rich</DIV><DIV><BR><DIV> <SPAN class="Apple-style-span" style="border-collapse: separate; border-spacing: 0px 0px; color: rgb(0, 0, 0); font-family: Courier; font-size: 12px; font-style: normal; font-variant: normal; font-weight: normal; letter-spacing: normal; line-height: normal; text-align: auto; -khtml-text-decorations-in-effect: none; text-indent: 0px; -apple-text-size-adjust: auto; text-transform: none; orphans: 2; white-space: normal; widows: 2; word-spacing: 0px; "><DIV style="margin-top: 0px; margin-right: 0px; margin-bottom: 0px; margin-left: 0px; ">Richard Foa Katz NSF IRFP Postdoc</DIV><DIV style="margin-top: 0px; margin-right: 0px; margin-bottom: 0px; margin-left: 0px; "><A href="http://www.damtp.cam.ac.uk/user/rfk22">http://www.damtp.cam.ac.uk/user/rfk22</A></DIV><BR class="Apple-interchange-newline"></SPAN> </DIV><BR></DIV></BODY></HTML>