File: lorentz.py

package info (click to toggle)
pyopengl 2.0.1.08-5.1
  • links: PTS
  • area: main
  • in suites: sarge
  • size: 19,484 kB
  • ctags: 9,036
  • sloc: pascal: 64,950; xml: 28,088; ansic: 20,696; python: 19,761; tcl: 668; makefile: 240; sh: 25
file content (82 lines) | stat: -rwxr-xr-x 1,869 bytes parent folder | download | duplicates (2)
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
#!

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


import string
__version__ = string.split('$Revision: 1.5 $')[1]
__date__ = string.join(string.split('$Date: 2001/08/30 19:17:17 $')[1:3], ' ')
__author__ = 'Tarn Weisner Burton <twburton@users.sourceforge.net>'

import sys
try:
    from Numeric import *
except:
    print "This Demo requires the Numeric extension, sorry."
    sys.exit()
from OpenGL.GL import *
from OpenGL.Tk import *

n, dt = 2000, 0.01
x, y, z = 0.01, 0.01, 0.01
frac = -1.0 * (8.0/3.0)
a = array((n,3), 'd')

def lorentz(o, x, y, z, n=2000, dt=0.01):
  """Generate Lorentz attractor.  Put graphic in a graphical object"""

  o.grob = glGenLists(10);
  glNewList(o.grob, GL_COMPILE);
  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
    glVertex3f(x, y, z)
  glEnd()
 
  glEnable(GL_LIGHTING)
  glEndList()

def redraw(o):
  """The main scene redraw function."""

  # Clear the background and depth buffer.
  glClearColor(1., 0., 1., 0.)
  glClear(GL_COLOR_BUFFER_BIT | GL_DEPTH_BUFFER_BIT)
  glColor3f(1., 1., 1.)
  glCallList(o.grob)

#
# Demo starts here really.

import Tkinter, sys

# Create the opengl widget here.

o = Opengl(None, width = 400, height = 400, double = 1)

# Register the redraw procedure for the widget.

o.redraw = redraw

o.pack(side = 'top', expand = 1, fill = 'both')
o.set_centerpoint(0., 0., 2000.)
o.set_eyepoint(13000.)

o.far = 15000.

lorentz(o, 0.01, 0.01, 0.01)

# Enter the tk mainloop.

Tkinter.mainloop()