1 |
# $Id$ |
# $Id$ |
2 |
from esys.finley import Rectangle |
from esys.escript import * |
3 |
from esys.escript.linearPDEs import Poisson |
from esys.escript.linearPDEs import Poisson |
4 |
from esys.escript import whereZero,saveVTK |
from esys.finley import Rectangle |
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 |
12 |
mypde.setValue(f=1,q=gammaD) |
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 |
saveVTK("u.xml",solution=u) |
saveVTK("u.xml",sol=u) |