Optimize speed of plotting mesh
- To: mathgroup at smc.vnet.net
- Subject: [mg99712] Optimize speed of plotting mesh
- From: "Gehricht at googlemail.com" <Gehricht at googlemail.com>
- Date: Tue, 12 May 2009 03:07:23 -0400 (EDT)
Dear MathGroup!
I have to bring coordinates with the corresponding mesh (at each
coordinate) together. The question is not, how do it but how to do it
more efficiently. The following steps are necessary
a) Generation of the coordinates from a dx in each spatial direction.
At the moment I am doing this with table
coordinates=Table[i*dx+j*dy+k*dz,{i,0,imax},{j,0,jmax},{k,0,kmax}]
which seems to be the true bottleneck (takes 24 sec's for 2299968
Points)
b) Bring coordinates and mesh together
(Flatten[#1]) & /@ Transpose[{coordinates, mesh}]
Takes about 6 secs
c) Delete Cases, in which the mesh data is lower than a certain value
choppedData=DeleteCases[data, {_, _, _, x_ /; x <= 10^-4}]
Also takes about 6 secs. Might be unnecessary, but I thought it is a
good idea to decrease the number of mesh points for the CounterPlot.
d) Plot the mesh
ListContourPlot3D[choppedData, Mesh -> None, Contours -> 5,
PlotRange -> All, MaxPlotPoints -> 50,
ContourStyle -> Yellow] // Timing
I guess, this can be improved by Interpolation of the data, though I
have not tried it yet. I have to use a large value for MaxPlotPoints,
because otherwise the mesh does not look smooth.
With thanks for your help
Yours Wolfgang