pyvisi/examples/escript_scalarData3D.py

#!/usr/bin/env python


# $Id: escript_scalarData3D.py,v 1.3 2006/01/03 08:46:20 paultcochrane Exp $

from esys.escript import *
from esys import bruce

import numarray

from pyvisi import *
from pyvisi.renderers.vtk import *

brickDomain = bruce.Brick(9,9,9,10,10,10)
brickFunctionSpace = escript.ContinuousFunction(brickDomain)

domainData = brickFunctionSpace.getX()

# plotting scalar data in a 3D array
scalarData3D = sin(domainData[0])

scene = Scene()
plot = IsosurfacePlot(scene)
plot.setData(scalarData3D)
scene.render(pause=True)

scene.save(fname="escript_scalarData3D_isosurfacePlot.png", format="png")

# vim: expandtab shiftwidth=4:
1	cochrane	417	#!/usr/bin/env python
2
3
4			# $Id: escript_scalarData3D.py,v 1.3 2006/01/03 08:46:20 paultcochrane Exp $
5
6			from esys.escript import *
7			from esys import bruce
8
9			import numarray
10
11			from pyvisi import *
12			from pyvisi.renderers.vtk import *
13
14			brickDomain = bruce.Brick(9,9,9,10,10,10)
15			brickFunctionSpace = escript.ContinuousFunction(brickDomain)
16
17			domainData = brickFunctionSpace.getX()
18
19			# plotting scalar data in a 3D array
20			scalarData3D = sin(domainData[0])
21
22			scene = Scene()
23			plot = IsosurfacePlot(scene)
24			plot.setData(scalarData3D)
25			scene.render(pause=True)
26
27			scene.save(fname="escript_scalarData3D_isosurfacePlot.png", format="png")
28
29			# vim: expandtab shiftwidth=4: