[yt-users] Hop in parallel?
Michael Kuhlen
mqk at astro.berkeley.edu
Fri Feb 24 14:12:32 PST 2012
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 *
* *
*********************************************************************
More information about the yt-users
mailing list