yt-users November 2014

yt-users@python.org

40 participants
41 discussions

Stars on volume rendering
by Elizabeth Harper-Clark 11 Dec '21

11 Dec '21

Hi all, Thanks for all your help over the last couple of days. One more question: - Can I plot particles on a volume rendered image? I have stars and I want to show where they are! Thanks, Libby -- Elizabeth Harper-Clark MA MSci PhD Candidate, Canadian Institute for Theoretical Astrophysics, UofT Sciences and Engineering Coordinator, Teaching Assistants' Training Program, UofT www.astro.utoronto.ca/~h-clark <http://www.astro.utoronto.ca/%7Eh-clark> h-clark(a)cita.utoronto.ca Astronomy office phone: +1-416-978-5759

8 14

Variance in phase/profiles
by david collins 11 Apr '20

11 Apr '20

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.

6 8

yt Community Survey
by Cameron Hummels 08 Dec '14

08 Dec '14

Fellow yt users: We at the yt development team want to make yt better for you, the user. We have prepared a short survey to gain your valuable feedback to help us understand what is working in yt and what is not. Could you please spare 2 minutes to fill out our survey below? We plan to collect responses and make available a selection of results by December 10, so please let us know what you think at the following URL: http://goo.gl/forms/hRNryOWTPO On behalf of the yt development team, Cameron -- Cameron Hummels Postdoctoral Researcher Steward Observatory University of Arizona http://chummels.org

2 2

domain-decomposition for volume rendering huge brick?
by Stuart Levy 03 Dec '14

03 Dec '14

Hello all, We're hoping to use yt parallel volume rendering on a very large generic brick - it's a simple rectangular unigrid slab, but containing something like 1.5e11 points, so much too large for load_uniform_grid() to load into memory in a single machine. I imagine it wouldn't be hard to do the domain decomposition by hand, loading a different chunk of grid into each MPI process. But then what? What would it take to invoke the volume renderer on each piece and composite them together? Would it help if the chunks were stored in a KDTree? Is there some example (one of the existing data loaders?) which I could follow?

4 10

Reading athena output
by almog yalin 02 Dec '14

02 Dec '14

Hi. I ran a 2D simulation with Athena and exported the data to vtk files. I was even able to load the file with yt and inspect some of its elements. What I would like to do is to compare the density from the simulation to an analytic solution (which I can evaluate at the grid points). The problem is I cannot figure out how to extract the relevant information from the yt.frontends.athena.data_structures.AthenaDataset object. What I expect are 2d arrays of the x coordinate, y coordinate and density. Thanks in Advance, Almog

3 2

Performance when using derived fields
by Daniel Fenn 02 Dec '14

02 Dec '14

Hi everyone. Is it normal for slices to take substantially longer to be created if they're made from a derived field rather than an intrinsic one? Specifically, I'm having an issue creating an axis-aligned slice using the divergence of the velocity field. It's taking around 6.5 minutes just to make the slice, whereas if I use temperature or density, it takes around 10 seconds or so for the same dataset. I also notice that the amount of time it takes is not dependent on the number of processors I'm using. I've used 1, 12, and 24 processors, with identical results, even though I'm calling enable_parallelism(), and I can see that all the processes are running. I read in the docs that slice operations aren't generally done in parallel, but in this case it seems that maybe it would be beneficial. A similar operation in VisIt completes much faster, so I'm wondering if I've misconfigured something, or if there is something I can do to speed things up. I'd appreciate any thoughts anyone has on the subject. Thanks, Dan

4 10

Volume Rendering and Field lines
by Dionysopoulou K. 01 Dec '14

01 Dec '14

Dear all, I would like to ask if there has been any progress on the volume rendering refactoring since late October. In particular, I would like to know if volume rendering and the new camera implementation allows for plotting field lines together with the rendering. If so, are there any examples on the web? Unfortunately, matplotlib contour plots with field lines are not good enough for me. Thanks in advance, Kiki Kyriaki Dionysopoulou ======================================================= Mathematical Sciences University of Southampton Southampton, SO17 1BJ, UK K.Dionysopoulou(a)soton.ac.uk<mailto:K.Dionysopoulou@soton.ac.uk>

2 1

Handling large files
by Adam Jacobs 30 Nov '14

30 Nov '14

Is it unreasonable for me to expect yt on a workstation with 4x3GHz cores and ~20 GB of RAM to be able to handle a ~100 GB dataset? I'm trying to select out a subset of the dataset using cut_region(), but still seem to run into hanging or running out of RAM. For example, when I try to do a write_out() on the cut region yt sucks up all 20 GB of RAM and I have to kill it. Is there a preferred method for loading in a subset of data? I don't need near the full 100 GB of data. I'm using yt dev and looking at BoxLib/Maestro data. -- Adam Jacobs Department of Physics and Astronomy, PhD Candidate Stony Brook University http://astro.sunysb.edu/amjacobs/

2 2

Is it possible to get grid size from an enzo AMR data?
by Reju Sam John 29 Nov '14

29 Nov '14

Dear all, Is it possible to take grid size from an enzo AMR data? Thanks -- Reju Sam John

2 6

Radial Velocity Projection
by John Regan 25 Nov '14

25 Nov '14

Hi All, I need to set the "center" field parameter for making a projection of the "radial_velocity". I see from the docs how to do that for a data container but its not clear what to do for a projection. The error I am getting at the moment when the "center" is not set is: File "/homeappl/home/regan/appl_taito/YT/Dev-3.0/yt/yt/data_objects/data_containers.py", line 249, in __getitem__ self.get_data(f) File "/homeappl/home/regan/appl_taito/YT/Dev-3.0/yt/yt/data_objects/data_containers.py", line 656, in get_data finfo.check_available(self) File "/homeappl/home/regan/appl_taito/YT/Dev-3.0/yt/yt/fields/derived_field.py", line 146, in check_available validator(data) File "/homeappl/home/regan/appl_taito/YT/Dev-3.0/yt/yt/fields/derived_field.py", line 235, in __call__ raise NeedsParameter(doesnt_have) yt.fields.field_exceptions.NeedsParameter: (['center']) Any tips much appreciated! Cheers, John

2 1

2022

2021

2020

2019

2018

2017

2016

2015

2014

2013

2012

2011

2010

2009

2008

2007

yt-users November 2014