pyvisi/examples/escript_twoByTwoTensorData2D.py

#!/usr/bin/env python


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

tensorDomain = bruce.Rectangle(9,9,10,10)
tensorFunctionSpace=escript.ContinuousFunction(tensorDomain)

tensorData2x2 = numarray.array([[1.0,2.0],[3.0,4.0]])

# plotting 2x2 tensors in a 2D array
twoByTwoTensorData2D = Data(tensorData2x2, tensorFunctionSpace, True)
scene = Scene()
plot = EllipsoidPlot(scene)
plot.setData(twoByTwoTensorData2D)
scene.render(pause=True)

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

# vim: expandtab shiftwidth=4:
1	cochrane	417	#!/usr/bin/env python
2
3
4			# $Id: escript_twoByTwoTensorData2D.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			tensorDomain = bruce.Rectangle(9,9,10,10)
15			tensorFunctionSpace=escript.ContinuousFunction(tensorDomain)
16
17			tensorData2x2 = numarray.array([[1.0,2.0],[3.0,4.0]])
18
19			# plotting 2x2 tensors in a 2D array
20			twoByTwoTensorData2D = Data(tensorData2x2, tensorFunctionSpace, True)
21			scene = Scene()
22			plot = EllipsoidPlot(scene)
23			plot.setData(twoByTwoTensorData2D)
24			scene.render(pause=True)
25
26			scene.save(fname="escript_twoByTwoTensorData2D_ellipsoidPlot.png",
27			format="png")
28
29			# vim: expandtab shiftwidth=4: