[yt-users] volume rendering in parallel

[Luigi Iapichino]

Dear all,

I attach to this e-mail a sample script for volume rendering. With YT 1.6 I 
used the old version of volume rendering (also in parallel) without problems, 
but now (I run on trunk, with an updated revision) I can run this script only 
on one CPU. If I try a parallel run, I get the following output:

--------------
mpirun -np 4 python2.6 ./sample-vr.py --parallel
yt         INFO       2010-08-20 16:53:08,959 Parallel computation enabled: 
0 / 4
yt         INFO       2010-08-20 16:53:08,959 Parallel computation enabled: 
2 / 4
yt         INFO       2010-08-20 16:53:08,959 Parallel computation enabled: 
1 / 4
yt         INFO       2010-08-20 16:53:08,959 Parallel computation enabled: 
3 / 4

... etc etc...

P001 yt         INFO       2010-08-20 16:54:51,234 Finishing 'Ray casting '
P002 yt         INFO       2010-08-20 16:54:51,383 Finishing 'Ray casting '
P003 yt         INFO       2010-08-20 16:55:02,526 Finishing 'Ray casting '
P000 yt         INFO       2010-08-20 16:55:02,791 Finishing 'Ray casting '
Traceback (most recent call last):
  File "./sample-vr.py", line 29, in <module>
    vp.ray_cast()
  
File "/home/hlrb2/pr95he/lu79qey3/remotefs/yt-trunk-svn/yt/extensions/volume_rendering/software_sampler.py", 
line 117, in ray_cast
    if finalize: self._finalize()
  
File "/home/hlrb2/pr95he/lu79qey3/remotefs/yt-trunk-svn/yt/extensions/volume_rendering/software_sampler.py", 
line 121, in _finalize
    im, f = self._mpi_catrgb((self.image, self.resolution))
  
File "/home/hlrb2/pr95he/lu79qey3/remotefs/yt-trunk-svn/yt/lagos/ParallelTools.py", 
line 165, in passage
    return func(self, data)
  
File "/home/hlrb2/pr95he/lu79qey3/remotefs/yt-trunk-svn/yt/lagos/ParallelTools.py", 
line 852, in _mpi_catrgb
    new_image[0:nsize[0],0:nsize[1],:] = data[:]
ValueError: shape mismatch: objects cannot be broadcast to a single shape

----------------

Any suggestion?
Thanks in advance,

  Luigi

-- 

---------------------------------------------------------------

Luigi Iapichino
Zentrum fuer Astronomie der Universitaet Heidelberg
Institut fuer Theoretische Astrophysik
Albert-Ueberle-Str. 2, D-69120 Heidelberg, Germany
Tel: +49 6221 548983, Fax: +49 6221 544221
e-mail: luigi at ita.uni-heidelberg.de
URL: http://www.ita.uni-heidelberg.de/~luigi/
-------------- next part --------------
A non-text attachment was scrubbed...
Name: sample-vr.py
Type: application/x-python
Size: 744 bytes
Desc: not available
URL: <http://lists.spacepope.org/pipermail/yt-users-spacepope.org/attachments/20100820/75029570/attachment.bin>