On Mon, Feb 7, 2011 at 5:02 AM, Paul Leopardi wrote:
Hi all,
On Sun, 6 Feb 2011 03:54:48 PM Paul Leopardi wrote:
> I'm having trouble using multiple figures with mplot3d.

I have appended an entire example script, below.

The script incrementally plots 3 curves, one in each of 3 figure windows. The
trouble is, once Figure 2 has finished plotting, the curve for Figure 1
disappears and is replaced by the curve for Figure 2, with the axes for Figure
1; once Figure 3 has finished plotting, the curves for Figures 1 and 2
disappear and are replaced by the curve for Figure 3, with the axes for Figure
1 and Figure 2, respectively.

The original code was written with incremental plotting because the points
took a long time to calculate. Without incremental plotting, the figures
stayed blank for a long time. The script below is very similar to my original
script, but does not depend on my GluCat library.

Best, Paul
---

# -*- coding: utf-8 -*-

# Imports needed for array calculation and plotting.
from numpy import array, floor, random, empty, cos, pi
from mpl_toolkits.mplot3d import Axes3D
import matplotlib.pyplot as plt

# Constants to control the plotting.
C=3 # Number of curves to plot.
P=1000 # Number of points overall.
R=2 # Scaling constant to use.
N=25 # Number of points in a curve segment.
M=P/N

# Array of points.
x=empty((3,P))
rgb=empty((3))

# Plot C curves.
for i in xrange(0,C):
# Initial point.
x0=random.randn(3)

# Plot a curve using a random bivector in R_{5,0}
# with appropriate scaling.
w=random.randn(3) * 2*pi*R/P

# Use a new figure for each curve.
fig=plt.figure(figsize=(15,12))
# ax=Axes3D(fig)
ax = fig.gca(projection='3d')
plt.show()

# Coordinate limits to determine the colour of the first curve segment.
minx=array([-x0[0],x0[1],-x0[2]])
maxx=minx.copy()

# Split the curve into M segments, each with an appropriate colour.
for j in range(0,M):

# Find N points forming a curve segment by
# exponentiating w*k for k from j*N to (j+1)*N-1.
abot=j*N
atop=abot+N
for k in xrange(abot,atop):
for h in range(0,3):
x[h,k]=x0[h]+cos(w[h]*k)

# Determine the colour of the curve segment.
amid=floor((abot+atop)/2)
for h in range(0,3):
sign=(-1)**(h+1)
minx[h]=min(minx[h],min(sign*x[h,abot:atop]))
maxx[h]=max(maxx[h],max(sign*x[h,abot:atop]))
rgb[h]=max(0.0,min((sign*x[h,amid]-minx[h])/(maxx[h]-minx[h]),1.0))

# Plot the curve segment using the chosen colour.
alow=(abot-1 if j>0 else abot)
ax.plot(x[0,alow:atop],x[1,alow:atop],x[2,alow:atop],c=rgb.tolist())
plt.draw()
plt.show()

Paul,

As formatted, the code would not run.  I presume that everything after "for j in range(0,M):" should be indented?  When I did that and ran it in ipython, I could not reproduce your problem.  What version of matplotlib are you running?

Ben Root