Hi Stephen, But is the derived quantity "IsBound" also not parallel in general? The code I'm using does not use clump object at all; all it's doing is picking out a sphere object centred around some point in the box (i.e. I have a separate code to find all the clumps and output all the relevant quantities such as the x,y,z of the densest point and the centre of mass, etc. onto a text file. In other words, the code I'm using is only retrieving the coordinates). I thought the reason why clump finding does not work in parallel was because parallel contour algorithm was not yet implemented, which, if I understand correctly, is separate from "IsBound". In any case, I'll try changing the opening angle. Thanks! Best Regards, Eve On 09/28/2011 11:19 AM, Stephen Skory wrote:
Hi Eve,
I'm trying to use mpi4py as this iteration takes a very long time on serial but I'm getting an error very similar to the error I received when I tried mpi4py with clump finding routine: Clump finding does not work in parallel - at least the expensive part of it (where the gravitational boundedness of the clump is tested) is definitely not parallel. You might want to experiment with changing your opening angle to see if it speeds things up enough to justify the loss in accuracy:
http://yt-project.org/doc/analysis_modules/clump_finding.html
Good luck!