Does anyone out there have a technique for getting the variance out of
a profile object? A profile object is good at getting <X> vs. B, I'd
then like to get < (X - <X>)^2 > vs B. Matt and I had spittballed the
possibility some time ago, but I was wondering if anyone out there had
successfully done it.
Sent from my computer.
I am very new to Ytini's mailing list, mailing lists as well as Houdini.
So please bear with me while I pick up the etiquette around this medium.
With that out of the way here's my first question:
After having succeeded installing Ytini for Houdini on a Mac OS X 10.14 I
would now like to achieve the same with my Houdini install on a Windows 10
I have tried to somewhat reproduce all steps described here
http://www.ytini.com/getstarted.html#collapseTwoOut (Houdini 17.5 Tested on
Mac OS X) on my Windows machine but wasn't surprised to learn that of
course Houdini could not find the ytini module when I tried to import it in
the python shell.
Is there somebody, that could point me . in the right direction of how to
install Ytini for Houdini on a Windows 10 OS environment?
Not afraid of reading so even hints a re more than welcome.
Many Thanks in advance
We've been working on updating our governance model for the yt-project to
help our project grow, to codify our decision-making process, and to
outline the governance structure of the project. I've submitted two
companion pull requests with drafts of this updated documentation:
(you can see the contents of the pull requests if you click on the "files
changed" tab; text highlighted in green is an addition and in red is
Our existing governance document is located at
The first pull request is to a new governance repository, where our
governance documentation will live separately from the YTEPs (yt
enhancement proposals). This will allow us to maintain and minorly update
our governance without having to update YTEP. The second pull request is to
the YTEP repository and generally outlines the core values and ideas we
want our governance structure to reflect. Hopefully all of the things I've
listed in the YTEP are reflected in the governance docs.
As members of the community, I'd like to solicit feedback from all of you
about these governance documents. Do these reflect our community values?
Should we add anything? Do you feel everything is clear? Is this too much
governance for our community right now? Is there something that's missing?
As users, do you feel our governance structure is supportive to you, or do
you think there is a way we could improve it to do so?
Feel free to reply here or comment on the pull requests! Our governance
will be better with your feedback.
Hi yt users!
I'm not sure if this is the place to ask this but I was following
The example to export to RADMC3D
Number of processors: 4
Number of threads in use: 1
Reading global frequencies/wavelengths...
Reading grid file and prepare grid tree...
At line 3196 of file ./amr_module.f90 (unit = 1, file = 'amr_grid.inp')
Fortran runtime error: Bad integer for item 5 in list input
I've run examples from radmc3dpy webpage and those worked. I'm not sure if
the problem is the format in which yt writes the amr input file or radmc3d
I am analyzing cosmological simulations with yt and I am trying to filter data by the radial velocity of a sphere with a defined center. I set the bulk velocity of the sphere using
v = sphere.quantities.bulk_velocity()
However, the numbers that I am having to input in order to filter out the data are extremely large even after setting the bulk velocity. Then I calculated the mean velocity of the sphere:
vx = sphere.mean('velocity_x',weight='density')
vy = sphere.mean('velocity_y',weight='density')
vz = sphere.mean('velocity_z',weight='density')
but even after setting the bulk velocity, the mean velocity of the sphere is not even close to being zero. The exact output is:
The bulk velocity is:
[ -3629508.65339038 -15948968.49831905 6561226.69040721] cm/s
The mean velocity after setting bulk velocity is:
[ -3495226.24285729 -16714324.95702405 5505530.59033089] cm/s
I'm curious if I'm doing something simple wrong or if there is a bigger issue. Any help would be much appreciated!
I ran both HOP and Rockstar halo finders on the same ART simulations,
and the mass functions match very well (attached hmf_m.pdf file).
However, when I compare virial radii (attached hmf_r.pdf file), they
don't match. Even if I try rescaling the Rockstar virial radius by 1.5,
the match is much worse than for mass functions.
How are the 'virial_radius' fields defined for both halo finders?
This is a z=7 simulations, so all mass definitions are the same, Mvir =
M200m = M200c.
I also attach the script to plot virial radii if you want to check it,
although it is trivial.
Many thanks (as usual),
I am currently attempting to compare a catalog with Mvir and Rvir units of
Mpc/h and Kpc/h comoving to a catalog with Mpc and Kpc proper. I am
currently using yt version 3.5.1, which I just checked for updates.
When doing this conversion of units within yt, there appears to be a
problem where an amount of h*h is not handeled properly. I have included a
screenshot of a jupyter session that outlines the problem as well as a .py
and jupyter script so it can be tested.
The problem arises when I try to calculate delta_vir, which is the density
of a galaxy/ critical density. The correct value for both catalogs I am
comparing is around 170. When I use yt to calculate the critical density at
z = 2, it get the correct value, which I convert to Mpc*h*h//kpc/kpc/kpc to
match with my /h catalog. Then, when I divide the densities, even thought
their units are identical, their division is different than if just the
float values are divided, by an exact factor of h*h. I wonder if there is
something wrong with how I am initializing my values with the ytquantity,
or if there is something wrong going on in the code.
[image: Screen Shot 2019-07-31 at 2.36.01 PM.png]
I'm running a halo finder using fof method on some arbitrary dataset which contains positions of all particles. The halo finder shows the number of groups it has found. but when saving halo data an error came up. I have debugged the code and see positions of all halos is set to zero and virial mass and radius is set to -1. I wonder lack of what thing and particle fields leads to this?