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.
Dear all,
I recently start to work with yt to analyze my FLASH simulation. The output of those simulation is an hdf5 file.
When I try to plot a slice of my simulation, some weird lines appear on it. Here is an example of the python script and the results it leads to.
"""
import yt
ds = yt.load("lasslab_hdf5_plt_cnt_0000.hdf5")
yt.SlicePlot(ds, 'z', "dens").show()
[cid:972219b9-edbb-4265-81cf-9a07e61cf360]
"""
Actually they also appear went plotting the results of a frb of a slice (even when I try to plot the x or y coordinate in 2d).
Do you have any idea of the reason for such lines, or a solution for my problem?
Best regards,
Gabriel Rigon
Hi all,
For some work, I had to use Millenium Simulation data. I have a catalog of
halos with their positions, and I wanted to get the halos in voids. Can yt
be of any use for such an analysis? The thing that makes it difficult is
periodic boundary conditions it doesn't allow you to use KDtree straight
away.
Please help,
Thanks,
With Regards,
Ankit
Hello yt-commuinty,
I was trying to use yt to analyse GIZMO data.
Is there a function which can translate GIZMO gas data,particle0, to
uniform grid with neutral hydrogen,velocity,temperature.etc?
I know the "deposit " field and covering_grid in yt, but "deposit" isn't
uniform and covering_grid is designed for grid-like data..
Thank you!
Edward
Hi all,
I'm experiencing a little trouble in reading some RAMSES files which
include radiative transfer.
In this (broadly used) implementation some of the fields do not follow the
order they have classical RAMSES files; e.g. the pressure field is later in
the fields list. Indeed when I read the files those fields do not have the
correct content.
Can someone please suggest a possible workaround? For instance a plugin for
this kind of files, if it exists, or suggest how I could swap the fields,
or rename them with the correct name?
I do not want simply to create new fields, without correcting the default
ones, as all derived fields (and possibly other valuable yt features) would
not work properly.
Thank you very much!
Best,
Salvatore
Hello Dear Users,
I'm recently trying to recreate some of the interesting 2-D radial
plots of stellar and gas surface density the ALMA folks are making (see
https://arxiv.org/abs/1703.07029 for instance) using my simulations. I've
been experimenting with using an FRB to do this, but I'm having trouble
resolving the innermost cells even though the FRB resolution matches my
grid's highest resolution.
I see that there was some discussion of including this in the code a
couple of years ago (as well as where I got the idea to try FRBs for this):
http://lists.spacepope.org/pipermail/yt-users-spacepope.org/2015-June/00620…
Did this ever make it in? I'd love to be able to use it if so.
Cordially,
Josh
--
Joshua Wall
Doctoral Candidate
Department of Physics
Drexel University
3141 Chestnut Street
Philadelphia, PA 19104
Hi People,
I was wondering if yt can produce smooth images of zoomed-in sliceplots,
where the pixelated AMR grids are not clearly visible. I was interested to
see the colors get mixed in smoothly (something like tetrahedralization),
is that possible anyway in yt?
Best
Tazkera
Hello guys,
our research group would like to visualize volume data, where the
transfer function is bivariate (see
http://lists.spacepope.org/pipermail/yt-users-spacepope.org/2017-April/0085…
). We have been struggling with this task, so I would like you to help
me understand some of the basic concepts I have come across and that I
find confusing.
1. I think that I somehow understand the concept of transfer functions
with multiple fields, but I have difficulties of creating a volume
source with multiple fields. I have somehow examined the source code,
and it seems to me that only the first field is used (judging by
https://github.com/yt-project/yt/blob/master/yt/visualization/volume_render…
). If I am wrong and it is meaningful and possible to construct a source
with multiple fields, could you please provide an example?
2. I have no doubts that transfer functions work, but I was unable to
find out how are they involved in the integration process. What I came
across were the integration routines and they have mentioned that the
transfer function is not used:
https://github.com/yt-project/yt/blob/master/yt/utilities/lib/ray_integrato…
. This is really interesting, how is it in reality?
3. You have mentioned the PlanckTransferFunction as a good example of
multivariate transfer function.
However, doing a Google search, hoping that we will find a code snippet
featuring a non-trivial example did not yield any results. Would you be
so kind and provide one?
Best regards,
Matej
https://github.com/matejak
Hi Matt,
Cheers for the info, I’ll take a look at your Unity code, I’m betting it’ll be better than my stab at it a few months back, I used someone else’s ray-marching code:
https://github.com/brianasu/unity-ray-marching/tree/volumetric-textures
Is it a mercurial repository? (I’ve just never used that before)
Hi Nathan,
My skillset doesn’t cover Open GL or general graphics either (yet, maybe in time), but I’m more than happy to help where I can and I have enough people here interested to get their data in VR that there’s a good amount of support. :)
I assume yt’s OpenGL volume rendering uses ray-marching?
Hi Suoqing,
Cheers, that’s cool, I’ll have to give that a try too, I’ll have to see if the data will look good on that.
Hi Jill,
That’s a good idea, I’d never heard of sketchfab, most of our datasets are more volume then surface, but it’s worth a try for some of my solar stuff.
Thanks guys, I’m getting loads of cool ideas from here.
Regards
Alex
________________________________________
From: yt-users-request(a)lists.spacepope.org [yt-users-request(a)lists.spacepope.org]
Sent: 23 May 2017 19:01
To: yt-users(a)lists.spacepope.org
Subject: yt-users Digest, Vol 111, Issue 44
Send yt-users mailing list submissions to
yt-users(a)lists.spacepope.org
To subscribe or unsubscribe via the World Wide Web, visit
http://lists.spacepope.org/listinfo.cgi/yt-users-spacepope.org
or, via email, send a message with subject or body 'help' to
yt-users-request(a)lists.spacepope.org
You can reach the person managing the list at
yt-users-owner(a)lists.spacepope.org
When replying, please edit your Subject line so it is more specific
than "Re: Contents of yt-users digest..."
Today's Topics:
1. Re: YT VR Visualisation (Naiman, Jill)
2. Re: YT VR Visualisation (Suoqing Ji)
3. Ramses MHD visualization (Daniel William Murray)
----------------------------------------------------------------------
Message: 1
Date: Tue, 23 May 2017 10:00:31 +1000
From: "Naiman, Jill" <jill.naiman(a)cfa.harvard.edu>
To: Discussion of the yt analysis package
<yt-users(a)lists.spacepope.org>
Subject: Re: [yt-users] YT VR Visualisation
Message-ID:
<CAFwWoV-K5rWnL_9LZVQ4ZVhvaH8znbMo0Vi1yrg0MeApeSxuXw(a)mail.gmail.com>
Content-Type: text/plain; charset="utf-8"
Hi All,
Just FYI - one can also generate points/surfaces with yt and upload them
somewhere like Sketchfab and use their VR mechanisms.
Looking forward to what comes of it!
Cheers,
-Jill
On Tue, May 23, 2017 at 2:14 AM, Nathan Goldbaum <nathan12343(a)gmail.com>
wrote:
>
>
> On Mon, May 22, 2017 at 11:07 AM, Matthew Turk <matthewturk(a)gmail.com>
> wrote:
>
>> Hi Alex,
>>
>> Yup, we have been, and with the vive and oculus specifically. Right
>> now though, it never got past much of a tech demo phase -- there's
>> code on my group's bitbucket account (bitbucket.org/data-exp-lab/ )
>> that runs in Unity, and Nathan Goldbaum got some of the stuff that's
>> in yt (the volume rendering) to work with the OpenVR python bindings
>> directly (not srue where that code is) but it's all still pretty early
>> on.
>
>
> This code isn't really usable yet. In order for this to actually work,
> we'll need a big overhaul of the OpenGL volume rendering in yt to use
> multiple passes.
>
> In addition there were significant performance issues that I never tracked
> down. Unfortunately I'm no expert on graphics programming and that's
> probably what's needed here.
>
>
>> Not sure it's quite shovel ready. There're some students coming
>> on soon that will take another look, but July might be a bit
>> optimistic... That being said, it'd be an awesome hack or
>> collaboration and if you are interested in trying to push it at all,
>> that'd be pretty rad.
>>
>> -Matt
>>
>> On Mon, May 22, 2017 at 10:10 AM, Alexander Hamilton
>> <A.Hamilton(a)2016.hull.ac.uk> wrote:
>> > Hi guys, I?m Alex from Hull University.
>> > We?re hosting the UK?s National Astronomy Meeting in July and I?ve been
>> > asked to run the hack day.
>> >
>> > I?ve managed to procure some VR headsets (HTC Vive, Oculus Rift and MS
>> > Hololens) and want to try and see if we could use them for visualising
>> some
>> > of our simulation data.
>> > I hear you guys may have been playing with VR, is that true? Is there
>> > anything that we could use/play-with at this point?
>> >
>> >
>> >
>> > Thanks
>> >
>> > Alex
>> >
>> > **************************************************
>> > To view the terms under which this email is
>> > distributed, please go to
>> > http://www2.hull.ac.uk/legal/disclaimer.aspx
>> > **************************************************
>> > _______________________________________________
>> > yt-users mailing list
>> > yt-users(a)lists.spacepope.org
>> > http://lists.spacepope.org/listinfo.cgi/yt-users-spacepope.org
>> >
>> _______________________________________________
>> yt-users mailing list
>> yt-users(a)lists.spacepope.org
>> http://lists.spacepope.org/listinfo.cgi/yt-users-spacepope.org
>>
>
>
> _______________________________________________
> yt-users mailing list
> yt-users(a)lists.spacepope.org
> http://lists.spacepope.org/listinfo.cgi/yt-users-spacepope.org
>
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.spacepope.org/pipermail/yt-users-spacepope.org/attachments/201…>
------------------------------
Message: 2
Date: Mon, 22 May 2017 17:18:20 -0700
From: Suoqing Ji <suoqing(a)physics.ucsb.edu>
To: Discussion of the yt analysis package
<yt-users(a)lists.spacepope.org>
Subject: Re: [yt-users] YT VR Visualisation
Message-ID: <ca3136c6-5fe9-4aac-92df-09d17d82bffc@Spark>
Content-Type: text/plain; charset="utf-8"
Hi Alex,
Just to add up one more point: yt is always able to generate pre-rendered virtual reality movies with stereo-spherical lens, which works perfectly with Google Cardboard.
To generate the movie, you could refer to this page: http://yt-project.org/docs/dev/visualizing/volume_rendering.html#spherical-…
An example move is at:?https://www.youtube.com/watch?v=ZYWY53X7UQE . You could also put your smartphone into Google Cardboard to watch the video with Youtube App.
Best wishes,
?
Suoqing Ji
Ph.D Candidate
Department of Physics
University of California, Santa Barbara
http://physics.ucsb.edu/~suoqing
On May 22, 2017, 5:00 PM -0700, Naiman, Jill <jill.naiman(a)cfa.harvard.edu>, wrote:
> Hi All,
>
> Just FYI - one can also generate points/surfaces with yt and upload them somewhere like Sketchfab and use their VR mechanisms.
>
> Looking forward to what comes of it!
>
> Cheers,
> -Jill
>
> > On Tue, May 23, 2017 at 2:14 AM, Nathan Goldbaum <nathan12343(a)gmail.com> wrote:
> > >
> > >
> > > > On Mon, May 22, 2017 at 11:07 AM, Matthew Turk <matthewturk(a)gmail.com> wrote:
> > > > > Hi Alex,
> > > > >
> > > > > Yup, we have been, and with the vive and oculus specifically.? Right
> > > > > now though, it never got past much of a tech demo phase -- there's
> > > > > code on my group's bitbucket account (bitbucket.org/data-exp-lab/ )
> > > > > that runs in Unity, and Nathan Goldbaum got some of the stuff that's
> > > > > in yt (the volume rendering) to work with the OpenVR python bindings
> > > > > directly (not srue where that code is) but it's all still pretty early
> > > > > on.
> > > >
> > > > This code isn't really usable yet. In order for this to actually work, we'll need a big overhaul of the OpenGL volume rendering in yt to use multiple passes.
> > > >
> > > > In addition there were significant performance issues that I never tracked down. Unfortunately I'm no expert on graphics programming and that's probably what's needed here.
> > > >
> > > > > Not sure it's quite shovel ready.? There're some students coming
> > > > > on soon that will take another look, but July might be a bit
> > > > > optimistic...? That being said, it'd be an awesome hack or
> > > > > collaboration and if you are interested in trying to push it at all,
> > > > > that'd be pretty rad.
> > > > >
> > > > > -Matt
> > > > >
> > > > > On Mon, May 22, 2017 at 10:10 AM, Alexander Hamilton
> > > > > <A.Hamilton(a)2016.hull.ac.uk> wrote:
> > > > > > Hi guys, I?m Alex from Hull University.
> > > > > > We?re hosting the UK?s National Astronomy Meeting in July and I?ve been
> > > > > > asked to run the hack day.
> > > > > >
> > > > > > I?ve managed to procure some VR headsets (HTC Vive, Oculus Rift and MS
> > > > > > Hololens) and want to try and see if we could use them for visualising some
> > > > > > of our simulation data.
> > > > > > I hear you guys may have been playing with VR, is that true? Is there
> > > > > > anything that we could use/play-with at this point?
> > > > > >
> > > > > >
> > > > > >
> > > > > > Thanks
> > > > > >
> > > > > >? ? ?Alex
> > > > > >
> > > > > > **************************************************
> > > > > > To view the terms under which this email is
> > > > > > distributed, please go to
> > > > > > http://www2.hull.ac.uk/legal/disclaimer.aspx
> > > > > > **************************************************
> > > > > > _______________________________________________
> > > > > > yt-users mailing list
> > > > > > yt-users(a)lists.spacepope.org
> > > > > > http://lists.spacepope.org/listinfo.cgi/yt-users-spacepope.org
> > > > > >
> > > > > _______________________________________________
> > > > > yt-users mailing list
> > > > > yt-users(a)lists.spacepope.org
> > > > > http://lists.spacepope.org/listinfo.cgi/yt-users-spacepope.org
> > >
> > >
> > > _______________________________________________
> > > yt-users mailing list
> > > yt-users(a)lists.spacepope.org
> > > http://lists.spacepope.org/listinfo.cgi/yt-users-spacepope.org
> > >
>
> _______________________________________________
> yt-users mailing list
> yt-users(a)lists.spacepope.org
> http://lists.spacepope.org/listinfo.cgi/yt-users-spacepope.org
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.spacepope.org/pipermail/yt-users-spacepope.org/attachments/201…>
------------------------------
Message: 3
Date: Tue, 23 May 2017 17:59:15 +0000
From: Daniel William Murray <dwmurray(a)uwm.edu>
To: "yt-users(a)lists.spacepope.org" <yt-users(a)lists.spacepope.org>
Subject: [yt-users] Ramses MHD visualization
Message-ID:
<MWHPR04MB1185DDA612C4C5D38349B6DDDCF90(a)MWHPR04MB1185.namprd04.prod.outlook.com>
Content-Type: text/plain; charset="iso-8859-1"
Dear yt folks,
I am trying to visualize some RAMSES MHD simulation data. The simulation is for gravitationally collapsing gas looking at star formation. I'm looking to obtain projection and slice plots and then annotate the magnetic fields onto the resulting image.
I'm on Mac OS X 10.9.5
I'm using conda 4.3.18
My version of yt is: 3.3.5-py27_0
and python 2.7.12
All of what I believe is the relevant code:
pf = yt.load('output_00284/info_00284.txt')
sp = pf.h.sphere(part_center, (zoom_width + 0.5, "pc"))
L = sp.quantities.angular_momentum_vector()
p = yt.OffAxisSlicePlot(pf, L, 'density', center=part_center, width =(zoom_width, "pc"))
p.annotate_magnetic_field()
p.save('test.png')
The error:
yt.utilities.exceptions.YTPlotCallbackError: annotate_magnetic_field callback failed with the following error: Could not find field '('all', 'magnetic_field_x')' in info_00284.
If I load the field list for the data set, I can see that the B field data is there, but under 'ramses', not 'all'
('ramses', 'x-Bfield-left'),
('ramses', 'x-Bfield-right'),
('ramses', 'x-velocity'),
('ramses', 'y-Bfield-left'),
('ramses', 'y-Bfield-right'),
('ramses', 'y-velocity'),
('ramses', 'z-Bfield-left'),
('ramses', 'z-Bfield-right'),
Is there a way to pass these fields to annotate_magnetic_field? Or a manner in which I can have this information load in the field ('all', 'magnetic_field_')
I should also note that I can create and save the plots with no issues, so long as I do not try to annotate the magnetic fields.
Cheers,
Daniel Murray
Graduate student at University of Wisconsin Milwaukee
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.spacepope.org/pipermail/yt-users-spacepope.org/attachments/201…>
------------------------------
Subject: Digest Footer
_______________________________________________
yt-users mailing list
yt-users(a)lists.spacepope.org
http://lists.spacepope.org/listinfo.cgi/yt-users-spacepope.org
------------------------------
End of yt-users Digest, Vol 111, Issue 44
*****************************************
**************************************************
To view the terms under which this email is
distributed, please go to
http://www2.hull.ac.uk/legal/disclaimer.aspx
**************************************************
Dear yt folks,
I am trying to visualize some RAMSES MHD simulation data. The simulation is for gravitationally collapsing gas looking at star formation. I'm looking to obtain projection and slice plots and then annotate the magnetic fields onto the resulting image.
I'm on Mac OS X 10.9.5
I'm using conda 4.3.18
My version of yt is: 3.3.5-py27_0
and python 2.7.12
All of what I believe is the relevant code:
pf = yt.load('output_00284/info_00284.txt')
sp = pf.h.sphere(part_center, (zoom_width + 0.5, "pc"))
L = sp.quantities.angular_momentum_vector()
p = yt.OffAxisSlicePlot(pf, L, 'density', center=part_center, width =(zoom_width, "pc"))
p.annotate_magnetic_field()
p.save('test.png')
The error:
yt.utilities.exceptions.YTPlotCallbackError: annotate_magnetic_field callback failed with the following error: Could not find field '('all', 'magnetic_field_x')' in info_00284.
If I load the field list for the data set, I can see that the B field data is there, but under 'ramses', not 'all'
('ramses', 'x-Bfield-left'),
('ramses', 'x-Bfield-right'),
('ramses', 'x-velocity'),
('ramses', 'y-Bfield-left'),
('ramses', 'y-Bfield-right'),
('ramses', 'y-velocity'),
('ramses', 'z-Bfield-left'),
('ramses', 'z-Bfield-right'),
Is there a way to pass these fields to annotate_magnetic_field? Or a manner in which I can have this information load in the field ('all', 'magnetic_field_')
I should also note that I can create and save the plots with no issues, so long as I do not try to annotate the magnetic fields.
Cheers,
Daniel Murray
Graduate student at University of Wisconsin Milwaukee