Richard.L February 2016
### Generating 3D graphic by PyOpenGL

I'm using PyOpenGL to generate the 3D sea surface according to "2D wave equation". The main purpose is to show the dynamic graphic of "2D wave equation".But it keeps telling me this error:

```
E:\WinPython-64bit-3.4.3.5\python-3.4.3.amd64\mytest>python seawave_2d_opengl.py
Traceback (most recent call last):
File "E:\WinPython-64bit-3.4.3.5\python-3.4.3.amd64\lib\site-packages\OpenGL\GLUT\
special.py", line 130, in safeCall
return function( *args, **named )
File "seawave_2d_opengl.py", line 106, in Draw
glVertex3f(x,y,z)
File "E:\WinPython-64bit-3.4.3.5\python-3.4.3.amd64\lib\site-packages\OpenGL\p
latform\baseplatform.py", line 402, in __call__
return self( *args, **named )
ctypes.ArgumentError: argument 3: <class 'TypeError'>: wrong type
GLUT Display callback <function Draw at 0x0000000004086950> with (),{} failed: r
eturning None argument 3: <class 'TypeError'>: wrong type
E:\WinPython-64bit-3.4.3.5\python-3.4.3.amd64\mytest>
```

Here is the code:

```
from numpy import linspace,zeros,sin,pi,exp,sqrt
from OpenGL.GL import *
from OpenGL.GLU import *
from OpenGL.GLUT import *
import sys
def solver0(I, f, c, bc, Lx, Ly, nx, ny, dt, tstop, user_action=None):
dx = Lx/float(nx)
dy = Ly/float(ny)
x = linspace(0, Lx, nx+1) #grid points in x dir
y = linspace(0, Ly, ny+1) #grid points in y dir
if dt <= 0: #max time step?
dt = (1/float(c))*(1/sqrt(1/dx**2 + 1/dy**2))
Cx2 = (c*dt/dx)**2
Cy2 = (c*dt/dy)**2 #help variables
dt2 = dt**2
up = zeros((nx+1,ny+1)) #solution array
u = up.copy() #solution at t-dt
um = up.copy() #solution at t-2*dt
#set initial condition:
t =0.0
for i in range(0,nx):
for j in range(0,ny):
u[i,j] = I(x[i], y[j])
for i in range(1,nx-1):
for j in range(1,ny-1):
um[i,j] = u[i,j] + \
0.5*Cx2*(u[i-1,j] -
```

```
```

```
```

```
```### Answers

jochen February 2016
As you can guess from the stacktrace, one of the arguments (the third?!?) in `glVertex3f(x,y,z)`

has a wrong type. The discussion in the comments made clear that `z`

is a two dimensional `ndarray`

while `glVertex3f()`

expects scalars. It looks like `solver0()`

computes an array of `z`

values instead of one z-value per call.

EDIT
I now sort of understand what `solver0()`

does. The function should be documented in the book where it was printed in. Although Stackoverflow is not meant to give interpretations of copy and paste code, I'll give a little overview of what I think the function does:

- Lx and Ly give the range of all x and y
- nx and ny give the number of x and y value between 0 and Lx, Ly that are used.
- The function computes an array of x and y values at wich the z-value (
`up`

) is computed.
- It computes
`up`

for several time values from 0 to `tstop`

with a step width `dt`

.
- If a user function
`user_action`

is given, it is called after `up`

was computed. The userfunction is called with `up, x, y, t`

as arguments.

To sum things up: One call of `solver0`

computes all x, y, and z values for a given range of x and y value and a given time span with a given resolution.

```
```#### Post Status

Asked in February 2016

Viewed 2,234 times

Voted 7

Answered 1 times
#### Search

## Leave an answer

```
```

```
```