def set_camera(self, width, height):
"""
Update the shown scene after the perspective has changed.
"""
rightt = self.mat[:3, 0]
upt = self.mat[:3, 1]
pt = self.mat[:3, 2] * self.d
t = self.mat[:3, 3]
self.proj_mat = create_perspective_projection_matrix(np.radians(self.fov), 1. * width / height, self.near, self.far)
gr3.setcameraprojectionparameters(self.fov, self.near, self.far)
self.lookat_mat = create_look_at_matrix(pt + t, t, upt)
gr3.cameralookat(pt[0] + t[0], pt[1] + t[1], pt[2] + t[2], t[0], t[1], t[2], upt[0], upt[1], upt[2])
gr.settextcolorind(2)
gr.textext(0.05, 0.08, 'y_{A}:%6.2f' % acceleration)
gr.updatews()
return
theta = 110.0 # initial angle
gamma = 0.1 # damping coefficient
L = 1 # pendulum length
t = 0
dt = 0.04
state = array([theta * pi / 180, 0])
gr3.init()
gr3.setcameraprojectionparameters(45, 1, 100)
gr3.cameralookat(0, -2, 6, 0, -2, 0, 0, 1, 0)
gr3.setbackgroundcolor(1, 1, 1, 1)
gr3.setlightdirection(1, 1, 10)
now = time.clock()
while t < 30:
start = now
t, state = rk4(t, dt, state, damped_pendulum_deriv)
theta, omega = state
acceleration = sqrt(2 * g * L * (1 - cos(theta)))
pendulum(t, theta, omega, acceleration)
now = time.clock()
def display():
global window_width, window_height, rx
gr3.setbackgroundcolor(1, 1, 1, 1)
# set up camera
gr3.setcameraprojectionparameters(45, 1, 200)
gr3.cameralookat(10 * math.cos(-rx * math.pi / 2),
10 * math.sin(-rx * math.pi / 2), 0, 0, 0, 0, 0, 0, 1)
gr3.drawimage(0, window_width, 0, window_height, window_width,
window_height, gr3.GR3_Drawable.GR3_DRAWABLE_OPENGL)
glViewport(0, 0, window_width, window_height)
glDisable(GL_LIGHTING)
glDisable(GL_DEPTH_TEST)
glMatrixMode(GL_MODELVIEW)
glPushMatrix()
glLoadIdentity()
glMatrixMode(GL_PROJECTION)
glPushMatrix()
glLoadIdentity()
glColor4f(1, 0, 0, 1)
x, y = 0, 0.22
glRasterPos2f(x * 2 - 1, y * 2 - 1)
for c in u"This is a GLUT window in which GR3 renders":
glutBitmapCharacter(GLUT_BITMAP_HELVETICA_18, ord(c))
y -= 0.05
glRasterPos2f(x * 2 - 1, y * 2 - 1)
for c in u"a scene. You can drag the mouse to rotate the":
glutBitmapCharacter(GLUT_BITMAP_HELVETICA_18, ord(c))
y -= 0.05
glRasterPos2f(x * 2 - 1, y * 2 - 1)
for c in u"molecule or right-click to open a context menu.":
glutBitmapCharacter(GLUT_BITMAP_HELVETICA_18, ord(c))
glEnable(GL_DEPTH_TEST)
glEnable(GL_LIGHTING)
glMatrixMode(GL_MODELVIEW)
glPopMatrix()
glMatrixMode(GL_PROJECTION)
glPopMatrix()
glutSwapBuffers()
gr.clearws()
selntran(0)
gr3.setquality(4)
gr3.drawimage(0, 0.5, 0.5, 1, 250, 250, gr3.GR3_Drawable.GR3_DRAWABLE_GKS)
gr.settextcolorind(1)
gr.settextfontprec(6, 0)
x = list(range(5))
y = list(range(5))
nominalWindowHeight = 500.0
pointSize = (8, 9, 10, 11, 12, 14, 18, 24, 36)
s = "i\\hbar\\frac{\\partial\\psi}{\\partial t} = \\frac{\\hbar^2}{2m}\\nabla^2\\psi + V(\\mathbf{r})\\psi"
x = 0.9
y = 0.9
gr.settextalign(3, 3)
for i in range(8):
gr.setcharheight(pointSize[i] / nominalWindowHeight)
gr.mathtex(x, y, s)
y -= 4 * pointSize[i] / nominalWindowHeight
gr.setcharheight(0.1)
gr.mathtex(0.9, 0.05, "Hello World!")
gr.settextcolorind(8)
gr.text(0.9, 0.05, "Hello World!")
gr.updatews()
gr3.setcameraprojectionparameters(45, 1, 200)
gr3.cameralookat(10 * math.cos(-rx * math.pi / 2),
10 * math.sin(-rx * math.pi / 2), 0, 0, 0, 0, 0, 0, 1)
def display():
global window_width, window_height, rx
gr3.setbackgroundcolor(1, 1, 1, 1)
# Kamera einstellen
gr3.setcameraprojectionparameters(45, 1, 200)
gr3.cameralookat(10 * math.cos(-rx * math.pi / 2),
10 * math.sin(-rx * math.pi / 2), 0, 0, 0, 0, 0, 0, 1)
gr3.drawimage(0, window_width, 0, window_height, window_width,
window_height, gr3.GR3_Drawable.GR3_DRAWABLE_OPENGL)
glViewport(0, 0, window_width, window_height)
glDisable(GL_LIGHTING)
glDisable(GL_DEPTH_TEST)
glMatrixMode(GL_MODELVIEW)
glPushMatrix()
glLoadIdentity()
glMatrixMode(GL_PROJECTION)
glPushMatrix()
glLoadIdentity()
glColor4f(1, 0, 0, 1)
x, y = 0, 0.22
glRasterPos2f(x * 2 - 1, y * 2 - 1)
for c in u"Dies ist ein GLUT-Fenster, in dem mit GR3 eine":
glutBitmapCharacter(GLUT_BITMAP_HELVETICA_18, ord(c))
y -= 0.05
glRasterPos2f(x * 2 - 1, y * 2 - 1)
for c in u"Szene gerendert wird. Mit der Maus kann man":
glutBitmapCharacter(GLUT_BITMAP_HELVETICA_18, ord(c))
y -= 0.05
glRasterPos2f(x * 2 - 1, y * 2 - 1)
for c in u"das dargestellte Molekül rotieren lassen.":
glutBitmapCharacter(GLUT_BITMAP_HELVETICA_18, ord(c))
y -= 0.04
glRasterPos2f(x * 2 - 1, y * 2 - 1)
for c in u"(Rechtsklick öffnet das Kontextmenü)":
glutBitmapCharacter(GLUT_BITMAP_HELVETICA_12, ord(c))
glEnable(GL_DEPTH_TEST)
glEnable(GL_LIGHTING)
glMatrixMode(GL_MODELVIEW)
glPopMatrix()
glMatrixMode(GL_PROJECTION)
glPopMatrix()
glutSwapBuffers()
gr.clearws()
selntran(0)
gr3.setquality(4)
gr3.drawimage(0, 0.5, 0.5, 1, 250, 250, gr3.GR3_Drawable.GR3_DRAWABLE_GKS)
gr.settextcolorind(1)
gr.settextfontprec(6, 0)
x = list(range(5))
y = list(range(5))
nominalWindowHeight = 500.0
pointSize = (8, 9, 10, 11, 12, 14, 18, 24, 36)
s = "i\\hbar\\frac{\\partial\\psi}{\\partial t} = \\frac{\\hbar^2}{2m}\\nabla^2\\psi + V(\\mathbf{r})\\psi"
x = 0.9
y = 0.9
gr.settextalign(3, 3)
for i in range(8):
gr.setcharheight(pointSize[i] / nominalWindowHeight)
gr.mathtex(x, y, s)
y -= 4 * pointSize[i] / nominalWindowHeight
gr.setcharheight(0.1)
gr.mathtex(0.9, 0.05, "Hello World!")
gr.settextcolorind(8)
gr.text(0.9, 0.05, "Hello World!")
gr.updatews()
gr3.setcameraprojectionparameters(45, 1, 200)
gr3.cameralookat(10 * math.cos(-rx * math.pi / 2),
10 * math.sin(-rx * math.pi / 2), 0, 0, 0, 0, 0, 0, 1)
def display():
global window_width, window_height, rx
gr3.setbackgroundcolor(1,1,1,1)
# Kamera einstellen
gr3.setcameraprojectionparameters(45, 1, 200)
gr3.cameralookat(10*math.cos(-rx*math.pi/2), 10*math.sin(-rx*math.pi/2), 0, 0, 0, 0, 0, 0, 1)
gr3.drawimage(0, window_width, 0, window_height, window_width, window_height, gr3.GR3_Drawable.GR3_DRAWABLE_OPENGL)
glViewport(0,0,window_width,window_height);
glDisable(GL_LIGHTING)
glDisable(GL_DEPTH_TEST)
glMatrixMode(GL_MODELVIEW)
glPushMatrix()
glLoadIdentity()
glMatrixMode(GL_PROJECTION)
glPushMatrix()
glLoadIdentity()
glColor4f(1,0,0,1)
x, y = 0, 0.22
glRasterPos2f(x*2-1,y*2-1)
for c in u"Dies ist ein GLUT-Fenster, in dem mit GR3 eine":
glutBitmapCharacter(GLUT_BITMAP_HELVETICA_18,ord(c))
y-=0.05
glRasterPos2f(x*2-1,y*2-1)
for c in u"Szene gerendert wird. Mit der Maus kann man":
glutBitmapCharacter(GLUT_BITMAP_HELVETICA_18,ord(c))
y-=0.05
glRasterPos2f(x*2-1,y*2-1)
for c in u"das dargestellte Molekül rotieren lassen.":
glutBitmapCharacter(GLUT_BITMAP_HELVETICA_18,ord(c))
y-=0.04
glRasterPos2f(x*2-1,y*2-1)
for c in u"(Rechtsklick öffnet das Kontextmenü)":
glutBitmapCharacter(GLUT_BITMAP_HELVETICA_12,ord(c))
glEnable(GL_DEPTH_TEST)
glEnable(GL_LIGHTING)
glMatrixMode(GL_MODELVIEW)
glPopMatrix()
glMatrixMode(GL_PROJECTION)
glPopMatrix()
glutSwapBuffers()
gr.clearws()
selntran(0)
gr3.setquality(4)
gr3.drawimage(0, 0.5, 0.5, 1, 250, 250, gr3.GR3_Drawable.GR3_DRAWABLE_GKS)
gr.settextcolorind(1)
gr.settextfontprec(6,0)
x = list(range(5))
y = list(range(5))
nominalWindowHeight = 500.0
pointSize = ( 8, 9, 10, 11, 12, 14, 18, 24, 36 )
s = "i\\hbar\\frac{\\partial\\psi}{\\partial t} = \\frac{\\hbar^2}{2m}\\nabla^2\\psi + V(\\mathbf{r})\\psi"
x = 0.9
y = 0.9;
gr.settextalign(3, 3)
for i in range(8):
gr.setcharheight(pointSize[i] / nominalWindowHeight)
gr.mathtex(x, y, s)
y -= 4 * pointSize[i] / nominalWindowHeight
gr.setcharheight(0.1)
gr.mathtex(0.9, 0.05, "Hello World!")
gr.settextcolorind(8)
gr.text(0.9, 0.05, "Hello World!")
gr.updatews()
gr3.setcameraprojectionparameters(45, 1, 200)
gr3.cameralookat(10*math.cos(-rx*math.pi/2), 10*math.sin(-rx*math.pi/2), 0, 0, 0, 0, 0, 0, 1)
gr.textext(0.05, 0.12, '\\omega:%7.2f' % omega)
gr.settextcolorind(2)
gr.textext(0.05, 0.08, 'y_{A}:%6.2f' % acceleration)
gr.updatews()
return
theta = 110.0 # initial angle
gamma = 0.1 # damping coefficient
L = 1 # pendulum length
t = 0
dt = 0.04
state = array([theta * pi / 180, 0])
gr3.init()
gr3.setcameraprojectionparameters(45, 1, 100)
gr3.cameralookat(0, -2, 6, 0, -2, 0, 0, 1, 0)
gr3.setbackgroundcolor(1, 1, 1, 1)
gr3.setlightdirection(1, 1, 10)
now = time.clock()
while t < 30:
start = now
t, state = rk4(t, dt, state, damped_pendulum_deriv)
theta, omega = state
acceleration = sqrt(2 * g * L * (1 - cos(theta)))
pendulum(t, theta, omega, acceleration)
now = time.clock()
