1 
# $Id$ 
# $Id$ 
2 
from esys.finley import Rectangle 
from esys.finley import Rectangle 
3 
from esys.linearPDEs import Poisson 
from esys.escript.linearPDEs import Poisson 
4 

from esys.escript import whereZero,saveVTK 
5 
# generate domain: 
# generate domain: 
6 
mydomain = Rectangle(l0=1.,l1=1.,n0=40, n1=20) 
mydomain = Rectangle(l0=1.,l1=1.,n0=40, n1=20) 
7 
# define characteristic function of Gamma^D 
# define characteristic function of Gamma^D 
8 
x = mydomain.getX() 
x = mydomain.getX() 
9 
gammaD = x[0].whereZero()+x[1].whereZero() 
gammaD = whereZero(x[0])+whereZero(x[1]) 
10 
# define PDE and get its solution u 
# define PDE and get its solution u 
11 
mypde = Poisson(domain=mydomain,f=1,q=gammaD) 
mypde = Poisson(domain=mydomain) 
12 

mypde.setValue(f=1,q=gammaD) 
13 
u = mypde.getSolution() 
u = mypde.getSolution() 
14 
# write u to an external file 
# write u to an external file 
15 
u.saveDX("u.dx") 
saveVTK("u.xml",solution=u) 