[yt-users] Hop in parallel?

[Michael Kuhlen]

Is running Hop in parallel currently broken? I get way fewer halos in
my catalog when I run Hop with 8 processors than with 1. Am I doing
something wrong?

This is my script:

$ cat find_halos.py
from yt.mods import *
pf = load('RD0004/RedshiftOutput0004')
halo_list = HaloFinder(pf, padding=0.02)
halo_list.write_out('RD0004/RedshiftOutput0004_Hop.out')

When I run 'python find_halos.py' I get 3716 halos, but when I run it
with 'mpiexec -np 8 python find_halos.py --parallel' I only get 244.
This is independent of the value of padding used. Note that the
parallel catalog is a subset of the serial one. The problem occurs
within HaloFinder, not during the write_out stage, since
len(halo_list) = 242 on all 8 processors.

Mike


-- 
*********************************************************************
*                                                                   *
*  Dr. Michael Kuhlen              Theoretical Astrophysics Center  *
*  email: mqk at astro.berkeley.edu   UC Berkeley                      *
*  cell phone: (831) 588-1468      B-116 Hearst Field Annex # 3411  *
*  skype username: mikekuhlen      Berkeley, CA 94720               *
*                                                                   *
*********************************************************************