[Yt-dev] RMS mass overdensity method

[Stephen Skory]

Yo,

to follow up, I realized that I can define a pf.h.sphere outside of the box, so I am finding particles using .sphere, but I am checking the spheres to see if they overlap with a box side. If they do, I reflect around the domain and make a sphere with a center outside the box, and get the particles that way.

http://paste.enzotools.org/show/64/

I've run this in serial and in parallel (2 procs) on a small dataset, and with 10,000 samples the two give answers different by six parts in one thousand. I think they would be identical if I pre-computed the random centers. Do we think I should do that?

Comments welcome! I'd like to eventually commit this, so I'd like to know what I should fix before I do that, if you can catch any bugs. Also, I haven't decided on a good name for this method. What I have is not very good.

Thanks!

 _______________________________________________________
sskory at physics.ucsd.edu           o__  Stephen Skory
http://physics.ucsd.edu/~sskory/ _.>/ _Graduate Student
________________________________(_)_\(_)_______________