for i in range(0,3*N):
rgb.append([q1[i],q2[i],q3[i]])
for i in range(0,N):
for j in range(0,M):
for k in range(0,3):
data[3*i+k,j] = rgb[i][k];
#####
duration = float(sys.argv[1]);
import time;
K = shape(data)
K = K[0]
t0 = time.time();
t = time.time();
while t < t0 + duration:
teh_display(data);
# time.sleep(0.5);
sv = data[0,:];
for j in range(0,K-1):
data[j,:] = data[j+1,:];
data[K-1,:] = sv;
t = time.time();
if len(sys.argv) == 5:
sleep(float(sys.argv[1]))
import time;
h = 0;
s = 0;
t0 = time.time();
t = time.time();
from time import sleep as wakeup;
viz = 1j*zeros([3*N,M]);
dataR = 0.01*random.rand(N,M);
dataG = 0.01*random.rand(N,M);
dataB = 0.01*random.rand(N,M);
while t < t0 + duration:
for i in range(0,10):
dataR = dynamics(dataR);
dataG = dynamics(dataG);
dataB = dynamics(dataB);
for i in range(0,N):
for j in range(0,M):
viz[3*i,j] = (1+real(dataR[i,j]))/2;
viz[3*i+1,j] = (1+real(dataG[i,j]))/2;
viz[3*i+2,j] = (1+real(dataB[i,j]))/2;
print(viz[0:3,0])
teh_display(viz);
t = time.time();
from time import sleep as wakeup
n=n+1;
if omega != 0:
h = sin(omega*2*pi*t);
s = sin(omega*sqrt(2.0)*pi*t);
rgb = colorsys.hsv_to_rgb((h+1)/2,(s+1)/2,1);
CM = colormap.solid_colormap(rgb[0],rgb[1],rgb[2]);
tw = zeros([3*N,M]);
for j in range(0,M):
tw[:,j] = colormap.i2c(sin(data[:,j]),colormap.WHITE_COLORMAP);
Q = zeros([3*N,M]);
for j in range(0,M):
Q[:,j] = colormap.i2c(dataQ[:,j],colormap.MATLAB_COLORMAP);
teh_display(tw*Q);
t = time.time();
c = poisson.rvs(mu,size=1);
for i in range(0,c):
x = floor(N*random.rand());
y = floor(M*random.rand());
d0[x][y] = d0[x][y] + 3.14159265;
theta = theta + 8*dt;
dataQ = dynamicsQ(theta);
data = dynamics(data);
wakeup(0.005)
