| 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) |
Parent Directory
| 
Revision Log
| 
Patch