Stephen,
Your suspicion was correct. For some reason, this simulation has a number
of particles with ids that are -1. I went all the way back to the very
first output of that simulation and found 126 out of 256^3 particles to have
ids of -1. I have absolutely no idea why this is. This simulation data is
more than a year old, so it is possible there was a bug in enzo at the
time. I just tested ParallelHOP out on a much newer dataset and it worked
without a problem.
Thanks a lot for your help and thanks again for writing this code.
Britton
On Tue, Dec 1, 2009 at 12:22 PM, Stephen Skory
Britton,
Ok, I got one of these:
Real index -1!!!
I didn't get a "Padded index -1!!!". You're suspicion seems to have been right.
So it seems to me that you probably have a particle with index == -1. There are a few things we can try. If you have just one with -1, we can change the default 'empty' value in those arrays to something else. I suspect you won't have just one, however, at which point you're kind of up a creek, because that's absolutely required for parallel HOP.
It may be worth your time to fix the -1 IDs to something unique in the enzo dataset.
I should add a check that screens for bad particle index values, like -1. I've learned something today.
Let me know if there's anything else I can help you with!
_______________________________________________________ sskory@physics.ucsd.edu o__ Stephen Skory http://physics.ucsd.edu/~sskory/ http://physics.ucsd.edu/%7Esskory/ _.>/ _Graduate Student ________________________________(_)_\(_)_______________
_______________________________________________ Yt-dev mailing list Yt-dev@lists.spacepope.org http://lists.spacepope.org/listinfo.cgi/yt-dev-spacepope.org