pyvisi/examples/escript_threeByThreeTensorData1D.py

#!/usr/bin/env python


# $Id: escript_threeByThreeTensorData1D.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 *

vectorDomain = bruce.Rectangle(9,1,10,1)
vectorFunctionSpace=escript.ContinuousFunction(vectorDomain)

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 1D array
threeByThreeTensorData1D = Data(tensorData3x3, vectorFunctionSpace, True)
scene = Scene()
plot = EllipsoidPlot(scene)
plot.setData(threeByThreeTensorData1D)
scene.render(pause=True)

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

# vim: expandtab shiftwidth=4:
1	#!/usr/bin/env python
2
3
4	# $Id: escript_threeByThreeTensorData1D.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	vectorDomain = bruce.Rectangle(9,1,10,1)
15	vectorFunctionSpace=escript.ContinuousFunction(vectorDomain)
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 1D array
20	threeByThreeTensorData1D = Data(tensorData3x3, vectorFunctionSpace, True)
21	scene = Scene()
22	plot = EllipsoidPlot(scene)
23	plot.setData(threeByThreeTensorData1D)
24	scene.render(pause=True)
25
26	scene.save(fname="escript_threeByThreeTensorData1D_ellipsoidPlot.png",
27	format="png")
28
29	# vim: expandtab shiftwidth=4: