pyvisi/examples/escript_threeByThreeTensorData3D.py

#!/usr/bin/env python


# $Id: escript_threeByThreeTensorData3D.py,v 1.3 2006/01/05 00:13:50 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)

tensorData3x3 = numarray.array([[1.0,2.0,3.0],[3.0,4.0,5.0],[5.0,6.0,7.0]])

# plotting 3x3 tensors in a 3D array
threeByThreeTensorData3D = Data(tensorData3x3, brickFunctionSpace, True)
scene = Scene()
plot = EllipsoidPlot(scene)
plot.setData(threeByThreeTensorData3D)
scene.render(pause=True)

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

# vim: expandtab shiftwidth=4:
1	#!/usr/bin/env python
2
3
4	# $Id: escript_threeByThreeTensorData3D.py,v 1.3 2006/01/05 00:13:50 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	tensorData3x3 = numarray.array([[1.0,2.0,3.0],[3.0,4.0,5.0],[5.0,6.0,7.0]])
18
19	# plotting 3x3 tensors in a 3D array
20	threeByThreeTensorData3D = Data(tensorData3x3, brickFunctionSpace, True)
21	scene = Scene()
22	plot = EllipsoidPlot(scene)
23	plot.setData(threeByThreeTensorData3D)
24	scene.render(pause=True)
25
26	scene.save(fname="escript_threeByThreeTensorData3D_ellipsoidPlot.png",
27	format="png")
28
29	# vim: expandtab shiftwidth=4: