Hi folks,
I've stumbled across across a research task that could really use yt's
Boolean data object feature, and according to the documentation it is not
implemented in yt-3. Have boolean data objects been replaced by something
else? I know there are various filtering methods in yt, but what I really
need to do is spatial filtering (i.e., data in a sphere with some chunks of
the sphere excised).
In principle I'd be happy to take a crack at re-implementing it, although I
might need some help. Could somebody give me a sense of why it wasn't
ported to yt-3 when the changeover was made?
Thanks,
Brian
New issue 1237: HaloMassFcn units (or coments) inconsistent
https://bitbucket.org/yt_analysis/yt/issues/1237/halomassfcn-units-or-com...
spthm:
Disclaimer: I'm not particularly experienced in HMFs. The below may not actually be a bug (just confusing/incorrect comments); if it is, I imagine it relates to issue #1228. It should possibly be a comment on that issue, but I can't confirm the two are related.
For reference, I'm going to link to the current stable yt, 3.2.3.
Based on comments in [`halo_mass_function.py`](https://bitbucket.org/yt_analysis/yt/src/bdea84d950999973d7b80cba34acc4c70b62a925/yt/analysis_modules/halo_mass_function/halo_mass_function.py?at=stable&fileviewer=file-view-default#halo_mass_function.py), I do not believe the stated output units are correct.
[`halo_mass_function.py:403`](https://bitbucket.org/yt_analysis/yt/src/bdea84d950999973d7b80cba34acc4c70b62a925/yt/analysis_modules/halo_mass_function/halo_mass_function.py?at=stable&fileviewer=file-view-default#halo_mass_function.py-403) multiplies `rho_crit_g_cm3_h2` by `hubble0**2`, resulting in units of `Msolar/Mpc^3` for `rho0` _if the comment at [`physical_ratios.py:103`](https://bitbucket.org/yt_analysis/yt/src/bdea84d950999973d7b80cba34acc4c70b62a925/yt/utilities/physical_ratios.py?at=yt&fileviewer=file-view-default#physical_ratios.py-103) is to be believed_. But line 402 directly above states `rho0 in units of h^2 Msolar/Mpc^3`.
This error carries through to line [`halo_mass_function.py:432`](https://bitbucket.org/yt_analysis/yt/src/bdea84d950999973d7b80cba34acc4c70b62a925/yt/analysis_modules/halo_mass_function/halo_mass_function.py?at=stable&fileviewer=file-view-default#halo_mass_function.py-432). It is claimed that `dndM_dM_analytic` has units of `h^3`, but it is actually `h`.
Finally, at [`halo_mass_function.py:436`](https://bitbucket.org/yt_analysis/yt/src/bdea84d950999973d7b80cba34acc4c70b62a925/yt/analysis_modules/halo_mass_function/halo_mass_function.py?at=stable&fileviewer=file-view-default#halo_mass_function.py-436) we multiply through by `h^3`, resulting units of `h^-2` for `dndM_dM_analytic` and `n_cumulative_analytic` (while it is claimed that all `h` are cancelled out).
Note that [`halo_mass_function.py:368`](https://bitbucket.org/yt_analysis/yt/src/bdea84d950999973d7b80cba34acc4c70b62a925/yt/analysis_modules/halo_mass_function/halo_mass_function.py?at=stable&fileviewer=file-view-default#halo_mass_function.py-368) has a similar issue, in that `h^2` units are claimed, but `rho0` has no `h` dependence. The proceeding usage of `rho0` at [`halo_mass_function.py:390`](https://bitbucket.org/yt_analysis/yt/src/bdea84d950999973d7b80cba34acc4c70b62a925/yt/analysis_modules/halo_mass_function/halo_mass_function.py?at=stable&fileviewer=file-view-default#halo_mass_function.py-390) implies that `rho0` indeed should have units `h^2` to achieve a `h^-1` for `thisradius`. I think this has a non-trivial impact on the values in `sigmaarray`, but I gave up trying to track it.
The masses do not appear to be affected, but both analytic output arrays seem to have erroneous factors of `h`, and I believe the sigma values are also computed with incorrect `R` values. Since I gave up trying to track how incorrect `R` impacts `sigma`, I don't know if this all agrees with the behaviour described by @ioodderskov in #1228.
All that said, if I make the obvious change of _not_ multiplying `rho` by `hubble0**2`, the analytic fits do a much worse job of fitting haloes in some of my Gadget-2 snapshots, whereas before agreement with P-S was very good. But I have no idea if this is a fluke, and the simulations are faulty, or whether all of the above is a misunderstanding on my part.
The multiplication by `hubble0**2` was introduced in ab0d4e0.
New issue 1236: While installing yt through script
https://bitbucket.org/yt_analysis/yt/issues/1236/while-installing-yt-thro...
Damodar Rajbhandari:
I have tried to install yt for the first time in my Ubuntu 16.04 LTS but i got an error "Import Error: cannot import name HTTPS Handler". Here are some lines from the yt_install.log
.....
using https://bitbucket.org/yt_analysis/yt-supplemental/
sending capabilities command
abort: No module named _ssl!
using https://bitbucket.org/yt_analysis/yt/
sending capabilities command
abort: No module named _ssl!
abort: repository /home/field/yt-x86_64/src/yt-hg/ not found!
2016-06-26 19:49:49 URL:https://bootstrap.pypa.io/get-pip.py [1524722/1524722] -> "get-pip.py" [1]
Traceback (most recent call last):
File "get-pip.py", line 19177, in <module>
main()
File "get-pip.py", line 194, in main
bootstrap(tmpdir=tmpdir)
File "get-pip.py", line 82, in bootstrap
import pip
File "/tmp/tmpzwM52k/pip.zip/pip/__init__.py", line 16, in <module>
File "/tmp/tmpzwM52k/pip.zip/pip/vcs/subversion.py", line 9, in <module>
File "/tmp/tmpzwM52k/pip.zip/pip/index.py", line 30, in <module>
File "/tmp/tmpzwM52k/pip.zip/pip/wheel.py", line 39, in <module>
File "/tmp/tmpzwM52k/pip.zip/pip/_vendor/distlib/scripts.py", line 14, in <module>
File "/tmp/tmpzwM52k/pip.zip/pip/_vendor/distlib/compat.py", line 31, in <module>
ImportError: cannot import name HTTPSHandler
using https://bitbucket.org/yt_analysis/yt-supplemental/
sending capabilities command
abort: No module named _ssl!
using https://bitbucket.org/yt_analysis/yt/
sending capabilities command
abort: No module named _ssl!
abort: repository /home/field/yt-x86_64/src/yt-hg/ not found!
2016-06-26 20:05:41 URL:https://bootstrap.pypa.io/get-pip.py [1524722/1524722] -> "get-pip.py.1" [1]
Traceback (most recent call last):
File "get-pip.py", line 19177, in <module>
main()
File "get-pip.py", line 194, in main
bootstrap(tmpdir=tmpdir)
File "get-pip.py", line 82, in bootstrap
import pip
File "/tmp/tmpJkD7Nt/pip.zip/pip/__init__.py", line 16, in <module>
File "/tmp/tmpJkD7Nt/pip.zip/pip/vcs/subversion.py", line 9, in <module>
File "/tmp/tmpJkD7Nt/pip.zip/pip/index.py", line 30, in <module>
File "/tmp/tmpJkD7Nt/pip.zip/pip/wheel.py", line 39, in <module>
File "/tmp/tmpJkD7Nt/pip.zip/pip/_vendor/distlib/scripts.py", line 14, in <module>
File "/tmp/tmpJkD7Nt/pip.zip/pip/_vendor/distlib/compat.py", line 31, in <module>
ImportError: cannot import name HTTPSHandler
Then, i try to install by typing sudo apt-get install python-yt. But i don't know what was going there!
So, again i tried, its showed me an error "make: *** [install-recursive] Error 1". Here are some lines from another yt_install.log
......
In file included from h5tools_dump.c:31:0:
../../src/H5private.h:1529:14: warning: ISO C90 does not support '__func__' predefined identifier [-Wpedantic]
#define FUNC __func__
^
h5tools_error.h:71:49: note: in expansion of macro 'FUNC'
H5Epush2(H5tools_ERR_STACK_g, __FILE__, FUNC, __LINE__, H5tools_ERR_CLS_g, maj_id, min_id, str); \
^
h5tools_error.h:101:5: note: in expansion of macro 'HERROR'
HERROR(H5E_tools_g, min_id, str); \
^
h5tools_dump.c:2688:9: note: in expansion of macro 'H5E_THROW'
H5E_THROW(FAIL, H5E_tools_min_id_g, "Could not close datatype's super class");
^
h5tools_dump.c: In function 'h5tools_dump_comment':
h5tools_dump.c:3460:59: warning: conversion to 'size_t {aka long unsigned int}' from 'ssize_t {aka long int}' may change the sign of the result [-Wsign-conversion]
cmt_bufsize = H5Oget_comment(obj_id, comment, cmt_bufsize);
^
h5tools_dump.c: In function 'h5tools_dump_attribute':
h5tools_dump.c:3489:39: warning: unused parameter 'oid' [-Wunused-parameter]
h5tools_context_t *ctx, hid_t oid, const char *attr_name, hid_t attr_id,
^
h5tools_dump.c: In function 'h5tools_print_packed_bits':
h5tools_dump.c:3628:48: warning: ISO C90 does not support 'long long' [-Wlong-long]
packed_bits_size = 8 * sizeof(long long);
^
h5tools_dump.c:3631:57: warning: ISO C90 does not support 'long long' [-Wlong-long]
packed_bits_size = 8 * sizeof(unsigned long long);
^
Makefile:553: recipe for target 'h5tools_dump.lo' failed
make[2]: *** [h5tools_dump.lo] Error 1
make[2]: Leaving directory '/home/field/yt-x86_64/src/hdf5-1.8.11/tools/lib'
Makefile:497: recipe for target 'install-recursive' failed
make[1]: *** [install-recursive] Error 1
make[1]: Leaving directory '/home/field/yt-x86_64/src/hdf5-1.8.11/tools'
Makefile:539: recipe for target 'install-recursive' failed
make: *** [install-recursive] Error 1
Can anyone can help me out from this issue.
Hi all,
This morning around 7am fido - our server located @ NCSA - was
compromised and used in DDoS attack on some random ip in China. As soon
as the secops noticed unusual behaviour they shut it down.
Fido serves as reverse proxy for various services related to yt.
As long as it's down anything that we were hosting at *.hub.yt and
*.use.yt won't be available.
This also affects: blog.yt-project.org and tests.yt-project.org.
At this time I don't have any ETA on when the services will be brought
back to life. I'll notify you as soon as there's any progress.
I'm sorry for all the inconvenience this situation creates.
Kacper
Hi all,
Does anyone have any pretty visuals (be it images, plots, or movies) that
they wouldn't mind me showing off at a yt tutorial? Of course you will be
duly credited and praised for your awesome visuals :)
I'll be giving the tutorial at McMaster this weekend as part of a python
for astronomers workshop before the Great Lakes Cosmology workshop.
Thanks!
-Nathan
New issue 1232: Absorption Spectrum Error (Gizmo file type)?
https://bitbucket.org/yt_analysis/yt/issues/1232/absorption-spectrum-erro...
Anonymous:
I tried using the cookbook code for creating an absorption spectrum (http://yt-project.org/docs/dev/cookbook/cosmological_analysis.html#creati...) for a Gizmo type cosmological simulation file. I created a light ray successfully, create and absorption spectrum via sp = AbsorptionSpectrum(900.0, 1400.0, 50000), and add a HI line to
the spectrum, following the cookbook. But whenever I get to the line:
wavelength, flux = sp.make_spectrum('lightray.h5',
output_file='spectrum.h5',
line_list_file='lines.txt',
use_peculiar_velocity=True)
I get the following error:
absorption_spectrum.pyc in _add_lines_to_spectrum(self, field_data, use_peculiar_velocity, output_absorbers_file, subgrid_resolution, observing_redshift, njobs):
444 # edges (ie the wings), then widen the wavelength
445 # window and repeat process.
--> 446 if (vtau[0] < min_tau and vtau[-1] < min_tau):
447 break
448 window_width_in_bins *= 2
IndexError: index 0 is out of bounds for axis 0 with size 0
I've looked at the source code but I can't figure out what's causing the problem. If I do not try to add any lines to my absorption spectrum at all, then it doesn't error. But since I'm following the cookbook code, I'm not sure what could be going wrong...
Help?
Thanks,
Kyle Stewart
