yt-users October 2019

yt-users@python.org

13 participants
14 discussions

by david collins

Hi, Everybody! 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. Thanks, d. -- Sent from my computer.

3 months

mass-weighted average magnetic field (QUIVERS) projected on to the plane plot

by nirmit sakre

I would like to annotate Mass weighted Magnetic field quivers on density projection plot so that i can display the magnetic field features. Do annotate_magnetic_field() , annotate_quiver() work correctly on projection plot or do they work only on slice plot correctly. If they work on projection plot then against which quantity magnetic field quivers are weighed or not weighted?

6 months, 1 week

Ytini installation for Houdini on Windows

by uran san

Hi guys, 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 OS. 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 Ur4n

7 months

How to pass value of the field corresponds to the index to the scalar function!!!

by Prateek Gupta

Dear Yt-users, I need to pass the value of the available/derived field corresponds to the index selected based on some criterion (using numpy.where()) to a scalar function. It should list out the indexes wheres the condition satisfies. It happens when this has been done outside the field but instead of getting 3D indexes it gives 1D YT_Array indexes I think. But I need to pass this value of field corresponds to the selected index to scalar function inside the field, which may give some problem (please find the attachment to see what exactly I am trying to do). The attached script can run on the ENZO test data set, enzo_tiny_cosmology. Many thanks and regards -Prateek

8 months, 1 week

Save the date: yt user/developer workshop in Edinburgh Summer 2020

by Britton Smith

Hi everyone, What: yt user/developer workshop Where: Edinburgh, UK When: June 29 to July 3, 2020 I'm very pleased to announce that there will be a yt user/developer workshop at the Higgs Centre for Theoretical Physics at University of Edinburgh from June 29 to July 3, 2020. The workshop will begin with a couple days of tutorials for new users and then transition into development activities. This will be a good opportunity to meet and join the yt community. More information to follow, but mark your calendars and diaries now. I hope to see you next summer! Britton

8 months, 2 weeks

Update on yt governance documents/structure

by Madicken Munk

Hi yt-users! 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: https://github.com/yt-project/governance/pull/1 and https://github.com/yt-project/ytep/pull/7 (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 removed) Our existing governance document is located at https://ytep.readthedocs.io/en/master/YTEPs/YTEP-1776.html 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. -Madicken

8 months, 3 weeks

Saving field data

by Edward Buie II

Hello hope you are doing well. I am trying to save some projection and slice field data within a sphere by doing the following : my_sphere = stuff.sphere([0.5,0.5,0.5],(214,"kpc")) prj1 = yt.SlicePlot(stuff,'z',('gas','radius'),data_source=my_sphere) rad = prj1.frb['gas','radius'] prj1.save() I then save rad to a pickle file to read in by another code. The full simulation domain is 800 kpc. When I load the pickle file and print out the values I am getting many 0's. And this is also the same when I make a projection of some field and save it in the same way to a pickle file. My question is, even though I am limiting the data to within a sphere, does rad also hold values that are not within that sphere and saving them as 0's? Is there a better way to go about this? Thank you. -- Edward N. Buie II, *Ph.D. Candidate* Astrophysics School of Earth and Space Exploration (SESE) Arizona State University

8 months, 3 weeks

create a derived field from a multi-dimensional array

by Yang Luo

Hi I would like to create a new filed, but I encounter one problem here. The dataset is in Arepo format, and there is one field ('PartType0', 'GFM_Metals'), the dimension of this field is (2716294, 10). For normal field, like ('PartType0', 'Masses'), the dimension is (2716294, ). So field ('PartType0', 'GFM_Metals') has more dimensions. Below is how I create the new field. As I only need one row from ('PartType0', 'GFM_Metals') , so that: def _nOxy(field, data): X_Oxygen = data['PartType0', 'GFM_Metals'][:,4] rho = data['PartType0', 'Density'] nOxy = X_Oxygen *rho/ ds.quan(mh,'g') return nOxy ds.add_field(('PartType0', 'nOxy'), function=_nOxy, units='cm**-3',particle_type=True) However, it shows error: IndexError: too many indices for array The error is from data['PartType0', 'GFM_Metals'][:,4], where yt does not support using [:,4]. So I am wondering if there is any other method to handle this. I am grateful for your suggestions. Thanks, Yang

8 months, 3 weeks

upcoming papers on yt performance

by Luigi Iapichino

Dear yt users and developers, we would like to draw your attention on two upcoming papers of ours on performance studies of yt. The first one will be on ArXiv on Friday ("Speeding simulation analysis up with yt and Intel Distribution for Python", authors: S. Cielo, L. Iapichino and F. Baruffa) and shows the performance gained when installing yt on top of a Conda environment based on the Intel Distribution for Python, rather than on the traditional Anaconda-based environment. The second paper is more general in scope (code optimisation of astrophysical codes on Intel Xeon Phi Knights Landing) but has a whole section on tests of different parallelisation strategies in yt (including Cython and OpenMP). This paper is still under review and we share a semi-final version with you under this link: https://syncandshare.lrz.de/dl/fiFjBeCCh63DWACCcf3cTjdm/?inline The aim of both works is to explore the potential offered by the software stack of modern HPC systems for applications like yt. Besides a custom installation procedure, no modification of the source code was required. We see these works as potential contributions to the yt user community, and therefore we would welcome comments, discussions and follow-up ideas about them. For those of you attending SC19 in November, one of us (Salvatore) will be there and is available for face-to-face meetings. Looking forward to your feedback, Best regards, Luigi and Salvatore -- ----------------------------------------- Dr. Luigi Iapichino Leibniz-Rechenzentrum High Performance Systems Division Boltzmannstr. 1, 85748 Garching, Germany Tel: +49 89 35831 7824 e-mail: luigi.iapichino(a)lrz.de www: http://iapichino.userweb.mwn.de/

9 months

extracting data from LinePlot or SlicePlot

by Tseng, Po-Hsun

Hi yt users, Is any ways to extract data from LinePlot or SlicePlot and put data into numpy array? thanks.

9 months

2020

2019

2018

2017

2016

2015

2014

2013

2012

2011

2010

2009

2008

2007

yt-users October 2019