19     def __init__(self, nx, ny, nz):  argument
20         self.n = np.array([nx, ny, nz], dtype='i')
40         nx, ny, nz = self.n
43         U = x.getArray(readonly=1).reshape(nx,ny,nz, order='f')
45         X, Y =  mgrid[0:1:1j*nx,0:1:1j*ny]
56         X, Y =  mgrid[0:1:1j*nx,0:1:1j*nz]
68         Z = U[nx//2,:,:]
79 def run_test(nx,ny,nz,samples,plot=False):  argument
90     x = PETSc.Vec().createSeq(nx*ny*nz)
92     app = AppCtx(nx, ny, nz)
93     f = PETSc.Vec().createSeq(nx*ny*nz)
126     nx = ny = nz = n+1  variable
127     wt = run_test(nx,ny,nz,samples,plot)
129                     % (nx,ny,nz,wt,samples))