Mailman 3 January 2017 - yt-dev

Issue #1308: profile variance vs standard deviation (yt_analysis/yt)
by Andrew Cunningham Jan. 10, 2017

Jan. 10, 2017

New issue 1308: profile variance vs standard deviation https://bitbucket.org/yt_analysis/yt/issues/1308/profile-variance-vs-standa… Andrew Cunningham: yt/data_objects/profiles.py :: _finalize_storage() sets: all_var = np.sqrt(all_var) which later is used to construct "self.variance". So, after this line, both "all_var" and the public-facing member "variance" is really a standard deviation. The cookbook example "profile_with_variance.py" correctly labels prof.variance as "standard … [View More]

1 0

Issue #1307: Add capability to visualize a single variable on multiple meshes (yt_analysis/yt)
by Alexander Lindsay Jan. 10, 2017

Jan. 10, 2017

New issue 1307: Add capability to visualize a single variable on multiple meshes https://bitbucket.org/yt_analysis/yt/issues/1307/add-capability-to-visualiz… Alexander Lindsay: - Starting with SlicePlot initially, it would be great to be able to visualize a variable that exists on multiple (most often geometrically and physically connected) mesh blocks. - This is probably unique to UnstructuredMesh (at least for now) Currently have a hacky solution...will refine over the next couple days.

1 0

Microns & Simulation-Specific Unit Base
by Axel Huebl Jan. 9, 2017

Jan. 9, 2017

Hi yt-devs! I have a quick user-question (sorry) and a devel question: First, how to convert to a unit like microns (10^-6 meters) or pico-seconds? I can't find an example here [1]. I tried .in_units("mum") .in_units("mu m) .in_units("1.e-6 m") .in_units("microns") .in_units("micro meter") .in_units("micro meters") .in_units("micro m") .in_units("μm") and only .in_units("1.e-6 * m") does not fail but looks ugly. An orthogonal question: in our domain (laser-plasma physics)… [View More]

3 6

Scheduled maintenance
by Kacper Kowalik Jan. 7, 2017

Jan. 7, 2017

Hi all, our cloud provider "will have a brief network outage starting at 10AM (CST/UTC-6) on Friday, January 13th in order to upgrade its network hardware. The upgrade will install new hardware and software and will double the bandwidth. Before the outage the new hardware will already be configured and connected to the network so the interruption should be brief, hopefully about 5 minutes." This affects: hub (only running notebooks), jenkins, blog, pastebin, slackbots (yt-fido, irc bridge). … [View More]

1 0

Issue #1306: Profile computes the standard deviation instead of variance (yt_analysis/yt)
by Greg Bryan Jan. 5, 2017

Jan. 5, 2017

New issue 1306: Profile computes the standard deviation instead of variance https://bitbucket.org/yt_analysis/yt/issues/1306/profile-computes-the-stand… Greg Bryan: Andrew Cunningham points out: In the following cookbook example: [http://yt-project.org/doc/cookbook/simple_plots.html#profiles-with-variance… URL) the variance along the profile is extracted: ``` #!python variance = prof.variance['gas', 'velocity_magnitude'].value ``` But this is actually the standard deviation. You can … [View More]

1 0

Issue #1305: Performance Issues with large BoxLib datasets (yt_analysis/yt)
by Chris Byrohl Jan. 4, 2017

Jan. 4, 2017

New issue 1305: Performance Issues with large BoxLib datasets https://bitbucket.org/yt_analysis/yt/issues/1305/performance-issues-with-la… Chris Byrohl: When loading a larger BoxLib dataset (~270Gb), yt is stuck for hours without making any progress. Interrupting points to _reconstruct_parent_child(self): ``` #!python --------------------------------------------------------------------------- KeyboardInterrupt Traceback (most recent call last) <ipython-input-14-… [View More]fa86699093be> in <module>() ----> 1 box['density'] /home/uni09/cosmo/cbyrohl/anaconda3/envs/py35/lib/python3.5/site-packages/yt/data_objects/data_containers.py in __getitem__(self, key) 263 Returns a single field. Will add if necessary. 264 """ --> 265 f = self._determine_fields([key])[0] 266 if f not in self.field_data and key not in self.field_data: 267 if f in self._container_fields: /home/uni09/cosmo/cbyrohl/anaconda3/envs/py35/lib/python3.5/site-packages/yt/data_objects/data_containers.py in _determine_fields(self, fields) 993 else: 994 fname = field --> 995 finfo = self.ds._get_field_info("unknown", fname) 996 if finfo.particle_type: 997 ftype = self._current_particle_type /home/uni09/cosmo/cbyrohl/anaconda3/envs/py35/lib/python3.5/site-packages/yt/data_objects/static_output.py in _get_field_info(self, ftype, fname) 622 _last_finfo = None 623 def _get_field_info(self, ftype, fname = None): --> 624 self.index 625 if fname is None: 626 if isinstance(ftype, DerivedField): /home/uni09/cosmo/cbyrohl/anaconda3/envs/py35/lib/python3.5/site-packages/yt/data_objects/static_output.py in index(self) 417 raise RuntimeError("You should not instantiate Dataset.") 418 self._instantiated_index = self._index_class( --> 419 self, dataset_type=self.dataset_type) 420 # Now we do things that we need an instantiated index for 421 # ...first off, we create our field_info now. /home/uni09/cosmo/cbyrohl/anaconda3/envs/py35/lib/python3.5/site-packages/yt/frontends/boxlib/data_structures.py in __init__(self, ds, dataset_type) 144 self.directory = ds.output_dir 145 --> 146 GridIndex.__init__(self, ds, dataset_type) 147 self._cache_endianness(self.grids[-1]) 148 /home/uni09/cosmo/cbyrohl/anaconda3/envs/py35/lib/python3.5/site-packages/yt/geometry/geometry_handler.py in __init__(self, ds, dataset_type) 48 49 mylog.debug("Setting up domain geometry.") ---> 50 self._setup_geometry() 51 52 mylog.debug("Initializing data grid data IO") /home/uni09/cosmo/cbyrohl/anaconda3/envs/py35/lib/python3.5/site-packages/yt/geometry/grid_geometry_handler.py in _setup_geometry(self) 52 53 mylog.debug("Constructing grid objects.") ---> 54 self._populate_grid_objects() 55 56 mylog.debug("Re-examining index") /home/uni09/cosmo/cbyrohl/anaconda3/envs/py35/lib/python3.5/site-packages/yt/frontends/boxlib/data_structures.py in _populate_grid_objects(self) 295 mylog.debug("Creating grid objects") 296 self.grids = np.array(self.grids, dtype='object') --> 297 self._reconstruct_parent_child() 298 for i, grid in enumerate(self.grids): 299 if (i % 1e4) == 0: mylog.debug("Prepared % 7i / % 7i grids", i, /home/uni09/cosmo/cbyrohl/anaconda3/envs/py35/lib/python3.5/site-packages/yt/frontends/boxlib/data_structures.py in _reconstruct_parent_child(self) 311 self.grid_levels[i] + 1, 312 self.grid_left_edge, self.grid_right_edge, --> 313 self.grid_levels, mask) 314 ids = np.where(mask.astype("bool")) # where is a tuple 315 grid._children_ids = ids[0] + grid._id_offset ``` The result of ``` #!python np.savez('data.npz', left_edge=self.grid_left_edge, right_edge=self.grid_right_edge, levels=self.grid_levels) ``` for the beginning of that routine can be found here: http://use.yt/upload/87b007b1 [View Less]

1 0

Testing infra is currently down
by Kacper Kowalik Jan. 4, 2017

Jan. 4, 2017

Hi all, we had a hdd failure on a machine hosting yt's jenkins and blog instances. PRs are not being tested at the moment. I don't have an ETA for a replacement yet. I will keep you posted. Cheers, Kacper

5 5