doc/examples/poisson.py

# $Id$
from esys.escript import *
from esys.escript.linearPDEs import Poisson
from esys.finley import Rectangle
# generate domain:
mydomain = Rectangle(l0=1.,l1=1.,n0=40, n1=20)
# define characteristic function of Gamma^D
x = mydomain.getX()
gammaD = whereZero(x[0])+whereZero(x[1])
# define PDE and get its solution u
mypde = Poisson(domain=mydomain)
mypde.setValue(f=1,q=gammaD)
u = mypde.getSolution()
# write u to an external file
saveVTK("u.xml",sol=u)
1	jgs	102	# $Id$
2	gross	565	from esys.escript import *
3	gross	567	from esys.escript.linearPDEs import Poisson
4	jgs	108	from esys.finley import Rectangle
5	jgs	102	# generate domain:
6	jgs	108	mydomain = Rectangle(l0=1.,l1=1.,n0=40, n1=20)
7	jgs	102	# define characteristic function of Gamma^D
8			x = mydomain.getX()
9	gross	327	gammaD = whereZero(x[0])+whereZero(x[1])
10	jgs	102	# define PDE and get its solution u
11	gross	327	mypde = Poisson(domain=mydomain)
12			mypde.setValue(f=1,q=gammaD)
13	jgs	102	u = mypde.getSolution()
14			# write u to an external file
15	gross	565	saveVTK("u.xml",sol=u)
Name	Value
svn:eol-style	native
svn:keywords	Author Date Id Revision