Streamlines in FLASH
Hi yt developers,
I'm trying to produce 3D streamlines using the AMR data from FLASH following the yt-project tutorial http://yt-project.org/doc/visualizing/streamlines.html. However, the execution was very slow (in particular, at the 'Streamlines' call). The code didn't complete one single streamline after an hour.
Here's the excerpt of my code: ————— c = ds.domain_center
streamlines = Streamlines(ds, c, 'velx', 'vely', 'velz', length=1.0*pc, get_magnitude=False) streamlines.integrate_through_volume()
fig=plt.figure() ax = Axes3D(fig) for stream in streamlines.streamlines: stream = stream[np.all(stream != 0.0, axis=1)] ax.plot3D(stream[:,0], stream[:,1], stream[:,2], alpha=0.1)
plt.savefig('streamlines.png') —————
I have two thoughts: 1. My FLASH file is quite large (~10GB). It may just be slow because of the large AMR grid size. Are there any ways to read only the velocities off the grid or only a subsection of the grid? 2. FLASH's octree data structure may be incompatible with the AMRKDTree construction (first step in streamlines creation). In the doc string of AMRKDTree, it says "Not applicable to particle or octree-based datasets.". If it's the case, I may have to find another way out.
Thank you very much for your time! Have a great weekend!
Best, Benny
Sorry everyone, you can ignore this last message. It is a duplicate that I mistakenly released from the spam filter.
On Fri, Jun 30, 2017 at 11:58 AM, Benny Tsang bthtsang@astro.as.utexas.edu wrote:
Hi yt developers,
I'm trying to produce 3D streamlines using the AMR data from FLASH following the yt-project tutorial http://yt-project.org/doc/visualizing/ streamlines.html. However, the execution was very slow (in particular, at the 'Streamlines' call). The code didn't complete one single streamline after an hour.
Here's the excerpt of my code: ————— c = ds.domain_center
streamlines = Streamlines(ds, c, 'velx', 'vely', 'velz', length=1.0*pc, get_magnitude=False) streamlines.integrate_through_volume()
fig=plt.figure() ax = Axes3D(fig) for stream in streamlines.streamlines: stream = stream[np.all(stream != 0.0, axis=1)] ax.plot3D(stream[:,0], stream[:,1], stream[:,2], alpha=0.1)
plt.savefig('streamlines.png') —————
I have two thoughts:
- My FLASH file is quite large (~10GB). It may just be slow because of
the large AMR grid size. Are there any ways to read only the velocities off the grid or only a subsection of the grid? 2. FLASH's octree data structure may be incompatible with the AMRKDTree construction (first step in streamlines creation). In the doc string of AMRKDTree, it says "Not applicable to particle or octree-based datasets.". If it's the case, I may have to find another way out.
Thank you very much for your time! Have a great weekend!
Best, Benny
yt-dev mailing list yt-dev@lists.spacepope.org http://lists.spacepope.org/listinfo.cgi/yt-dev-spacepope.org
participants (2)
-
Benny Tsang -
Britton Smith