Mailman 3 - yt-users

Error in build of kdtree
by Eric Hallman 26 Apr '10

26 Apr '10

I am getting this error after Forthon install and make in the kdtree dir: login4% make Forthon -F gfortran --compile_first fKD_source --no2underscores --with- numpy --fopt "-O3" fKD fKD_source.f90 IN FORTHON_BUILDER Traceback (most recent call last): File "/work/01059/tg803229/raft-0.1a/local/bin/Forthon", line 2, in <module> import Forthon.Forthon_builder File "/work/01059/tg803229/raft-0.1a/local/lib/python2.6/site- packages/Forthon/Forthon_builder.py", line 142, in <module> twounderscores=twounderscores) File "/work/01059/tg803229/raft-0.1a/local/lib/python2.6/site- packages/Forthon/compilers.py", line 66, in __init__ if self.linux_pathscale() is not None: break File "/work/01059/tg803229/raft-0.1a/local/lib/python2.6/site- packages/Forthon/compilers.py", line 352, in linux_pathscale if usecompiler('pathscale','pathf95'): NameError: global name 'usecompiler' is not defined make: *** [fKD] Error 1 Any ideas would help. Thanks! Eric Hallman Google Voice: (774) 469-0278 hallman13(a)gmail.com

3 4

libpng and new volume rendering methods
by Matthew Turk 25 Apr '10

25 Apr '10

Hi everyone, (if you're on yt-1.6, feel free to ignore this message!) I've backported the multivariate volume rendering (for instance: http://yt.enzotools.org/attachment/wiki/Screenshots/hiiregion_rgb.png http://yt.enzotools.org/attachment/wiki/Screenshots/Reion512_600px.png http://yt.enzotools.org/attachment/wiki/Screenshots/SimB-BH.png ) to the yt trunk; this also comes with the ability to project through a simulation along a non-orthogonal vector, so you can create images that, for instance, peer down the L vector of a gas cloud but include contributions from all the intervening gas -- it's much faster and higher quality than stacking cutting planes, which was the old method for doing this. (Example: http://yt.enzotools.org/attachment/wiki/Screenshots/halo_0029.png ) The new backend for volume rendering should be API compatible with the old, but if you update your installation you will have to rebuild the C extensions either by doing: $ yt instinfo -u or by re-installing. The new multivariate volume rendering as well as the off-axis projecting have not yet been documented. This also comes with a direct PNG writer, which is much faster than using matplotlib for outputting raw images; it also comes with a couple colormaps. This requires that libpng be installed. yt will do its best to locate libpng, but if it cannot find it, you will have to supply it in png.cfg, much like for hdf5 -- however, the installation script has been updating to install libpng if requested, so if you run into any problems, rerun a new version of the installation script. (Running it over top of an old installation should be fine, as it will just conduct the necessary upgrades.) For the next release, both the new volume rendering mechanism and the image writer will be documented in detail; the release is probably about a week or two away. Let me know if you're eager to get your hands on the API docs for either of these things and I'll shoot you some sample scripts. If you run into any problems, let us know! Best, Matt

1 0

'dx' field from ortho_ray
by Michael Kuhlen 24 Apr '10

24 Apr '10

It seems that ortho_ray produces a 'dx' field (also 'CellLength') that is not the same order as the other data fields. This is again for a 1D simulation, this time with a few levels of AMR on. ray = pf.h.ortho_ray(0, [0.5, 0.5]) print ray['x'] yt INFO 2010-04-23 16:54:42,181 Getting field x from 6 [ 0.005 0.015 0.025 0.035 0.045 0.055 0.065 0.075 0.085 0.095 0.105 0.115 0.125 0.135 0.145 0.155 0.165 0.175 0.185 0.195 0.205 0.215 0.225 0.235 0.245 0.255 0.265 0.275 0.285 0.295 0.305 0.315 0.325 0.335 0.345 0.355 0.365 0.375 0.385 0.395 0.405 0.415 0.425 0.435 0.445 0.455 0.465 0.475 0.485 0.495 0.505 0.515 0.525 0.535 0.545 0.555 0.5625 0.5675 0.5725 0.5775 0.5825 0.5875 0.595 0.605 0.615 0.625 0.635 0.645 0.655 0.665 0.675 0.685 0.695 0.705 0.7125 0.7175 0.72125 0.72375 0.725625 0.726875 0.728125 0.7290625 0.7296875 0.7303125 0.7309375 0.7315625 0.7321875 0.7328125 0.7334375 0.7340625 0.7346875 0.7353125 0.7359375 0.736875 0.738125 0.739375 0.74125 0.74375 0.7475 0.7525 0.7575 0.765 0.775 0.785 0.795 0.805 0.815 0.825 0.835 0.845 0.855 0.865 0.875 0.885 0.895 0.905 0.915 0.925 0.935 0.945 0.955 0.965 0.975 0.985 0.995 ] print ray['dx'] yt INFO 2010-04-23 16:54:42,232 Getting field dx from 6 [ 0.01 0.01 0.01 0.01 0.01 0.01 0.01 0.01 0.01 0.01 0.01 0.01 0.01 0.01 0.01 0.01 0.01 0.01 0.01 0.01 0.01 0.01 0.01 0.01 0.01 0.01 0.01 0.01 0.01 0.01 0.01 0.01 0.01 0.01 0.01 0.01 0.01 0.01 0.01 0.01 0.01 0.01 0.01 0.01 0.01 0.01 0.01 0.01 0.01 0.01 0.01 0.01 0.01 0.01 0.01 0.01 0.01 0.01 0.01 0.01 0.01 0.01 0.01 0.01 0.01 0.01 0.01 0.01 0.01 0.01 0.01 0.01 0.01 0.01 0.01 0.01 0.01 0.01 0.01 0.01 0.01 0.01 0.01 0.01 0.01 0.01 0.01 0.01 0.01 0.01 0.01 0.01 0.005 0.005 0.005 0.005 0.005 0.005 0.005 0.005 0.005 0.005 0.005 0.0025 0.0025 0.0025 0.0025 0.00125 0.00125 0.00125 0.00125 0.00125 0.00125 0.000625 0.000625 0.000625 0.000625 0.000625 0.000625 0.000625 0.000625 0.000625 0.000625 0.000625 0.000625] Any ideas how to fix this bug? Or some other way to access the ray's 'dx' values in the correct order? I suppose I could simply calculate the dx's myself from the 'x' field... Thanks, Mike -- ********************************************************************* * * * Dr. Michael Kuhlen Theoretical Astrophysics Center * * email: mqk(a)astro.berkeley.edu UC Berkeley * * cell phone: (831) 588-1468 601 Campbell Hall * * skype username: mikekuhlen Berkeley, CA 94720 * * * *********************************************************************

2 1

Volume rendering - filenames
by Mike Butler 22 Apr '10

22 Apr '10

Hello, This is probably a simple thing to fix, but I've been trying to use the volume rendering package on my data and am running into an error. If I follow the steps at http://paste.enzotools.org/show/404/ when I get to mi,ma = dd.quantities["Extrema"]("Density")[0] I get the error message at the bottom of this email. It seems that it is looking for files ending in .cpu000# rather than .grid.cpu000# like the default enzo output. I've tried doing a rename but then I get the same error but the other way around. Making a copy of the file with each filename works if I'm just doing a single image, but I'm trying to make a movie so doing that for 100 files for 8 cpus each would be a bit tedious. If anyone knows a way around this or I'm doing something wrong, I'd appreciate any suggestions. Thanks -Mike In [8]: mi,ma=dd.quantities["Extrema"]("Density")[0] ERROR: An unexpected error occurred while tokenizing input The following traceback may be corrupted or invalid The error message is: ('EOF in multi-line statement', (105, 0)) ERROR: An unexpected error occurred while tokenizing input The following traceback may be corrupted or invalid The error message is: ('EOF in multi-line statement', (72, 0)) ERROR: An unexpected error occurred while tokenizing input The following traceback may be corrupted or invalid The error message is: ('EOF in multi-line statement', (119, 0)) --------------------------------------------------------------------------- ReadingError Traceback (most recent call last) /home/butler85/yt/src/yt-trunk-svn/scripts/iyt in <module>() ----> 1 2 3 4 5 /home/butler85/yt/src/yt-trunk-svn/yt/lagos/DerivedQuantities.pyc in __call__(self, *args, **kwargs) 70 self._data_source.pf.h.io) 71 if lazy_reader and not self.force_unlazy: ---> 72 return self._call_func_lazy(args, kwargs) 73 else: 74 return self._call_func_unlazy(args, kwargs) /home/butler85/yt/src/yt-trunk-svn/yt/lagos/DerivedQuantities.pyc in _call_func_lazy(self, args, kwargs) 77 self.retvals = [ [] for i in range(self.n_ret)] 78 for gi,g in enumerate(self._get_grids()): ---> 79 rv = self.func(GridChildMaskWrapper(g, self._data_source), *args, **kwargs) 80 for i in range(self.n_ret): self.retvals[i].append(rv[i]) 81 g.clear_data() /home/butler85/yt/src/yt-trunk-svn/yt/lagos/DerivedQuantities.pyc in _Extrema(data, fields, filter) 385 mins, maxs = [], [] 386 for field in fields: --> 387 if data[field].size < 1: 388 mins.append(1e90) 389 maxs.append(-1e90) /home/butler85/yt/src/yt-trunk-svn/yt/lagos/DerivedQuantities.pyc in __getitem__(self, item) 41 return getattr(self.grid, attr) 42 def __getitem__(self, item): ---> 43 return self.data_source._get_data_from_grid(self.grid, item) 44 45 class DerivedQuantity(ParallelAnalysisInterface): /home/butler85/yt/src/yt-trunk-svn/yt/lagos/BaseDataTypes.pyc in save_state(self, grid, field) 41 old_keys = grid.data.keys() 42 grid.field_parameters = self.field_parameters ---> 43 tr = func(self, grid, field) 44 grid.field_parameters = old_params 45 grid.data = dict( [(k, grid.data[k]) for k in old_keys] ) /home/butler85/yt/src/yt-trunk-svn/yt/lagos/BaseDataTypes.pyc in _get_data_from_grid(self, grid, field) 1660 else: 1661 pointI = self._get_point_indices(grid) -> 1662 if grid[field].size == 1: # dx, dy, dz, cellvolume 1663 t = grid[field] * na.ones(grid.ActiveDimensions, dtype='float64') 1664 return t[pointI].ravel() /home/butler85/yt/src/yt-trunk-svn/yt/lagos/BaseGridType.pyc in __getitem__(self, key) 133 """ 134 if not self.data.has_key(key): --> 135 self.get_data(key) 136 return self.data[key] 137 /home/butler85/yt/src/yt-trunk-svn/yt/lagos/BaseGridType.pyc in get_data(self, field) 166 return self.data[field] 167 try: --> 168 temp = self.hierarchy.io.pop(self, field) 169 self[field] = na.multiply(temp, conv_factor, temp) 170 except self.hierarchy.io._read_exception, exc: /home/butler85/yt/src/yt-trunk-svn/yt/lagos/DataReadingFuncs.pyc in pop(self, grid, field) 55 else: 56 # We only read the one set and do not store it if it isn't pre-loaded ---> 57 return self._read_data_set(grid, field) 58 59 def peek(self, grid, field): /home/butler85/yt/src/yt-trunk-svn/yt/lagos/DataReadingFuncs.pyc in _read_data_set(self, grid, field) 232 def _read_data_set(self, grid, field): 233 return HDF5LightReader.ReadData(grid.filename, --> 234 "/Grid%08i/%s" % (grid.id, field)).swapaxes(0,2) 235 236 def _read_data_slice(self, grid, field, axis, coord): ReadingError: ReadHDF5DataSet: /home/butler85/sim/starsfuvamrrec/Cloud_128AMR4StarsFUV_0101.cpu0000 does not exist, or no read permissions

2 1

bug in ortho_ray and derived fields?
by Michael Kuhlen 21 Apr '10

21 Apr '10

Hi I'm experiencing some weird behavior with ortho_ray() when used with derived fields. I suspect this is a bug, but I'm new to yt, so maybe I'm doing something incorrectly. This is for a 1D Enzo simulation, and the only fields in the simulation dataset are Density, x-velocity, and TotalEnergy. I defined two derived fields, InternalEnergy and Pressure: ### define InternalEnergy field def _InternalEnergy(field, data): return data['TotalEnergy'] - 0.5*data['x-velocity']*data['x-velocity'] add_field('InternalEnergy', function=_InternalEnergy, units=r'\rm{erg}/\rm{g}') ### define Pressure field def _Pressure(field, data): return (data.pf['Gamma'] - 1.0) * data['Density'] * data['InternalEnergy'] add_field('Pressure', function=_Pressure, units=r'\rm{dyne}/\rm{cm}^{2}') Next I want to extract a ray of Density, x-velocity, Pressure, and InternalEnergy along the only dimension. I figured ortho_ray is the right tool. ray = pf.h.ortho_ray(0, [0.5, 0.5], ['x-velocity', 'Density', 'InternalEnergy', 'Pressure']) print ray.keys() This produces: ['Pressure', 'x', 'x-velocity', 'InternalEnergy', 'Density'] Perfect, everything as desired. However, if I change the order of the input fields, some disappear from the ray. Examples: (1) ray = pf.h.ortho_ray(0, [0.5, 0.5], ['Density', 'x-velocity', 'Pressure', 'InternalEnergy']) print ray.keys() ['x', 'x-velocity', 'Pressure', 'Density'] -- 'InternalEnergy' is missing. (2) ray = pf.h.ortho_ray(0, [0.5, 0.5], ['Pressure', 'Density', 'x-velocity']) print ray.keys() ['Pressure'] -- 'Density' and 'x-velocity' are missing. (3) ray = pf.h.ortho_ray(0, [0.5, 0.5], ['InternalEnergy', 'Density', 'x-velocity']) print ray.keys() ['InternalEnergy', 'Density'] -- 'x-velocity' is missing. The pattern I'm seeing here is that whenever 'InternalEnergy' is specified as an input only one field after it makes it into the ray. When 'Pressure' is specified, no further fields make it. I don't imagine this behavior is intended, but maybe it's due to some user error on my part? Thanks in advance for any clarifications. Mike -- ********************************************************************* * * * Dr. Michael Kuhlen Theoretical Astrophysics Center * * email: mqk(a)astro.berkeley.edu UC Berkeley * * cell phone: (831) 588-1468 601 Campbell Hall * * skype username: mikekuhlen Berkeley, CA 94720 * * * *********************************************************************

2 4

Problem with yt and h5py
by Brian O'Shea 18 Apr '10

18 Apr '10

Dear yt-users, I'm having a strange problem with yt (svn trunk, r1693). If I run this script: from yt.mods import * fn = "DD0050/DD0050" pf = load(fn) dd = pf.h.all_data() I get the following error: yt INFO 2010-04-18 13:41:48,183 Getting the binary hierarchy Traceback (most recent call last): File "<stdin>", line 1, in <module> File "/mnt/home/oshea/yt-x86_64/src/yt-trunk-svn/yt/lagos/OutputTypes.py", line 150, in _get_hierarchy self.__hierarchy = self._hierarchy_class(self, data_style=self.data_style) File "/mnt/home/oshea/yt-x86_64/src/yt-trunk-svn/yt/lagos/HierarchyType.py", line 374, in __init__ AMRHierarchy.__init__(self, pf, data_style) File "/mnt/home/oshea/yt-x86_64/src/yt-trunk-svn/yt/lagos/HierarchyType.py", line 55, in __init__ self._parse_hierarchy() File "/mnt/home/oshea/yt-x86_64/src/yt-trunk-svn/yt/lagos/HierarchyType.py", line 435, in _parse_hierarchy if self._parse_binary_hierarchy(): return File "/mnt/home/oshea/yt-x86_64/src/yt-trunk-svn/yt/lagos/HierarchyType.py", line 506, in _parse_binary_hierarchy self.grid_dimensions[:] = f["/ActiveDimensions"][:] File "/mnt/home/oshea/yt-x86_64/lib/python2.6/site-packages/h5py-1.2.0-py2.6-linux-x86_64.egg/h5py/highlevel.py", line 961, in __getitem__ basetype = self.id.dtype File "h5d.pyx", line 109, in h5py.h5d.DatasetID.dtype.__get__ (h5py/h5d.c:1423) File "h5t.pyx", line 364, in h5py.h5t.TypeID.dtype.__get__ (h5py/h5t.c:2919) File "h5t.pyx", line 367, in h5py.h5t.TypeID.py_dtype (h5py/h5t.c:2976) TypeError: No NumPy equivalent for TypeIntegerID exists I installed this version of yt on MSU's supercomputer using the install script at http://svn.enzotools.org/yt/trunk/doc/install_script.sh (the supercomputer uses gcc v4.1.2), and as far as I can tell everything relevant, including hdf5 and h5py, installed successfully. The same simple set of commands works fine on my laptop with the same dataset. I'm pointing my LD_LIBRARY_PATH environmental variable to the correct directory, so I don't think it's that. Any ideas about what's going on here? Some more system details: The output of uname -a for this system is is "Linux dev-amd09 2.6.16.60-0.21-smp #1 SMP Tue May 6 12:41:02 UTC 2008 x86_64 x86_64 x86_64 GNU/Linux". The system uses gcc v4.1.2 and I'm installing yt _without_ wxPython. Thanks, Brian

1 1

2d radial profiles
by david collins 17 Apr '10

17 Apr '10

Hi, all-- Is it possible to make a 2d radial profile from a projection with the existing profile tool? I think a better question is how do I make an extraction from a projection that I can then hand to BinnedProfile1d? Thanks, d. -- Sent from my Stone Tablet and carried by my Pterodactyl.

4 5

If statement for particles existing?
by Elizabeth Harper-Clark 17 Apr '10

17 Apr '10

Hi guys, I sometimes have particles in my sims and when i do I want to plot them. At the moment I have: if sp["ParticleMassMsun"].size > 0: p.modify["particles"](1.0, p_size=10.0, col='k') but if there are no particles I get this error message: yt INFO 2010-04-16 11:58:29,408 Getting ParticleMassMsun using ParticleIO Traceback (most recent call last): File "Density_movie.py", line 29, in <module> if sp["ParticleMassMsun"].size > 0: File "/scratch/YT-x86_64/lib/python2.5/site-packages/yt-2.0dev-py2.5-linux-x86_64.egg/yt/lagos/BaseDataTypes.py", line 249, in __getitem__ self.get_data(key) File "/scratch/YT-x86_64/lib/python2.5/site-packages/yt-2.0dev-py2.5-linux-x86_64.egg/yt/lagos/BaseDataTypes.py", line 1636, in get_data self[field] = self.particles[field] File "/scratch/YT-x86_64/lib/python2.5/site-packages/yt-2.0dev-py2.5-linux-x86_64.egg/yt/lagos/ParticleIO.py", line 45, in __getitem__ return self.get_data(key) File "/scratch/YT-x86_64/lib/python2.5/site-packages/yt-2.0dev-py2.5-linux-x86_64.egg/yt/lagos/ParticleIO.py", line 91, in get_data conv_factors) File "/scratch/YT-x86_64/lib/python2.5/site-packages/yt-2.0dev-py2.5-linux-x86_64.egg/yt/lagos/DataReadingFuncs.py", line 206, in _read_particles filenames, ids, conv_factors = zip(*sorted(zip(filenames, ids, conv_factors))) ValueError: need more than 0 values to unpack Can anyone think of a better if statement to check if i have particles that doesn't crash YT if there isn't? Many thanks, Libby -- Elizabeth Harper-Clark MA MSci PhD Candidate, Astrophysics, UofT www.astro.utoronto.ca/~h-clark h-clark(a)cita.utoronto.ca AIM: edphc1 MSN: edphc1(a)hotmail.com Skype: eharperclark Office phone: 416-978-5759

3 2

solution to wxPython installation problem
by Michael Kuhlen 15 Apr '10

15 Apr '10

I just installed yt 1.6 (stable branch) for the first time. I ran into some trouble with the wxPython installation. I realize that it's an optional component, but since I found a solution to my problem, I thought I'd share it here. I'm installing yt on an x86_64 machine running Ubuntu 9.04, 2.6.28-18. The error I got during the wxPython installation was: src/gtk/_gdi_wrap.cpp:31407: error: ‘class wxGraphicsRenderer’ has no member named ‘CreateBitmap’ Some googling revealed this: http://groups.google.com/group/wxPython-users/browse_thread/thread/129ba27e… Applying the patch discussed at the end of that thread allowed the compilation of wxPython to complete successfully. To be more specific, you need to add this one line wxGraphicsBitmap CreateBitmap( const wxBitmap &bitmap ) const { return wxNullGraphicsBitmap; } after line 4197 in the file $YT_DIR/src/wxPython-src-2.8.10.1/wxPython/src/gtk/_gdi_wrap.cpp. Might be worth including this patch in the version of wxPython that is downloaded from yt.enzotools.org. Mike -- ********************************************************************* * * * Dr. Michael Kuhlen Theoretical Astrophysics Center * * email: mqk(a)astro.berkeley.edu UC Berkeley * * cell phone: (831) 588-1468 601 Campbell Hall * * skype username: mikekuhlen Berkeley, CA 94720 * * * *********************************************************************

2 1

projection depth
by Luigi Iapichino 13 Apr '10

13 Apr '10

Dear all, by comparing volume renderings and projections in my simulation I realised that I don't have control on the projection depth, in the latter case. This is the relevant part of my script: ... L = 0.1 # projection width pf = load(data_dir+fn) # load data reg = pf.h.region([0.5, 0.5, 0.5],[0.45,0.45,0.45],[0.55,0.55,0.55]) # a small region #reg = pf.h.region([0.5, 0.5, 0.5],[0,0,0],[1,1,1]) # the whole box pc = PlotCollection(reg, center=[0.5,0.5,0.5]) p = pc.add_projection("Density", 0, weight_field="Density",fig_size=(7.5,6)) p.modify["contour"]("Density", ncont=5,take_log=True) pc.set_width(L, '1') pc.save("%s_test-proj" % pf,override=True) If I set L = 0.2, I don't see in the central region of the projection more structures than with L=0.1 (only the plot width is, of course, larger). Moreover, the result do not depend on the region that I select (either third or fourth line). It seems to me that, in all cases, the projection is actually performed along the whole simulation box, and not on the region that I select. In other words, I can control the plot width, but not the projection depth. In the volume rendering, on the other hand, the plot width sets the depth, too. Is anything wrong or missing in my script? If needed, I can send you a few explanatory plots and the whole sample script for reproducing my issue. Looking forward to your suggestions, 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(a)ita.uni-heidelberg.de URL: http://www.ita.uni-heidelberg.de/~luigi/

4 8