[yt-svn] commit/yt: 7 new changesets
commits-noreply at bitbucket.org
commits-noreply at bitbucket.org
Sat Jul 19 13:50:05 PDT 2014
7 new commits in yt:
https://bitbucket.org/yt_analysis/yt/commits/57e178ad7d1d/
Changeset: 57e178ad7d1d
Branch: yt-3.0
User: ngoldbaum
Date: 2014-07-18 20:31:40
Summary: Removing references to removed analysis modules in the API docs.
Affected #: 1 file
diff -r 894d44f837d4d259550bb7b741091424226375b7 -r 57e178ad7d1d45efc886593325c87ac3a0442bf5 doc/source/reference/api/api.rst
--- a/doc/source/reference/api/api.rst
+++ b/doc/source/reference/api/api.rst
@@ -437,17 +437,7 @@
~yt.analysis_modules.halo_finding.halo_objects.LoadTextHaloes
~yt.analysis_modules.halo_finding.halo_objects.LoadRockstarHalos
-We have several methods that work to create merger trees:
-
-.. autosummary::
- :toctree: generated/
-
- ~yt.analysis_modules.halo_merger_tree.merger_tree.MergerTree
- ~yt.analysis_modules.halo_merger_tree.merger_tree.MergerTreeConnect
- ~yt.analysis_modules.halo_merger_tree.merger_tree.MergerTreeDotOutput
- ~yt.analysis_modules.halo_merger_tree.merger_tree.MergerTreeTextOutput
-
-You can use Halo catalogs generatedl externally as well:
+You can use Halo catalogs generated externally as well:
.. autosummary::
:toctree: generated/
@@ -456,19 +446,6 @@
~yt.analysis_modules.halo_merger_tree.enzofof_merger_tree.EnzoFOFMergerTree
~yt.analysis_modules.halo_merger_tree.enzofof_merger_tree.plot_halo_evolution
-Halo Profiling
-^^^^^^^^^^^^^^
-
-yt provides a comprehensive halo profiler that can filter, center, and analyze
-halos en masse.
-
-.. autosummary::
- :toctree: generated/
-
- ~yt.analysis_modules.halo_profiler.multi_halo_profiler.HaloProfiler
- ~yt.analysis_modules.halo_profiler.multi_halo_profiler.VirialFilter
-
-
Two Point Functions
^^^^^^^^^^^^^^^^^^^
@@ -512,16 +489,6 @@
Extension Types
---------------
-Coordinate Transformations
-^^^^^^^^^^^^^^^^^^^^^^^^^^
-
-
-.. autosummary::
- :toctree: generated/
-
- ~yt.analysis_modules.coordinate_transformation.transforms.arbitrary_regrid
- ~yt.analysis_modules.coordinate_transformation.transforms.spherical_regrid
-
Cosmology, Star Particle Analysis, and Simulated Observations
^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
@@ -576,17 +543,6 @@
~yt.analysis_modules.radmc3d_export.RadMC3DInterface.RadMC3DLayer
~yt.analysis_modules.radmc3d_export.RadMC3DInterface.RadMC3DWriter
-Radial Column Density
-^^^^^^^^^^^^^^^^^^^^^
-
-If you'd like to calculate the column density out to a given point, from a
-specified center, yt can provide that information.
-
-.. autosummary::
- :toctree: generated/
-
- ~yt.analysis_modules.radial_column_density.radial_column_density.RadialColumnDensity
-
Volume Rendering
^^^^^^^^^^^^^^^^
https://bitbucket.org/yt_analysis/yt/commits/436b57a86065/
Changeset: 436b57a86065
Branch: yt-3.0
User: ngoldbaum
Date: 2014-07-18 21:22:59
Summary: Fixes for doc build errors, some style fixes.
Affected #: 6 files
diff -r 894d44f837d4d259550bb7b741091424226375b7 -r 436b57a86065d7352259ad4ef2aa830c821c83a3 doc/source/bootcamp/2)_Data_Inspection.ipynb
--- a/doc/source/bootcamp/2)_Data_Inspection.ipynb
+++ b/doc/source/bootcamp/2)_Data_Inspection.ipynb
@@ -1,7 +1,7 @@
{
"metadata": {
"name": "",
- "signature": "sha256:9d67e9e4ca5ce92dcd0658025dbfbd28be47b47ca8d4531fdac16cc2c2fa038b"
+ "signature": "sha256:a8fe78715c1f3900c37c675d84320fe65f0ba8734abba60fd12e74d957e5d8ee"
},
"nbformat": 3,
"nbformat_minor": 0,
@@ -180,7 +180,7 @@
"cell_type": "code",
"collapsed": false,
"input": [
- "print ds.index.grids[0]"
+ "print ds.index.grids[1]"
],
"language": "python",
"metadata": {},
@@ -190,7 +190,7 @@
"cell_type": "code",
"collapsed": false,
"input": [
- "g = ds.index.grids[0]\n",
+ "g = ds.index.grids[1]\n",
"print g"
],
"language": "python",
diff -r 894d44f837d4d259550bb7b741091424226375b7 -r 436b57a86065d7352259ad4ef2aa830c821c83a3 doc/source/examining/Loading_Generic_Particle_Data.ipynb
--- a/doc/source/examining/Loading_Generic_Particle_Data.ipynb
+++ b/doc/source/examining/Loading_Generic_Particle_Data.ipynb
@@ -1,7 +1,7 @@
{
"metadata": {
"name": "",
- "signature": "sha256:6da8ec00f414307f27544fbdbc6b4fa476e5e96809003426279b2a1c898b4546"
+ "signature": "sha256:66184ef8ba7b430812fef4db39a425ea1b8d07ebb0f7609823951b55262dffbc"
},
"nbformat": 3,
"nbformat_minor": 0,
@@ -31,7 +31,8 @@
],
"language": "python",
"metadata": {},
- "outputs": []
+ "outputs": [],
+ "prompt_number": 1
},
{
"cell_type": "markdown",
@@ -51,7 +52,8 @@
],
"language": "python",
"metadata": {},
- "outputs": []
+ "outputs": [],
+ "prompt_number": 2
},
{
"cell_type": "markdown",
@@ -80,7 +82,8 @@
],
"language": "python",
"metadata": {},
- "outputs": []
+ "outputs": [],
+ "prompt_number": 3
},
{
"cell_type": "markdown",
@@ -116,7 +119,8 @@
],
"language": "python",
"metadata": {},
- "outputs": []
+ "outputs": [],
+ "prompt_number": 4
},
{
"cell_type": "code",
@@ -126,7 +130,24 @@
],
"language": "python",
"metadata": {},
- "outputs": []
+ "outputs": [
+ {
+ "metadata": {},
+ "output_type": "pyout",
+ "prompt_number": 5,
+ "text": [
+ "[('io', 'particle_position_y'),\n",
+ " ('io', 'particle_mass'),\n",
+ " ('io', 'particle_position_x'),\n",
+ " ('io', 'particle_position_z'),\n",
+ " ('all', 'particle_position_z'),\n",
+ " ('all', 'particle_position_x'),\n",
+ " ('all', 'particle_position_y'),\n",
+ " ('all', 'particle_mass')]"
+ ]
+ }
+ ],
+ "prompt_number": 5
},
{
"cell_type": "code",
@@ -136,7 +157,68 @@
],
"language": "python",
"metadata": {},
- "outputs": []
+ "outputs": [
+ {
+ "metadata": {},
+ "output_type": "pyout",
+ "prompt_number": 6,
+ "text": [
+ "[('all', 'mesh_id'),\n",
+ " ('all', 'particle_mass'),\n",
+ " ('all', 'particle_ones'),\n",
+ " ('all', 'particle_phi_spherical'),\n",
+ " ('all', 'particle_position'),\n",
+ " ('all', 'particle_position_x'),\n",
+ " ('all', 'particle_position_y'),\n",
+ " ('all', 'particle_position_z'),\n",
+ " ('all', 'particle_radius'),\n",
+ " ('all', 'particle_radius_spherical'),\n",
+ " ('all', 'particle_theta_spherical'),\n",
+ " ('deposit', 'all_cic'),\n",
+ " ('deposit', 'all_count'),\n",
+ " ('deposit', 'all_density'),\n",
+ " ('deposit', 'all_mass'),\n",
+ " ('deposit', 'io_cic'),\n",
+ " ('deposit', 'io_count'),\n",
+ " ('deposit', 'io_density'),\n",
+ " ('deposit', 'io_mass'),\n",
+ " ('index', 'cell_volume'),\n",
+ " ('index', 'contours'),\n",
+ " ('index', 'cylindrical_r'),\n",
+ " ('index', 'cylindrical_theta'),\n",
+ " ('index', 'cylindrical_z'),\n",
+ " ('index', 'disk_angle'),\n",
+ " ('index', 'dx'),\n",
+ " ('index', 'dy'),\n",
+ " ('index', 'dz'),\n",
+ " ('index', 'grid_indices'),\n",
+ " ('index', 'grid_level'),\n",
+ " ('index', 'height'),\n",
+ " ('index', 'ones'),\n",
+ " ('index', 'ones_over_dx'),\n",
+ " ('index', 'radius'),\n",
+ " ('index', 'spherical_phi'),\n",
+ " ('index', 'spherical_r'),\n",
+ " ('index', 'spherical_theta'),\n",
+ " ('index', 'x'),\n",
+ " ('index', 'y'),\n",
+ " ('index', 'z'),\n",
+ " ('index', 'zeros'),\n",
+ " ('io', 'mesh_id'),\n",
+ " ('io', 'particle_mass'),\n",
+ " ('io', 'particle_ones'),\n",
+ " ('io', 'particle_phi_spherical'),\n",
+ " ('io', 'particle_position'),\n",
+ " ('io', 'particle_position_x'),\n",
+ " ('io', 'particle_position_y'),\n",
+ " ('io', 'particle_position_z'),\n",
+ " ('io', 'particle_radius'),\n",
+ " ('io', 'particle_radius_spherical'),\n",
+ " ('io', 'particle_theta_spherical')]"
+ ]
+ }
+ ],
+ "prompt_number": 6
},
{
"cell_type": "code",
@@ -147,7 +229,38 @@
],
"language": "python",
"metadata": {},
- "outputs": []
+ "outputs": [
+ {
+ "ename": "ValueError",
+ "evalue": "could not broadcast input array from shape (5000000) into shape (15872)",
+ "output_type": "pyerr",
+ "traceback": [
+ "\u001b[0;31m---------------------------------------------------------------------------\u001b[0m\n\u001b[0;31mValueError\u001b[0m Traceback (most recent call last)",
+ "\u001b[0;32m<ipython-input-7-196e90bd4e01>\u001b[0m in \u001b[0;36m<module>\u001b[0;34m()\u001b[0m\n\u001b[0;32m----> 1\u001b[0;31m \u001b[0mslc\u001b[0m \u001b[0;34m=\u001b[0m \u001b[0myt\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0mSlicePlot\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0mds\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0;36m2\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0;34m(\u001b[0m\u001b[0;34m'deposit'\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0;34m'all_cic'\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[0m\u001b[1;32m 2\u001b[0m \u001b[0mslc\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0mset_width\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0;36m8\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0;34m'Mpc'\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n",
+ "\u001b[0;32m/Users/goldbaum/Documents/yt-hg/yt/visualization/plot_window.pyc\u001b[0m in \u001b[0;36mSlicePlot\u001b[0;34m(pf, normal, fields, axis, *args, **kwargs)\u001b[0m\n\u001b[1;32m 1882\u001b[0m \u001b[0;32mdel\u001b[0m \u001b[0mkwargs\u001b[0m\u001b[0;34m[\u001b[0m\u001b[0;34m'north_vector'\u001b[0m\u001b[0;34m]\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[1;32m 1883\u001b[0m \u001b[0;34m\u001b[0m\u001b[0m\n\u001b[0;32m-> 1884\u001b[0;31m \u001b[0;32mreturn\u001b[0m \u001b[0mAxisAlignedSlicePlot\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0mpf\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0mnormal\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0mfields\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0;34m*\u001b[0m\u001b[0margs\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0;34m**\u001b[0m\u001b[0mkwargs\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[0m",
+ "\u001b[0;32m/Users/goldbaum/Documents/yt-hg/yt/visualization/plot_window.pyc\u001b[0m in \u001b[0;36m__init__\u001b[0;34m(self, pf, axis, fields, center, width, axes_unit, origin, fontsize, field_parameters, window_size, aspect)\u001b[0m\n\u001b[1;32m 1089\u001b[0m slc = pf.slice(axis, center[axis],\n\u001b[1;32m 1090\u001b[0m field_parameters = field_parameters, center=center)\n\u001b[0;32m-> 1091\u001b[0;31m \u001b[0mslc\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0mget_data\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0mfields\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[0m\u001b[1;32m 1092\u001b[0m PWViewerMPL.__init__(self, slc, bounds, origin=origin,\n\u001b[1;32m 1093\u001b[0m \u001b[0mfontsize\u001b[0m\u001b[0;34m=\u001b[0m\u001b[0mfontsize\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0mfields\u001b[0m\u001b[0;34m=\u001b[0m\u001b[0mfields\u001b[0m\u001b[0;34m,\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n",
+ "\u001b[0;32m/Users/goldbaum/Documents/yt-hg/yt/data_objects/data_containers.pyc\u001b[0m in \u001b[0;36mget_data\u001b[0;34m(self, fields)\u001b[0m\n\u001b[1;32m 633\u001b[0m \u001b[0;34m\u001b[0m\u001b[0m\n\u001b[1;32m 634\u001b[0m \u001b[0mfields_to_generate\u001b[0m \u001b[0;34m+=\u001b[0m \u001b[0mgen_fluids\u001b[0m \u001b[0;34m+\u001b[0m \u001b[0mgen_particles\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[0;32m--> 635\u001b[0;31m \u001b[0mself\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0m_generate_fields\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0mfields_to_generate\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[0m\u001b[1;32m 636\u001b[0m \u001b[0;32mfor\u001b[0m \u001b[0mfield\u001b[0m \u001b[0;32min\u001b[0m \u001b[0mself\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0mfield_data\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0mkeys\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m:\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[1;32m 637\u001b[0m \u001b[0;32mif\u001b[0m \u001b[0mfield\u001b[0m \u001b[0;32mnot\u001b[0m \u001b[0;32min\u001b[0m \u001b[0mofields\u001b[0m\u001b[0;34m:\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n",
+ "\u001b[0;32m/Users/goldbaum/Documents/yt-hg/yt/data_objects/data_containers.pyc\u001b[0m in \u001b[0;36m_generate_fields\u001b[0;34m(self, fields_to_generate)\u001b[0m\n\u001b[1;32m 653\u001b[0m \u001b[0mfi\u001b[0m \u001b[0;34m=\u001b[0m \u001b[0mself\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0mpf\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0m_get_field_info\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0;34m*\u001b[0m\u001b[0mfield\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[1;32m 654\u001b[0m \u001b[0;32mtry\u001b[0m\u001b[0;34m:\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[0;32m--> 655\u001b[0;31m \u001b[0mfd\u001b[0m \u001b[0;34m=\u001b[0m \u001b[0mself\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0m_generate_field\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0mfield\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[0m\u001b[1;32m 656\u001b[0m \u001b[0;32mif\u001b[0m \u001b[0mtype\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0mfd\u001b[0m\u001b[0;34m)\u001b[0m \u001b[0;34m==\u001b[0m \u001b[0mnp\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0mndarray\u001b[0m\u001b[0;34m:\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[1;32m 657\u001b[0m \u001b[0mfd\u001b[0m \u001b[0;34m=\u001b[0m \u001b[0mself\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0mpf\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0marr\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0mfd\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0mfi\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0munits\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n",
+ "\u001b[0;32m/Users/goldbaum/Documents/yt-hg/yt/data_objects/data_containers.pyc\u001b[0m in \u001b[0;36m_generate_field\u001b[0;34m(self, field)\u001b[0m\n\u001b[1;32m 255\u001b[0m \u001b[0mtr\u001b[0m \u001b[0;34m=\u001b[0m \u001b[0mself\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0m_generate_particle_field\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0mfield\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[1;32m 256\u001b[0m \u001b[0;32melse\u001b[0m\u001b[0;34m:\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[0;32m--> 257\u001b[0;31m \u001b[0mtr\u001b[0m \u001b[0;34m=\u001b[0m \u001b[0mself\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0m_generate_fluid_field\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0mfield\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[0m\u001b[1;32m 258\u001b[0m \u001b[0;32mif\u001b[0m \u001b[0mtr\u001b[0m \u001b[0;32mis\u001b[0m \u001b[0mNone\u001b[0m\u001b[0;34m:\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[1;32m 259\u001b[0m \u001b[0;32mraise\u001b[0m \u001b[0mYTCouldNotGenerateField\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0mfield\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0mself\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0mpf\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n",
+ "\u001b[0;32m/Users/goldbaum/Documents/yt-hg/yt/data_objects/data_containers.pyc\u001b[0m in \u001b[0;36m_generate_fluid_field\u001b[0;34m(self, field)\u001b[0m\n\u001b[1;32m 273\u001b[0m \u001b[0mfinfo\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0mcheck_available\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0mgen_obj\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[1;32m 274\u001b[0m \u001b[0;32mexcept\u001b[0m \u001b[0mNeedsGridType\u001b[0m \u001b[0;32mas\u001b[0m \u001b[0mngt_exception\u001b[0m\u001b[0;34m:\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[0;32m--> 275\u001b[0;31m \u001b[0mrv\u001b[0m \u001b[0;34m=\u001b[0m \u001b[0mself\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0m_generate_spatial_fluid\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0mfield\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0mngt_exception\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0mghost_zones\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[0m\u001b[1;32m 276\u001b[0m \u001b[0;32melse\u001b[0m\u001b[0;34m:\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[1;32m 277\u001b[0m \u001b[0mrv\u001b[0m \u001b[0;34m=\u001b[0m \u001b[0mfinfo\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0mgen_obj\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n",
+ "\u001b[0;32m/Users/goldbaum/Documents/yt-hg/yt/data_objects/data_containers.pyc\u001b[0m in \u001b[0;36m_generate_spatial_fluid\u001b[0;34m(self, field, ngz)\u001b[0m\n\u001b[1;32m 289\u001b[0m \u001b[0mo\u001b[0m \u001b[0;34m=\u001b[0m \u001b[0mself\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0m_current_chunk\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0mobjs\u001b[0m\u001b[0;34m[\u001b[0m\u001b[0;36m0\u001b[0m\u001b[0;34m]\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[1;32m 290\u001b[0m \u001b[0;32mwith\u001b[0m \u001b[0mo\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0m_activate_cache\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m:\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[0;32m--> 291\u001b[0;31m \u001b[0mind\u001b[0m \u001b[0;34m+=\u001b[0m \u001b[0mo\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0mselect\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0mself\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0mselector\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0mself\u001b[0m\u001b[0;34m[\u001b[0m\u001b[0mfield\u001b[0m\u001b[0;34m]\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0mrv\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0mind\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[0m\u001b[1;32m 292\u001b[0m \u001b[0;32melse\u001b[0m\u001b[0;34m:\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[1;32m 293\u001b[0m \u001b[0mchunks\u001b[0m \u001b[0;34m=\u001b[0m \u001b[0mself\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0mindex\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0m_chunk\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0mself\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0;34m\"spatial\"\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0mngz\u001b[0m \u001b[0;34m=\u001b[0m \u001b[0mngz\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n",
+ "\u001b[0;32m/Users/goldbaum/Documents/yt-hg/yt/data_objects/data_containers.pyc\u001b[0m in \u001b[0;36m__getitem__\u001b[0;34m(self, key)\u001b[0m\n\u001b[1;32m 218\u001b[0m \u001b[0;32mreturn\u001b[0m \u001b[0mself\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0mfield_data\u001b[0m\u001b[0;34m[\u001b[0m\u001b[0mf\u001b[0m\u001b[0;34m]\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[1;32m 219\u001b[0m \u001b[0;32melse\u001b[0m\u001b[0;34m:\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[0;32m--> 220\u001b[0;31m \u001b[0mself\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0mget_data\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0mf\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[0m\u001b[1;32m 221\u001b[0m \u001b[0;31m# fi.units is the unit expression string. We depend on the registry\u001b[0m\u001b[0;34m\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[1;32m 222\u001b[0m \u001b[0;31m# hanging off the dataset to define this unit object.\u001b[0m\u001b[0;34m\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n",
+ "\u001b[0;32m/Users/goldbaum/Documents/yt-hg/yt/data_objects/data_containers.pyc\u001b[0m in \u001b[0;36mget_data\u001b[0;34m(self, fields)\u001b[0m\n\u001b[1;32m 633\u001b[0m \u001b[0;34m\u001b[0m\u001b[0m\n\u001b[1;32m 634\u001b[0m \u001b[0mfields_to_generate\u001b[0m \u001b[0;34m+=\u001b[0m \u001b[0mgen_fluids\u001b[0m \u001b[0;34m+\u001b[0m \u001b[0mgen_particles\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[0;32m--> 635\u001b[0;31m \u001b[0mself\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0m_generate_fields\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0mfields_to_generate\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[0m\u001b[1;32m 636\u001b[0m \u001b[0;32mfor\u001b[0m \u001b[0mfield\u001b[0m \u001b[0;32min\u001b[0m \u001b[0mself\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0mfield_data\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0mkeys\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m:\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[1;32m 637\u001b[0m \u001b[0;32mif\u001b[0m \u001b[0mfield\u001b[0m \u001b[0;32mnot\u001b[0m \u001b[0;32min\u001b[0m \u001b[0mofields\u001b[0m\u001b[0;34m:\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n",
+ "\u001b[0;32m/Users/goldbaum/Documents/yt-hg/yt/data_objects/data_containers.pyc\u001b[0m in \u001b[0;36m_generate_fields\u001b[0;34m(self, fields_to_generate)\u001b[0m\n\u001b[1;32m 653\u001b[0m \u001b[0mfi\u001b[0m \u001b[0;34m=\u001b[0m \u001b[0mself\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0mpf\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0m_get_field_info\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0;34m*\u001b[0m\u001b[0mfield\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[1;32m 654\u001b[0m \u001b[0;32mtry\u001b[0m\u001b[0;34m:\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[0;32m--> 655\u001b[0;31m \u001b[0mfd\u001b[0m \u001b[0;34m=\u001b[0m \u001b[0mself\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0m_generate_field\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0mfield\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[0m\u001b[1;32m 656\u001b[0m \u001b[0;32mif\u001b[0m \u001b[0mtype\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0mfd\u001b[0m\u001b[0;34m)\u001b[0m \u001b[0;34m==\u001b[0m \u001b[0mnp\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0mndarray\u001b[0m\u001b[0;34m:\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[1;32m 657\u001b[0m \u001b[0mfd\u001b[0m \u001b[0;34m=\u001b[0m \u001b[0mself\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0mpf\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0marr\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0mfd\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0mfi\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0munits\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n",
+ "\u001b[0;32m/Users/goldbaum/Documents/yt-hg/yt/data_objects/data_containers.pyc\u001b[0m in \u001b[0;36m_generate_field\u001b[0;34m(self, field)\u001b[0m\n\u001b[1;32m 255\u001b[0m \u001b[0mtr\u001b[0m \u001b[0;34m=\u001b[0m \u001b[0mself\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0m_generate_particle_field\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0mfield\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[1;32m 256\u001b[0m \u001b[0;32melse\u001b[0m\u001b[0;34m:\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[0;32m--> 257\u001b[0;31m \u001b[0mtr\u001b[0m \u001b[0;34m=\u001b[0m \u001b[0mself\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0m_generate_fluid_field\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0mfield\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[0m\u001b[1;32m 258\u001b[0m \u001b[0;32mif\u001b[0m \u001b[0mtr\u001b[0m \u001b[0;32mis\u001b[0m \u001b[0mNone\u001b[0m\u001b[0;34m:\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[1;32m 259\u001b[0m \u001b[0;32mraise\u001b[0m \u001b[0mYTCouldNotGenerateField\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0mfield\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0mself\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0mpf\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n",
+ "\u001b[0;32m/Users/goldbaum/Documents/yt-hg/yt/data_objects/data_containers.pyc\u001b[0m in \u001b[0;36m_generate_fluid_field\u001b[0;34m(self, field)\u001b[0m\n\u001b[1;32m 275\u001b[0m \u001b[0mrv\u001b[0m \u001b[0;34m=\u001b[0m \u001b[0mself\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0m_generate_spatial_fluid\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0mfield\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0mngt_exception\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0mghost_zones\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[1;32m 276\u001b[0m \u001b[0;32melse\u001b[0m\u001b[0;34m:\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[0;32m--> 277\u001b[0;31m \u001b[0mrv\u001b[0m \u001b[0;34m=\u001b[0m \u001b[0mfinfo\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0mgen_obj\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[0m\u001b[1;32m 278\u001b[0m \u001b[0;32mreturn\u001b[0m \u001b[0mrv\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[1;32m 279\u001b[0m \u001b[0;34m\u001b[0m\u001b[0m\n",
+ "\u001b[0;32m/Users/goldbaum/Documents/yt-hg/yt/fields/derived_field.pyc\u001b[0m in \u001b[0;36m__call__\u001b[0;34m(self, data)\u001b[0m\n\u001b[1;32m 182\u001b[0m \"for %s\" % (self.name,))\n\u001b[1;32m 183\u001b[0m \u001b[0;32mwith\u001b[0m \u001b[0mself\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0munit_registry\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0mdata\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m:\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[0;32m--> 184\u001b[0;31m \u001b[0mdd\u001b[0m \u001b[0;34m=\u001b[0m \u001b[0mself\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0m_function\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0mself\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0mdata\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[0m\u001b[1;32m 185\u001b[0m \u001b[0;32mfor\u001b[0m \u001b[0mfield_name\u001b[0m \u001b[0;32min\u001b[0m \u001b[0mdata\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0mkeys\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m:\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[1;32m 186\u001b[0m \u001b[0;32mif\u001b[0m \u001b[0mfield_name\u001b[0m \u001b[0;32mnot\u001b[0m \u001b[0;32min\u001b[0m \u001b[0moriginal_fields\u001b[0m\u001b[0;34m:\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n",
+ "\u001b[0;32m/Users/goldbaum/Documents/yt-hg/yt/fields/particle_fields.pyc\u001b[0m in \u001b[0;36mparticle_cic\u001b[0;34m(field, data)\u001b[0m\n\u001b[1;32m 114\u001b[0m \u001b[0;34m\u001b[0m\u001b[0m\n\u001b[1;32m 115\u001b[0m \u001b[0;32mdef\u001b[0m \u001b[0mparticle_cic\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0mfield\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0mdata\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m:\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[0;32m--> 116\u001b[0;31m \u001b[0mpos\u001b[0m \u001b[0;34m=\u001b[0m \u001b[0mdata\u001b[0m\u001b[0;34m[\u001b[0m\u001b[0mptype\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0mcoord_name\u001b[0m\u001b[0;34m]\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[0m\u001b[1;32m 117\u001b[0m \u001b[0md\u001b[0m \u001b[0;34m=\u001b[0m \u001b[0mdata\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0mdeposit\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0mpos\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0;34m[\u001b[0m\u001b[0mdata\u001b[0m\u001b[0;34m[\u001b[0m\u001b[0mptype\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0mmass_name\u001b[0m\u001b[0;34m]\u001b[0m\u001b[0;34m]\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0mmethod\u001b[0m \u001b[0;34m=\u001b[0m \u001b[0;34m\"cic\"\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[1;32m 118\u001b[0m \u001b[0md\u001b[0m \u001b[0;34m=\u001b[0m \u001b[0mdata\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0mapply_units\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0md\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0mdata\u001b[0m\u001b[0;34m[\u001b[0m\u001b[0mptype\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0mmass_name\u001b[0m\u001b[0;34m]\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0munits\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n",
+ "\u001b[0;32m/Users/goldbaum/Documents/yt-hg/yt/data_objects/octree_subset.pyc\u001b[0m in \u001b[0;36m__getitem__\u001b[0;34m(self, key)\u001b[0m\n\u001b[1;32m 74\u001b[0m \u001b[0;34m\u001b[0m\u001b[0m\n\u001b[1;32m 75\u001b[0m \u001b[0;32mdef\u001b[0m \u001b[0m__getitem__\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0mself\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0mkey\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m:\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[0;32m---> 76\u001b[0;31m \u001b[0mtr\u001b[0m \u001b[0;34m=\u001b[0m \u001b[0msuper\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0mOctreeSubset\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0mself\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0m__getitem__\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0mkey\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[0m\u001b[1;32m 77\u001b[0m \u001b[0;32mtry\u001b[0m\u001b[0;34m:\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[1;32m 78\u001b[0m \u001b[0mfields\u001b[0m \u001b[0;34m=\u001b[0m \u001b[0mself\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0m_determine_fields\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0mkey\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n",
+ "\u001b[0;32m/Users/goldbaum/Documents/yt-hg/yt/data_objects/data_containers.pyc\u001b[0m in \u001b[0;36m__getitem__\u001b[0;34m(self, key)\u001b[0m\n\u001b[1;32m 218\u001b[0m \u001b[0;32mreturn\u001b[0m \u001b[0mself\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0mfield_data\u001b[0m\u001b[0;34m[\u001b[0m\u001b[0mf\u001b[0m\u001b[0;34m]\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[1;32m 219\u001b[0m \u001b[0;32melse\u001b[0m\u001b[0;34m:\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[0;32m--> 220\u001b[0;31m \u001b[0mself\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0mget_data\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0mf\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[0m\u001b[1;32m 221\u001b[0m \u001b[0;31m# fi.units is the unit expression string. We depend on the registry\u001b[0m\u001b[0;34m\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[1;32m 222\u001b[0m \u001b[0;31m# hanging off the dataset to define this unit object.\u001b[0m\u001b[0;34m\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n",
+ "\u001b[0;32m/Users/goldbaum/Documents/yt-hg/yt/data_objects/data_containers.pyc\u001b[0m in \u001b[0;36mget_data\u001b[0;34m(self, fields)\u001b[0m\n\u001b[1;32m 627\u001b[0m \u001b[0;34m\u001b[0m\u001b[0m\n\u001b[1;32m 628\u001b[0m read_particles, gen_particles = self.index._read_particle_fields(\n\u001b[0;32m--> 629\u001b[0;31m particles, self, self._current_chunk)\n\u001b[0m\u001b[1;32m 630\u001b[0m \u001b[0;32mfor\u001b[0m \u001b[0mf\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0mv\u001b[0m \u001b[0;32min\u001b[0m \u001b[0mread_particles\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0mitems\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m:\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[1;32m 631\u001b[0m \u001b[0mself\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0mfield_data\u001b[0m\u001b[0;34m[\u001b[0m\u001b[0mf\u001b[0m\u001b[0;34m]\u001b[0m \u001b[0;34m=\u001b[0m \u001b[0mself\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0mpf\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0marr\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0mv\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0minput_units\u001b[0m \u001b[0;34m=\u001b[0m \u001b[0mfinfos\u001b[0m\u001b[0;34m[\u001b[0m\u001b[0mf\u001b[0m\u001b[0;34m]\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0munits\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n",
+ "\u001b[0;32m/Users/goldbaum/Documents/yt-hg/yt/geometry/geometry_handler.pyc\u001b[0m in \u001b[0;36m_read_particle_fields\u001b[0;34m(self, fields, dobj, chunk)\u001b[0m\n\u001b[1;32m 231\u001b[0m \u001b[0mself\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0m_chunk_io\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0mdobj\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0mcache\u001b[0m \u001b[0;34m=\u001b[0m \u001b[0mFalse\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m,\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[1;32m 232\u001b[0m \u001b[0mselector\u001b[0m\u001b[0;34m,\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[0;32m--> 233\u001b[0;31m fields_to_read)\n\u001b[0m\u001b[1;32m 234\u001b[0m \u001b[0;32mfor\u001b[0m \u001b[0mfield\u001b[0m \u001b[0;32min\u001b[0m \u001b[0mfields_to_read\u001b[0m\u001b[0;34m:\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[1;32m 235\u001b[0m \u001b[0mftype\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0mfname\u001b[0m \u001b[0;34m=\u001b[0m \u001b[0mfield\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n",
+ "\u001b[0;32m/Users/goldbaum/Documents/yt-hg/yt/utilities/io_handler.pyc\u001b[0m in \u001b[0;36m_read_particle_selection\u001b[0;34m(self, chunks, selector, fields)\u001b[0m\n\u001b[1;32m 172\u001b[0m \u001b[0;31m#mylog.debug(\"Filling %s from %s to %s with %s\",\u001b[0m\u001b[0;34m\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[1;32m 173\u001b[0m \u001b[0;31m# field_f, my_ind, my_ind+vals.shape[0], field_r)\u001b[0m\u001b[0;34m\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[0;32m--> 174\u001b[0;31m \u001b[0mrv\u001b[0m\u001b[0;34m[\u001b[0m\u001b[0mfield_f\u001b[0m\u001b[0;34m]\u001b[0m\u001b[0;34m[\u001b[0m\u001b[0mmy_ind\u001b[0m\u001b[0;34m:\u001b[0m\u001b[0mmy_ind\u001b[0m \u001b[0;34m+\u001b[0m \u001b[0mvals\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0mshape\u001b[0m\u001b[0;34m[\u001b[0m\u001b[0;36m0\u001b[0m\u001b[0;34m]\u001b[0m\u001b[0;34m,\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0;34m]\u001b[0m \u001b[0;34m=\u001b[0m \u001b[0mvals\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[0m\u001b[1;32m 175\u001b[0m \u001b[0mind\u001b[0m\u001b[0;34m[\u001b[0m\u001b[0mfield_f\u001b[0m\u001b[0;34m]\u001b[0m \u001b[0;34m+=\u001b[0m \u001b[0mvals\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0mshape\u001b[0m\u001b[0;34m[\u001b[0m\u001b[0;36m0\u001b[0m\u001b[0;34m]\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[1;32m 176\u001b[0m \u001b[0;31m# Now we need to truncate all our fields, since we allow for\u001b[0m\u001b[0;34m\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n",
+ "\u001b[0;31mValueError\u001b[0m: could not broadcast input array from shape (5000000) into shape (15872)"
+ ]
+ }
+ ],
+ "prompt_number": 7
}
],
"metadata": {}
diff -r 894d44f837d4d259550bb7b741091424226375b7 -r 436b57a86065d7352259ad4ef2aa830c821c83a3 doc/source/examining/tipsy_and_yt.ipynb
--- a/doc/source/examining/tipsy_and_yt.ipynb
+++ b/doc/source/examining/tipsy_and_yt.ipynb
@@ -1,7 +1,7 @@
{
"metadata": {
"name": "",
- "signature": "sha256:a80c1b224c121c67e57acfa9183c5660a332a37556a492e230476b424827885f"
+ "signature": "sha256:2ae8b1599fa35495fa1bb8deb1c67094e3529e70093b30e20354122cd9403d9d"
},
"nbformat": 3,
"nbformat_minor": 0,
@@ -35,7 +35,7 @@
"cell_type": "code",
"collapsed": false,
"input": [
- "from yt.mods import *"
+ "import yt"
],
"language": "python",
"metadata": {},
@@ -59,7 +59,7 @@
"cell_type": "code",
"collapsed": false,
"input": [
- "ds = load('TipsyGalaxy/galaxy.00300', n_ref=8)"
+ "ds = yt.load('TipsyGalaxy/galaxy.00300', n_ref=8)"
],
"language": "python",
"metadata": {},
@@ -111,7 +111,8 @@
"collapsed": false,
"input": [
"%matplotlib inline\n",
- "import matplotlib.pyplot as plt"
+ "import matplotlib.pyplot as plt\n",
+ "import numpy as np"
],
"language": "python",
"metadata": {},
@@ -155,7 +156,7 @@
"cell_type": "code",
"collapsed": false,
"input": [
- "SlicePlot(ds, 'z', ('gas','density'), width=(40, 'kpc'), center='m')"
+ "yt.SlicePlot(ds, 'z', ('gas','density'), width=(40, 'kpc'), center='m')"
],
"language": "python",
"metadata": {},
@@ -165,7 +166,7 @@
"cell_type": "code",
"collapsed": false,
"input": [
- "ProjectionPlot(ds, 'z', ('gas','density'), width=(40, 'kpc'), center='m')"
+ "yt.ProjectionPlot(ds, 'z', ('gas','density'), width=(40, 'kpc'), center='m')"
],
"language": "python",
"metadata": {},
@@ -182,7 +183,7 @@
"cell_type": "code",
"collapsed": false,
"input": [
- "SlicePlot(ds, 'z', ('gas', 'FeMassFrac'), width=(40, 'kpc'), center='m')"
+ "yt.SlicePlot(ds, 'z', ('gas', 'FeMassFrac'), width=(40, 'kpc'), center='m')"
],
"language": "python",
"metadata": {},
diff -r 894d44f837d4d259550bb7b741091424226375b7 -r 436b57a86065d7352259ad4ef2aa830c821c83a3 doc/source/visualizing/_cb_docstrings.inc
--- a/doc/source/visualizing/_cb_docstrings.inc
+++ b/doc/source/visualizing/_cb_docstrings.inc
@@ -441,5 +441,5 @@
import yt
ds = yt.load("IsolatedGalaxy/galaxy0030/galaxy0030")
p = yt.SlicePlot(ds, 'z', 'density', center='c', width=(20, 'kpc'))
- p.annotate_timestamp()
+ p.annotate_timestamp(-9, -9)
p.save()
diff -r 894d44f837d4d259550bb7b741091424226375b7 -r 436b57a86065d7352259ad4ef2aa830c821c83a3 doc/source/visualizing/plots.rst
--- a/doc/source/visualizing/plots.rst
+++ b/doc/source/visualizing/plots.rst
@@ -847,7 +847,7 @@
import yt
ds = yt.load("IsolatedGalaxy/galaxy0030/galaxy0030")
- p = ProjectionPlot(ds, "x", "density", center='m', width=(10,'kpc'),
+ p = yt.ProjectionPlot(ds, "x", "density", center='m', width=(10,'kpc'),
weight_field='density')
p.show()
diff -r 894d44f837d4d259550bb7b741091424226375b7 -r 436b57a86065d7352259ad4ef2aa830c821c83a3 yt/data_objects/profiles.py
--- a/yt/data_objects/profiles.py
+++ b/yt/data_objects/profiles.py
@@ -1339,6 +1339,8 @@
except KeyError:
field_ex = list(extrema[bin_field])
if units is not None and bin_field in units:
+ if isinstance(field_ex[0], tuple):
+ field_ex = [data_source.pf.quan(*f) for f in field_ex]
fe = data_source.pf.arr(field_ex, units[bin_field])
fe.convert_to_units(bf_units)
field_ex = [fe[0].v, fe[1].v]
https://bitbucket.org/yt_analysis/yt/commits/3180321c28a6/
Changeset: 3180321c28a6
Branch: yt-3.0
User: ngoldbaum
Date: 2014-07-18 21:23:40
Summary: Merging with other docs PR.
Affected #: 1 file
diff -r 436b57a86065d7352259ad4ef2aa830c821c83a3 -r 3180321c28a61079b4e32676d02bb233046a4023 doc/source/reference/api/api.rst
--- a/doc/source/reference/api/api.rst
+++ b/doc/source/reference/api/api.rst
@@ -437,17 +437,7 @@
~yt.analysis_modules.halo_finding.halo_objects.LoadTextHaloes
~yt.analysis_modules.halo_finding.halo_objects.LoadRockstarHalos
-We have several methods that work to create merger trees:
-
-.. autosummary::
- :toctree: generated/
-
- ~yt.analysis_modules.halo_merger_tree.merger_tree.MergerTree
- ~yt.analysis_modules.halo_merger_tree.merger_tree.MergerTreeConnect
- ~yt.analysis_modules.halo_merger_tree.merger_tree.MergerTreeDotOutput
- ~yt.analysis_modules.halo_merger_tree.merger_tree.MergerTreeTextOutput
-
-You can use Halo catalogs generatedl externally as well:
+You can use Halo catalogs generated externally as well:
.. autosummary::
:toctree: generated/
@@ -456,19 +446,6 @@
~yt.analysis_modules.halo_merger_tree.enzofof_merger_tree.EnzoFOFMergerTree
~yt.analysis_modules.halo_merger_tree.enzofof_merger_tree.plot_halo_evolution
-Halo Profiling
-^^^^^^^^^^^^^^
-
-yt provides a comprehensive halo profiler that can filter, center, and analyze
-halos en masse.
-
-.. autosummary::
- :toctree: generated/
-
- ~yt.analysis_modules.halo_profiler.multi_halo_profiler.HaloProfiler
- ~yt.analysis_modules.halo_profiler.multi_halo_profiler.VirialFilter
-
-
Two Point Functions
^^^^^^^^^^^^^^^^^^^
@@ -512,16 +489,6 @@
Extension Types
---------------
-Coordinate Transformations
-^^^^^^^^^^^^^^^^^^^^^^^^^^
-
-
-.. autosummary::
- :toctree: generated/
-
- ~yt.analysis_modules.coordinate_transformation.transforms.arbitrary_regrid
- ~yt.analysis_modules.coordinate_transformation.transforms.spherical_regrid
-
Cosmology, Star Particle Analysis, and Simulated Observations
^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
@@ -576,17 +543,6 @@
~yt.analysis_modules.radmc3d_export.RadMC3DInterface.RadMC3DLayer
~yt.analysis_modules.radmc3d_export.RadMC3DInterface.RadMC3DWriter
-Radial Column Density
-^^^^^^^^^^^^^^^^^^^^^
-
-If you'd like to calculate the column density out to a given point, from a
-specified center, yt can provide that information.
-
-.. autosummary::
- :toctree: generated/
-
- ~yt.analysis_modules.radial_column_density.radial_column_density.RadialColumnDensity
-
Volume Rendering
^^^^^^^^^^^^^^^^
https://bitbucket.org/yt_analysis/yt/commits/2e1effc71929/
Changeset: 2e1effc71929
Branch: yt-3.0
User: ngoldbaum
Date: 2014-07-18 21:27:00
Summary: Removing output that should not have been saved.
Affected #: 1 file
diff -r 3180321c28a61079b4e32676d02bb233046a4023 -r 2e1effc719298feb546edf3222fabb0e3d1e8ed5 doc/source/examining/Loading_Generic_Particle_Data.ipynb
--- a/doc/source/examining/Loading_Generic_Particle_Data.ipynb
+++ b/doc/source/examining/Loading_Generic_Particle_Data.ipynb
@@ -1,7 +1,7 @@
{
"metadata": {
"name": "",
- "signature": "sha256:66184ef8ba7b430812fef4db39a425ea1b8d07ebb0f7609823951b55262dffbc"
+ "signature": "sha256:6da8ec00f414307f27544fbdbc6b4fa476e5e96809003426279b2a1c898b4546"
},
"nbformat": 3,
"nbformat_minor": 0,
@@ -31,8 +31,7 @@
],
"language": "python",
"metadata": {},
- "outputs": [],
- "prompt_number": 1
+ "outputs": []
},
{
"cell_type": "markdown",
@@ -52,8 +51,7 @@
],
"language": "python",
"metadata": {},
- "outputs": [],
- "prompt_number": 2
+ "outputs": []
},
{
"cell_type": "markdown",
@@ -82,8 +80,7 @@
],
"language": "python",
"metadata": {},
- "outputs": [],
- "prompt_number": 3
+ "outputs": []
},
{
"cell_type": "markdown",
@@ -119,8 +116,7 @@
],
"language": "python",
"metadata": {},
- "outputs": [],
- "prompt_number": 4
+ "outputs": []
},
{
"cell_type": "code",
@@ -130,24 +126,7 @@
],
"language": "python",
"metadata": {},
- "outputs": [
- {
- "metadata": {},
- "output_type": "pyout",
- "prompt_number": 5,
- "text": [
- "[('io', 'particle_position_y'),\n",
- " ('io', 'particle_mass'),\n",
- " ('io', 'particle_position_x'),\n",
- " ('io', 'particle_position_z'),\n",
- " ('all', 'particle_position_z'),\n",
- " ('all', 'particle_position_x'),\n",
- " ('all', 'particle_position_y'),\n",
- " ('all', 'particle_mass')]"
- ]
- }
- ],
- "prompt_number": 5
+ "outputs": []
},
{
"cell_type": "code",
@@ -157,68 +136,7 @@
],
"language": "python",
"metadata": {},
- "outputs": [
- {
- "metadata": {},
- "output_type": "pyout",
- "prompt_number": 6,
- "text": [
- "[('all', 'mesh_id'),\n",
- " ('all', 'particle_mass'),\n",
- " ('all', 'particle_ones'),\n",
- " ('all', 'particle_phi_spherical'),\n",
- " ('all', 'particle_position'),\n",
- " ('all', 'particle_position_x'),\n",
- " ('all', 'particle_position_y'),\n",
- " ('all', 'particle_position_z'),\n",
- " ('all', 'particle_radius'),\n",
- " ('all', 'particle_radius_spherical'),\n",
- " ('all', 'particle_theta_spherical'),\n",
- " ('deposit', 'all_cic'),\n",
- " ('deposit', 'all_count'),\n",
- " ('deposit', 'all_density'),\n",
- " ('deposit', 'all_mass'),\n",
- " ('deposit', 'io_cic'),\n",
- " ('deposit', 'io_count'),\n",
- " ('deposit', 'io_density'),\n",
- " ('deposit', 'io_mass'),\n",
- " ('index', 'cell_volume'),\n",
- " ('index', 'contours'),\n",
- " ('index', 'cylindrical_r'),\n",
- " ('index', 'cylindrical_theta'),\n",
- " ('index', 'cylindrical_z'),\n",
- " ('index', 'disk_angle'),\n",
- " ('index', 'dx'),\n",
- " ('index', 'dy'),\n",
- " ('index', 'dz'),\n",
- " ('index', 'grid_indices'),\n",
- " ('index', 'grid_level'),\n",
- " ('index', 'height'),\n",
- " ('index', 'ones'),\n",
- " ('index', 'ones_over_dx'),\n",
- " ('index', 'radius'),\n",
- " ('index', 'spherical_phi'),\n",
- " ('index', 'spherical_r'),\n",
- " ('index', 'spherical_theta'),\n",
- " ('index', 'x'),\n",
- " ('index', 'y'),\n",
- " ('index', 'z'),\n",
- " ('index', 'zeros'),\n",
- " ('io', 'mesh_id'),\n",
- " ('io', 'particle_mass'),\n",
- " ('io', 'particle_ones'),\n",
- " ('io', 'particle_phi_spherical'),\n",
- " ('io', 'particle_position'),\n",
- " ('io', 'particle_position_x'),\n",
- " ('io', 'particle_position_y'),\n",
- " ('io', 'particle_position_z'),\n",
- " ('io', 'particle_radius'),\n",
- " ('io', 'particle_radius_spherical'),\n",
- " ('io', 'particle_theta_spherical')]"
- ]
- }
- ],
- "prompt_number": 6
+ "outputs": []
},
{
"cell_type": "code",
@@ -229,38 +147,7 @@
],
"language": "python",
"metadata": {},
- "outputs": [
- {
- "ename": "ValueError",
- "evalue": "could not broadcast input array from shape (5000000) into shape (15872)",
- "output_type": "pyerr",
- "traceback": [
- "\u001b[0;31m---------------------------------------------------------------------------\u001b[0m\n\u001b[0;31mValueError\u001b[0m Traceback (most recent call last)",
- "\u001b[0;32m<ipython-input-7-196e90bd4e01>\u001b[0m in \u001b[0;36m<module>\u001b[0;34m()\u001b[0m\n\u001b[0;32m----> 1\u001b[0;31m \u001b[0mslc\u001b[0m \u001b[0;34m=\u001b[0m \u001b[0myt\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0mSlicePlot\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0mds\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0;36m2\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0;34m(\u001b[0m\u001b[0;34m'deposit'\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0;34m'all_cic'\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[0m\u001b[1;32m 2\u001b[0m \u001b[0mslc\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0mset_width\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0;36m8\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0;34m'Mpc'\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n",
- "\u001b[0;32m/Users/goldbaum/Documents/yt-hg/yt/visualization/plot_window.pyc\u001b[0m in \u001b[0;36mSlicePlot\u001b[0;34m(pf, normal, fields, axis, *args, **kwargs)\u001b[0m\n\u001b[1;32m 1882\u001b[0m \u001b[0;32mdel\u001b[0m \u001b[0mkwargs\u001b[0m\u001b[0;34m[\u001b[0m\u001b[0;34m'north_vector'\u001b[0m\u001b[0;34m]\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[1;32m 1883\u001b[0m \u001b[0;34m\u001b[0m\u001b[0m\n\u001b[0;32m-> 1884\u001b[0;31m \u001b[0;32mreturn\u001b[0m \u001b[0mAxisAlignedSlicePlot\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0mpf\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0mnormal\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0mfields\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0;34m*\u001b[0m\u001b[0margs\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0;34m**\u001b[0m\u001b[0mkwargs\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[0m",
- "\u001b[0;32m/Users/goldbaum/Documents/yt-hg/yt/visualization/plot_window.pyc\u001b[0m in \u001b[0;36m__init__\u001b[0;34m(self, pf, axis, fields, center, width, axes_unit, origin, fontsize, field_parameters, window_size, aspect)\u001b[0m\n\u001b[1;32m 1089\u001b[0m slc = pf.slice(axis, center[axis],\n\u001b[1;32m 1090\u001b[0m field_parameters = field_parameters, center=center)\n\u001b[0;32m-> 1091\u001b[0;31m \u001b[0mslc\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0mget_data\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0mfields\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[0m\u001b[1;32m 1092\u001b[0m PWViewerMPL.__init__(self, slc, bounds, origin=origin,\n\u001b[1;32m 1093\u001b[0m \u001b[0mfontsize\u001b[0m\u001b[0;34m=\u001b[0m\u001b[0mfontsize\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0mfields\u001b[0m\u001b[0;34m=\u001b[0m\u001b[0mfields\u001b[0m\u001b[0;34m,\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n",
- "\u001b[0;32m/Users/goldbaum/Documents/yt-hg/yt/data_objects/data_containers.pyc\u001b[0m in \u001b[0;36mget_data\u001b[0;34m(self, fields)\u001b[0m\n\u001b[1;32m 633\u001b[0m \u001b[0;34m\u001b[0m\u001b[0m\n\u001b[1;32m 634\u001b[0m \u001b[0mfields_to_generate\u001b[0m \u001b[0;34m+=\u001b[0m \u001b[0mgen_fluids\u001b[0m \u001b[0;34m+\u001b[0m \u001b[0mgen_particles\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[0;32m--> 635\u001b[0;31m \u001b[0mself\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0m_generate_fields\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0mfields_to_generate\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[0m\u001b[1;32m 636\u001b[0m \u001b[0;32mfor\u001b[0m \u001b[0mfield\u001b[0m \u001b[0;32min\u001b[0m \u001b[0mself\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0mfield_data\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0mkeys\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m:\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[1;32m 637\u001b[0m \u001b[0;32mif\u001b[0m \u001b[0mfield\u001b[0m \u001b[0;32mnot\u001b[0m \u001b[0;32min\u001b[0m \u001b[0mofields\u001b[0m\u001b[0;34m:\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n",
- "\u001b[0;32m/Users/goldbaum/Documents/yt-hg/yt/data_objects/data_containers.pyc\u001b[0m in \u001b[0;36m_generate_fields\u001b[0;34m(self, fields_to_generate)\u001b[0m\n\u001b[1;32m 653\u001b[0m \u001b[0mfi\u001b[0m \u001b[0;34m=\u001b[0m \u001b[0mself\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0mpf\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0m_get_field_info\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0;34m*\u001b[0m\u001b[0mfield\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[1;32m 654\u001b[0m \u001b[0;32mtry\u001b[0m\u001b[0;34m:\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[0;32m--> 655\u001b[0;31m \u001b[0mfd\u001b[0m \u001b[0;34m=\u001b[0m \u001b[0mself\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0m_generate_field\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0mfield\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[0m\u001b[1;32m 656\u001b[0m \u001b[0;32mif\u001b[0m \u001b[0mtype\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0mfd\u001b[0m\u001b[0;34m)\u001b[0m \u001b[0;34m==\u001b[0m \u001b[0mnp\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0mndarray\u001b[0m\u001b[0;34m:\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[1;32m 657\u001b[0m \u001b[0mfd\u001b[0m \u001b[0;34m=\u001b[0m \u001b[0mself\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0mpf\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0marr\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0mfd\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0mfi\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0munits\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n",
- "\u001b[0;32m/Users/goldbaum/Documents/yt-hg/yt/data_objects/data_containers.pyc\u001b[0m in \u001b[0;36m_generate_field\u001b[0;34m(self, field)\u001b[0m\n\u001b[1;32m 255\u001b[0m \u001b[0mtr\u001b[0m \u001b[0;34m=\u001b[0m \u001b[0mself\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0m_generate_particle_field\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0mfield\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[1;32m 256\u001b[0m \u001b[0;32melse\u001b[0m\u001b[0;34m:\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[0;32m--> 257\u001b[0;31m \u001b[0mtr\u001b[0m \u001b[0;34m=\u001b[0m \u001b[0mself\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0m_generate_fluid_field\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0mfield\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[0m\u001b[1;32m 258\u001b[0m \u001b[0;32mif\u001b[0m \u001b[0mtr\u001b[0m \u001b[0;32mis\u001b[0m \u001b[0mNone\u001b[0m\u001b[0;34m:\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[1;32m 259\u001b[0m \u001b[0;32mraise\u001b[0m \u001b[0mYTCouldNotGenerateField\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0mfield\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0mself\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0mpf\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n",
- "\u001b[0;32m/Users/goldbaum/Documents/yt-hg/yt/data_objects/data_containers.pyc\u001b[0m in \u001b[0;36m_generate_fluid_field\u001b[0;34m(self, field)\u001b[0m\n\u001b[1;32m 273\u001b[0m \u001b[0mfinfo\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0mcheck_available\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0mgen_obj\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[1;32m 274\u001b[0m \u001b[0;32mexcept\u001b[0m \u001b[0mNeedsGridType\u001b[0m \u001b[0;32mas\u001b[0m \u001b[0mngt_exception\u001b[0m\u001b[0;34m:\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[0;32m--> 275\u001b[0;31m \u001b[0mrv\u001b[0m \u001b[0;34m=\u001b[0m \u001b[0mself\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0m_generate_spatial_fluid\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0mfield\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0mngt_exception\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0mghost_zones\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[0m\u001b[1;32m 276\u001b[0m \u001b[0;32melse\u001b[0m\u001b[0;34m:\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[1;32m 277\u001b[0m \u001b[0mrv\u001b[0m \u001b[0;34m=\u001b[0m \u001b[0mfinfo\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0mgen_obj\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n",
- "\u001b[0;32m/Users/goldbaum/Documents/yt-hg/yt/data_objects/data_containers.pyc\u001b[0m in \u001b[0;36m_generate_spatial_fluid\u001b[0;34m(self, field, ngz)\u001b[0m\n\u001b[1;32m 289\u001b[0m \u001b[0mo\u001b[0m \u001b[0;34m=\u001b[0m \u001b[0mself\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0m_current_chunk\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0mobjs\u001b[0m\u001b[0;34m[\u001b[0m\u001b[0;36m0\u001b[0m\u001b[0;34m]\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[1;32m 290\u001b[0m \u001b[0;32mwith\u001b[0m \u001b[0mo\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0m_activate_cache\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m:\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[0;32m--> 291\u001b[0;31m \u001b[0mind\u001b[0m \u001b[0;34m+=\u001b[0m \u001b[0mo\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0mselect\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0mself\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0mselector\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0mself\u001b[0m\u001b[0;34m[\u001b[0m\u001b[0mfield\u001b[0m\u001b[0;34m]\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0mrv\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0mind\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[0m\u001b[1;32m 292\u001b[0m \u001b[0;32melse\u001b[0m\u001b[0;34m:\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[1;32m 293\u001b[0m \u001b[0mchunks\u001b[0m \u001b[0;34m=\u001b[0m \u001b[0mself\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0mindex\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0m_chunk\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0mself\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0;34m\"spatial\"\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0mngz\u001b[0m \u001b[0;34m=\u001b[0m \u001b[0mngz\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n",
- "\u001b[0;32m/Users/goldbaum/Documents/yt-hg/yt/data_objects/data_containers.pyc\u001b[0m in \u001b[0;36m__getitem__\u001b[0;34m(self, key)\u001b[0m\n\u001b[1;32m 218\u001b[0m \u001b[0;32mreturn\u001b[0m \u001b[0mself\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0mfield_data\u001b[0m\u001b[0;34m[\u001b[0m\u001b[0mf\u001b[0m\u001b[0;34m]\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[1;32m 219\u001b[0m \u001b[0;32melse\u001b[0m\u001b[0;34m:\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[0;32m--> 220\u001b[0;31m \u001b[0mself\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0mget_data\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0mf\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[0m\u001b[1;32m 221\u001b[0m \u001b[0;31m# fi.units is the unit expression string. We depend on the registry\u001b[0m\u001b[0;34m\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[1;32m 222\u001b[0m \u001b[0;31m# hanging off the dataset to define this unit object.\u001b[0m\u001b[0;34m\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n",
- "\u001b[0;32m/Users/goldbaum/Documents/yt-hg/yt/data_objects/data_containers.pyc\u001b[0m in \u001b[0;36mget_data\u001b[0;34m(self, fields)\u001b[0m\n\u001b[1;32m 633\u001b[0m \u001b[0;34m\u001b[0m\u001b[0m\n\u001b[1;32m 634\u001b[0m \u001b[0mfields_to_generate\u001b[0m \u001b[0;34m+=\u001b[0m \u001b[0mgen_fluids\u001b[0m \u001b[0;34m+\u001b[0m \u001b[0mgen_particles\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[0;32m--> 635\u001b[0;31m \u001b[0mself\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0m_generate_fields\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0mfields_to_generate\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[0m\u001b[1;32m 636\u001b[0m \u001b[0;32mfor\u001b[0m \u001b[0mfield\u001b[0m \u001b[0;32min\u001b[0m \u001b[0mself\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0mfield_data\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0mkeys\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m:\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[1;32m 637\u001b[0m \u001b[0;32mif\u001b[0m \u001b[0mfield\u001b[0m \u001b[0;32mnot\u001b[0m \u001b[0;32min\u001b[0m \u001b[0mofields\u001b[0m\u001b[0;34m:\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n",
- "\u001b[0;32m/Users/goldbaum/Documents/yt-hg/yt/data_objects/data_containers.pyc\u001b[0m in \u001b[0;36m_generate_fields\u001b[0;34m(self, fields_to_generate)\u001b[0m\n\u001b[1;32m 653\u001b[0m \u001b[0mfi\u001b[0m \u001b[0;34m=\u001b[0m \u001b[0mself\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0mpf\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0m_get_field_info\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0;34m*\u001b[0m\u001b[0mfield\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[1;32m 654\u001b[0m \u001b[0;32mtry\u001b[0m\u001b[0;34m:\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[0;32m--> 655\u001b[0;31m \u001b[0mfd\u001b[0m \u001b[0;34m=\u001b[0m \u001b[0mself\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0m_generate_field\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0mfield\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[0m\u001b[1;32m 656\u001b[0m \u001b[0;32mif\u001b[0m \u001b[0mtype\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0mfd\u001b[0m\u001b[0;34m)\u001b[0m \u001b[0;34m==\u001b[0m \u001b[0mnp\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0mndarray\u001b[0m\u001b[0;34m:\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[1;32m 657\u001b[0m \u001b[0mfd\u001b[0m \u001b[0;34m=\u001b[0m \u001b[0mself\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0mpf\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0marr\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0mfd\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0mfi\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0munits\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n",
- "\u001b[0;32m/Users/goldbaum/Documents/yt-hg/yt/data_objects/data_containers.pyc\u001b[0m in \u001b[0;36m_generate_field\u001b[0;34m(self, field)\u001b[0m\n\u001b[1;32m 255\u001b[0m \u001b[0mtr\u001b[0m \u001b[0;34m=\u001b[0m \u001b[0mself\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0m_generate_particle_field\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0mfield\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[1;32m 256\u001b[0m \u001b[0;32melse\u001b[0m\u001b[0;34m:\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[0;32m--> 257\u001b[0;31m \u001b[0mtr\u001b[0m \u001b[0;34m=\u001b[0m \u001b[0mself\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0m_generate_fluid_field\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0mfield\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[0m\u001b[1;32m 258\u001b[0m \u001b[0;32mif\u001b[0m \u001b[0mtr\u001b[0m \u001b[0;32mis\u001b[0m \u001b[0mNone\u001b[0m\u001b[0;34m:\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[1;32m 259\u001b[0m \u001b[0;32mraise\u001b[0m \u001b[0mYTCouldNotGenerateField\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0mfield\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0mself\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0mpf\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n",
- "\u001b[0;32m/Users/goldbaum/Documents/yt-hg/yt/data_objects/data_containers.pyc\u001b[0m in \u001b[0;36m_generate_fluid_field\u001b[0;34m(self, field)\u001b[0m\n\u001b[1;32m 275\u001b[0m \u001b[0mrv\u001b[0m \u001b[0;34m=\u001b[0m \u001b[0mself\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0m_generate_spatial_fluid\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0mfield\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0mngt_exception\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0mghost_zones\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[1;32m 276\u001b[0m \u001b[0;32melse\u001b[0m\u001b[0;34m:\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[0;32m--> 277\u001b[0;31m \u001b[0mrv\u001b[0m \u001b[0;34m=\u001b[0m \u001b[0mfinfo\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0mgen_obj\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[0m\u001b[1;32m 278\u001b[0m \u001b[0;32mreturn\u001b[0m \u001b[0mrv\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[1;32m 279\u001b[0m \u001b[0;34m\u001b[0m\u001b[0m\n",
- "\u001b[0;32m/Users/goldbaum/Documents/yt-hg/yt/fields/derived_field.pyc\u001b[0m in \u001b[0;36m__call__\u001b[0;34m(self, data)\u001b[0m\n\u001b[1;32m 182\u001b[0m \"for %s\" % (self.name,))\n\u001b[1;32m 183\u001b[0m \u001b[0;32mwith\u001b[0m \u001b[0mself\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0munit_registry\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0mdata\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m:\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[0;32m--> 184\u001b[0;31m \u001b[0mdd\u001b[0m \u001b[0;34m=\u001b[0m \u001b[0mself\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0m_function\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0mself\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0mdata\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[0m\u001b[1;32m 185\u001b[0m \u001b[0;32mfor\u001b[0m \u001b[0mfield_name\u001b[0m \u001b[0;32min\u001b[0m \u001b[0mdata\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0mkeys\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m:\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[1;32m 186\u001b[0m \u001b[0;32mif\u001b[0m \u001b[0mfield_name\u001b[0m \u001b[0;32mnot\u001b[0m \u001b[0;32min\u001b[0m \u001b[0moriginal_fields\u001b[0m\u001b[0;34m:\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n",
- "\u001b[0;32m/Users/goldbaum/Documents/yt-hg/yt/fields/particle_fields.pyc\u001b[0m in \u001b[0;36mparticle_cic\u001b[0;34m(field, data)\u001b[0m\n\u001b[1;32m 114\u001b[0m \u001b[0;34m\u001b[0m\u001b[0m\n\u001b[1;32m 115\u001b[0m \u001b[0;32mdef\u001b[0m \u001b[0mparticle_cic\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0mfield\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0mdata\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m:\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[0;32m--> 116\u001b[0;31m \u001b[0mpos\u001b[0m \u001b[0;34m=\u001b[0m \u001b[0mdata\u001b[0m\u001b[0;34m[\u001b[0m\u001b[0mptype\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0mcoord_name\u001b[0m\u001b[0;34m]\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[0m\u001b[1;32m 117\u001b[0m \u001b[0md\u001b[0m \u001b[0;34m=\u001b[0m \u001b[0mdata\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0mdeposit\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0mpos\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0;34m[\u001b[0m\u001b[0mdata\u001b[0m\u001b[0;34m[\u001b[0m\u001b[0mptype\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0mmass_name\u001b[0m\u001b[0;34m]\u001b[0m\u001b[0;34m]\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0mmethod\u001b[0m \u001b[0;34m=\u001b[0m \u001b[0;34m\"cic\"\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[1;32m 118\u001b[0m \u001b[0md\u001b[0m \u001b[0;34m=\u001b[0m \u001b[0mdata\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0mapply_units\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0md\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0mdata\u001b[0m\u001b[0;34m[\u001b[0m\u001b[0mptype\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0mmass_name\u001b[0m\u001b[0;34m]\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0munits\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n",
- "\u001b[0;32m/Users/goldbaum/Documents/yt-hg/yt/data_objects/octree_subset.pyc\u001b[0m in \u001b[0;36m__getitem__\u001b[0;34m(self, key)\u001b[0m\n\u001b[1;32m 74\u001b[0m \u001b[0;34m\u001b[0m\u001b[0m\n\u001b[1;32m 75\u001b[0m \u001b[0;32mdef\u001b[0m \u001b[0m__getitem__\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0mself\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0mkey\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m:\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[0;32m---> 76\u001b[0;31m \u001b[0mtr\u001b[0m \u001b[0;34m=\u001b[0m \u001b[0msuper\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0mOctreeSubset\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0mself\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0m__getitem__\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0mkey\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[0m\u001b[1;32m 77\u001b[0m \u001b[0;32mtry\u001b[0m\u001b[0;34m:\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[1;32m 78\u001b[0m \u001b[0mfields\u001b[0m \u001b[0;34m=\u001b[0m \u001b[0mself\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0m_determine_fields\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0mkey\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n",
- "\u001b[0;32m/Users/goldbaum/Documents/yt-hg/yt/data_objects/data_containers.pyc\u001b[0m in \u001b[0;36m__getitem__\u001b[0;34m(self, key)\u001b[0m\n\u001b[1;32m 218\u001b[0m \u001b[0;32mreturn\u001b[0m \u001b[0mself\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0mfield_data\u001b[0m\u001b[0;34m[\u001b[0m\u001b[0mf\u001b[0m\u001b[0;34m]\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[1;32m 219\u001b[0m \u001b[0;32melse\u001b[0m\u001b[0;34m:\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[0;32m--> 220\u001b[0;31m \u001b[0mself\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0mget_data\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0mf\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[0m\u001b[1;32m 221\u001b[0m \u001b[0;31m# fi.units is the unit expression string. We depend on the registry\u001b[0m\u001b[0;34m\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[1;32m 222\u001b[0m \u001b[0;31m# hanging off the dataset to define this unit object.\u001b[0m\u001b[0;34m\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n",
- "\u001b[0;32m/Users/goldbaum/Documents/yt-hg/yt/data_objects/data_containers.pyc\u001b[0m in \u001b[0;36mget_data\u001b[0;34m(self, fields)\u001b[0m\n\u001b[1;32m 627\u001b[0m \u001b[0;34m\u001b[0m\u001b[0m\n\u001b[1;32m 628\u001b[0m read_particles, gen_particles = self.index._read_particle_fields(\n\u001b[0;32m--> 629\u001b[0;31m particles, self, self._current_chunk)\n\u001b[0m\u001b[1;32m 630\u001b[0m \u001b[0;32mfor\u001b[0m \u001b[0mf\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0mv\u001b[0m \u001b[0;32min\u001b[0m \u001b[0mread_particles\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0mitems\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m:\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[1;32m 631\u001b[0m \u001b[0mself\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0mfield_data\u001b[0m\u001b[0;34m[\u001b[0m\u001b[0mf\u001b[0m\u001b[0;34m]\u001b[0m \u001b[0;34m=\u001b[0m \u001b[0mself\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0mpf\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0marr\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0mv\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0minput_units\u001b[0m \u001b[0;34m=\u001b[0m \u001b[0mfinfos\u001b[0m\u001b[0;34m[\u001b[0m\u001b[0mf\u001b[0m\u001b[0;34m]\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0munits\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n",
- "\u001b[0;32m/Users/goldbaum/Documents/yt-hg/yt/geometry/geometry_handler.pyc\u001b[0m in \u001b[0;36m_read_particle_fields\u001b[0;34m(self, fields, dobj, chunk)\u001b[0m\n\u001b[1;32m 231\u001b[0m \u001b[0mself\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0m_chunk_io\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0mdobj\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0mcache\u001b[0m \u001b[0;34m=\u001b[0m \u001b[0mFalse\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m,\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[1;32m 232\u001b[0m \u001b[0mselector\u001b[0m\u001b[0;34m,\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[0;32m--> 233\u001b[0;31m fields_to_read)\n\u001b[0m\u001b[1;32m 234\u001b[0m \u001b[0;32mfor\u001b[0m \u001b[0mfield\u001b[0m \u001b[0;32min\u001b[0m \u001b[0mfields_to_read\u001b[0m\u001b[0;34m:\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[1;32m 235\u001b[0m \u001b[0mftype\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0mfname\u001b[0m \u001b[0;34m=\u001b[0m \u001b[0mfield\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n",
- "\u001b[0;32m/Users/goldbaum/Documents/yt-hg/yt/utilities/io_handler.pyc\u001b[0m in \u001b[0;36m_read_particle_selection\u001b[0;34m(self, chunks, selector, fields)\u001b[0m\n\u001b[1;32m 172\u001b[0m \u001b[0;31m#mylog.debug(\"Filling %s from %s to %s with %s\",\u001b[0m\u001b[0;34m\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[1;32m 173\u001b[0m \u001b[0;31m# field_f, my_ind, my_ind+vals.shape[0], field_r)\u001b[0m\u001b[0;34m\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[0;32m--> 174\u001b[0;31m \u001b[0mrv\u001b[0m\u001b[0;34m[\u001b[0m\u001b[0mfield_f\u001b[0m\u001b[0;34m]\u001b[0m\u001b[0;34m[\u001b[0m\u001b[0mmy_ind\u001b[0m\u001b[0;34m:\u001b[0m\u001b[0mmy_ind\u001b[0m \u001b[0;34m+\u001b[0m \u001b[0mvals\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0mshape\u001b[0m\u001b[0;34m[\u001b[0m\u001b[0;36m0\u001b[0m\u001b[0;34m]\u001b[0m\u001b[0;34m,\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0;34m]\u001b[0m \u001b[0;34m=\u001b[0m \u001b[0mvals\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[0m\u001b[1;32m 175\u001b[0m \u001b[0mind\u001b[0m\u001b[0;34m[\u001b[0m\u001b[0mfield_f\u001b[0m\u001b[0;34m]\u001b[0m \u001b[0;34m+=\u001b[0m \u001b[0mvals\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0mshape\u001b[0m\u001b[0;34m[\u001b[0m\u001b[0;36m0\u001b[0m\u001b[0;34m]\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[1;32m 176\u001b[0m \u001b[0;31m# Now we need to truncate all our fields, since we allow for\u001b[0m\u001b[0;34m\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n",
- "\u001b[0;31mValueError\u001b[0m: could not broadcast input array from shape (5000000) into shape (15872)"
- ]
- }
- ],
- "prompt_number": 7
+ "outputs": []
}
],
"metadata": {}
https://bitbucket.org/yt_analysis/yt/commits/a3eb6511c590/
Changeset: a3eb6511c590
Branch: yt-3.0
User: ngoldbaum
Date: 2014-07-18 22:01:37
Summary: Updating to bootcamp to no longer reference the yt-doc repo.
Affected #: 1 file
diff -r 2e1effc719298feb546edf3222fabb0e3d1e8ed5 -r a3eb6511c59060f0d288068b339bc09beb012c62 doc/source/bootcamp/index.rst
--- a/doc/source/bootcamp/index.rst
+++ b/doc/source/bootcamp/index.rst
@@ -12,19 +12,20 @@
on time, you can non-interactively go through the following pages and view the
worked examples.
-To execute the bootcamp interactively, you need to download the repository
-and start the IPython notebook. The easiest way to get the repository is to
-use your already-installed mercurial program to grab it:
+To execute the bootcamp interactively, you need to download the repository and
+start the IPython notebook. If you do not already have the yt repository, the
+easiest way to get the repository is to use your already-installed mercurial
+program to grab it:
.. code-block:: bash
- hg clone https://bitbucket.org/yt_analysis/yt-doc
+ hg clone https://bitbucket.org/yt_analysis/yt
Now start the IPython notebook from within the repository:
.. code-block:: bash
- cd yt-doc/source/bootcamp
+ cd yt/doc/source/bootcamp
yt notebook
This command will give you information about the notebook server and how to
https://bitbucket.org/yt_analysis/yt/commits/f244da87178a/
Changeset: f244da87178a
Branch: yt-3.0
User: ngoldbaum
Date: 2014-07-18 23:30:43
Summary: Merging to clear conflicts.
Affected #: 335 files
diff -r a3eb6511c59060f0d288068b339bc09beb012c62 -r f244da87178a1bb4396b3e3959004ac38f14c6ac doc/cheatsheet.tex
--- a/doc/cheatsheet.tex
+++ b/doc/cheatsheet.tex
@@ -208,38 +208,38 @@
After that, simulation data is generally accessed in yt using {\it Data Containers} which are Python objects
that define a region of simulation space from which data should be selected.
\settowidth{\MyLen}{\texttt{multicol} }
-\texttt{pf = load(}{\it dataset}\texttt{)} \textemdash\ Reference a single snapshot.\\
-\texttt{dd = pf.h.all\_data()} \textemdash\ Select the entire volume.\\
+\texttt{ds = load(}{\it dataset}\texttt{)} \textemdash\ Reference a single snapshot.\\
+\texttt{dd = ds.all\_data()} \textemdash\ Select the entire volume.\\
\texttt{a = dd[}{\it field\_name}\texttt{]} \textemdash\ Saves the contents of {\it field} into the
numpy array \texttt{a}. Similarly for other data containers.\\
-\texttt{pf.h.field\_list} \textemdash\ A list of available fields in the snapshot. \\
-\texttt{pf.h.derived\_field\_list} \textemdash\ A list of available derived fields
+\texttt{ds.field\_list} \textemdash\ A list of available fields in the snapshot. \\
+\texttt{ds.derived\_field\_list} \textemdash\ A list of available derived fields
in the snapshot. \\
-\texttt{val, loc = pf.h.find\_max("Density")} \textemdash\ Find the \texttt{val}ue of
+\texttt{val, loc = ds.find\_max("Density")} \textemdash\ Find the \texttt{val}ue of
the maximum of the field \texttt{Density} and its \texttt{loc}ation. \\
-\texttt{sp = pf.sphere(}{\it cen}\texttt{,}{\it radius}\texttt{)} \textemdash\ Create a spherical data
+\texttt{sp = ds.sphere(}{\it cen}\texttt{,}{\it radius}\texttt{)} \textemdash\ Create a spherical data
container. {\it cen} may be a coordinate, or ``max'' which
centers on the max density point. {\it radius} may be a float in
code units or a tuple of ({\it length, unit}).\\
-\texttt{re = pf.region({\it cen}, {\it left edge}, {\it right edge})} \textemdash\ Create a
+\texttt{re = ds.region({\it cen}, {\it left edge}, {\it right edge})} \textemdash\ Create a
rectilinear data container. {\it cen} is required but not used.
{\it left} and {\it right edge} are coordinate values that define the region.
-\texttt{di = pf.disk({\it cen}, {\it normal}, {\it radius}, {\it height})} \textemdash\
+\texttt{di = ds.disk({\it cen}, {\it normal}, {\it radius}, {\it height})} \textemdash\
Create a cylindrical data container centered at {\it cen} along the
direction set by {\it normal},with total length
2$\times${\it height} and with radius {\it radius}. \\
- \texttt{bl = pf.boolean({\it constructor})} \textemdash\ Create a boolean data
+ \texttt{bl = ds.boolean({\it constructor})} \textemdash\ Create a boolean data
container. {\it constructor} is a list of pre-defined non-boolean
data containers with nested boolean logic using the
``AND'', ``NOT'', or ``OR'' operators. E.g. {\it constructor=}
{\it [sp, ``NOT'', (di, ``OR'', re)]} gives a volume defined
by {\it sp} minus the patches covered by {\it di} and {\it re}.\\
-\texttt{pf.h.save\_object(sp, {\it ``sp\_for\_later''})} \textemdash\ Save an object (\texttt{sp}) for later use.\\
-\texttt{sp = pf.h.load\_object({\it ``sp\_for\_later''})} \textemdash\ Recover a saved object.\\
+\texttt{ds.save\_object(sp, {\it ``sp\_for\_later''})} \textemdash\ Save an object (\texttt{sp}) for later use.\\
+\texttt{sp = ds.load\_object({\it ``sp\_for\_later''})} \textemdash\ Recover a saved object.\\
\subsection{Defining New Fields \& Quantities}
@@ -261,15 +261,15 @@
\subsection{Slices and Projections}
\settowidth{\MyLen}{\texttt{multicol} }
-\texttt{slc = SlicePlot(pf, {\it axis}, {\it field}, {\it center=}, {\it width=}, {\it weight\_field=}, {\it additional parameters})} \textemdash\ Make a slice plot
+\texttt{slc = SlicePlot(ds, {\it axis}, {\it field}, {\it center=}, {\it width=}, {\it weight\_field=}, {\it additional parameters})} \textemdash\ Make a slice plot
perpendicular to {\it axis} of {\it field} weighted by {\it weight\_field} at (code-units) {\it center} with
{\it width} in code units or a (value, unit) tuple. Hint: try {\it SlicePlot?} in IPython to see additional parameters.\\
\texttt{slc.save({\it file\_prefix})} \textemdash\ Save the slice to a png with name prefix {\it file\_prefix}.
\texttt{.save()} works similarly for the commands below.\\
-\texttt{prj = ProjectionPlot(pf, {\it axis}, {\it field}, {\it addit. params})} \textemdash\ Make a projection. \\
-\texttt{prj = OffAxisSlicePlot(pf, {\it normal}, {\it fields}, {\it center=}, {\it width=}, {\it depth=},{\it north\_vector=},{\it weight\_field=})} \textemdash Make an off-axis slice. Note this takes an array of fields. \\
-\texttt{prj = OffAxisProjectionPlot(pf, {\it normal}, {\it fields}, {\it center=}, {\it width=}, {\it depth=},{\it north\_vector=},{\it weight\_field=})} \textemdash Make an off axis projection. Note this takes an array of fields. \\
+\texttt{prj = ProjectionPlot(ds, {\it axis}, {\it field}, {\it addit. params})} \textemdash\ Make a projection. \\
+\texttt{prj = OffAxisSlicePlot(ds, {\it normal}, {\it fields}, {\it center=}, {\it width=}, {\it depth=},{\it north\_vector=},{\it weight\_field=})} \textemdash Make an off-axis slice. Note this takes an array of fields. \\
+\texttt{prj = OffAxisProjectionPlot(ds, {\it normal}, {\it fields}, {\it center=}, {\it width=}, {\it depth=},{\it north\_vector=},{\it weight\_field=})} \textemdash Make an off axis projection. Note this takes an array of fields. \\
\subsection{Plot Annotations}
\settowidth{\MyLen}{\texttt{multicol} }
@@ -365,8 +365,8 @@
\subsection{FAQ}
\settowidth{\MyLen}{\texttt{multicol}}
-\texttt{pf.field\_info[`field'].take\_log = False} \textemdash\ When plotting \texttt{field}, do not take log.
-Must enter \texttt{pf.h} before this command. \\
+\texttt{ds.field\_info[`field'].take\_log = False} \textemdash\ When plotting \texttt{field}, do not take log.
+Must enter \texttt{ds.index} before this command. \\
%\rule{0.3\linewidth}{0.25pt}
diff -r a3eb6511c59060f0d288068b339bc09beb012c62 -r f244da87178a1bb4396b3e3959004ac38f14c6ac doc/coding_styleguide.txt
--- a/doc/coding_styleguide.txt
+++ b/doc/coding_styleguide.txt
@@ -49,7 +49,7 @@
* Don't create a new class to replicate the functionality of an old class --
replace the old class. Too many options makes for a confusing user
experience.
- * Parameter files are a last resort.
+ * Parameter files external to yt are a last resort.
* The usage of the **kwargs construction should be avoided. If they cannot
be avoided, they must be explained, even if they are only to be passed on to
a nested function.
@@ -61,7 +61,7 @@
* Hard-coding parameter names that are the same as those in Enzo. The
following translation table should be of some help. Note that the
parameters are now properties on a Dataset subclass: you access them
- like pf.refine_by .
+ like ds.refine_by .
* RefineBy => refine_by
* TopGridRank => dimensionality
* TopGridDimensions => domain_dimensions
diff -r a3eb6511c59060f0d288068b339bc09beb012c62 -r f244da87178a1bb4396b3e3959004ac38f14c6ac doc/docstring_example.txt
--- a/doc/docstring_example.txt
+++ b/doc/docstring_example.txt
@@ -73,7 +73,7 @@
Examples
--------
These are written in doctest format, and should illustrate how to
- use the function. Use the variables 'pf' for the parameter file, 'pc' for
+ use the function. Use the variables 'ds' for the dataset, 'pc' for
a plot collection, 'c' for a center, and 'L' for a vector.
>>> a=[1,2,3]
diff -r a3eb6511c59060f0d288068b339bc09beb012c62 -r f244da87178a1bb4396b3e3959004ac38f14c6ac doc/docstring_idioms.txt
--- a/doc/docstring_idioms.txt
+++ b/doc/docstring_idioms.txt
@@ -19,7 +19,7 @@
useful variable names that correspond to specific instances that the user is
presupposed to have created.
- * `pf`: a parameter file, loaded successfully
+ * `ds`: a dataset, loaded successfully
* `sp`: a sphere
* `c`: a 3-component "center"
* `L`: a 3-component vector that corresponds to either angular momentum or a
diff -r a3eb6511c59060f0d288068b339bc09beb012c62 -r f244da87178a1bb4396b3e3959004ac38f14c6ac doc/helper_scripts/parse_cb_list.py
--- a/doc/helper_scripts/parse_cb_list.py
+++ b/doc/helper_scripts/parse_cb_list.py
@@ -2,7 +2,7 @@
import inspect
from textwrap import TextWrapper
-pf = load("RD0005-mine/RedshiftOutput0005")
+ds = load("RD0005-mine/RedshiftOutput0005")
output = open("source/visualizing/_cb_docstrings.inc", "w")
diff -r a3eb6511c59060f0d288068b339bc09beb012c62 -r f244da87178a1bb4396b3e3959004ac38f14c6ac doc/helper_scripts/parse_dq_list.py
--- a/doc/helper_scripts/parse_dq_list.py
+++ b/doc/helper_scripts/parse_dq_list.py
@@ -2,7 +2,7 @@
import inspect
from textwrap import TextWrapper
-pf = load("RD0005-mine/RedshiftOutput0005")
+ds = load("RD0005-mine/RedshiftOutput0005")
output = open("source/analyzing/_dq_docstrings.inc", "w")
@@ -29,7 +29,7 @@
docstring = docstring))
#docstring = "\n".join(tw.wrap(docstring))))
-dd = pf.h.all_data()
+dd = ds.all_data()
for n,func in sorted(dd.quantities.functions.items()):
print n, func
write_docstring(output, n, func[1])
diff -r a3eb6511c59060f0d288068b339bc09beb012c62 -r f244da87178a1bb4396b3e3959004ac38f14c6ac doc/helper_scripts/parse_object_list.py
--- a/doc/helper_scripts/parse_object_list.py
+++ b/doc/helper_scripts/parse_object_list.py
@@ -2,7 +2,7 @@
import inspect
from textwrap import TextWrapper
-pf = load("RD0005-mine/RedshiftOutput0005")
+ds = load("RD0005-mine/RedshiftOutput0005")
output = open("source/analyzing/_obj_docstrings.inc", "w")
@@ -27,7 +27,7 @@
f.write(template % dict(clsname = clsname, sig = sig, clsproxy=clsproxy,
docstring = 'physical-object-api'))
-for n,c in sorted(pf.h.__dict__.items()):
+for n,c in sorted(ds.__dict__.items()):
if hasattr(c, '_con_args'):
print n
write_docstring(output, n, c)
diff -r a3eb6511c59060f0d288068b339bc09beb012c62 -r f244da87178a1bb4396b3e3959004ac38f14c6ac doc/helper_scripts/show_fields.py
--- a/doc/helper_scripts/show_fields.py
+++ b/doc/helper_scripts/show_fields.py
@@ -17,15 +17,15 @@
everywhere, "Enzo" fields in Enzo datasets, "Orion" fields in Orion datasets,
and so on.
-Try using the ``pf.field_list`` and ``pf.derived_field_list`` to view the
+Try using the ``ds.field_list`` and ``ds.derived_field_list`` to view the
native and derived fields available for your dataset respectively. For example
to display the native fields in alphabetical order:
.. notebook-cell::
from yt.mods import *
- pf = load("Enzo_64/DD0043/data0043")
- for i in sorted(pf.field_list):
+ ds = load("Enzo_64/DD0043/data0043")
+ for i in sorted(ds.field_list):
print i
.. note:: Universal fields will be overridden by a code-specific field.
diff -r a3eb6511c59060f0d288068b339bc09beb012c62 -r f244da87178a1bb4396b3e3959004ac38f14c6ac doc/source/analyzing/_obj_docstrings.inc
--- a/doc/source/analyzing/_obj_docstrings.inc
+++ b/doc/source/analyzing/_obj_docstrings.inc
@@ -1,12 +1,12 @@
-.. class:: boolean(self, regions, fields=None, pf=None, **field_parameters):
+.. class:: boolean(self, regions, fields=None, ds=None, **field_parameters):
For more information, see :ref:`physical-object-api`
(This is a proxy for :class:`~yt.data_objects.data_containers.AMRBooleanRegionBase`.)
-.. class:: covering_grid(self, level, left_edge, dims, fields=None, pf=None, num_ghost_zones=0, use_pbar=True, **field_parameters):
+.. class:: covering_grid(self, level, left_edge, dims, fields=None, ds=None, num_ghost_zones=0, use_pbar=True, **field_parameters):
For more information, see :ref:`physical-object-api`
(This is a proxy for :class:`~yt.data_objects.data_containers.AMRCoveringGridBase`.)
@@ -24,13 +24,13 @@
(This is a proxy for :class:`~yt.data_objects.data_containers.AMRCuttingPlaneBase`.)
-.. class:: disk(self, center, normal, radius, height, fields=None, pf=None, **field_parameters):
+.. class:: disk(self, center, normal, radius, height, fields=None, ds=None, **field_parameters):
For more information, see :ref:`physical-object-api`
(This is a proxy for :class:`~yt.data_objects.data_containers.AMRCylinderBase`.)
-.. class:: ellipsoid(self, center, A, B, C, e0, tilt, fields=None, pf=None, **field_parameters):
+.. class:: ellipsoid(self, center, A, B, C, e0, tilt, fields=None, ds=None, **field_parameters):
For more information, see :ref:`physical-object-api`
(This is a proxy for :class:`~yt.data_objects.data_containers.AMREllipsoidBase`.)
@@ -48,79 +48,79 @@
(This is a proxy for :class:`~yt.data_objects.data_containers.AMRFixedResCuttingPlaneBase`.)
-.. class:: fixed_res_proj(self, axis, level, left_edge, dims, fields=None, pf=None, **field_parameters):
+.. class:: fixed_res_proj(self, axis, level, left_edge, dims, fields=None, ds=None, **field_parameters):
For more information, see :ref:`physical-object-api`
(This is a proxy for :class:`~yt.data_objects.data_containers.AMRFixedResProjectionBase`.)
-.. class:: grid_collection(self, center, grid_list, fields=None, pf=None, **field_parameters):
+.. class:: grid_collection(self, center, grid_list, fields=None, ds=None, **field_parameters):
For more information, see :ref:`physical-object-api`
(This is a proxy for :class:`~yt.data_objects.data_containers.AMRGridCollectionBase`.)
-.. class:: grid_collection_max_level(self, center, max_level, fields=None, pf=None, **field_parameters):
+.. class:: grid_collection_max_level(self, center, max_level, fields=None, ds=None, **field_parameters):
For more information, see :ref:`physical-object-api`
(This is a proxy for :class:`~yt.data_objects.data_containers.AMRMaxLevelCollectionBase`.)
-.. class:: inclined_box(self, origin, box_vectors, fields=None, pf=None, **field_parameters):
+.. class:: inclined_box(self, origin, box_vectors, fields=None, ds=None, **field_parameters):
For more information, see :ref:`physical-object-api`
(This is a proxy for :class:`~yt.data_objects.data_containers.AMRInclinedBoxBase`.)
-.. class:: ortho_ray(self, axis, coords, fields=None, pf=None, **field_parameters):
+.. class:: ortho_ray(self, axis, coords, fields=None, ds=None, **field_parameters):
For more information, see :ref:`physical-object-api`
(This is a proxy for :class:`~yt.data_objects.data_containers.AMROrthoRayBase`.)
-.. class:: overlap_proj(self, axis, field, weight_field=None, max_level=None, center=None, pf=None, source=None, node_name=None, field_cuts=None, preload_style='level', serialize=True, **field_parameters):
+.. class:: overlap_proj(self, axis, field, weight_field=None, max_level=None, center=None, ds=None, source=None, node_name=None, field_cuts=None, preload_style='level', serialize=True, **field_parameters):
For more information, see :ref:`physical-object-api`
(This is a proxy for :class:`~yt.data_objects.data_containers.AMRProjBase`.)
-.. class:: periodic_region(self, center, left_edge, right_edge, fields=None, pf=None, **field_parameters):
+.. class:: periodic_region(self, center, left_edge, right_edge, fields=None, ds=None, **field_parameters):
For more information, see :ref:`physical-object-api`
(This is a proxy for :class:`~yt.data_objects.data_containers.AMRPeriodicRegionBase`.)
-.. class:: periodic_region_strict(self, center, left_edge, right_edge, fields=None, pf=None, **field_parameters):
+.. class:: periodic_region_strict(self, center, left_edge, right_edge, fields=None, ds=None, **field_parameters):
For more information, see :ref:`physical-object-api`
(This is a proxy for :class:`~yt.data_objects.data_containers.AMRPeriodicRegionStrictBase`.)
-.. class:: proj(self, axis, field, weight_field=None, max_level=None, center=None, pf=None, source=None, node_name=None, field_cuts=None, preload_style=None, serialize=True, style='integrate', **field_parameters):
+.. class:: proj(self, axis, field, weight_field=None, max_level=None, center=None, ds=None, source=None, node_name=None, field_cuts=None, preload_style=None, serialize=True, style='integrate', **field_parameters):
For more information, see :ref:`physical-object-api`
(This is a proxy for :class:`~yt.data_objects.data_containers.AMRQuadTreeProjBase`.)
-.. class:: ray(self, start_point, end_point, fields=None, pf=None, **field_parameters):
+.. class:: ray(self, start_point, end_point, fields=None, ds=None, **field_parameters):
For more information, see :ref:`physical-object-api`
(This is a proxy for :class:`~yt.data_objects.data_containers.AMRRayBase`.)
-.. class:: region(self, center, left_edge, right_edge, fields=None, pf=None, **field_parameters):
+.. class:: region(self, center, left_edge, right_edge, fields=None, ds=None, **field_parameters):
For more information, see :ref:`physical-object-api`
(This is a proxy for :class:`~yt.data_objects.data_containers.AMRRegionBase`.)
-.. class:: region_strict(self, center, left_edge, right_edge, fields=None, pf=None, **field_parameters):
+.. class:: region_strict(self, center, left_edge, right_edge, fields=None, ds=None, **field_parameters):
For more information, see :ref:`physical-object-api`
(This is a proxy for :class:`~yt.data_objects.data_containers.AMRRegionStrictBase`.)
-.. class:: slice(self, axis, coord, fields=None, center=None, pf=None, node_name=False, **field_parameters):
+.. class:: slice(self, axis, coord, fields=None, center=None, ds=None, node_name=False, **field_parameters):
For more information, see :ref:`physical-object-api`
(This is a proxy for :class:`~yt.data_objects.data_containers.AMRSliceBase`.)
@@ -132,13 +132,13 @@
(This is a proxy for :class:`~yt.data_objects.data_containers.AMRSmoothedCoveringGridBase`.)
-.. class:: sphere(self, center, radius, fields=None, pf=None, **field_parameters):
+.. class:: sphere(self, center, radius, fields=None, ds=None, **field_parameters):
For more information, see :ref:`physical-object-api`
(This is a proxy for :class:`~yt.data_objects.data_containers.AMRSphereBase`.)
-.. class:: streamline(self, positions, length=1.0, fields=None, pf=None, **field_parameters):
+.. class:: streamline(self, positions, length=1.0, fields=None, ds=None, **field_parameters):
For more information, see :ref:`physical-object-api`
(This is a proxy for :class:`~yt.data_objects.data_containers.AMRStreamlineBase`.)
diff -r a3eb6511c59060f0d288068b339bc09beb012c62 -r f244da87178a1bb4396b3e3959004ac38f14c6ac doc/source/analyzing/analysis_modules/Halo_Analysis.ipynb
--- a/doc/source/analyzing/analysis_modules/Halo_Analysis.ipynb
+++ b/doc/source/analyzing/analysis_modules/Halo_Analysis.ipynb
@@ -44,7 +44,7 @@
"tmpdir = tempfile.mkdtemp()\n",
"\n",
"# Load the data set with the full simulation information\n",
- "data_pf = load('Enzo_64/RD0006/RedshiftOutput0006')"
+ "data_ds = load('Enzo_64/RD0006/RedshiftOutput0006')"
],
"language": "python",
"metadata": {},
@@ -62,7 +62,7 @@
"collapsed": false,
"input": [
"# Load the rockstar data files\n",
- "halos_pf = load('rockstar_halos/halos_0.0.bin')"
+ "halos_ds = load('rockstar_halos/halos_0.0.bin')"
],
"language": "python",
"metadata": {},
@@ -80,7 +80,7 @@
"collapsed": false,
"input": [
"# Instantiate a catalog using those two paramter files\n",
- "hc = HaloCatalog(data_pf=data_pf, halos_pf=halos_pf, \n",
+ "hc = HaloCatalog(data_ds=data_ds, halos_ds=halos_ds, \n",
" output_dir=os.path.join(tmpdir, 'halo_catalog'))"
],
"language": "python",
@@ -295,9 +295,9 @@
"cell_type": "code",
"collapsed": false,
"input": [
- "halos_pf = load(os.path.join(tmpdir, 'halo_catalog/halo_catalog.0.h5'))\n",
+ "halos_ds = load(os.path.join(tmpdir, 'halo_catalog/halo_catalog.0.h5'))\n",
"\n",
- "hc_reloaded = HaloCatalog(halos_pf=halos_pf,\n",
+ "hc_reloaded = HaloCatalog(halos_ds=halos_ds,\n",
" output_dir=os.path.join(tmpdir, 'halo_catalog'))"
],
"language": "python",
@@ -407,4 +407,4 @@
"metadata": {}
}
]
-}
\ No newline at end of file
+}
diff -r a3eb6511c59060f0d288068b339bc09beb012c62 -r f244da87178a1bb4396b3e3959004ac38f14c6ac doc/source/analyzing/analysis_modules/PPVCube.ipynb
--- a/doc/source/analyzing/analysis_modules/PPVCube.ipynb
+++ b/doc/source/analyzing/analysis_modules/PPVCube.ipynb
@@ -222,7 +222,7 @@
"cell_type": "code",
"collapsed": false,
"input": [
- "pf = load(\"cube.fits\")"
+ "ds = load(\"cube.fits\")"
],
"language": "python",
"metadata": {},
@@ -233,7 +233,7 @@
"collapsed": false,
"input": [
"# Specifying no center gives us the center slice\n",
- "slc = SlicePlot(pf, \"z\", [\"density\"])\n",
+ "slc = SlicePlot(ds, \"z\", [\"density\"])\n",
"slc.show()"
],
"language": "python",
@@ -246,9 +246,9 @@
"input": [
"import yt.units as u\n",
"# Picking different velocities for the slices\n",
- "new_center = pf.domain_center\n",
- "new_center[2] = pf.spec2pixel(-1.0*u.km/u.s)\n",
- "slc = SlicePlot(pf, \"z\", [\"density\"], center=new_center)\n",
+ "new_center = ds.domain_center\n",
+ "new_center[2] = ds.spec2pixel(-1.0*u.km/u.s)\n",
+ "slc = SlicePlot(ds, \"z\", [\"density\"], center=new_center)\n",
"slc.show()"
],
"language": "python",
@@ -259,8 +259,8 @@
"cell_type": "code",
"collapsed": false,
"input": [
- "new_center[2] = pf.spec2pixel(0.7*u.km/u.s)\n",
- "slc = SlicePlot(pf, \"z\", [\"density\"], center=new_center)\n",
+ "new_center[2] = ds.spec2pixel(0.7*u.km/u.s)\n",
+ "slc = SlicePlot(ds, \"z\", [\"density\"], center=new_center)\n",
"slc.show()"
],
"language": "python",
@@ -271,8 +271,8 @@
"cell_type": "code",
"collapsed": false,
"input": [
- "new_center[2] = pf.spec2pixel(-0.3*u.km/u.s)\n",
- "slc = SlicePlot(pf, \"z\", [\"density\"], center=new_center)\n",
+ "new_center[2] = ds.spec2pixel(-0.3*u.km/u.s)\n",
+ "slc = SlicePlot(ds, \"z\", [\"density\"], center=new_center)\n",
"slc.show()"
],
"language": "python",
@@ -290,7 +290,7 @@
"cell_type": "code",
"collapsed": false,
"input": [
- "prj = ProjectionPlot(pf, \"z\", [\"density\"], proj_style=\"sum\")\n",
+ "prj = ProjectionPlot(ds, \"z\", [\"density\"], proj_style=\"sum\")\n",
"prj.set_log(\"density\", True)\n",
"prj.set_zlim(\"density\", 1.0e-3, 0.2)\n",
"prj.show()"
@@ -303,4 +303,4 @@
"metadata": {}
}
]
-}
\ No newline at end of file
+}
diff -r a3eb6511c59060f0d288068b339bc09beb012c62 -r f244da87178a1bb4396b3e3959004ac38f14c6ac doc/source/analyzing/analysis_modules/clump_finding.rst
--- a/doc/source/analyzing/analysis_modules/clump_finding.rst
+++ b/doc/source/analyzing/analysis_modules/clump_finding.rst
@@ -84,8 +84,8 @@
from yt.mods import *
- pf = load("DD0000")
- sp = pf.sphere([0.5, 0.5, 0.5], radius=0.1)
+ ds = load("DD0000")
+ sp = ds.sphere([0.5, 0.5, 0.5], radius=0.1)
ratio = sp.quantities["IsBound"](truncate=False, include_thermal_energy=True,
treecode=True, opening_angle=2.0)
@@ -97,8 +97,8 @@
from yt.mods import *
- pf = load("DD0000")
- sp = pf.sphere([0.5, 0.5, 0.5], radius=0.1)
+ ds = load("DD0000")
+ sp = ds.sphere([0.5, 0.5, 0.5], radius=0.1)
ratio = sp.quantities["IsBound"](truncate=False, include_thermal_energy=True,
treecode=False)
diff -r a3eb6511c59060f0d288068b339bc09beb012c62 -r f244da87178a1bb4396b3e3959004ac38f14c6ac doc/source/analyzing/analysis_modules/ellipsoid_analysis.rst
--- a/doc/source/analyzing/analysis_modules/ellipsoid_analysis.rst
+++ b/doc/source/analyzing/analysis_modules/ellipsoid_analysis.rst
@@ -58,8 +58,8 @@
from yt.mods import *
from yt.analysis_modules.halo_finding.api import *
- pf=load('Enzo_64/RD0006/RedshiftOutput0006')
- halo_list = parallelHF(pf)
+ ds=load('Enzo_64/RD0006/RedshiftOutput0006')
+ halo_list = parallelHF(ds)
halo_list.dump('MyHaloList')
Ellipsoid Parameters
@@ -69,8 +69,8 @@
from yt.mods import *
from yt.analysis_modules.halo_finding.api import *
- pf=load('Enzo_64/RD0006/RedshiftOutput0006')
- haloes = LoadHaloes(pf, 'MyHaloList')
+ ds=load('Enzo_64/RD0006/RedshiftOutput0006')
+ haloes = LoadHaloes(ds, 'MyHaloList')
Once the halo information is saved you can load it into the data
object "haloes", you can get loop over the list of haloes and do
@@ -107,7 +107,7 @@
.. code-block:: python
- ell = pf.ellipsoid(ell_param[0],
+ ell = ds.ellipsoid(ell_param[0],
ell_param[1],
ell_param[2],
ell_param[3],
diff -r a3eb6511c59060f0d288068b339bc09beb012c62 -r f244da87178a1bb4396b3e3959004ac38f14c6ac doc/source/analyzing/analysis_modules/halo_analysis.rst
--- a/doc/source/analyzing/analysis_modules/halo_analysis.rst
+++ b/doc/source/analyzing/analysis_modules/halo_analysis.rst
@@ -8,6 +8,7 @@
:maxdepth: 1
halo_catalogs
+ halo_transition
halo_finding
halo_mass_function
halo_analysis_example
diff -r a3eb6511c59060f0d288068b339bc09beb012c62 -r f244da87178a1bb4396b3e3959004ac38f14c6ac doc/source/analyzing/analysis_modules/halo_catalogs.rst
--- a/doc/source/analyzing/analysis_modules/halo_catalogs.rst
+++ b/doc/source/analyzing/analysis_modules/halo_catalogs.rst
@@ -7,9 +7,11 @@
together into a single framework. This framework is substantially
different from the limited framework included in yt-2.x and is only
backwards compatible in that output from old halo finders may be loaded.
+For a direct translation of various halo analysis tasks using yt-2.x
+to yt-3.0 please see :ref:`halo_transition`.
A catalog of halos can be created from any initial dataset given to halo
-catalog through data_pf. These halos can be found using friends-of-friends,
+catalog through data_ds. These halos can be found using friends-of-friends,
HOP, and Rockstar. The finder_method keyword dictates which halo finder to
use. The available arguments are 'fof', 'hop', and'rockstar'. For more
details on the relative differences between these halo finders see
@@ -19,32 +21,32 @@
from yt.mods import *
from yt.analysis_modules.halo_analysis.api import HaloCatalog
- data_pf = load('Enzo_64/RD0006/RedshiftOutput0006')
- hc = HaloCatalog(data_pf=data_pf, finder_method='hop')
+ data_ds = load('Enzo_64/RD0006/RedshiftOutput0006')
+ hc = HaloCatalog(data_ds=data_ds, finder_method='hop')
A halo catalog may also be created from already run rockstar outputs.
This method is not implemented for previously run friends-of-friends or
HOP finders. Even though rockstar creates one file per processor,
specifying any one file allows the full catalog to be loaded. Here we
only specify the file output by the processor with ID 0. Note that the
-argument for supplying a rockstar output is `halos_pf`, not `data_pf`.
+argument for supplying a rockstar output is `halos_ds`, not `data_ds`.
.. code-block:: python
- halos_pf = load(path+'rockstar_halos/halos_0.0.bin')
- hc = HaloCatalog(halos_pf=halos_pf)
+ halos_ds = load(path+'rockstar_halos/halos_0.0.bin')
+ hc = HaloCatalog(halos_ds=halos_ds)
Although supplying only the binary output of the rockstar halo finder
is sufficient for creating a halo catalog, it is not possible to find
any new information about the identified halos. To associate the halos
with the dataset from which they were found, supply arguments to both
-halos_pf and data_pf.
+halos_ds and data_ds.
.. code-block:: python
- halos_pf = load(path+'rockstar_halos/halos_0.0.bin')
- data_pf = load('Enzo_64/RD0006/RedshiftOutput0006')
- hc = HaloCatalog(data_pf=data_pf, halos_pf=halos_pf)
+ halos_ds = load(path+'rockstar_halos/halos_0.0.bin')
+ data_ds = load('Enzo_64/RD0006/RedshiftOutput0006')
+ hc = HaloCatalog(data_ds=data_ds, halos_ds=halos_ds)
A data container can also be supplied via keyword data_source,
associated with either dataset, to control the spatial region in
@@ -215,8 +217,8 @@
.. code-block:: python
- hpf = load(path+"halo_catalogs/catalog_0046/catalog_0046.0.h5")
- hc = HaloCatalog(halos_pf=hpf,
+ hds = load(path+"halo_catalogs/catalog_0046/catalog_0046.0.h5")
+ hc = HaloCatalog(halos_ds=hds,
output_dir="halo_catalogs/catalog_0046")
hc.add_callback("load_profiles", output_dir="profiles",
filename="virial_profiles")
diff -r a3eb6511c59060f0d288068b339bc09beb012c62 -r f244da87178a1bb4396b3e3959004ac38f14c6ac doc/source/analyzing/analysis_modules/halo_finders.rst
--- /dev/null
+++ b/doc/source/analyzing/analysis_modules/halo_finders.rst
@@ -0,0 +1,192 @@
+.. _halo_finding:
+
+Halo Finding
+============
+
+There are four methods of finding particle haloes in yt. The
+recommended and default method is called HOP, a method described
+in `Eisenstein and Hut (1998)
+<http://adsabs.harvard.edu/abs/1998ApJ...498..137E>`_. A basic
+friends-of-friends (e.g. `Efstathiou et al. (1985)
+<http://adsabs.harvard.edu/abs/1985ApJS...57..241E>`_) halo
+finder is also implemented. Finally Rockstar (`Behroozi et a.
+(2011) <http://adsabs.harvard.edu/abs/2011arXiv1110.4372B>`_) is
+a 6D-phase space halo finder developed by Peter Behroozi that
+excels in finding subhalos and substrcture, but does not allow
+multiple particle masses.
+
+HOP
+---
+
+The version of HOP used in yt is an upgraded version of the
+`publicly available HOP code
+<http://cmb.as.arizona.edu/~eisenste/hop/hop.html>`_. Support
+for 64-bit floats and integers has been added, as well as
+parallel analysis through spatial decomposition. HOP builds
+groups in this fashion:
+
+ 1. Estimates the local density at each particle using a
+ smoothing kernel.
+ 2. Builds chains of linked particles by 'hopping' from one
+ particle to its densest neighbor. A particle which is
+ its own densest neighbor is the end of the chain.
+ 3. All chains that share the same densest particle are
+ grouped together.
+ 4. Groups are included, linked together, or discarded
+ depending on the user-supplied over density
+ threshold parameter. The default is 160.0.
+
+Please see the `HOP method paper
+<http://adsabs.harvard.edu/abs/1998ApJ...498..137E>`_ for
+full details.
+
+.. warning:: The FoF halo finder in yt is not thoroughly tested!
+ It is probably fine to use, but you are strongly encouraged
+ to check your results against the data for errors.
+
+Rockstar Halo Finding
+---------------------
+
+Rockstar uses an adaptive hierarchical refinement of friends-of-friends
+groups in six phase-space dimensions and one time dimension, which
+allows for robust (grid-independent, shape-independent, and noise-
+resilient) tracking of substructure. The code is prepackaged with yt,
+but also `separately available <http://code.google.com/p/rockstar>`_. The lead
+developer is Peter Behroozi, and the methods are described in `Behroozi
+et al. 2011 <http://rockstar.googlecode.com/files/rockstar_ap101911.pdf>`_.
+
+.. note:: At the moment, Rockstar does not support multiple particle masses,
+ instead using a fixed particle mass. This will not affect most dark matter
+ simulations, but does make it less useful for finding halos from the stellar
+ mass. In simulations where the highest-resolution particles all have the
+ same mass (ie: zoom-in grid based simulations), one can set up a particle
+ filter to select the lowest mass particles and perform the halo finding
+ only on those.
+
+To run the Rockstar Halo finding, you must launch python with MPI and
+parallelization enabled. While Rockstar itself does not require MPI to run,
+the MPI libraries allow yt to distribute particle information across multiple
+nodes.
+
+.. warning:: At the moment, running Rockstar inside of yt on multiple compute nodes
+ connected by an Infiniband network can be problematic. Therefore, for now
+ we recommend forcing the use of the non-Infiniband network (e.g. Ethernet)
+ using this flag: ``--mca btl ^openib``.
+ For example, here is how Rockstar might be called using 24 cores:
+ ``mpirun -n 24 --mca btl ^openib python ./run_rockstar.py --parallel``.
+
+The script above configures the Halo finder, launches a server process which
+disseminates run information and coordinates writer-reader processes.
+Afterwards, it launches reader and writer tasks, filling the available MPI
+slots, which alternately read particle information and analyze for halo
+content.
+
+The RockstarHaloFinder class has these options that can be supplied to the
+halo catalog through the ``finder_kwargs`` argument:
+
+ * ``dm_type``, the index of the dark matter particle. Default is 1.
+ * ``outbase``, This is where the out*list files that Rockstar makes should be
+ placed. Default is 'rockstar_halos'.
+ * ``num_readers``, the number of reader tasks (which are idle most of the
+ time.) Default is 1.
+ * ``num_writers``, the number of writer tasks (which are fed particles and
+ do most of the analysis). Default is MPI_TASKS-num_readers-1.
+ If left undefined, the above options are automatically
+ configured from the number of available MPI tasks.
+ * ``force_res``, the resolution that Rockstar uses for various calculations
+ and smoothing lengths. This is in units of Mpc/h.
+ If no value is provided, this parameter is automatically set to
+ the width of the smallest grid element in the simulation from the
+ last data snapshot (i.e. the one where time has evolved the
+ longest) in the time series:
+ ``ds_last.index.get_smallest_dx() * ds_last['mpch']``.
+ * ``total_particles``, if supplied, this is a pre-calculated
+ total number of dark matter
+ particles present in the simulation. For example, this is useful
+ when analyzing a series of snapshots where the number of dark
+ matter particles should not change and this will save some disk
+ access time. If left unspecified, it will
+ be calculated automatically. Default: ``None``.
+ * ``dm_only``, if set to ``True``, it will be assumed that there are
+ only dark matter particles present in the simulation.
+ This option does not modify the halos found by Rockstar, however
+ this option can save disk access time if there are no star particles
+ (or other non-dark matter particles) in the simulation. Default: ``False``.
+
+Rockstar dumps halo information in a series of text (halo*list and
+out*list) and binary (halo*bin) files inside the ``outbase`` directory.
+We use the halo list classes to recover the information.
+
+Inside the ``outbase`` directory there is a text file named ``datasets.txt``
+that records the connection between ds names and the Rockstar file names.
+
+Parallel HOP and FOF
+--------------------
+
+Both the HOP and FoF halo finders can run in parallel using simple
+spatial decomposition. In order to run them in parallel it is helpful
+to understand how it works. Below in the first plot (i) is a simplified
+depiction of three haloes labeled 1,2 and 3:
+
+.. image:: _images/ParallelHaloFinder.png
+ :width: 500
+
+Halo 3 is twice reflected around the periodic boundary conditions.
+
+In (ii), the volume has been sub-divided into four equal subregions,
+A,B,C and D, shown with dotted lines. Notice that halo 2 is now in
+two different subregions, C and D, and that halo 3 is now in three,
+A, B and D. If the halo finder is run on these four separate subregions,
+halo 1 is be identified as a single halo, but haloes 2 and 3 are split
+up into multiple haloes, which is incorrect. The solution is to give
+each subregion padding to oversample into neighboring regions.
+
+In (iii), subregion C has oversampled into the other three regions,
+with the periodic boundary conditions taken into account, shown by
+dot-dashed lines. The other subregions oversample in a similar way.
+
+The halo finder is then run on each padded subregion independently
+and simultaneously. By oversampling like this, haloes 2 and 3 will
+both be enclosed fully in at least one subregion and identified
+completely.
+
+Haloes identified with centers of mass inside the padded part of a
+subregion are thrown out, eliminating the problem of halo duplication.
+The centers for the three haloes are shown with stars. Halo 1 will
+belong to subregion A, 2 to C and 3 to B.
+
+To run with parallel halo finding, you must supply a value for
+padding in the finder_kwargs argument. The ``padding`` parameter
+is in simulation units and defaults to 0.02. This parameter is how
+much padding is added to each of the six sides of a subregion.
+This value should be 2x-3x larger than the largest expected halo
+in the simulation. It is unlikely, of course, that the largest
+object in the simulation will be on a subregion boundary, but there
+is no way of knowing before the halo finder is run.
+
+.. code-block:: python
+
+ from yt.mods import *
+ from yt.analysis_modules.halo_analysis.api import *
+ ds = load("data0001")
+ hc= HaloCatalog(data_ds =ds,finder_method='hop'
+ finder_kwargs={'padding':0.02})
+ # --or--
+ hc= HaloCatalog(data_ds =ds,finder_method='fof'
+ finder_kwargs={'padding':0.02})
+
+
+In general, a little bit of padding goes a long way, and too much
+just slows down the analysis and doesn't improve the answer (but
+doesn't change it). It may be worth your time to run the parallel
+halo finder at a few paddings to find the right amount, especially
+if you're analyzing many similar datasets.
+
+Rockstar Installation
+=====================
+
+The Rockstar is slightly patched and modified to run as a library inside of
+yt. By default it will be built with yt using the ``install_script.sh``.
+If it wasn't installed, please make sure that the installation setting
+``INST_ROCKSTAR=1`` is defined in the ``install_script.sh`` and re-run
+the installation script.
diff -r a3eb6511c59060f0d288068b339bc09beb012c62 -r f244da87178a1bb4396b3e3959004ac38f14c6ac doc/source/analyzing/analysis_modules/halo_mass_function.rst
--- a/doc/source/analyzing/analysis_modules/halo_mass_function.rst
+++ b/doc/source/analyzing/analysis_modules/halo_mass_function.rst
@@ -60,8 +60,8 @@
from yt.mods import *
from yt.analysis_modules.halo_mass_function.api import *
- pf = load("data0030")
- hmf = HaloMassFcn(pf, halo_file="FilteredQuantities.out", num_sigma_bins=200,
+ ds = load("data0030")
+ hmf = HaloMassFcn(ds, halo_file="FilteredQuantities.out", num_sigma_bins=200,
mass_column=5)
Attached to ``hmf`` is the convenience function ``write_out``, which saves
@@ -102,8 +102,8 @@
from yt.mods import *
from yt.analysis_modules.halo_mass_function.api import *
- pf = load("data0030")
- hmf = HaloMassFcn(pf, halo_file="FilteredQuantities.out",
+ ds = load("data0030")
+ hmf = HaloMassFcn(ds, halo_file="FilteredQuantities.out",
sigma8input=0.9, primordial_index=1., omega_baryon0=0.06,
fitting_function=4)
hmf.write_out(prefix='hmf')
diff -r a3eb6511c59060f0d288068b339bc09beb012c62 -r f244da87178a1bb4396b3e3959004ac38f14c6ac doc/source/analyzing/analysis_modules/halo_profiling.rst
--- a/doc/source/analyzing/analysis_modules/halo_profiling.rst
+++ /dev/null
@@ -1,451 +0,0 @@
-.. _halo_profiling:
-
-Halo Profiling
-==============
-.. sectionauthor:: Britton Smith <brittonsmith at gmail.com>,
- Stephen Skory <s at skory.us>
-
-The ``HaloProfiler`` provides a means of performing analysis on multiple halos
-in a parallel-safe way.
-
-The halo profiler performs three primary functions: radial profiles,
-projections, and custom analysis. See the cookbook for a recipe demonstrating
-all of these features.
-
-Configuring the Halo Profiler
------------------------------
-
-The only argument required to create a ``HaloProfiler`` object is the path
-to the dataset.
-
-.. code-block:: python
-
- from yt.analysis_modules.halo_profiler.api import *
- hp = HaloProfiler("enzo_tiny_cosmology/DD0046/DD0046")
-
-Most of the halo profiler's options are configured with additional keyword
-arguments:
-
- * **output_dir** (*str*): if specified, all output will be put into this path
- instead of in the dataset directories. Default: None.
-
- * **halos** (*str*): "multiple" for profiling more than one halo. In this mode
- halos are read in from a list or identified with a
- `halo finder <../cookbook/running_halofinder.html>`_. In "single" mode, the
- one and only halo center is identified automatically as the location of the
- peak in the density field. Default: "multiple".
-
- * **halo_list_file** (*str*): name of file containing the list of halos.
- The halo profiler will look for this file in the data directory.
- Default: "HopAnalysis.out".
-
- * **halo_list_format** (*str* or *dict*): the format of the halo list file.
- "yt_hop" for the format given by yt's halo finders. "enzo_hop" for the
- format written by enzo_hop. This keyword can also be given in the form of a
- dictionary specifying the column in which various properties can be found.
- For example, {"id": 0, "center": [1, 2, 3], "mass": 4, "radius": 5}.
- Default: "yt_hop".
-
- * **halo_finder_function** (*function*): If halos is set to multiple and the
- file given by halo_list_file does not exit, the halo finding function
- specified here will be called. Default: HaloFinder (yt_hop).
-
- * **halo_finder_args** (*tuple*): args given with call to halo finder function.
- Default: None.
-
- * **halo_finder_kwargs** (*dict*): kwargs given with call to halo finder
- function. Default: None.
-
- * **recenter** (*string* or function name): The name of a function
- that will be used to move the center of the halo for the purposes of
- analysis. See explanation and examples, below. Default: None, which
- is equivalent to the center of mass of the halo as output by the halo
- finder.
-
- * **halo_radius** (*float*): if no halo radii are provided in the halo list
- file, this parameter is used to specify the radius out to which radial
- profiles will be made. This keyword is also used when halos is set to
- single. Default: 0.1.
-
- * **radius_units** (*str*): the units of **halo_radius**.
- Default: "1" (code units).
-
- * **n_profile_bins** (*int*): the number of bins in the radial profiles.
- Default: 50.
-
- * **profile_output_dir** (*str*): the subdirectory, inside the data directory,
- in which radial profile output files will be created. The directory will be
- created if it does not exist. Default: "radial_profiles".
-
- * **projection_output_dir** (*str*): the subdirectory, inside the data
- directory, in which projection output files will be created. The directory
- will be created if it does not exist. Default: "projections".
-
- * **projection_width** (*float*): the width of halo projections.
- Default: 8.0.
-
- * **projection_width_units** (*str*): the units of projection_width.
- Default: "mpc".
-
- * **project_at_level** (*int* or "max"): the maximum refinement level to be
- included in projections. Default: "max" (maximum level within the dataset).
-
- * **velocity_center** (*list*): the method in which the halo bulk velocity is
- calculated (used for calculation of radial and tangential velocities. Valid
- options are:
- - ["bulk", "halo"] (Default): the velocity provided in the halo list
- - ["bulk", "sphere"]: the bulk velocity of the sphere centered on the halo center.
- - ["max", field]: the velocity of the cell that is the location of the maximum of the field specified.
-
- * **filter_quantities** (*list*): quantities from the original halo list
- file to be written out in the filtered list file. Default: ['id','center'].
-
- * **use_critical_density** (*bool*): if True, the definition of overdensity
- for virial quantities is calculated with respect to the critical
- density. If False, overdensity is with respect to mean matter density,
- which is lower by a factor of Omega_M. Default: False.
-
-Profiles
---------
-
-Once the halo profiler object has been instantiated, fields can be added for
-profiling with the :meth:`add_profile` method:
-
-.. code-block:: python
-
- hp.add_profile('cell_volume', weight_field=None, accumulation=True)
- hp.add_profile('TotalMassMsun', weight_field=None, accumulation=True)
- hp.add_profile('density', weight_field=None, accumulation=False)
- hp.add_profile('temperature', weight_field='cell_mass', accumulation=False)
- hp.make_profiles(njobs=-1, prefilters=["halo['mass'] > 1e13"],
- filename='VirialQuantities.h5')
-
-The :meth:`make_profiles` method will begin the profiling. Use the
-**njobs** keyword to control the number of jobs over which the
-profiling is divided. Setting to -1 results in a single processor per
-halo. Setting to 1 results in all available processors working on the
-same halo. The prefilters keyword tells the profiler to skip all halos with
-masses (as loaded from the halo finder) less than a given amount. See below
-for more information. Additional keyword arguments are:
-
- * **filename** (*str*): If set, a file will be written with all of the
- filtered halos and the quantities returned by the filter functions.
- Default: None.
-
- * **prefilters** (*list*): A single dataset can contain thousands or tens of
- thousands of halos. Significant time can be saved by not profiling halos
- that are certain to not pass any filter functions in place. Simple filters
- based on quantities provided in the initial halo list can be used to filter
- out unwanted halos using this parameter. Default: None.
-
- * **njobs** (*int*): The number of jobs over which to split the profiling.
- Set to -1 so that each halo is done by a single processor. Default: -1.
-
- * **dynamic** (*bool*): If True, distribute halos using a task queue. If
- False, distribute halos evenly over all jobs. Default: False.
-
- * **profile_format** (*str*): The file format for the radial profiles,
- 'ascii' or 'hdf5'. Default: 'ascii'.
-
-.. image:: _images/profiles.png
- :width: 500
-
-Radial profiles of Overdensity (left) and Temperature (right) for five halos.
-
-Projections
------------
-
-The process of making projections is similar to that of profiles:
-
-.. code-block:: python
-
- hp.add_projection('density', weight_field=None)
- hp.add_projection('temperature', weight_field='density')
- hp.add_projection('metallicity', weight_field='density')
- hp.make_projections(axes=[0, 1, 2], save_cube=True, save_images=True,
- halo_list="filtered", njobs=-1)
-
-If **save_cube** is set to True, the projection data
-will be written to a set of hdf5 files
-in the directory given by **projection_output_dir**.
-The keyword, **halo_list**, can be
-used to select between the full list of halos ("all"),
-the filtered list ("filtered"), or
-an entirely new list given in the form of a file name.
-See :ref:`filter_functions` for a
-discussion of filtering halos. Use the **njobs** keyword to control
-the number of jobs over which the profiling is divided. Setting to -1
-results in a single processor per halo. Setting to 1 results in all
-available processors working on the same halo. The keyword arguments are:
-
- * **axes** (*list*): A list of the axes to project along, using the usual
- 0,1,2 convention. Default=[0,1,2].
-
- * **halo_list** (*str*) {'filtered', 'all'}: Which set of halos to make
- profiles of, either ones passed by the halo filters (if enabled/added), or
- all halos. Default='filtered'.
-
- * **save_images** (*bool*): Whether or not to save images of the projections.
- Default=False.
-
- * **save_cube** (*bool*): Whether or not to save the HDF5 files of the halo
- projections. Default=True.
-
- * **njobs** (*int*): The number of jobs over which to split the projections.
- Set to -1 so that each halo is done by a single processor. Default: -1.
-
- * **dynamic** (*bool*): If True, distribute halos using a task queue. If
- False, distribute halos evenly over all jobs. Default: False.
-
-.. image:: _images/projections.png
- :width: 500
-
-Projections of Density (top) and Temperature,
-weighted by Density (bottom), in the x (left),
-y (middle), and z (right) directions for a single halo with a width of 8 Mpc.
-
-Halo Filters
-------------
-
-Filters can be added to create a refined list of
-halos based on their profiles or to avoid
-profiling halos altogether based on information
-given in the halo list file.
-
-.. _filter_functions:
-
-Filter Functions
-^^^^^^^^^^^^^^^^
-
-It is often the case that one is looking to
-identify halos with a specific set of
-properties. This can be accomplished through the creation
-of filter functions. A filter
-function can take as many args and kwargs as you like,
-as long as the first argument is a
-profile object, or at least a dictionary which contains
-the profile arrays for each field.
-Filter functions must return a list of two things.
-The first is a True or False indicating
-whether the halo passed the filter.
-The second is a dictionary containing quantities
-calculated for that halo that will be written to a
-file if the halo passes the filter.
-A sample filter function based on virial quantities can be found in
-``yt/analysis_modules/halo_profiler/halo_filters.py``.
-
-Halo filtering takes place during the call to :meth:`make_profiles`.
-The :meth:`add_halo_filter` method is used to add a filter to be used
-during the profiling:
-
-.. code-block:: python
-
- hp.add_halo_filter(HP.VirialFilter, must_be_virialized=True,
- overdensity_field='ActualOverdensity',
- virial_overdensity=200,
- virial_filters=[['TotalMassMsun','>=','1e14']],
- virial_quantities=['TotalMassMsun','RadiusMpc'],
- use_log=True)
-
-The addition above will calculate and return virial quantities,
-mass and radius, for an
-overdensity of 200. In order to pass the filter, at least one
-point in the profile must be
-above the specified overdensity and the virial mass must be at
-least 1e14 solar masses. The **use_log** keyword indicates that interpolation
-should be done in log space. If
-the VirialFilter function has been added to the filter list,
-the halo profiler will make
-sure that the fields necessary for calculating virial quantities are added.
-As many filters as desired can be added. If filters have been added,
-the next call to :meth:`make_profiles` will filter by all of
-the added filter functions:
-
-.. code-block:: python
-
- hp.make_profiles(filename="FilteredQuantities.out")
-
-If the **filename** keyword is set, a file will be written with all of the
-filtered halos and the quantities returned by the filter functions.
-
-.. note:: If the profiles have already been run, the halo profiler will read
- in the previously created output files instead of re-running the profiles.
- The halo profiler will check to make sure the output file contains all of
- the requested halo fields. If not, the profile will be made again from
- scratch.
-
-.. _halo_profiler_pre_filters:
-
-Pre-filters
-^^^^^^^^^^^
-
-A single dataset can contain thousands or tens of thousands of halos.
-Significant time can
-be saved by not profiling halos that are certain to not pass any filter
-functions in place.
-Simple filters based on quantities provided in the initial halo list
-can be used to filter
-out unwanted halos using the **prefilters** keyword:
-
-.. code-block:: python
-
- hp.make_profiles(filename="FilteredQuantities.out",
- prefilters=["halo['mass'] > 1e13"])
-
-Arguments provided with the **prefilters** keyword should be given
-as a list of strings.
-Each string in the list will be evaluated with an *eval*.
-
-.. note:: If a VirialFilter function has been added with a filter based
- on mass (as in the example above), a prefilter will be automatically
- added to filter out halos with masses greater or less than (depending
- on the conditional of the filter) a factor of ten of the specified
- virial mass.
-
-Recentering the Halo For Analysis
----------------------------------
-
-It is possible to move the center of the halo to a new point using an
-arbitrary function for making profiles.
-By default, the center is provided by the halo finder,
-which outputs the center of mass of the particles. For the purposes of
-analysis, it may be important to recenter onto a gas density maximum,
-or a temperature minimum.
-
-There are a number of built-in functions to do this, listed below.
-Each of the functions uses mass-weighted fields for the calculations
-of new center points.
-To use
-them, supply the HaloProfiler with the ``recenter`` option and
-the name of the function, as in the example below.
-
-.. code-block:: python
-
- hp = HaloProfiler("enzo_tiny_cosmology/DD0046/DD0046",
- recenter="Max_Dark_Matter_Density")
-
-Additional options are:
-
- * *Min_Dark_Matter_Density* - Recenter on the point of minimum dark matter
- density in the halo.
-
- * *Max_Dark_Matter_Density* - Recenter on the point of maximum dark matter
- density in the halo.
-
- * *CoM_Dark_Matter_Density* - Recenter on the center of mass of the dark
- matter density field. This will be very similar to what the halo finder
- provides, but not precisely similar.
-
- * *Min_Gas_Density* - Recenter on the point of minimum gas density in the
- halo.
-
- * *Max_Gas_Density* - Recenter on the point of maximum gas density in the
- halo.
-
- * *CoM_Gas_Density* - Recenter on the center of mass of the gas density field
- in the halo.
-
- * *Min_Total_Density* - Recenter on the point of minimum total (gas + dark
- matter) density in the halo.
-
- * *Max_Total_Density* - Recenter on the point of maximum total density in the
- halo.
-
- * *CoM_Total_Density* - Recenter on the center of mass for the total density
- in the halo.
-
- * *Min_Temperature* - Recenter on the point of minimum temperature in the
- halo.
-
- * *Max_Temperature* - Recenter on the point of maximum temperature in the
- halo.
-
-It is also possible to supply a user-defined function to the HaloProfiler.
-This can be used if the pre-defined functions above are not sufficient.
-The function takes a single argument, a data container for the halo,
-which is a sphere. The function returns a 3-list with the new center.
-
-In this example below, a function is used such that the halos will be
-re-centered on the point of absolute minimum temperature, that is not
-mass weighted.
-
-.. code-block:: python
-
- from yt.mods import *
-
- def find_min_temp(sphere):
- ma, mini, mx, my, mz, mg = sphere.quantities['MinLocation']('temperature')
- return [mx,my,mz]
-
- hp = HaloProfiler("enzo_tiny_cosmology/DD0046/DD0046", recenter=find_min_temp)
-
-It is possible to make more complicated functions. This example below extends
-the example above to include a distance control that prevents the center from
-being moved too far. If the recenter moves too far, ``[-1, -1, -1]`` is
-returned which will prevent the halo from being profiled.
-Any triplet of values less than the ``domain_left_edge`` will suffice.
-There will be a note made in the output (stderr) showing which halos were
-skipped.
-
-.. code-block:: python
-
- from yt.mods import *
- from yt.utilities.math_utils import periodic_dist
-
- def find_min_temp_dist(sphere):
- old = sphere.center
- ma, mini, mx, my, mz, mg = sphere.quantities['MinLocation']('temperature')
- d = sphere.pf['kpc'] * periodic_dist(old, [mx, my, mz],
- sphere.pf.domain_right_edge - sphere.pf.domain_left_edge)
- # If new center farther than 5 kpc away, don't recenter
- if d > 5.: return [-1, -1, -1]
- return [mx,my,mz]
-
- hp = HaloProfiler("enzo_tiny_cosmology/DD0046/DD0046",
- recenter=find_min_temp_dist)
-
-Custom Halo Analysis
---------------------
-
-Besides radial profiles and projections, the halo profiler has the
-ability to run custom analysis functions on each halo. Custom halo
-analysis functions take two arguments: a halo dictionary containing
-the id, center, etc; and a sphere object. The example function shown
-below creates a 2D profile of the total mass in bins of density and
-temperature for a given halo.
-
-.. code-block:: python
-
- from yt.mods import *
- from yt.data_objects.profiles import BinnedProfile2D
-
- def halo_2D_profile(halo, sphere):
- "Make a 2D profile for a halo."
- my_profile = BinnedProfile2D(sphere,
- 128, 'density', 1e-30, 1e-24, True,
- 128, 'temperature', 1e2, 1e7, True,
- end_collect=False)
- my_profile.add_fields('cell_mass', weight=None, fractional=False)
- my_filename = os.path.join(sphere.pf.fullpath, '2D_profiles',
- 'Halo_%04d.h5' % halo['id'])
- my_profile.write_out_h5(my_filename)
-
-Using the :meth:`analyze_halo_spheres` function, the halo profiler
-will create a sphere centered on each halo, and perform the analysis
-from the custom routine.
-
-.. code-block:: python
-
- hp.analyze_halo_sphere(halo_2D_profile, halo_list='filtered',
- analysis_output_dir='2D_profiles',
- njobs=-1, dynamic=False)
-
-Just like with the :meth:`make_projections` function, the keyword,
-**halo_list**, can be used to select between the full list of halos
-("all"), the filtered list ("filtered"), or an entirely new list given
-in the form of a file name. If the **analysis_output_dir** keyword is
-set, the halo profiler will make sure the desired directory exists in
-a parallel-safe manner. Use the **njobs** keyword to control the
-number of jobs over which the profiling is divided. Setting to -1
-results in a single processor per halo. Setting to 1 results in all
-available processors working on the same halo.
diff -r a3eb6511c59060f0d288068b339bc09beb012c62 -r f244da87178a1bb4396b3e3959004ac38f14c6ac doc/source/analyzing/analysis_modules/halo_transition.rst
--- /dev/null
+++ b/doc/source/analyzing/analysis_modules/halo_transition.rst
@@ -0,0 +1,106 @@
+
+Getting up to Speed with Halo Analysis in yt-3.0
+================================================
+
+If you're used to halo analysis in yt-2.x, heres a guide to
+how to update your analysis pipeline to take advantage of
+the new halo catalog infrastructure.
+
+Finding Halos
+-------------
+
+Previously, halos were found using calls to ``HaloFinder``,
+``FOFHaloFinder`` and ``RockstarHaloFinder``. Now it is
+encouraged that you find the halos upon creation of the halo catalog
+by supplying a value to the ``finder_method`` keyword when calling
+``HaloCatalog``. Currently, only halos found using rockstar or a
+previous instance of a halo catalog are able to be loaded
+using the ``halos_ds`` keyword.
+
+To pass additional arguments to the halo finders
+themselves, supply a dictionary to ``finder_kwargs`` where
+each key in the dictionary is a keyword of the halo finder
+and the corresponding value is the value to be passed for
+that keyword.
+
+Getting Halo Information
+------------------------
+All quantities that used to be present in a ``halo_list`` are
+still able to be found but are not necessarily included by default.
+Every halo will by default have the following properties:
+
+* particle_position_i (where i can be x,y,z)
+* particle_mass
+* virial_radius
+* particle_identifier
+
+If other quantities are desired, they can be included by adding
+the corresponding quantity before the catalog is created. See
+the full halo catalog documentation for further information about
+how to add these quantities and what quantities are available.
+
+You no longer have to iteratre over halos in the ``halo_list``.
+Now a halo dataset can be treated as a regular dataset and
+all quantities are available by accessing ``all_data``.
+Specifically, all quantities can be accessed as shown:
+
+.. code-block:: python
+ from yt.mods import *
+ from yt.analysis_modules.halo_analysis.api import HaloCatalog
+ data_ds = load('Enzo_64/RD0006/RedshiftOutput0006')
+ hc = HaloCatalog(data_ds=data_ds, finder_method='hop')
+ hc.create()
+ ad = hc.all_data()
+ masses = ad['particle_mass'][:]
+
+
+Prefiltering Halos
+------------------
+
+Prefiltering halos before analysis takes place is now done
+by adding a filter before the call to create. An example
+is shown below
+
+.. code-block:: python
+ from yt.mods import *
+ from yt.analysis_modules.halo_analysis.api import HaloCatalog
+ data_ds = load('Enzo_64/RD0006/RedshiftOutput0006')
+ hc = HaloCatalog(data_ds=data_ds, finder_method='hop')
+ hc.add_filter("quantity_value", "particle_mass", ">", 1e13, "Msun")
+ hc.create()
+
+Profiling Halos
+---------------
+
+The halo profiler available in yt-2.x has been removed, and
+profiling functionality is now completely contained within the
+halo catalog. A complete example of how to profile halos by
+radius using the new infrastructure is given in
+:ref:`halo_analysis_example`.
+
+Plotting Halos
+--------------
+
+Annotating halo locations onto a slice or projection works in
+the same way as in yt-2.x, but now a halo catalog must be
+passed to the annotate halo call rather than a halo list.
+
+.. code-block:: python
+ from yt.mods import *
+ from yt.analysis_modules.halo_analysis.api import HaloCatalog
+
+ data_ds = load('Enzo_64/RD0006/RedshiftOutput0006')
+ hc = HaloCatalog(data_ds=data_ds, finder_method='hop')
+ hc.create()
+
+ prj = ProjectionPlot(data_ds, 'z', 'density')
+ prj.annotate_halos(hc)
+ prj.save()
+
+Written Data
+------------
+
+Data is now written out in the form of h5 files rather than
+text files. The directory they are written out to is
+controlled by the keyword ``output_dir``. Each quantity
+is a field in the file.
diff -r a3eb6511c59060f0d288068b339bc09beb012c62 -r f244da87178a1bb4396b3e3959004ac38f14c6ac doc/source/analyzing/analysis_modules/hmf_howto.rst
--- a/doc/source/analyzing/analysis_modules/hmf_howto.rst
+++ b/doc/source/analyzing/analysis_modules/hmf_howto.rst
@@ -27,8 +27,8 @@
.. code-block:: python
from yt.mods import *
- pf = load("data0001")
- halo_list = HaloFinder(pf)
+ ds = load("data0001")
+ halo_list = HaloFinder(ds)
halo_list.write_out("HopAnalysis.out")
The only important columns of data in the text file ``HopAnalysis.out``
@@ -79,8 +79,8 @@
from yt.mods import *
from yt.analysis_modules.halo_mass_function.api import *
- pf = load("data0001")
- hmf = HaloMassFcn(pf, halo_file="VirialHaloes.out",
+ ds = load("data0001")
+ hmf = HaloMassFcn(ds, halo_file="VirialHaloes.out",
sigma8input=0.9, primordial_index=1., omega_baryon0=0.06,
fitting_function=4, mass_column=5, num_sigma_bins=200)
hmf.write_out(prefix='hmf')
@@ -107,9 +107,9 @@
from yt.analysis_modules.halo_mass_function.api import *
# If desired, start loop here.
- pf = load("data0001")
+ ds = load("data0001")
- halo_list = HaloFinder(pf)
+ halo_list = HaloFinder(ds)
halo_list.write_out("HopAnalysis.out")
hp = HP.HaloProfiler("data0001", halo_list_file='HopAnalysis.out')
@@ -120,7 +120,7 @@
virial_quantities=['TotalMassMsun','RadiusMpc'])
hp.make_profiles(filename="VirialHaloes.out")
- hmf = HaloMassFcn(pf, halo_file="VirialHaloes.out",
+ hmf = HaloMassFcn(ds, halo_file="VirialHaloes.out",
sigma8input=0.9, primordial_index=1., omega_baryon0=0.06,
fitting_function=4, mass_column=5, num_sigma_bins=200)
hmf.write_out(prefix='hmf')
diff -r a3eb6511c59060f0d288068b339bc09beb012c62 -r f244da87178a1bb4396b3e3959004ac38f14c6ac doc/source/analyzing/analysis_modules/light_cone_generator.rst
--- a/doc/source/analyzing/analysis_modules/light_cone_generator.rst
+++ b/doc/source/analyzing/analysis_modules/light_cone_generator.rst
@@ -60,7 +60,7 @@
when gathering datasets for time series. Default: True.
* **set_parameters** (*dict*): Dictionary of parameters to attach to
- pf.parameters. Default: None.
+ ds.parameters. Default: None.
* **output_dir** (*string*): The directory in which images and data files
will be written. Default: 'LC'.
diff -r a3eb6511c59060f0d288068b339bc09beb012c62 -r f244da87178a1bb4396b3e3959004ac38f14c6ac doc/source/analyzing/analysis_modules/merger_tree.rst
--- a/doc/source/analyzing/analysis_modules/merger_tree.rst
+++ b/doc/source/analyzing/analysis_modules/merger_tree.rst
@@ -2,8 +2,9 @@
Halo Merger Tree
================
-.. sectionauthor:: Stephen Skory <sskory at physics.ucsd.edu>
-.. versionadded:: 1.7
+
+.. note:: At the moment the merger tree is not yet implemented using new
+ halo catalog functionality.
The Halo Merger Tree extension is capable of building a database of halo mergers
over a set of time-ordered Enzo datasets. The fractional contribution of older
diff -r a3eb6511c59060f0d288068b339bc09beb012c62 -r f244da87178a1bb4396b3e3959004ac38f14c6ac doc/source/analyzing/analysis_modules/photon_simulator.rst
--- a/doc/source/analyzing/analysis_modules/photon_simulator.rst
+++ b/doc/source/analyzing/analysis_modules/photon_simulator.rst
@@ -48,7 +48,7 @@
.. code:: python
- pf = load("MHDSloshing/virgo_low_res.0054.vtk",
+ ds = load("MHDSloshing/virgo_low_res.0054.vtk",
parameters={"time_unit":(1.0,"Myr"),
"length_unit":(1.0,"Mpc"),
"mass_unit":(1.0e14,"Msun")})
@@ -423,7 +423,7 @@
evacuated two "bubbles" of radius 30 kpc at a distance of 50 kpc from
the center.
-Now, we create a parameter file out of this dataset:
+Now, we create a yt Dataset object out of this dataset:
.. code:: python
@@ -445,7 +445,7 @@
.. code:: python
- sphere = ds.sphere(pf.domain_center, (1.0,"Mpc"))
+ sphere = ds.sphere(ds.domain_center, (1.0,"Mpc"))
A = 6000.
exp_time = 2.0e5
diff -r a3eb6511c59060f0d288068b339bc09beb012c62 -r f244da87178a1bb4396b3e3959004ac38f14c6ac doc/source/analyzing/analysis_modules/radial_column_density.rst
--- a/doc/source/analyzing/analysis_modules/radial_column_density.rst
+++ b/doc/source/analyzing/analysis_modules/radial_column_density.rst
@@ -41,15 +41,15 @@
from yt.mods import *
from yt.analysis_modules.radial_column_density.api import *
- pf = load("data0030")
+ ds = load("data0030")
- rcdnumdens = RadialColumnDensity(pf, 'NumberDensity', [0.5, 0.5, 0.5],
+ rcdnumdens = RadialColumnDensity(ds, 'NumberDensity', [0.5, 0.5, 0.5],
max_radius = 0.5)
def _RCDNumberDensity(field, data, rcd = rcdnumdens):
return rcd._build_derived_field(data)
add_field('RCDNumberDensity', _RCDNumberDensity, units=r'1/\rm{cm}^2')
- dd = pf.h.all_data()
+ dd = ds.all_data()
print dd['RCDNumberDensity']
The field ``RCDNumberDensity`` can be used just like any other derived field
diff -r a3eb6511c59060f0d288068b339bc09beb012c62 -r f244da87178a1bb4396b3e3959004ac38f14c6ac doc/source/analyzing/analysis_modules/radmc3d_export.rst
--- a/doc/source/analyzing/analysis_modules/radmc3d_export.rst
+++ b/doc/source/analyzing/analysis_modules/radmc3d_export.rst
@@ -41,8 +41,8 @@
.. code-block:: python
- pf = load("galaxy0030/galaxy0030")
- writer = RadMC3DWriter(pf)
+ ds = load("galaxy0030/galaxy0030")
+ writer = RadMC3DWriter(ds)
writer.write_amr_grid()
writer.write_dust_file("DustDensity", "dust_density.inp")
@@ -87,8 +87,8 @@
return (x_co/mu_h)*data["density"]
add_field("NumberDensityCO", function=_NumberDensityCO)
- pf = load("galaxy0030/galaxy0030")
- writer = RadMC3DWriter(pf)
+ ds = load("galaxy0030/galaxy0030")
+ writer = RadMC3DWriter(ds)
writer.write_amr_grid()
writer.write_line_file("NumberDensityCO", "numberdens_co.inp")
This diff is so big that we needed to truncate the remainder.
https://bitbucket.org/yt_analysis/yt/commits/c44a150175c9/
Changeset: c44a150175c9
Branch: yt-3.0
User: ngoldbaum
Date: 2014-07-19 18:51:31
Summary: Merging with tip to resolve conflict.
Affected #: 13 files
diff -r f244da87178a1bb4396b3e3959004ac38f14c6ac -r c44a150175c9fb3ee59ad73bda43af21dd2e6c21 doc/source/analyzing/fields.rst
--- /dev/null
+++ b/doc/source/analyzing/fields.rst
@@ -0,0 +1,19 @@
+Particle Fields
+====================================
+Naturally, particle fields contain properties of particles rather than
+grid cells. Many of these fields have corresponding grid fields that
+can be populated by "depositing" the particle values onto a yt grid.
+
+General Particle Fields
+------------------------------------
+Every particle will contain both a ``particle_position`` and ``particle_velocity``
+that tracks the position and velocity (respectively) in code units.
+
+
+SPH Fields
+------------------------------------
+For gas particles from SPH simulations, each particle will typically carry
+a field for the smoothing length `h`, which is roughly equivalent to
+`(m/\rho)^{1/3}`, where `m` and `rho` are the particle mass and density
+respectively. This can be useful for doing neighbour finding.
+
diff -r f244da87178a1bb4396b3e3959004ac38f14c6ac -r c44a150175c9fb3ee59ad73bda43af21dd2e6c21 doc/source/analyzing/objects.rst
--- a/doc/source/analyzing/objects.rst
+++ b/doc/source/analyzing/objects.rst
@@ -152,6 +152,11 @@
Combining Objects: Boolean Data Objects
---------------------------------------
+.. note:: Boolean Data Objects have not yet been ported to yt 3.0 from
+ yt 2.x. If you are interested in aiding in this port, please contact
+ the yt-dev mailing list. Until it is ported, this functionality below
+ will not work.
+
A special type of data object is the *boolean* data object.
It works only on three-dimensional objects.
It is built by relating already existing data objects with boolean operators.
diff -r f244da87178a1bb4396b3e3959004ac38f14c6ac -r c44a150175c9fb3ee59ad73bda43af21dd2e6c21 doc/source/examining/loading_data.rst
--- a/doc/source/examining/loading_data.rst
+++ b/doc/source/examining/loading_data.rst
@@ -156,7 +156,7 @@
yt has support for reading Gadget data in both raw binary and HDF5 formats. It
is able to access the particles as it would any other particle dataset, and it
can apply smoothing kernels to the data to produce both quantitative analysis
-and visualization.
+and visualization. See also the section :ref:`loading-sph-data`
Gadget data in HDF5 format can be loaded with the ``load`` command:
@@ -367,7 +367,8 @@
yt also supports loading Tipsy data. Many of its characteristics are similar
to how Gadget data is loaded; specifically, it shares its definition of
indexing and mesh-identification with that described in
-:ref:`particle-indexing-criteria`.
+:ref:`particle-indexing-criteria`. Like with gadget, see
+:ref:`loading-sph-data for more details`.
.. code-block:: python
@@ -903,3 +904,20 @@
---------------------
.. notebook:: Loading_Generic_Particle_Data.ipynb
+
+.. _loading_sph_data:
+
+SPH Particle Data
+-----------------
+For all of the SPH frontends, yt uses a cython-based SPH to created deposit
+mesh fields from individual particle fields. This uses a standard M4 smoothing
+kernel and the ``SmoothingLength`` field to calculate SPH sums, filling in the
+mesh fields. This gives you the ability to both track individual particles
+(useful for tasks like following contiguous clouds of gas that would be require
+a clump finder in grid data) as well as doing standard grid-based analysis.
+The ``SmoothingLength`` variable is also useful for determining which particles
+can interact with each other, since particles more distant than twice the
+smoothing length do not typically see each other in SPH simulations. By
+changing the value of the ``SmoothingLength`` and then re-depositing particles
+onto the grid, you can also effectively mimic what your data would look like at
+lower resolution.
diff -r f244da87178a1bb4396b3e3959004ac38f14c6ac -r c44a150175c9fb3ee59ad73bda43af21dd2e6c21 doc/source/visualizing/_cb_docstrings.inc
--- a/doc/source/visualizing/_cb_docstrings.inc
+++ b/doc/source/visualizing/_cb_docstrings.inc
@@ -388,23 +388,6 @@
ds = yt.load("IsolatedGalaxy/galaxy0030/galaxy0030")
p = yt.ProjectionPlot(ds, 'z', 'density', center='c', width=(20, 'kpc'))
p.annotate_title('Density Plot')
- s.save()
-
--------------
-
-.. function:: annotate_title(self, title='Plot'):
-
- (This is a proxy for :class:`~yt.visualization.plot_modifications.TitleCallback`.)
-
- Accepts a *title* and adds it to the plot
-
-.. python-script::
-
- from yt.mods import *
- ds = load("IsolatedGalaxy/galaxy0030/galaxy0030")
- p = ProjectionPlot(ds, 'z', 'density', center='c', width=(20, 'kpc'))
- p.annotate_title('Density plot')
->>>>>>> other
p.save()
Overplot quivers for the velocity field
diff -r f244da87178a1bb4396b3e3959004ac38f14c6ac -r c44a150175c9fb3ee59ad73bda43af21dd2e6c21 yt/analysis_modules/cosmological_observation/light_cone/light_cone.py
--- a/yt/analysis_modules/cosmological_observation/light_cone/light_cone.py
+++ b/yt/analysis_modules/cosmological_observation/light_cone/light_cone.py
@@ -389,169 +389,6 @@
attrs={"field_of_view": str(field_of_view),
"image_resolution": str(image_resolution)})
- def rerandomize_light_cone_solution(self, new_seed, recycle=True, filename=None):
- """
- When making a projection for a light cone, only randomizations along the
- line of sight make any given projection unique, since the lateral shifting
- and tiling is done after the projection is made. Therefore, multiple light
- cones can be made from a single set of projections by introducing different
- lateral random shifts and keeping all the original shifts along the line of
- sight. This routine will take in a new random seed and rerandomize the
- parts of the light cone that do not contribute to creating a unique
- projection object. Additionally, this routine is built such that if the
- same random seed is given for the rerandomizing, the solution will be
- identical to the original.
-
- This routine has now been updated to be a general solution rescrambler.
- If the keyword recycle is set to True, then it will recycle. Otherwise, it
- will create a completely new solution.
-
- new_sed : float
- The new random seed.
- recycle : bool
- If True, the new solution will have the same shift in the line of
- sight as the original solution. Since the projections of each
- slice are serialized and stored for the entire width of the box
- (even if the width used is left than the total box), the projection
- data can be deserialized instead of being remade from scratch.
- This can greatly speed up the creation of a large number of light
- cone projections.
- Default: True.
- filename : string
- If given, a text file detailing the solution will be written out.
- Default: None.
-
- """
-
- # Get rid of old halo mask, if one was there.
- self.halo_mask = []
-
- # Clean ds objects out of light cone solution.
- for my_slice in self.light_cone_solution:
- if my_slice.has_key('object'):
- del my_slice['object']
-
- if recycle:
- mylog.debug("Recycling solution made with %s with new seed %s." %
- (self.original_random_seed, new_seed))
- self.recycle_random_seed = int(new_seed)
- else:
- mylog.debug("Creating new solution with random seed %s." % new_seed)
- self.original_random_seed = int(new_seed)
- self.recycle_random_seed = 0
-
- self.recycle_solution = recycle
-
- # Keep track of fraction of volume in common between the original and
- # recycled solution.
- my_volume = 0.0
- total_volume = 0.0
-
- # For box coherence, keep track of effective depth travelled.
- box_fraction_used = 0.0
-
- # Seed random number generator with new seed.
- np.random.seed(int(new_seed))
-
- for q, output in enumerate(self.light_cone_solution):
- # It is necessary to make the same number of calls to the random
- # number generator so the original solution willbe produced if the
- # same seed is given.
-
- # Get projection axis and center.
- # If using box coherence, only get random axis and center if enough
- # of the box has been used, or if box_fraction_used will be greater
- # than 1 after this slice.
- if (q == 0) or (self.minimum_coherent_box_fraction == 0) or \
- (box_fraction_used > self.minimum_coherent_box_fraction) or \
- (box_fraction_used + self.light_cone_solution[q]['box_depth_fraction'] > 1.0):
- # Get random projection axis and center.
- # If recycling, axis will get thrown away since it is used in
- # creating a unique projection object.
- newAxis = np.random.randint(0, 3)
-
- newCenter = [np.random.random() for i in range(3)]
- box_fraction_used = 0.0
- else:
- # Same axis and center as previous slice, but with depth center shifted.
- newAxis = self.light_cone_solution[q-1]['projection_axis']
- newCenter = copy.deepcopy(self.light_cone_solution[q-1]['projection_center'])
- newCenter[newAxis] += \
- 0.5 * (self.light_cone_solution[q]['box_depth_fraction'] +
- self.light_cone_solution[q-1]['box_depth_fraction'])
- if newCenter[newAxis] >= 1.0:
- newCenter[newAxis] -= 1.0
-
- if recycle:
- output['projection_axis'] = self.master_solution[q]['projection_axis']
- else:
- output['projection_axis'] = newAxis
-
- box_fraction_used += self.light_cone_solution[q]['box_depth_fraction']
-
- # Make list of rectangle corners to calculate common volume.
- newCube = np.zeros(shape=(len(newCenter), 2))
- oldCube = np.zeros(shape=(len(newCenter), 2))
- for w in range(len(newCenter)):
- if (w == self.master_solution[q]['projection_axis']):
- oldCube[w] = [self.master_solution[q]['projection_center'][w] -
- 0.5 * self.master_solution[q]['box_depth_fraction'],
- self.master_solution[q]['projection_center'][w] +
- 0.5 * self.master_solution[q]['box_depth_fraction']]
- else:
- oldCube[w] = [self.master_solution[q]['projection_center'][w] -
- 0.5 * self.master_solution[q]['box_width_fraction'],
- self.master_solution[q]['projection_center'][w] +
- 0.5 * self.master_solution[q]['box_width_fraction']]
-
- if (w == output['projection_axis']):
- if recycle:
- newCube[w] = oldCube[w]
- else:
- newCube[w] = \
- [newCenter[w] -
- 0.5 * self.master_solution[q]['box_depth_fraction'],
- newCenter[w] +
- 0.5 * self.master_solution[q]['box_depth_fraction']]
- else:
- newCube[w] = [newCenter[w] -
- 0.5 * self.master_solution[q]['box_width_fraction'],
- newCenter[w] +
- 0.5 * self.master_solution[q]['box_width_fraction']]
-
- my_volume += common_volume(oldCube, newCube,
- periodic=np.array([[0, 1],
- [0, 1],
- [0, 1]]))
- total_volume += output['box_depth_fraction'] * \
- output['box_width_fraction']**2
-
- # Replace centers for every axis except the line of sight axis.
- for w in range(len(newCenter)):
- if not(recycle and
- (w == self.light_cone_solution[q]['projection_axis'])):
- self.light_cone_solution[q]['projection_center'][w] = \
- newCenter[w]
-
- if recycle:
- mylog.debug("Fractional common volume between master and recycled solution is %.2e" % \
- (my_volume / total_volume))
- else:
- mylog.debug("Fraction of total volume in common with old solution is %.2e." % \
- (my_volume / total_volume))
- self.master_solution = [copy.deepcopy(q) \
- for q in self.light_cone_solution]
-
- # Write solution to a file.
- if filename is not None:
- self._save_light_cone_solution(filename=filename)
-
- def restore_master_solution(self):
- "Reset the active light cone solution to the master solution."
- self.light_cone_solution = [copy.deepcopy(q) \
- for q in self.master_solution]
->>>>>>> other
-
@parallel_root_only
def _save_light_cone_solution(self, filename="light_cone.dat"):
"Write out a text file with information on light cone solution."
diff -r f244da87178a1bb4396b3e3959004ac38f14c6ac -r c44a150175c9fb3ee59ad73bda43af21dd2e6c21 yt/analysis_modules/halo_finding/halo_objects.py
--- a/yt/analysis_modules/halo_finding/halo_objects.py
+++ b/yt/analysis_modules/halo_finding/halo_objects.py
@@ -1813,11 +1813,12 @@
class GenericHaloFinder(HaloList, ParallelAnalysisInterface):
- def __init__(self, ds, ds, dm_only=True, padding=0.0):
+ def __init__(self, ds, data_source, dm_only=True, padding=0.0):
ParallelAnalysisInterface.__init__(self)
self.ds = ds
self.index = ds.index
- self.center = (np.array(ds.right_edge) + np.array(ds.left_edge)) / 2.0
+ self.center = (np.array(data_source.right_edge) +
+ np.array(data_source.left_edge)) / 2.0
def _parse_halolist(self, threshold_adjustment):
groups = []
diff -r f244da87178a1bb4396b3e3959004ac38f14c6ac -r c44a150175c9fb3ee59ad73bda43af21dd2e6c21 yt/data_objects/data_containers.py
--- a/yt/data_objects/data_containers.py
+++ b/yt/data_objects/data_containers.py
@@ -114,6 +114,10 @@
self.set_field_parameter(key, val)
@property
+ def pf(self):
+ return getattr(self, 'ds', None)
+
+ @property
def index(self):
if self._index is not None:
return self._index
diff -r f244da87178a1bb4396b3e3959004ac38f14c6ac -r c44a150175c9fb3ee59ad73bda43af21dd2e6c21 yt/utilities/fits_image.py
--- a/yt/utilities/fits_image.py
+++ b/yt/utilities/fits_image.py
@@ -13,6 +13,7 @@
import numpy as np
from yt.funcs import mylog, iterable, fix_axis, ensure_list
from yt.visualization.fixed_resolution import FixedResolutionBuffer
+from yt.visualization.plot_window import get_sanitized_center
from yt.data_objects.construction_data_containers import YTCoveringGridBase
from yt.utilities.on_demand_imports import _astropy
from yt.units.yt_array import YTQuantity
@@ -293,19 +294,20 @@
The axis of the slice. One of "x","y","z", or 0,1,2.
fields : string or list of strings
The fields to slice
- coord : float, tuple, or YTQuantity
- The coordinate of the slice along *axis*. Can be a (value,
- unit) tuple, a YTQuantity, or a float. If a float, it will be
- interpreted as in units of code_length.
+ center : A sequence floats, a string, or a tuple.
+ The coordinate of the center of the image. If set to 'c', 'center' or
+ left blank, the plot is centered on the middle of the domain. If set to
+ 'max' or 'm', the center will be located at the maximum of the
+ ('gas', 'density') field. Units can be specified by passing in center
+ as a tuple containing a coordinate and string unit name or by passing
+ in a YTArray. If a list or unitless array is supplied, code units are
+ assumed.
"""
- def __init__(self, ds, axis, fields, coord, **kwargs):
+ def __init__(self, ds, axis, fields, center="c", **kwargs):
fields = ensure_list(fields)
axis = fix_axis(axis, ds)
- if isinstance(coord, tuple):
- coord = ds.quan(coord[0], coord[1]).in_units("code_length").value
- elif isinstance(coord, YTQuantity):
- coord = coord.in_units("code_length").value
- slc = ds.slice(axis, coord, **kwargs)
+ center = get_sanitized_center(center, ds)
+ slc = ds.slice(axis, center[axis], **kwargs)
w, frb = construct_image(slc)
super(FITSSlice, self).__init__(frb, fields=fields, wcs=w)
for i, field in enumerate(fields):
diff -r f244da87178a1bb4396b3e3959004ac38f14c6ac -r c44a150175c9fb3ee59ad73bda43af21dd2e6c21 yt/utilities/grid_data_format/writer.py
--- a/yt/utilities/grid_data_format/writer.py
+++ b/yt/utilities/grid_data_format/writer.py
@@ -227,7 +227,7 @@
if unit_name in dataset_units:
value, units = dataset_units[unit_name]
else:
- attr = getattr(pf, unit_name)
+ attr = getattr(ds, unit_name)
value = float(attr)
units = str(attr.units)
d = g.create_dataset(unit_name, data=value)
diff -r f244da87178a1bb4396b3e3959004ac38f14c6ac -r c44a150175c9fb3ee59ad73bda43af21dd2e6c21 yt/visualization/plot_window.py
--- a/yt/visualization/plot_window.py
+++ b/yt/visualization/plot_window.py
@@ -807,10 +807,10 @@
pass
elif np.nanmax(image) == np.nanmin(image):
msg = "Plot image for field %s has zero dynamic " \
- "range. Min = Max = %d." % (f, np.nanmax(image))
+ "range. Min = Max = %f." % (f, np.nanmax(image))
elif np.nanmax(image) <= 0:
msg = "Plot image for field %s has no positive " \
- "values. Max = %d." % (f, np.nanmax(image))
+ "values. Max = %f." % (f, np.nanmax(image))
elif not np.any(np.isfinite(image)):
msg = "Plot image for field %s is filled with NaNs." % (f,)
if msg is not None:
diff -r f244da87178a1bb4396b3e3959004ac38f14c6ac -r c44a150175c9fb3ee59ad73bda43af21dd2e6c21 yt/visualization/profile_plotter.py
--- a/yt/visualization/profile_plotter.py
+++ b/yt/visualization/profile_plotter.py
@@ -797,19 +797,19 @@
positive_values = data[data > 0.0]
if len(positive_values) == 0:
mylog.warning("Profiled field %s has no positive "
- "values. Max = %d." %
+ "values. Max = %f." %
(f, np.nanmax(data)))
mylog.warning("Switching to linear colorbar scaling.")
- zmin = data.min()
+ zmin = np.nanmin(data)
z_scale = 'linear'
self._field_transform[f] = linear_transform
else:
zmin = positive_values.min()
self._field_transform[f] = log_transform
else:
- zmin = data.min()
+ zmin = np.nanmin(data)
self._field_transform[f] = linear_transform
- zlim = [zmin, data.max()]
+ zlim = [zmin, np.nanmax(data)]
fp = self._font_properties
f = self.profile.data_source._determine_fields(f)[0]
diff -r f244da87178a1bb4396b3e3959004ac38f14c6ac -r c44a150175c9fb3ee59ad73bda43af21dd2e6c21 yt/visualization/volume_rendering/transfer_function_helper.py
--- a/yt/visualization/volume_rendering/transfer_function_helper.py
+++ b/yt/visualization/volume_rendering/transfer_function_helper.py
@@ -19,7 +19,6 @@
from yt.visualization.volume_rendering.api import ColorTransferFunction
from yt.visualization._mpl_imports import FigureCanvasAgg
from matplotlib.figure import Figure
-from IPython.core.display import Image
from yt.extern.six.moves import StringIO
import numpy as np
@@ -190,6 +189,7 @@
ax.set_ylim(y.max()*1.0e-3, y.max()*2)
if fn is None:
+ from IPython.core.display import Image
f = StringIO()
canvas.print_figure(f)
f.seek(0)
Repository URL: https://bitbucket.org/yt_analysis/yt/
--
This is a commit notification from bitbucket.org. You are receiving
this because you have the service enabled, addressing the recipient of
this email.
More information about the yt-svn
mailing list