Mailman 3 August 2018 - yt-users

Porting .ply files to Blender with colors?
by Jill P. Naiman Aug. 15, 2022

Aug. 15, 2022

Hi All, I was curious if anyone has had any luck porting .ply files generated with yt into Blender with colors? It seems to port in the surfaces just fine, but they are a uniform color. From what I can gather, Blender imports .ply files with face color instead of vertex so I tried the surf.export_ply routine with "face" instead of "vertex" for sample_type, but no luck ("vertex" doesn't work either). Let me know if I can provide any more useful info. Thanks for any help, -J

3 2

Stars on volume rendering
by Elizabeth Harper-Clark Dec. 11, 2021

Dec. 11, 2021

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 April 11, 2020

April 11, 2020

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

3-D Mass power spectra
by 梁家豪 June 3, 2019

June 3, 2019

Dear yt Can current yt calculate 3-D Mass power spectra? I checked the website but I didn't find any information. I think calculating 3-D Mass power spectra is a very useful for cosmological simulations. So I guess maybe yt supports this function now....? Thanks in advance

4 7

Internal memory usage
by Nick Gnedin Sept. 7, 2018

Sept. 7, 2018

Folks, I would like to test the integrity of a large dataset, for example by finding the max value of some data field. The dataset is too large to be fully loaded into memory, so I would like to read it slice by slice and compute max in each slice, and then combine them together. To that end I wrote the following script: import yt, yt.funcs import sys import numpy as np f = sys.argv[1] ds = yt.load(f) n = ds.domain_dimensions[0] dmax = np.zeros((n,)) for i in range(n): d = ds.r[i:(i+1),:,:][('artio', 'HVAR_GAS_DENSITY')] dmax[i] = np.amax(d) print(i,dmax[i],yt.funcs.get_memory_usage()) del d print(np.amax(dmax)) However, when I run it, I get: yt : [INFO ] 2018-08-28 10:52:58,075 Created 4945 chunks for ARTIO 0 16945762.0 202.5546875 yt : [INFO ] 2018-08-28 10:52:59,388 Created 1232 chunks for ARTIO 1 2416576.25 221.36328125 yt : [INFO ] 2018-08-28 10:53:01,599 Created 9635 chunks for ARTIO 2 10419311.0 269.06640625 ... yt : [INFO ] 2018-08-28 10:53:23,397 Created 1474 chunks for ARTIO 13 2395590.5 698.91015625 yt : [INFO ] 2018-08-28 10:53:25,594 Created 9747 chunks for ARTIO 14 11139424.0 739.16015625 The number of chunks created varies in each iteration of the loop, but the total memory usage is still steadily climbing up. Could you advise what I am doing wrong? Many thanks, Nick

2 1

Large Memory Usage
by Jared Coughlin Aug. 31, 2018

Aug. 31, 2018

Hello! I'm using the sph-viz branch of yt and I'm trying to make a density projection plot from a gadget-2 snapshot. This works like a charm when I use a snapshot with 428^3 particles in a 15 Mpc box. However, when I try to do it on a snapshot with 1024^3 particles in a 40 Mpc box, my job keeps getting killed because it is trying to exceed the available system memory. I'm submitting it to a node with 256 GB of RAM, and I'm requesting all of it, which means that more than this is trying to be used. The snapshot itself is ~ 40 GB. The memory usage appears to spike when allocating for the KDTree. Before that, it appears to generate the .ewah file without issue. I was just wondering if anyone had any thoughts on why this large spike in memory might be occurring and how I might go about fixing it? If not, no worries, I can try and investigate it myself, but it's a pain to do with pdb. My gadget simulation used a maximum of ~ 173 GB of RAM, so I feel like the kdtree shouldn't be using as much memory as it's trying to. Thanks, and sorry for the bother! Sincerely, -Jared

3 10

clump finder
by Stephanie Tonnesen Aug. 31, 2018

Aug. 31, 2018

Dear yt-users, I am having some trouble with the yt clump finder. Mostly, my issue is that I don't understand what I am getting when I annotate the projectionplot. Instead of the nice lines in the example, I get strange filled-in circles? My figure is attached, and here is my code: ################## ds = yt.load("/mnt/ceph/users/stonnesen/runs/JO204wind/DD"+loop[i]+"/JO204cw"+loop[i]) ad = ds.all_data() data_source = ad.cut_region(["obj['z'] > 0.27"]) #ds.disk([0.5, 0.5, 0.5], [0., 0., 1.],(80, 'kpc'), (145, 'kpc')) # the field to be used for contouring field = ("gas", "density") # This is the multiplicative interval between contours. step = 2.0 # Now we set some sane min/max values between which we want to find contours # This is how we tell the clump finder what to look for -- it won't look for # contours connected below or above these threshold values. c_min = 2e-25 #10**np.floor(np.log10(data_source[field]).min()+3 ) c_max = 10**np.floor(np.log10(data_source[field]).max()+1) # Now find get our 'base' clump -- this one just covers the whole domain. master_clump = Clump(data_source, 'density') # Add a "validator" to weed out clumps with less than 16 cells. # As many validators can be added as you want. master_clump.add_validator("min_cells", 16) # Calculate center of mass for all clumps. master_clump.add_info_item("center_of_mass") # Begin clump finding. find_clumps(master_clump, c_min, c_max, step) # Save the clump tree as a reloadable dataset fn = master_clump.save_as_dataset(fields=["density"]) # Traverse the clump hierarchy to get a list of all of the 'leaf' clumps leaf_clumps = get_lowest_clumps(master_clump) # To visualize these clumps, a list of clumps can be supplied to # the "clumps" callback on a plot. First, we create a projection plot: prj = yt.ProjectionPlot(ds, 0, 'density', center=[0.5,0.5,0.5],width=(160,'kpc')) prj.set_buff_size(4000) # Next we annotate the plot with contours on the borders of the clumps prj.annotate_clumps(leaf_clumps) # Save the plot to disk. prj.save('JO204cw'+loop[i]+'_clumps_16pix') ########################## Also, when I have > print (leaf_clumps[0]["gas","density"]) I get 22 outputs, which seems like a lot less than what is getting put on the figure? Anyway, I am lost and any help would be greatly appreciated. Best, Stephanie -- Dr. Stephanie Tonnesen Associate Research Scientist CCA, Flatiron Institute New York, NY stonnes(a)gmail.com

3 6

set_buff_size
by Stephanie Tonnesen Aug. 30, 2018

Aug. 30, 2018

Hi yt-users! So I have a question about set_buff_size. I would like my pixels to be the same size as my more refined cell in a projectionplot, so have a couple lines of code like this: plot6z=yt.ProjectionPlot(ds,"x","density",center=[0.5,0.5,0.33],width=(76,'kpc')) plot6z.set_buff_size(2000) plot6z.save('JO201scc_76kpc_'+loop[i]+'t'+times[i]) When I take a look at the png file that is made, there are 800 pixels across instead of the 2000 that I expect. I just checked the docs here: http://yt-project.org/doc/visualizing/plots.html and based on the SlicePlot example I feel like it should work! Then again, it looks like in the example the number of pixels also doesn't change? I would appreciate any/all advice! Thanks!! Stephanie -- Dr. Stephanie Tonnesen Associate Research Scientist CCA, Flatiron Institute New York, NY stonnes(a)gmail.com

2 4

ANN: ytree 2.2 released
by Britton Smith Aug. 28, 2018

Aug. 28, 2018

Hi everyone, ytree version 2.2 is released! ytree is an extension of yt for working with merger-tree data from multiple sources. This release adds a couple new features as well as improvements for i/o and memory usage. - New feature: vector fields return x/y/z components of a vector as a single array with automatically created magnitude fields. - New feature: a select_halos function to search all trees for halos matching specified criteria. - Enhancement: significant refactors of the field and i/o systems to reduce memory, speed up field access, and simplify implementing new frontends. - Enhancement: upgrades to the testing infrastructure. ytree can be installed with pip: pip install ytree or from source, by going to https://github.com/brittonsmith/ytree git clone https://github.com/brittonsmith/ytree cd ytree pip install -e . Full documentation can be found at http://ytree.readthedocs.io/. As always, ytree would be nothing without the wonderful yt community. Thanks to all those amazing people! Britton

1 0

ARTIO question
by Nick Gnedin Aug. 27, 2018

Aug. 27, 2018

Folks, I need to make a small change in the ARTIO backend - I added an option to compress grid data, and that required a small change in the meaning of some disk variables. In order to test these changes, I made them inside yt/frontends/artio/artio_headers on my local anaconda install, but I do not know how to compile them into yt - sorry for a such a dummy question. Any advice? Nick Gnedin

2 1

Re: yt-users Digest, Vol 126, Issue 10
by Slavin, Jonathan Aug. 24, 2018

Aug. 24, 2018

Alex, Not sure why, but your message acquired some odd extra characters. In any case if you were saying that you can do something like: x = 43*yt.units.pc print('{:.4f}'.format(x)) then we must be using different versions of yt or something. If I do that I get: --------------------------------------------------------------------------- TypeError Traceback (most recent call last) <ipython-input-6-2df4b250bc03> in <module>() ----> 1 print('{:.4f}'.format(x)) TypeError: unsupported format string passed to YTQuantity.__format__ I can do things like print('{:.4f}'.format(x.v.item())) though that's a bit cumbersome. I'm using yt version 3.4.1 and python 3.6.4 Jon On Fri, Aug 24, 2018 at 12:00 PM, <yt-users-request(a)python.org> wrote: > Send yt-users mailing list submissions to > yt-users(a)python.org > > To subscribe or unsubscribe via the World Wide Web, visit > https://mail.python.org/mm3/mailman3/lists/yt-users.python.org/ > or, via email, send a message with subject or body 'help' to > yt-users-request(a)python.org > > You can reach the person managing the list at > yt-users-owner(a)python.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: formatting of YTQuantities (Hill, Alex) > > > ---------------------------------------------------------------------- > > Date: Fri, 24 Aug 2018 14:36:23 +0000 > From: "Hill, Alex" <ashill(a)astro.ubc.ca> > Subject: [yt-users] Re: formatting of YTQuantities > To: Discussion of the yt analysis package <yt-users(a)python.org> > Message-ID: <0064EDE6-8538-4BA4-9A35-88BFCD4A29B8(a)astro.ubc.ca> > Content-Type: multipart/signed; > boundary="Apple-Mail=_D985D977-7CFB-4A1E-809B-69B20CE817D5"; > protocol="application/pkcs7-signature"; micalg=sha1 > > --Apple-Mail=_D985D977-7CFB-4A1E-809B-69B20CE817D5 > Content-Transfer-Encoding: quoted-printable > Content-Type: text/plain; > charset=utf-8 > > Hi Jon, > > For me, > > asdf =3D 43*yt.units.pc > print( '{:.4f}'.format(asdf) ) > > outputs > > 43.0000 > > and > > print( '{:.4f} {}'.format(asdf, asdf.units) ) > > outputs > > 43.0000 pc > > Is that what you want? It=E2=80=99s not quite clean, but not too bad. > > Cheers, > Alex > > --------- > Alex Hill > University of British Columbia and Space Science Institute > email: ashill(a)astro.ubc.ca > > Based at > Dominion Radio Astrophysical Observatory, National Research Council = > Canada > Mailing: PO Box 248, Penticton, BC V2A 6K3, Canada > Physical: 717 White Lake Road, Kaleden, BC V0H 1K0, Canada > > phone: +1 250-497-2356 > > > On Aug 23, 2018, at 06:21, Slavin, Jonathan <jslavin(a)cfa.harvard.edu> = > wrote: > >=20 > > Hi, > >=20 > > I'm wondering if there's some way to format the way a YTQuantity = > prints, for example the number of decimal places. Simply using print() = > leads to many decimal places and includes the units. It's nice to have = > the units included automatically, but would also be good to have control = > over the number of decimal places. =20 > >=20 > > Actually, from looking at the source I see how the __str__ method for = > YTArray and YTQuantity objects is defined, which would indicate that no = > such formatting is possible, as built in -- though I suppose one could = > use the value and units parts to do the outputting. It might be nice to = > have a simple method built in like say valunit(x) which would return a = > tuple of the value and unit so then you could do > > print('{:.4f} {}'.format(x.valunit())) > >=20 > > Jon > >=20 > > --=20 > > ________________________________________________________ > > Jonathan D. Slavin Harvard-Smithsonian CfA > > jslavin(a)cfa.harvard.edu 60 Garden Street, MS 83 > > phone: (617) 496-7981 Cambridge, MA 02138-1516 > > cell: (781) 363-0035 USA > > ________________________________________________________ > >=20 > > _______________________________________________ > > yt-users mailing list -- yt-users(a)python.org > > To unsubscribe send an email to yt-users-leave(a)python.org > > > --Apple-Mail=_D985D977-7CFB-4A1E-809B-69B20CE817D5 > Content-Disposition: attachment; filename="smime.p7s" > Content-Type: application/pkcs7-signature; name="smime.p7s" > Content-Transfer-Encoding: base64 > > MIAGCSqGSIb3DQEHAqCAMIACAQExCzAJBgUrDgMCGgUAMIAGCSqGSIb3DQEH > AQAAoIILIDCCBTIw > ggQaoAMCAQICEQDFlowlhsLICe0pbW5w+5UkMA0GCSqGSIb3DQEBCwUAMIGXMQs > wCQYDVQQGEwJH > QjEbMBkGA1UECBMSR3JlYXRlciBNYW5jaGVzdGVyMRAwDgYDVQQHEwdTYWxm > b3JkMRowGAYDVQQK > ExFDT01PRE8gQ0EgTGltaXRlZDE9MDsGA1UEAxM0Q09NT0RPIFJTQSBDbGll > bnQgQXV0aGVudGlj > YXRpb24gYW5kIFNlY3VyZSBFbWFpbCBDQTAeFw0xODA3MDQwMDAwMDBaFw0x > OTA3MDQyMzU5NTla > MCQxIjAgBgkqhkiG9w0BCQEWE2FzaGlsbEBhc3Ryby51YmMuY2EwggEiMA0G > CSqGSIb3DQEBAQUA > A4IBDwAwggEKAoIBAQDQwbnx4Ficrq2UzJZ/ptxVG5vHYghfYDpZnQ0RwgyEi+ > seF7fW+FT7cK/c > 3dCLBMtKrX3/1emA6xMkIamVakeQw2Est9AFu0vSHJQvs6X0GOiEoVnD+bWYhzcD// > DJNY17c52p > CPz05xqzkYQFLoRrzzol8UaHsKcl10UKI0iKUA2sbWV2AiF4QhD6t4tADQRw > RYdTE+5RLU/zZzRj > 0lwN2mHmzkNrg0HjymNtdN7zz4IYnVXHImPqJlR3VgPXQgzpjX/Yo7ZJ+ > w4GpE7mvqehcZjKe3+N > qKP5X8+Zp7eF3A2tsk2+Wyzrf86hbE4LZOruB2YzbaMSy5LSEH > zjOAvVAgMBAAGjggHpMIIB5TAf > BgNVHSMEGDAWgBSCr2yM+MX+lmF86B89K3FIXsSLwDAdBgNVHQ4EFg > QUgGjEUtyOKsQBtWkxscuN > eVndV78wDgYDVR0PAQH/BAQDAgWgMAwGA1UdEwEB/wQCMAAwIAYDVR0lBBkwFwYIKwYBBQU > HAwQG > CysGAQQBsjEBAwUCMBEGCWCGSAGG+EIBAQQEAwIFIDBGBgNVHSAEPzA9MDs > GDCsGAQQBsjEBAgEB > ATArMCkGCCsGAQUFBwIBFh1odHRwczovL3NlY3VyZS5jb21vZG8ubmV0L0NQ > UzBaBgNVHR8EUzBR > ME+gTaBLhklodHRwOi8vY3JsLmNvbW9kb2NhLmNvbS9DT01PRE9SU0FDbGllbnR > BdXRoZW50aWNh > dGlvbmFuZFNlY3VyZUVtYWlsQ0EuY3JsMIGLBggrBgEFBQcBAQR/ > MH0wVQYIKwYBBQUHMAKGSWh0 > dHA6Ly9jcnQuY29tb2RvY2EuY29tL0NPTU9ET1JTQUNsaWVudEF1dGhlbnRp > Y2F0aW9uYW5kU2Vj > dXJlRW1haWxDQS5jcnQwJAYIKwYBBQUHMAGGGGh0dHA6Ly9vY3NwLmNvbW9k > b2NhLmNvbTAeBgNV > HREEFzAVgRNhc2hpbGxAYXN0cm8udWJjLmNhMA0GCSqGSIb3DQEBCwUAA4IB > AQAPwkdRlOCH90lM > /Hp1wbNLbCvTaPgmv9K7dy2M6EPhyxOCHCWOQpD8hWLQWDw24qMrdYUc/Sv+ > UKZkvgjQXU2JKtTu > aIq505XLY0CX2ApOGCejOQromBdLjbKeg0j+Z4Nii2qud4TXjGrsQeFcHFQ0Nm3iXq > aDcCgKO9W3 > JE8fnQynERZUn395GzGILUHFPEG4n8XJgn5z5IPH6ZSnUSzhZ4+ > jUPCL19LrG67UCjyf2aTjO46U > +JxRn6zQPxRAnO8Fwi8Fc0weEWLh5T8NGPiQaVFK9ojkwdWLWnKqu7/O/ > h0FlHZPKl0VYAEfi/RA > qF5hJFsXI4Z8v9p3ZhbxgxdpMIIF5jCCA86gAwIBAgIQapvhODv/ > K2ufAdXZuKdSVjANBgkqhkiG > 9w0BAQwFADCBhTELMAkGA1UEBhMCR0IxGzAZBgNVBAgTEkdyZWF0ZXIgTWFu > Y2hlc3RlcjEQMA4G > A1UEBxMHU2FsZm9yZDEaMBgGA1UEChMRQ09NT0RPIENBIExpbWl0ZWQxKzAp > BgNVBAMTIkNPTU9E > TyBSU0EgQ2VydGlmaWNhdGlvbiBBdXRob3JpdHkwHhcNMTMwMTEwMDAwMDAw > WhcNMjgwMTA5MjM1 > OTU5WjCBlzELMAkGA1UEBhMCR0IxGzAZBgNVBAgTEkdyZWF0ZXIgTWFuY2hl > c3RlcjEQMA4GA1UE > BxMHU2FsZm9yZDEaMBgGA1UEChMRQ09NT0RPIENBIExpbWl0ZWQxPTA7BgNV > BAMTNENPTU9ETyBS > U0EgQ2xpZW50IEF1dGhlbnRpY2F0aW9uIGFuZCBTZWN1cmUgRW1haWwgQ0Ew > ggEiMA0GCSqGSIb3 > DQEBAQUAA4IBDwAwggEKAoIBAQC+s55XrCh2dUAWxzgDmNPGGHYhUPMleQ > tMtaDRfTpYPpynMS6n > 9jR22YRq2tA9NEjk6vW7rN/5sYFLIP1of3l0NKZ6fLWfF2VgJ5cij > KYy/qlAckY1wgOkUMgzKlWl > VJGyK+UlNEQ1/5ErCsHq9x9aU/x1KwTdF/LCrT03Rl/FwFrf1XTCwa2QZYL55AqLPikFlgqOt > zk0 > 6kb2qvGlnHJvijjI03BOrNpo+kZGpcHsgyO1/u1OZTaOo8wvEU17VVeP1cHWse9tGKT > DyUGg2hJZ > jrqck39UIm/nKbpDSZ0JsMoIw/JtOOg0JC56VzQgBo7ictReTQE5LFLG > 3yQK+xS1AgMBAAGjggE8 > MIIBODAfBgNVHSMEGDAWgBS7r34CPfqm8TyEjq3uOJjs2TIy1DAdBgNVHQ4E > FgQUgq9sjPjF/pZh > fOgfPStxSF7Ei8AwDgYDVR0PAQH/BAQDAgGGMBIGA1UdEwEB/ > wQIMAYBAf8CAQAwEQYDVR0gBAow > CDAGBgRVHSAAMEwGA1UdHwRFMEMwQaA/oD2GO2h0dHA6Ly9jcmwuY29tb2RvY2 > EuY29tL0NPTU9E > T1JTQUNlcnRpZmljYXRpb25BdXRob3JpdHkuY3JsMHEGCCsGAQUFBwEBBGUw > YzA7BggrBgEFBQcw > AoYvaHR0cDovL2NydC5jb21vZG9jYS5jb20vQ09NT0RPUlNBQWRkVHJ1c3RD > QS5jcnQwJAYIKwYB > BQUHMAGGGGh0dHA6Ly9vY3NwLmNvbW9kb2NhLmNvbTANBgkqhkiG9w0BAQwF > AAOCAgEAeFyygSg0 > TzzuX1bOn5dW7I+iaxf28/ZJCAbU2C81zd9A/tNx4+jsQgwRGiHjZrAYayZrrm78hOx7aEpk > fNPQ > IHGG6Fvq3EzWf/Lvx7/hk6zSPwIal9v5IkDcZoFD7f3iT7Pdk > HJY9B51csvU50rxpEg1OyOT8fk2 > zvvPBuM4qQNqbGWlnhMpIMwpWZT89RY0wpJO+2V6eXEGGHsROs3njeP9DqqqAJaBa4w > BeKOdGCWn > 1/Jp2oY6dyNmNppI4ZNMUH4Tam85S1j6E95u4+1Nuru84OrMIzqvISE2HN/ > 56ebTOWlcrurffade > 2022O/tUU1gb4jfWCcyvB8czm12FgX/y/lRjmDbEA08QJNB2729Y+ > io1IYO3ztveBdvUCIYZojTq > /OCR6MvnzS6X72HP0PRLRTiOSEmIDsS5N5w/8IW1Hva5hEFy6fDAfd9yI+O+ > IMMAj1KcL/Zo9jzJ > 16HO5m60ttl1Enk8MQkz/W3JlHaeI5iKFn4UJu1/cP2YHXYPiWf2JyBzsLBrGk1II+ > 3yL8aorYew > 6CQvdVifC3HtwlSam9V1niiCfOBe2C12TdKGu05LWIA3ZkFcWJGaNXOZ6Ggy > h/TqvXG5v7zmEVDN > XFnHn9tFpMpOUvxhcsjycBtH0dZ0WrNw6gH+HF8TIhCnH3+ > zzWuDN0Rk6h9KVkfKehIxggO6MIID > tgIBATCBrTCBlzELMAkGA1UEBhMCR0IxGzAZBgNVBAgTEkdyZWF0ZXIgTWFu > Y2hlc3RlcjEQMA4G > A1UEBxMHU2FsZm9yZDEaMBgGA1UEChMRQ09NT0RPIENBIExpbWl0ZWQxPTA7 > BgNVBAMTNENPTU9E > TyBSU0EgQ2xpZW50IEF1dGhlbnRpY2F0aW9uIGFuZCBTZWN1cmUgRW1haWwg > Q0ECEQDFlowlhsLI > Ce0pbW5w+5UkMAkGBSsOAwIaBQCgggHhMBgGCSqGSIb3DQEJAzELBgkqhkiG9w0BBwEwH > AYJKoZI > hvcNAQkFMQ8XDTE4MDgyNDE0MzYyM1owIwYJKoZIhvcNAQkEMRYEFEUDj2AB > 1bAIz/IWzePCRtat > u3cUMIG+BgkrBgEEAYI3EAQxgbAwga0wgZcxCzAJBgNVBAYTAkdCMRswGQYDVQQIExJH > cmVhdGVy > IE1hbmNoZXN0ZXIxEDAOBgNVBAcTB1NhbGZvcmQxGjAYBgNVBAoTEUNPTU9E > TyBDQSBMaW1pdGVk > MT0wOwYDVQQDEzRDT01PRE8gUlNBIENsaWVudCBBdXRoZW50aWNhdGlvbiBh > bmQgU2VjdXJlIEVt > YWlsIENBAhEAxZaMJYbCyAntKW1ucPuVJDCBwAYLKoZIhvcNAQkQAgsxgbCg > ga0wgZcxCzAJBgNV > BAYTAkdCMRswGQYDVQQIExJHcmVhdGVyIE1hbmNoZXN0ZXIxEDAOBgNVBAcT > B1NhbGZvcmQxGjAY > BgNVBAoTEUNPTU9ETyBDQSBMaW1pdGVkMT0wOwYDVQQDEzRDT01PRE8gUlNB > IENsaWVudCBBdXRo > ZW50aWNhdGlvbiBhbmQgU2VjdXJlIEVtYWlsIENBAhEAxZaMJYbCyAntKW1u > cPuVJDANBgkqhkiG > 9w0BAQEFAASCAQAf1DYeDzYVAY9jqjLFnxUdieQxHX8ADKJ767oWMrAX63Tu > /5uGGVOoE866FNGU > DW2tJUpk3irkFoTuOSruKAQLg2j+F7ZGyCHB0ChmWj1Zfs1yykPNWBDfhE > jifl9rd3bOzdtqHicX > 04W4WpEjnuRmo83tPhtYUN0oZoT0IPaJm65k8EsGLQkjY1z// > wvKNsZC1vwwBdQojZvrC1kyjoWv > rni/LLuvAXJePx8qroiyr7G4gMgUEpC+YdEo04SN5rsBw9QlHxjeNTwNXvoBzr > Zg6eZxBcCahOLt > QJPRap/MX13Ntabw135oBDtTlygJb98xzzoNxajTDj9WNI6c+xcOAAAAAAAA > > --Apple-Mail=_D985D977-7CFB-4A1E-809B-69B20CE817D5-- > > ------------------------------ > > Subject: Digest Footer > > _______________________________________________ > yt-users mailing list -- yt-users(a)python.org > To unsubscribe send an email to yt-users-leave(a)python.org > > > ------------------------------ > > End of yt-users Digest, Vol 126, Issue 10 > ***************************************** > -- ________________________________________________________ Jonathan D. Slavin Harvard-Smithsonian CfA jslavin(a)cfa.harvard.edu 60 Garden Street, MS 83 phone: (617) 496-7981 Cambridge, MA 02138-1516 cell: (781) 363-0035 USA ________________________________________________________

2 1

formatting of YTQuantities
by Slavin, Jonathan Aug. 24, 2018

Aug. 24, 2018

Hi, I'm wondering if there's some way to format the way a YTQuantity prints, for example the number of decimal places. Simply using print() leads to many decimal places and includes the units. It's nice to have the units included automatically, but would also be good to have control over the number of decimal places. Actually, from looking at the source I see how the __str__ method for YTArray and YTQuantity objects is defined, which would indicate that no such formatting is possible, as built in -- though I suppose one could use the value and units parts to do the outputting. It might be nice to have a simple method built in like say valunit(x) which would return a tuple of the value and unit so then you could do print('{:.4f} {}'.format(x.valunit())) Jon -- ________________________________________________________ Jonathan D. Slavin Harvard-Smithsonian CfA jslavin(a)cfa.harvard.edu 60 Garden Street, MS 83 phone: (617) 496-7981 Cambridge, MA 02138-1516 cell: (781) 363-0035 USA ________________________________________________________

3 3

yt build error, 'numpy' has no '__version__'
by David Collins Aug. 8, 2018

Aug. 8, 2018

Hi, Everybody-- I'm trying to build yt, and I keep getting the following error (full stack at the end) File "setup.py", line 317, in finalize_options if LooseVersion(numpy.__version__) < LooseVersion('1.10.4'): AttributeError: module 'numpy' has no attribute '__version__' Has anyone else seen this, and have a fix? I'm installing this from the install_script.sh, and it doesn't seem to matter what options I select for building. This is on an Ubuntu 16.04 box, I recently did a large update, after which I started getting this error. Thanks! d. EXECUTING: ln -s /home/dcollins4096/local-yt-2018-08-06/yt-conda/envs/py27/bin/hg /home/dcollins4096/local-yt-2018-08-06/yt-conda/bin EXECUTING: /home/dcollins4096/local-yt-2018-08-06/yt-conda/bin/git clone https://github.com/yt-project/yt_conda /home/dcollins4096/local-yt-2018-08-06/yt-conda/src/yt_conda Cloning into '/home/dcollins4096/local-yt-2018-08-06/yt-conda/src/yt_conda'... EXECUTING: /home/dcollins4096/local-yt-2018-08-06/yt-conda/bin/git clone https://github.com/yt-project/yt /home/dcollins4096/local-yt-2018-08-06/yt-conda/src/yt-git Cloning into '/home/dcollins4096/local-yt-2018-08-06/yt-conda/src/yt-git'... EXECUTING: /home/dcollins4096/local-yt-2018-08-06/yt-conda/bin/git -C /home/dcollins4096/local-yt-2018-08-06/yt-conda/src/yt-git checkout master Already on 'master' Your branch is up to date with 'origin/master'. Using OpenMP to compile parallel extensions running develop running egg_info creating yt.egg-info writing yt.egg-info/PKG-INFO writing dependency_links to yt.egg-info/dependency_links.txt writing entry points to yt.egg-info/entry_points.txt writing requirements to yt.egg-info/requires.txt writing top-level names to yt.egg-info/top_level.txt writing manifest file 'yt.egg-info/SOURCES.txt' Traceback (most recent call last): File "setup.py", line 399, in <module> python_requires='>=2.7,!=3.0.*,!=3.1.*,!=3.2.*,!=3.3.*,!=3.4.*' File "/home/dcollins4096/local-yt-2018-08-06/yt-conda/lib/python3.6/site-packages/setuptools/__init__.py", line 131, in setup return distutils.core.setup(**attrs) File "/home/dcollins4096/local-yt-2018-08-06/yt-conda/lib/python3.6/distutils/core.py", line 148, in setup dist.run_commands() File "/home/dcollins4096/local-yt-2018-08-06/yt-conda/lib/python3.6/distutils/dist.py", line 955, in run_commands self.run_command(cmd) File "/home/dcollins4096/local-yt-2018-08-06/yt-conda/lib/python3.6/distutils/dist.py", line 974, in run_command cmd_obj.run() File "/home/dcollins4096/local-yt-2018-08-06/yt-conda/lib/python3.6/site-packages/setuptools/command/develop.py", line 38, in run self.install_for_development() File "/home/dcollins4096/local-yt-2018-08-06/yt-conda/lib/python3.6/site-packages/setuptools/command/develop.py", line 134, in install_for_development self.run_command('egg_info') File "/home/dcollins4096/local-yt-2018-08-06/yt-conda/lib/python3.6/distutils/cmd.py", line 313, in run_command self.distribution.run_command(command) File "/home/dcollins4096/local-yt-2018-08-06/yt-conda/lib/python3.6/distutils/dist.py", line 974, in run_command cmd_obj.run() File "/home/dcollins4096/local-yt-2018-08-06/yt-conda/lib/python3.6/site-packages/setuptools/command/egg_info.py", line 278, in run self.find_sources() File "/home/dcollins4096/local-yt-2018-08-06/yt-conda/lib/python3.6/site-packages/setuptools/command/egg_info.py", line 293, in find_sources mm.run() File "/home/dcollins4096/local-yt-2018-08-06/yt-conda/lib/python3.6/site-packages/setuptools/command/egg_info.py", line 524, in run self.add_defaults() File "/home/dcollins4096/local-yt-2018-08-06/yt-conda/lib/python3.6/site-packages/setuptools/command/egg_info.py", line 560, in add_defaults sdist.add_defaults(self) File "/home/dcollins4096/local-yt-2018-08-06/yt-conda/lib/python3.6/site-packages/setuptools/command/py36compat.py", line 36, in add_defaults self._add_defaults_ext() File "/home/dcollins4096/local-yt-2018-08-06/yt-conda/lib/python3.6/site-packages/setuptools/command/py36compat.py", line 119, in _add_defaults_ext build_ext = self.get_finalized_command('build_ext') File "/home/dcollins4096/local-yt-2018-08-06/yt-conda/lib/python3.6/distutils/cmd.py", line 299, in get_finalized_command cmd_obj.ensure_finalized() File "/home/dcollins4096/local-yt-2018-08-06/yt-conda/lib/python3.6/distutils/cmd.py", line 107, in ensure_finalized self.finalize_options() File "setup.py", line 317, in finalize_options if LooseVersion(numpy.__version__) < LooseVersion('1.10.4'): AttributeError: module 'numpy' has no attribute '__version__' -- -- Sent from a computer. -- -- Sent from a computer.

3 5

Annotate sphere
by Sushilkumar Aug. 2, 2018

Aug. 2, 2018

Dear Nathan: Thank you On Wed, Aug 1, 2018 at 11:00 AM <yt-users-request(a)python.org> wrote: > Send yt-users mailing list submissions to > yt-users(a)python.org > > To subscribe or unsubscribe via the World Wide Web, visit > https://mail.python.org/mm3/mailman3/lists/yt-users.python.org/ > or, via email, send a message with subject or body 'help' to > yt-users-request(a)python.org > > You can reach the person managing the list at > yt-users-owner(a)python.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: Annotate sphere (Nathan Goldbaum) > > > ---------------------------------------------------------------------- > > Date: Tue, 31 Jul 2018 15:26:44 -0500 > From: Nathan Goldbaum <nathan12343(a)gmail.com> > Subject: [yt-users] Re: Annotate sphere > To: Discussion of the yt analysis package <yt-users(a)python.org> > Message-ID: > <CAJXewOmdDDEkDHYGxgLSvO7JSj5tkca1WeMRLpb0E= > xueYxzGA(a)mail.gmail.com> > Content-Type: multipart/alternative; > boundary="000000000000c1fdcc0572516664" > > --000000000000c1fdcc0572516664 > Content-Type: text/plain; charset="UTF-8" > > On Tue, Jul 31, 2018 at 3:03 PM, Sushilkumar <sushil.sush19us(a)gmail.com> > wrote: > > > Dear yt: > > > > Is there a way to annotate sphere to a plot? I used annotate_sphere but > it > > plots an outline of a circle but I would like a "filled circle" to > > represent a star? > > > > I have a python script (see below) to create a sphere, but not sure how > to > > integrate it with yt > > > > You can use the sphere callback as-is, you just need to pass in the > necessary circle_args keyword argument to tell matplotlib to fill in the > circle. > > Take a look at this notebook: > https://gist.github.com/ngoldbaum/c6e0a4ce4c399e2f5bace6033ced2146 > > circle_args should take any keyword arguments accepted by the matplotlib > circle patch object: > > https://matplotlib.org/api/_as_gen/matplotlib.patches.Circle.html > > Hope that helps, > > Nathan > > > > > > Any suggestions > > Thank you in advance > > > > #----------------------------------------------------------- > > --------------- > > > > def create_sphere_coords(radius=10): > > """ > > function just returns a set of x,y,z coordinates for generating > > a sphere. > > """ > > r = radius > > pi = np.pi > > cos = np.cos > > sin = np.sin > > phi, theta = np.mgrid[0:pi:101j, 0:2 * pi:101j] > > > > x = r * sin(phi) * cos(theta) > > y = r * sin(phi) * sin(theta) > > z = r * cos(phi) > > > > return (x,y,z) > > > > #----------------------------------------------------------- > > --------------- > > > > --000000000000c1fdcc0572516664 > Content-Type: text/html; charset="UTF-8" > Content-Transfer-Encoding: quoted-printable > > <div dir=3D"ltr"><br><div class=3D"gmail_extra"><br><div > class=3D"gmail_quo= > te">On Tue, Jul 31, 2018 at 3:03 PM, Sushilkumar <span dir=3D"ltr"><<a > h= > ref=3D"mailto:sushil.sush19us@gmail.com" > target=3D"_blank">sushil.sush19us@= > gmail.com</a>></span> wrote:<br><blockquote class=3D"gmail_quote" > style= > =3D"margin:0px 0px 0px 0.8ex;border-left:1px solid > rgb(204,204,204);padding= > -left:1ex"><div dir=3D"ltr"><div>Dear yt:</div><div><br></div><div>Is > there= > a way to annotate sphere to a plot? I used annotate_sphere but it plots > an= > outline of a circle but I would like a "filled circle" to > repres= > ent a star? <br></div><div><br></div><div>I have a python script (see > below= > ) to create a sphere, but not sure how to integrate it with > yt<br></div></d= > iv></blockquote><div><br></div><div>You can use the sphere callback as-is, > = > you just need to pass in the necessary circle_args keyword argument to > tell= > matplotlib to fill in the circle.</div><div><br></div><div>Take a look at > = > this notebook: <a href=3D" > https://gist.github.com/ngoldbaum/c6e0a4ce4c399e2= > f5bace6033ced2146 > <https://gist.github.com/ngoldbaum/c6e0a4ce4c399e2=f5bace6033ced2146>"> > https://gist.github.com/ngoldbaum/c6e0a4ce4c399e2f5bace6= > 033ced2146 > <https://gist.github.com/ngoldbaum/c6e0a4ce4c399e2f5bace6=033ced2146></a></div><div><br></div><div>circle_args > should take any keyword= > arguments accepted by the matplotlib circle patch > object:</div><div><br></= > div><div><a href=3D" > https://matplotlib.org/api/_as_gen/matplotlib.patches.C= > ircle.html > <https://matplotlib.org/api/_as_gen/matplotlib.patches.C=ircle.html>"> > https://matplotlib.org/api/_as_gen/matplotlib.patches.Circle.ht= > ml</a></div><div><br></div><div>Hope that > helps,</div><div><br></div><div>N= > athan<br></div><div>=C2=A0</div><blockquote class=3D"gmail_quote" > style=3D"= > margin:0px 0px 0px 0.8ex;border-left:1px solid > rgb(204,204,204);padding-lef= > t:1ex"><div dir=3D"ltr"><div></div><div><br></div><div>Any suggestions > <br>= > </div><div>Thank you in > advance</div><div><br></div><div>#-----------------= > > ------------<wbr>------------------------------<wbr>---------------<br><br>= > def create_sphere_coords(radius=3D<wbr>10):<br>=C2=A0=C2=A0=C2=A0 > "&qu= > ot;"<br>=C2=A0=C2=A0=C2=A0 function just returns a set of x,y,z > coordi= > nates for generating <br>=C2=A0=C2=A0=C2=A0 a > sphere.<br>=C2=A0=C2=A0=C2=A0= > """<br>=C2=A0=C2=A0=C2=A0 r =3D radius<br>=C2=A0=C2=A0=C2= > =A0 pi =3D np.pi<br>=C2=A0=C2=A0=C2=A0 cos =3D > np.cos<br>=C2=A0=C2=A0=C2=A0= > sin =3D np.sin<br>=C2=A0=C2=A0=C2=A0 phi, theta =3D np.mgrid[0:pi:101j, > 0:= > 2 * pi:101j]<br>=C2=A0=C2=A0=C2=A0 <br>=C2=A0=C2=A0=C2=A0 x =3D r * > sin(phi= > ) * cos(theta)<br>=C2=A0=C2=A0=C2=A0 y =3D r * sin(phi) * > sin(theta)<br>=C2= > =A0=C2=A0=C2=A0 z =3D r * cos(phi)<br><br>=C2=A0=C2=A0=C2=A0 return > (x,y,z)= > > <br><br>#-----------------------------<wbr>------------------------------<w= > br>---------------<br></div></div></blockquote></div></div></div> > > --000000000000c1fdcc0572516664-- > > ------------------------------ > > Subject: Digest Footer > > _______________________________________________ > yt-users mailing list -- yt-users(a)python.org > To unsubscribe send an email to yt-users-leave(a)python.org > > > ------------------------------ > > End of yt-users Digest, Vol 125, Issue 19 > ***************************************** > -- *SK2* *"**Claiming that something can move faster than light is a good conversation-stopper in physics. People edge away from you in cocktail parties; friends never return phone calls. You just don’t mess with Albert Einstein.**"*

1 0