[Yt-dev] RMS mass overdensity method

[Stephen Skory]

Matt,

I set logging to DEBUG, and there aren't any hanging barriers. There is a barrier as part of the _mpi_allsum call, and I see all the processors reach it, but after that I don't see any barrier calls at all. I've removed or commented out all of my manual MPI.COMM_WORLD.Barrier()s.

I've been working on trying to narrow down the focus, and it seems that cutting a sphere out many times causes problems. Below, between lines 67 and 73 you can see that I've de-randomized the center, so the sphere is the same every time. This eliminates potential problems of reading data outside each sub-box. I am also not retrieving ParticleMassMsun at all. It still crashes at variable points, one time similar to a standard run, one time a bit further.

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

The error messages with this aren't any more helpful, unfortunately.

Do you have ideas of what I should try next? Is there anywhere I can put some prints that would help us debug this? Should I switch tatics and instead extract the particle data out and use a KDTree?

http://www.scipy.org/Cookbook/KDTree

Thanks!

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