particle tracking on a DA

[Richard Katz]

Hi All,

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.

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.

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?

Cheers
Rich

Richard Foa Katz     NSF IRFP Postdoc
http://www.damtp.cam.ac.uk/user/rfk22


-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/petsc-dev/attachments/20070926/aa510a40/attachment.html>