"""GLUT replacement for the original lorentz.py demonstration code

This is the original lorentz.py demo, with a few changes to
run under GLUT and not require the TK examination widget
extensions.

Code is also broken up somewhat more, with the actual drawing
of the curve in its own function which can be called instead of
the display-list if desired.
"""

# This is statement is required by the build system to query build info
if __name__ == '__build__':
	raise Exception

__version__='$Revision: 1.1.2.1 $'[11:-2]
__date__ = '$Date: 2004/11/15 01:10:45 $'[6:-2]


from OpenGL.GL import *
from OpenGL.GLU import *
from OpenGL.GLUT import *
import time

n, dt = 2000, 0.01
x, y, z = 0.01, 0.01, 0.01
frac = -1.0 * (8.0/3.0)

def lorentz(x, y, z, n=2000, dt=0.01):
	"""Generate Lorentz attractor as a Display-list"""
	target = glGenLists(1);
	glNewList(target, GL_COMPILE);
	try:
		drawLorentz( x,y,z,n,dt )
	finally:
		glEndList()
	return target
def drawLorentz( x, y, z, n=2000, dt=0.01):
	"""Do the actual drawing & calculation of lorentz"""
	glDisable(GL_LIGHTING)
	glBegin(GL_LINE_STRIP)

	glVertex3f(x, y, z)
	frac = -1.0 * (8.0/3.0)
	for i in range(0, n):
		xp = x + (-10.0 * x * dt + 10.0 * y * dt)
		yp = y + ( 28.0 * x * dt - y * dt - x * dt * z *dt)
		zp = z + ( frac * z * dt + x * dt * y * dt)
		x=xp
		y=yp
		z=zp
		glColor3f( i/float(n),i/float(n), 0 )
		glVertex3f(x, y, z)
	glEnd()

	glEnable(GL_LIGHTING)

LORENTZ_LIST = lorentz( 0.01, 0.01, 0.01 )

def display( ):
	"""Callback function for displaying the scene

	This defines a unit-square environment in which to draw,
	i.e. width is one drawing unit, as is height
	"""
	glClearColor(0.5, 0.5, 0.5, 0)
	glClear(GL_COLOR_BUFFER_BIT | GL_DEPTH_BUFFER_BIT)

	# establish the projection matrix (perspective)
	glMatrixMode(GL_PROJECTION)
	glLoadIdentity()
	gluPerspective(
		45, # field of view in degrees
		glutGet(GLUT_WINDOW_WIDTH)/float(glutGet(GLUT_WINDOW_HEIGHT) or 1), # aspect ratio
		1, # near clipping plane
		30000, # far clipping plane
	)
	# and then the model view matrix
	glMatrixMode(GL_MODELVIEW)
	glLoadIdentity()
	gluLookAt(
		0,0,13000, # eyepoint
		0,0,2000, # center-of-view
		0,1,0, # up-vector
	)
	rotation()
	drawLorentz( 0.01, 0.01, 0.01 )
	#glCallList( LORENTZ_LIST )
	glutSwapBuffers()

def idle( ):
	glutPostRedisplay()

starttime = time.time()

def rotation( period = 10):
	"""Do rotation of the scene at given rate"""
	angle = (((time.time()-starttime)%period)/period)* 360
	glRotate( angle, 0,1,0)
	return angle

if __name__ == "__main__":
	print """You should see a curved line representing a Lorentz
attractor, rendered in 3D, rotation about the origin."""
	import sys
	glutInit(sys.argv)
	glutInitDisplayMode(GLUT_RGBA | GLUT_DOUBLE | GLUT_DEPTH)
	glutCreateWindow('Lorentz Attractor Demo')
	glutDisplayFunc(display)
	glutIdleFunc(display)
	glutMainLoop()