Re: [yt-users] Problem in simple volume rendering (yt-users Digest, Vol 36, Issue 18)
Dear Matt,
Thank you very much for your quick reply.
1. The balck_center I found is
back_center = [-0.28306926 -0.1132277 -0.39629696]
2. The physical size I used is from -8.2605e11 cm to 8.2605e11 cm (same for x,y,z directions)
The domain is correctly printing when running my yt script
yt INFO 2011-02-24 13:08:20,186 Parameters: domain_left_edge = [ -8.26050000e+11 -8.26050000e+11 -8.26050000e+11]
yt INFO 2011-02-24 13:08:20,186 Parameters: domain_right_edge = [ 8.26050000e+11 8.26050000e+11 8.26050000e+11]
3. I set the width = 1 as used in the example code.
(all value are the same as the example code, except center, c=[0.,0.,0.] and the datafile is the Flash output file)
Thank you very much for your kindly assistance.
Best,
Kuo-Chuan
I have got some warning from yt, but I think that may not relate to my problem.
The whole output are listed as bellow:
[kpan2@taranaki yt-scripts]$ ./ShowPlots.py
yt WARNING 2011-02-24 13:21:00,791 Setting 1.0 in code units to be 1.0 cm
yt WARNING 2011-02-24 13:21:00,793 No time units. Setting 1.0 = 1 second.
yt INFO 2011-02-24 13:21:00,794 Parameters: current_time = 5635.73778789
yt ERROR 2011-02-24 13:21:00,794 Missing domain_dimensions in parameter file definition!
yt INFO 2011-02-24 13:21:00,794 Parameters: domain_left_edge = [ -8.26050000e+11 -8.26050000e+11 -8.26050000e+11]
yt INFO 2011-02-24 13:21:00,794 Parameters: domain_right_edge = [ 8.26050000e+11 8.26050000e+11 8.26050000e+11]
yt INFO 2011-02-24 13:21:00,795 Parameters: cosmological_simulation = 0.0
yt WARNING 2011-02-24 13:21:00,795 Setting 1.0 in code units to be 1.0 cm
yt WARNING 2011-02-24 13:21:00,795 No time units. Setting 1.0 = 1 second.
Warning: invalid value encountered in sqrt
Warning: invalid value encountered in sqrt
Adding gaussian at -5.8560788466 with width 0.02 and colors (0.2867691866508944, 0.015686274509803921, 0.031372549019607843, 0.01)
Adding gaussian at -4.32198464433 with width 0.02 and colors (0.23702881439314105, 0.20784313725490194, 0.41568627450980389, 0.025118864315095794)
Adding gaussian at -2.78789044206 with width 0.02 and colors (0.40392156862745099, 0.40392156862745099, 0.80784313725490198, 0.063095734448019331)
Adding gaussian at -1.25379623979 with width 0.02 and colors (0.59607843137254901, 0.59607843137254901, 0.59115561118064242, 0.15848931924611143)
Adding gaussian at 0.280297962486 with width 0.02 and colors (0.792156862745098, 0.792156862745098, 0.21568627450980413, 0.39810717055349731)
Adding gaussian at 1.81439216476 with width 0.02 and colors (0.98431372549019602, 0.98431372549019602, 0.94080651128375881, 1.0)
THIS MANY GRIDS! 0| | ETA: --:--:--
Partitioning 100% |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||| Time: 00:00:00
KC>>> back_center = [ 0.21693074 0.3867723 0.10370304]
KC>>> bricks = []
KC>>> LE = []
Traceback (most recent call last):
File "./ShowPlots.py", line 161, in <module>
test.plotVR()
File "./ShowPlots.py", line 126, in plotVR
vp.ray_cast()
File "/home/kpan2/softwares/yt-x86_64/src/yt-hg/yt/visualization/volume_rendering/software_sampler.py", line 111, in ray_cast
LE = na.array(LE) - self.back_center
ValueError: shape mismatch: objects cannot be broadcast to a single shape
Exception AttributeError: "'NoneType' object has no attribute 'chdir'" in
Send yt-users mailing list submissions to yt-users@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@lists.spacepope.org
You can reach the person managing the list at yt-users-owner@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. Problem in simple volume rendering (Kuo-Chuan Pan) 2. Re: Problem in simple volume rendering (Matthew Turk) 3. Forthon and yt installation problems. (Andrew Davis) 4. Re: Forthon and yt installation problems. (Matthew Turk)
----------------------------------------------------------------------
Message: 1 Date: Thu, 24 Feb 2011 04:12:18 -0600 From: Kuo-Chuan Pan
To: yt-users@lists.spacepope.org Subject: [yt-users] Problem in simple volume rendering Message-ID: <7301A4B8-2601-4AFB-960D-D6CAC55F6D29@astro.illinois.edu> Content-Type: text/plain; charset=us-ascii Dear Sir/Madam,
I am new in yt and trying to visualize my FLASH simulation data (3D) with volume rendering. I followed the cookbook in:
http://hg.enzotools.org/cookbook/raw-file/tip/recipes/simple_volume_renderin...
However I got an error message:
File "/home/kpan2/softwares/yt-x86_64/src/yt-hg/yt/visualization/volume_rendering/software_sampler.py", line 110, in ray_cast LE = na.array(LE) - self.back_center ValueError: shape mismatch: objects cannot be broadcast to a single shape
I looked the software_sampler.py file and found all LE, RE, self.bricks are empty lists ( = []), but self.back_center is not and correct. I guess there are some problems when setting self.bricks.
If anyone has any suggestions, I would be very happy to hear. Thank you.
Best regards, Kuo-Chuan
------------------------------
Message: 2 Date: Thu, 24 Feb 2011 07:49:45 -0500 From: Matthew Turk
To: Discussion of the yt analysis package Cc: Kuo-Chuan Pan Subject: Re: [yt-users] Problem in simple volume rendering Message-ID: Content-Type: text/plain; charset=ISO-8859-1 Hi Kuo-Chuan,
Sorry to hear it didn't work for you! Do you think you could help us out with a bit more information, and then we'll see if we can get it working for you?
* What is the exact value of back_center that you find? * What is the physical size of your domain? * What is the width you supplied to the volume_rendering object?
It might also be easier if you could try it with the "Camera" interface as described here:
http://yt.enzotools.org/doc/visualizing/volume_rendering.html#volume-renderi...
That should allow for a bit more flexibility in describing the scene, and may help us narrow down the error.
Best wishes,
Matt
On Thu, Feb 24, 2011 at 5:12 AM, Kuo-Chuan Pan
wrote: Dear Sir/Madam,
I am new in yt and trying to visualize my FLASH simulation data (3D) with volume rendering. I followed the cookbook in:
? ?http://hg.enzotools.org/cookbook/raw-file/tip/recipes/simple_volume_renderin...
However I got an error message:
?File "/home/kpan2/softwares/yt-x86_64/src/yt-hg/yt/visualization/volume_rendering/software_sampler.py", line 110, in ray_cast ? ?LE = na.array(LE) - self.back_center ValueError: shape mismatch: objects cannot be broadcast to a single shape
I looked the software_sampler.py file and found all LE, RE, self.bricks are empty lists ( = []), but self.back_center is not and correct. I guess there are some problems when setting self.bricks.
If anyone has any suggestions, I would be very happy to hear. Thank you.
Best regards, Kuo-Chuan _______________________________________________ yt-users mailing list yt-users@lists.spacepope.org http://lists.spacepope.org/listinfo.cgi/yt-users-spacepope.org
------------------------------
Message: 3 Date: Thu, 24 Feb 2011 11:44:37 -0500 From: Andrew Davis
To: yt-users@lists.spacepope.org Subject: [yt-users] Forthon and yt installation problems. Message-ID: <20110224114437.djff25pj404cc8sg@www.mail.yale.edu> Content-Type: text/plain; charset=ISO-8859-1 Hi folks,
I'm updating my yt to the latest version, and so I'm just starting from scratch in installing everything on a local cluster. I am running into errors, however with the install script (line 479). When running the script, I get the following output regarding Forthon:
Installing Forthon-0.8.4 (arguments: 'Forthon-0.8.4') Doing yt update, wiping local changes and updating to branch stable Installing yt install_script.sh: line 479: cd: yt/utilities/kdtree: No such file or directory Failure. Check /home/ajd53/YT_NEW/yt-x86_64/yt_install.log.
I haven't been able to find any directory called yt/utilities, which is probably why it can't find it. At the end of the yt_install.log file, I see the following error, which doesn't seem to be related to the install error at first glance:
abort: data/doc/install_script.sh.i@908327e16660: no match found! /home/ajd53/YT_NEW/yt-x86_64/bin/python2.6: can't open file 'setup.py': [Errno 2] No such file or directory
Any thoughts?
Thanks, andrew
************************************************************* ** Andrew J. Davis andrew.davis@yale.edu ** ** Dept. of Astronomy 203-432-5119 ** ** Yale University www.astro.yale.edu/adavis ** *************************************************************
------------------------------
Message: 4 Date: Thu, 24 Feb 2011 11:51:38 -0500 From: Matthew Turk
To: Discussion of the yt analysis package Subject: Re: [yt-users] Forthon and yt installation problems. Message-ID: Content-Type: text/plain; charset=ISO-8859-1 Hi Andrew,
Sorry about the problem -- this looks like an issue with our hg server; I don't believe it's directly with Forthon. There's something up with the webserver, so you should be able to do this:
cd /home/ajd53/YT_NEW/yt-x86_64/src rm -rf yt-hg /home/ajd53/YT_NEW/yt-x86_64/bin/hg clone --pull --uncompressed http://hg.enzotools.org/yt ./yt-hg
and that *should* get you the full set of yt, and then re-run the install script with no modifications. Please let me know if that doesn't work.
-Matt
PS The webserver for serving hg has been an issue for a while now, and I'm embarrassed by the problem. I had hoped to move us to a different host several weeks ago, but we've been trying to get the username "yt." I think I need to give up that particular goal and settle for something different. I'll have an update on this soon.
On Thu, Feb 24, 2011 at 11:44 AM, Andrew Davis
wrote: Hi folks,
I'm updating my yt to the latest version, and so I'm just starting from scratch in installing everything on a local cluster. ?I am running into errors, however with the install script (line 479). ?When running the script, I get the following output regarding Forthon:
Installing Forthon-0.8.4 (arguments: 'Forthon-0.8.4') Doing yt update, wiping local changes and updating to branch stable Installing yt install_script.sh: line 479: cd: yt/utilities/kdtree: No such file or directory Failure. ?Check /home/ajd53/YT_NEW/yt-x86_64/yt_install.log.
I haven't been able to find any directory called yt/utilities, which is probably why it can't find it. At the end of the yt_install.log file, I see the following error, which doesn't seem to be related to the install error at first glance:
abort: data/doc/install_script.sh.i@908327e16660: no match found! /home/ajd53/YT_NEW/yt-x86_64/bin/python2.6: can't open file 'setup.py': [Errno 2] No such file or directory
Any thoughts?
Thanks, andrew
************************************************************* ** ?Andrew J. Davis ? ? ? ? ? ? andrew.davis@yale.edu ? ? ?** ** ?Dept. of Astronomy ? ? ? ? ?203-432-5119 ? ? ? ? ? ? ? ** ** ?Yale University ? ? ? ? ? ? www.astro.yale.edu/adavis ?** *************************************************************
_______________________________________________ yt-users mailing list yt-users@lists.spacepope.org http://lists.spacepope.org/listinfo.cgi/yt-users-spacepope.org
------------------------------
_______________________________________________ yt-users mailing list yt-users@lists.spacepope.org http://lists.spacepope.org/listinfo.cgi/yt-users-spacepope.org
End of yt-users Digest, Vol 36, Issue 18 ****************************************
Hi Kuo-Chuan,
3. I set the width = 1 as used in the example code.
I think this is your problem. Flash uses cgs units, and so you are trying to volume render a volume of width 1 cm. This is far tinier than a single cell, and so no data (NoneType) is returned to the render. Instead, try setting width = pf.domain_right_edge - pf.domain_left_edge, assuming your parameter file is called pf. This should set the width to be the whole domain. I hope this helps, Jeff
Hi Kuo-Chuan,
You can also get the full domain by using the unit 'unitary' in the
same way you would use other units. For example,
1.0/pf['unitary']
should get you the largest box side length, so it should be the whole domain.
When I update the recipes and documentation next, I will be more clear
that "1.0" refers to a particular method of describing the domain, and
I will move it to using "unitary" everywhere instead of the implicit
domain size.
Let us know how it goes!
Best,
Matt
On Thu, Feb 24, 2011 at 2:33 PM, j s oishi
Hi Kuo-Chuan,
3. I set the width = 1 as used in the example code.
I think this is your problem. Flash uses cgs units, and so you are trying to volume render a volume of width 1 cm. This is far tinier than a single cell, and so no data (NoneType) is returned to the render.
Instead, try setting width = pf.domain_right_edge - pf.domain_left_edge, assuming your parameter file is called pf. This should set the width to be the whole domain.
I hope this helps,
Jeff _______________________________________________ yt-users mailing list yt-users@lists.spacepope.org http://lists.spacepope.org/listinfo.cgi/yt-users-spacepope.org
Dear Matt, Thank you so much for your help! Changing the width to pf.domain_right_edge -pf.domain_left_edge seems works now! but it is still calculating (> 2hrs). I thought width=1 refer to the whole domain. After the calculation finished, I will let you know whether it works or not. Thank you so much. Best, Kuo-Chuan On Feb 24, 2011, at 3:46 PM, Matthew Turk wrote:
Hi Kuo-Chuan,
You can also get the full domain by using the unit 'unitary' in the same way you would use other units. For example,
1.0/pf['unitary']
should get you the largest box side length, so it should be the whole domain.
When I update the recipes and documentation next, I will be more clear that "1.0" refers to a particular method of describing the domain, and I will move it to using "unitary" everywhere instead of the implicit domain size.
Let us know how it goes!
Best,
Matt
On Thu, Feb 24, 2011 at 2:33 PM, j s oishi
wrote: Hi Kuo-Chuan,
3. I set the width = 1 as used in the example code.
I think this is your problem. Flash uses cgs units, and so you are trying to volume render a volume of width 1 cm. This is far tinier than a single cell, and so no data (NoneType) is returned to the render.
Instead, try setting width = pf.domain_right_edge - pf.domain_left_edge, assuming your parameter file is called pf. This should set the width to be the whole domain.
I hope this helps,
Jeff _______________________________________________ yt-users mailing list yt-users@lists.spacepope.org http://lists.spacepope.org/listinfo.cgi/yt-users-spacepope.org
_______________________________________________ yt-users mailing list yt-users@lists.spacepope.org http://lists.spacepope.org/listinfo.cgi/yt-users-spacepope.org
Hi Kuo-Chuan,
We have seen really long calculation times in the past. If you use
the camera interface and specify:
no_ghost=True
when calling .camera( ... ) it should speed up the rendering time
considerably. Try the camera interface, and give a shot at the
no_ghost option, and see if that improves it?
-Matt
On Thu, Feb 24, 2011 at 5:05 PM, Kuo-Chuan Pan
Dear Matt,
Thank you so much for your help! Changing the width to pf.domain_right_edge -pf.domain_left_edge seems works now! but it is still calculating (> 2hrs).
I thought width=1 refer to the whole domain.
After the calculation finished, I will let you know whether it works or not. Thank you so much.
Best, Kuo-Chuan
On Feb 24, 2011, at 3:46 PM, Matthew Turk wrote:
Hi Kuo-Chuan,
You can also get the full domain by using the unit 'unitary' in the same way you would use other units. For example,
1.0/pf['unitary']
should get you the largest box side length, so it should be the whole domain.
When I update the recipes and documentation next, I will be more clear that "1.0" refers to a particular method of describing the domain, and I will move it to using "unitary" everywhere instead of the implicit domain size.
Let us know how it goes!
Best,
Matt
On Thu, Feb 24, 2011 at 2:33 PM, j s oishi
wrote: Hi Kuo-Chuan,
3. I set the width = 1 as used in the example code.
I think this is your problem. Flash uses cgs units, and so you are trying to volume render a volume of width 1 cm. This is far tinier than a single cell, and so no data (NoneType) is returned to the render.
Instead, try setting width = pf.domain_right_edge - pf.domain_left_edge, assuming your parameter file is called pf. This should set the width to be the whole domain.
I hope this helps,
Jeff _______________________________________________ yt-users mailing list yt-users@lists.spacepope.org http://lists.spacepope.org/listinfo.cgi/yt-users-spacepope.org
_______________________________________________ yt-users mailing list yt-users@lists.spacepope.org http://lists.spacepope.org/listinfo.cgi/yt-users-spacepope.org
_______________________________________________ yt-users mailing list yt-users@lists.spacepope.org http://lists.spacepope.org/listinfo.cgi/yt-users-spacepope.org
Dear Matt, Thanks very much for your suggestion. Using the camera interface with no_ghost=True improve the calculation time to <5min for N=512! (9 level AMR with each block contains 8^3 zones) I found for my case, I followed the cooked http://yt.enzotools.org/doc/visualizing/volume_rendering.html#volume-renderi... there is a line: cam = vr.Camera(c, L, W, (N,N), tf) but I have to put pf in the Camera class, otherwise it will not work. cam = Camera(c, L, W, (N,N), tf, pf=pf, no_ghost=True) Not sure, that is a typo in the cookbook or not. Sincerely yours, Kuo-Chuan ================================== Kuo-Chuan Pan Graduate Student Department of Astronomy University of Illinois at Urbana-Champaign Email: astroodie[at]gmail.com kpan2[at]astro.uiuc.edu On Feb 24, 2011, at 4:10 PM, Matthew Turk wrote:
Hi Kuo-Chuan,
We have seen really long calculation times in the past. If you use the camera interface and specify:
no_ghost=True
when calling .camera( ... ) it should speed up the rendering time considerably. Try the camera interface, and give a shot at the no_ghost option, and see if that improves it?
-Matt
On Thu, Feb 24, 2011 at 5:05 PM, Kuo-Chuan Pan
wrote: Dear Matt,
Thank you so much for your help! Changing the width to pf.domain_right_edge -pf.domain_left_edge seems works now! but it is still calculating (> 2hrs).
I thought width=1 refer to the whole domain.
After the calculation finished, I will let you know whether it works or not. Thank you so much.
Best, Kuo-Chuan
On Feb 24, 2011, at 3:46 PM, Matthew Turk wrote:
Hi Kuo-Chuan,
You can also get the full domain by using the unit 'unitary' in the same way you would use other units. For example,
1.0/pf['unitary']
should get you the largest box side length, so it should be the whole domain.
When I update the recipes and documentation next, I will be more clear that "1.0" refers to a particular method of describing the domain, and I will move it to using "unitary" everywhere instead of the implicit domain size.
Let us know how it goes!
Best,
Matt
On Thu, Feb 24, 2011 at 2:33 PM, j s oishi
wrote: Hi Kuo-Chuan,
3. I set the width = 1 as used in the example code.
I think this is your problem. Flash uses cgs units, and so you are trying to volume render a volume of width 1 cm. This is far tinier than a single cell, and so no data (NoneType) is returned to the render.
Instead, try setting width = pf.domain_right_edge - pf.domain_left_edge, assuming your parameter file is called pf. This should set the width to be the whole domain.
I hope this helps,
Jeff _______________________________________________ yt-users mailing list yt-users@lists.spacepope.org http://lists.spacepope.org/listinfo.cgi/yt-users-spacepope.org
_______________________________________________ yt-users mailing list yt-users@lists.spacepope.org http://lists.spacepope.org/listinfo.cgi/yt-users-spacepope.org
_______________________________________________ yt-users mailing list yt-users@lists.spacepope.org http://lists.spacepope.org/listinfo.cgi/yt-users-spacepope.org
_______________________________________________ yt-users mailing list yt-users@lists.spacepope.org http://lists.spacepope.org/listinfo.cgi/yt-users-spacepope.org
Hi Kuo-Chuan,
Thanks for the note, and I'm glad it is much faster now. Sam
Skillman, who implemented the no_ghost option (over my objection,
which I freely confess was a mistake) has gone ahead and turned it on
by default, but added a warning about it. If you end up with
artifacts, turning on no_ghost should fix them. I also think the
performance hit is worse for FLASH than for Enzo, so we don't see it
quite as often.
I've fixed the documentation as well, but not rebuilt/reuploaded the
new docs yet.
-Matt
On Thu, Feb 24, 2011 at 6:05 PM, Kuo-Chuan Pan
Dear Matt, Thanks very much for your suggestion. Using the camera interface with no_ghost=True improve the calculation time to <5min for N=512! (9 level AMR with each block contains 8^3 zones) I found for my case, I followed the cooked http://yt.enzotools.org/doc/visualizing/volume_rendering.html#volume-renderi... there is a line:
cam = vr.Camera(c, L, W, (N,N), tf)
but I have to put pf in the Camera class, otherwise it will not work.
cam = Camera(c, L, W, (N,N), tf, pf=pf, no_ghost=True)
Not sure, that is a typo in the cookbook or not.
Sincerely yours, Kuo-Chuan
================================== Kuo-Chuan Pan Graduate Student Department of Astronomy University of Illinois at Urbana-Champaign Email: astroodie[at]gmail.com kpan2[at]astro.uiuc.edu On Feb 24, 2011, at 4:10 PM, Matthew Turk wrote:
Hi Kuo-Chuan,
We have seen really long calculation times in the past. If you use the camera interface and specify:
no_ghost=True
when calling .camera( ... ) it should speed up the rendering time considerably. Try the camera interface, and give a shot at the no_ghost option, and see if that improves it?
-Matt
On Thu, Feb 24, 2011 at 5:05 PM, Kuo-Chuan Pan
wrote: Dear Matt,
Thank you so much for your help!
Changing the width to pf.domain_right_edge -pf.domain_left_edge seems works now!
but it is still calculating (> 2hrs).
I thought width=1 refer to the whole domain.
After the calculation finished, I will let you know whether it works or not.
Thank you so much.
Best,
Kuo-Chuan
On Feb 24, 2011, at 3:46 PM, Matthew Turk wrote:
Hi Kuo-Chuan,
You can also get the full domain by using the unit 'unitary' in the
same way you would use other units. For example,
1.0/pf['unitary']
should get you the largest box side length, so it should be the whole domain.
When I update the recipes and documentation next, I will be more clear
that "1.0" refers to a particular method of describing the domain, and
I will move it to using "unitary" everywhere instead of the implicit
domain size.
Let us know how it goes!
Best,
Matt
On Thu, Feb 24, 2011 at 2:33 PM, j s oishi
wrote: Hi Kuo-Chuan,
3. I set the width = 1 as used in the example code.
I think this is your problem. Flash uses cgs units, and so you are
trying to volume render a volume of width 1 cm. This is far tinier
than a single cell, and so no data (NoneType) is returned to the
render.
Instead, try setting width = pf.domain_right_edge -
pf.domain_left_edge, assuming your parameter file is called pf. This
should set the width to be the whole domain.
I hope this helps,
Jeff
_______________________________________________
yt-users mailing list
yt-users@lists.spacepope.org
http://lists.spacepope.org/listinfo.cgi/yt-users-spacepope.org
_______________________________________________
yt-users mailing list
yt-users@lists.spacepope.org
http://lists.spacepope.org/listinfo.cgi/yt-users-spacepope.org
_______________________________________________
yt-users mailing list
yt-users@lists.spacepope.org
http://lists.spacepope.org/listinfo.cgi/yt-users-spacepope.org
_______________________________________________ yt-users mailing list yt-users@lists.spacepope.org http://lists.spacepope.org/listinfo.cgi/yt-users-spacepope.org
_______________________________________________ yt-users mailing list yt-users@lists.spacepope.org http://lists.spacepope.org/listinfo.cgi/yt-users-spacepope.org
Hi Kuo-Chuan, Matt, If you end up with artifacts, turning on no_ghost should fix them. Rather, setting no_ghost=False should fix the artifacts. As Matt said, it is now True by default. Sam
I've fixed the documentation as well, but not rebuilt/reuploaded the new docs yet.
-Matt
On Thu, Feb 24, 2011 at 6:05 PM, Kuo-Chuan Pan
wrote: Dear Matt, Thanks very much for your suggestion. Using the camera interface with no_ghost=True improve the calculation time to <5min for N=512! (9 level AMR with each block contains 8^3 zones) I found for my case, I followed the cooked
http://yt.enzotools.org/doc/visualizing/volume_rendering.html#volume-renderi...
there is a line:
cam = vr.Camera(c, L, W, (N,N), tf)
but I have to put pf in the Camera class, otherwise it will not work.
cam = Camera(c, L, W, (N,N), tf, pf=pf, no_ghost=True)
Not sure, that is a typo in the cookbook or not.
Sincerely yours, Kuo-Chuan
================================== Kuo-Chuan Pan Graduate Student Department of Astronomy University of Illinois at Urbana-Champaign Email: astroodie[at]gmail.com kpan2[at]astro.uiuc.edu On Feb 24, 2011, at 4:10 PM, Matthew Turk wrote:
Hi Kuo-Chuan,
We have seen really long calculation times in the past. If you use the camera interface and specify:
no_ghost=True
when calling .camera( ... ) it should speed up the rendering time considerably. Try the camera interface, and give a shot at the no_ghost option, and see if that improves it?
-Matt
On Thu, Feb 24, 2011 at 5:05 PM, Kuo-Chuan Pan
wrote:
Dear Matt,
Thank you so much for your help!
Changing the width to pf.domain_right_edge -pf.domain_left_edge seems works now!
but it is still calculating (> 2hrs).
I thought width=1 refer to the whole domain.
After the calculation finished, I will let you know whether it works or not.
Thank you so much.
Best,
Kuo-Chuan
On Feb 24, 2011, at 3:46 PM, Matthew Turk wrote:
Hi Kuo-Chuan,
You can also get the full domain by using the unit 'unitary' in the
same way you would use other units. For example,
1.0/pf['unitary']
should get you the largest box side length, so it should be the whole domain.
When I update the recipes and documentation next, I will be more clear
that "1.0" refers to a particular method of describing the domain, and
I will move it to using "unitary" everywhere instead of the implicit
domain size.
Let us know how it goes!
Best,
Matt
On Thu, Feb 24, 2011 at 2:33 PM, j s oishi
wrote: Hi Kuo-Chuan,
3. I set the width = 1 as used in the example code.
I think this is your problem. Flash uses cgs units, and so you are
trying to volume render a volume of width 1 cm. This is far tinier
than a single cell, and so no data (NoneType) is returned to the
render.
Instead, try setting width = pf.domain_right_edge -
pf.domain_left_edge, assuming your parameter file is called pf. This
should set the width to be the whole domain.
I hope this helps,
Jeff
_______________________________________________
yt-users mailing list
yt-users@lists.spacepope.org
http://lists.spacepope.org/listinfo.cgi/yt-users-spacepope.org
_______________________________________________
yt-users mailing list
yt-users@lists.spacepope.org
http://lists.spacepope.org/listinfo.cgi/yt-users-spacepope.org
_______________________________________________
yt-users mailing list
yt-users@lists.spacepope.org
http://lists.spacepope.org/listinfo.cgi/yt-users-spacepope.org
_______________________________________________ yt-users mailing list yt-users@lists.spacepope.org http://lists.spacepope.org/listinfo.cgi/yt-users-spacepope.org
_______________________________________________ yt-users mailing list yt-users@lists.spacepope.org http://lists.spacepope.org/listinfo.cgi/yt-users-spacepope.org
_______________________________________________ yt-users mailing list yt-users@lists.spacepope.org http://lists.spacepope.org/listinfo.cgi/yt-users-spacepope.org
Dear Matt and Sam, Thanks very much for your useful suggestion. Somehow no_ghost =False will take a much longer time in my data, but the image with no_ghost=True looks fine so far. Thank you, Kuo-Chuan On Feb 24, 2011, at 6:10 PM, Sam Skillman wrote:
Hi Kuo-Chuan, Matt,
If you end up with artifacts, turning on no_ghost should fix them.
Rather, setting no_ghost=False should fix the artifacts. As Matt said, it is now True by default.
Sam
I've fixed the documentation as well, but not rebuilt/reuploaded the new docs yet.
-Matt
On Thu, Feb 24, 2011 at 6:05 PM, Kuo-Chuan Pan
wrote: Dear Matt, Thanks very much for your suggestion. Using the camera interface with no_ghost=True improve the calculation time to <5min for N=512! (9 level AMR with each block contains 8^3 zones) I found for my case, I followed the cooked http://yt.enzotools.org/doc/visualizing/volume_rendering.html#volume-renderi... there is a line:
cam = vr.Camera(c, L, W, (N,N), tf)
but I have to put pf in the Camera class, otherwise it will not work.
cam = Camera(c, L, W, (N,N), tf, pf=pf, no_ghost=True)
Not sure, that is a typo in the cookbook or not.
Sincerely yours, Kuo-Chuan
================================== Kuo-Chuan Pan Graduate Student Department of Astronomy University of Illinois at Urbana-Champaign Email: astroodie[at]gmail.com kpan2[at]astro.uiuc.edu On Feb 24, 2011, at 4:10 PM, Matthew Turk wrote:
Hi Kuo-Chuan,
We have seen really long calculation times in the past. If you use the camera interface and specify:
no_ghost=True
when calling .camera( ... ) it should speed up the rendering time considerably. Try the camera interface, and give a shot at the no_ghost option, and see if that improves it?
-Matt
On Thu, Feb 24, 2011 at 5:05 PM, Kuo-Chuan Pan
wrote: Dear Matt,
Thank you so much for your help!
Changing the width to pf.domain_right_edge -pf.domain_left_edge seems works now!
but it is still calculating (> 2hrs).
I thought width=1 refer to the whole domain.
After the calculation finished, I will let you know whether it works or not.
Thank you so much.
Best,
Kuo-Chuan
On Feb 24, 2011, at 3:46 PM, Matthew Turk wrote:
Hi Kuo-Chuan,
You can also get the full domain by using the unit 'unitary' in the
same way you would use other units. For example,
1.0/pf['unitary']
should get you the largest box side length, so it should be the whole domain.
When I update the recipes and documentation next, I will be more clear
that "1.0" refers to a particular method of describing the domain, and
I will move it to using "unitary" everywhere instead of the implicit
domain size.
Let us know how it goes!
Best,
Matt
On Thu, Feb 24, 2011 at 2:33 PM, j s oishi
wrote: Hi Kuo-Chuan,
3. I set the width = 1 as used in the example code.
I think this is your problem. Flash uses cgs units, and so you are
trying to volume render a volume of width 1 cm. This is far tinier
than a single cell, and so no data (NoneType) is returned to the
render.
Instead, try setting width = pf.domain_right_edge -
pf.domain_left_edge, assuming your parameter file is called pf. This
should set the width to be the whole domain.
I hope this helps,
Jeff
_______________________________________________
yt-users mailing list
yt-users@lists.spacepope.org
http://lists.spacepope.org/listinfo.cgi/yt-users-spacepope.org
_______________________________________________
yt-users mailing list
yt-users@lists.spacepope.org
http://lists.spacepope.org/listinfo.cgi/yt-users-spacepope.org
_______________________________________________
yt-users mailing list
yt-users@lists.spacepope.org
http://lists.spacepope.org/listinfo.cgi/yt-users-spacepope.org
_______________________________________________ yt-users mailing list yt-users@lists.spacepope.org http://lists.spacepope.org/listinfo.cgi/yt-users-spacepope.org
_______________________________________________ yt-users mailing list yt-users@lists.spacepope.org http://lists.spacepope.org/listinfo.cgi/yt-users-spacepope.org
_______________________________________________ yt-users mailing list yt-users@lists.spacepope.org http://lists.spacepope.org/listinfo.cgi/yt-users-spacepope.org
_______________________________________________ yt-users mailing list yt-users@lists.spacepope.org http://lists.spacepope.org/listinfo.cgi/yt-users-spacepope.org
Following up on this discussion, Is it possible to use the camera interface to make an off-axis projection? Thanks, John On Feb 24, 2011, at 8:18 PM, Kuo-Chuan Pan wrote:
Dear Matt and Sam,
Thanks very much for your useful suggestion. Somehow no_ghost =False will take a much longer time in my data, but the image with no_ghost=True looks fine so far.
Thank you, Kuo-Chuan
On Feb 24, 2011, at 6:10 PM, Sam Skillman wrote:
Hi Kuo-Chuan, Matt,
If you end up with artifacts, turning on no_ghost should fix them.
Rather, setting no_ghost=False should fix the artifacts. As Matt said, it is now True by default.
Sam
I've fixed the documentation as well, but not rebuilt/reuploaded the new docs yet.
-Matt
On Thu, Feb 24, 2011 at 6:05 PM, Kuo-Chuan Pan
wrote: Dear Matt, Thanks very much for your suggestion. Using the camera interface with no_ghost=True improve the calculation time to <5min for N=512! (9 level AMR with each block contains 8^3 zones) I found for my case, I followed the cooked http://yt.enzotools.org/doc/visualizing/volume_rendering.html#volume-renderi... there is a line:
cam = vr.Camera(c, L, W, (N,N), tf)
but I have to put pf in the Camera class, otherwise it will not work.
cam = Camera(c, L, W, (N,N), tf, pf=pf, no_ghost=True)
Not sure, that is a typo in the cookbook or not.
Sincerely yours, Kuo-Chuan
================================== Kuo-Chuan Pan Graduate Student Department of Astronomy University of Illinois at Urbana-Champaign Email: astroodie[at]gmail.com kpan2[at]astro.uiuc.edu On Feb 24, 2011, at 4:10 PM, Matthew Turk wrote:
Hi Kuo-Chuan,
We have seen really long calculation times in the past. If you use the camera interface and specify:
no_ghost=True
when calling .camera( ... ) it should speed up the rendering time considerably. Try the camera interface, and give a shot at the no_ghost option, and see if that improves it?
-Matt
On Thu, Feb 24, 2011 at 5:05 PM, Kuo-Chuan Pan
wrote: Dear Matt,
Thank you so much for your help!
Changing the width to pf.domain_right_edge -pf.domain_left_edge seems works now!
but it is still calculating (> 2hrs).
I thought width=1 refer to the whole domain.
After the calculation finished, I will let you know whether it works or not.
Thank you so much.
Best,
Kuo-Chuan
On Feb 24, 2011, at 3:46 PM, Matthew Turk wrote:
Hi Kuo-Chuan,
You can also get the full domain by using the unit 'unitary' in the
same way you would use other units. For example,
1.0/pf['unitary']
should get you the largest box side length, so it should be the whole domain.
When I update the recipes and documentation next, I will be more clear
that "1.0" refers to a particular method of describing the domain, and
I will move it to using "unitary" everywhere instead of the implicit
domain size.
Let us know how it goes!
Best,
Matt
On Thu, Feb 24, 2011 at 2:33 PM, j s oishi
wrote: Hi Kuo-Chuan,
3. I set the width = 1 as used in the example code.
I think this is your problem. Flash uses cgs units, and so you are
trying to volume render a volume of width 1 cm. This is far tinier
than a single cell, and so no data (NoneType) is returned to the
render.
Instead, try setting width = pf.domain_right_edge -
pf.domain_left_edge, assuming your parameter file is called pf. This
should set the width to be the whole domain.
I hope this helps,
Jeff
_______________________________________________
yt-users mailing list
yt-users@lists.spacepope.org
http://lists.spacepope.org/listinfo.cgi/yt-users-spacepope.org
_______________________________________________
yt-users mailing list
yt-users@lists.spacepope.org
http://lists.spacepope.org/listinfo.cgi/yt-users-spacepope.org
_______________________________________________
yt-users mailing list
yt-users@lists.spacepope.org
http://lists.spacepope.org/listinfo.cgi/yt-users-spacepope.org
_______________________________________________ yt-users mailing list yt-users@lists.spacepope.org http://lists.spacepope.org/listinfo.cgi/yt-users-spacepope.org
_______________________________________________ yt-users mailing list yt-users@lists.spacepope.org http://lists.spacepope.org/listinfo.cgi/yt-users-spacepope.org
_______________________________________________ yt-users mailing list yt-users@lists.spacepope.org http://lists.spacepope.org/listinfo.cgi/yt-users-spacepope.org
_______________________________________________ yt-users mailing list yt-users@lists.spacepope.org http://lists.spacepope.org/listinfo.cgi/yt-users-spacepope.org
_______________________________________________ yt-users mailing list yt-users@lists.spacepope.org http://lists.spacepope.org/listinfo.cgi/yt-users-spacepope.org
participants (5)
-
j s oishi
-
John ZuHone
-
Kuo-Chuan Pan
-
Matthew Turk
-
Sam Skillman