1 |
# $Id$ |
# $Id$ |
2 |
import esys.finley |
from esys.finley import Rectangle |
3 |
from esys.linearPDEs import Poisson |
from esys.linearPDEs import Poisson |
4 |
# generate domain: |
# generate domain: |
5 |
mydomain = esys.finley.Rectangle(l0=1.,l1=1.,n0=40, n1=20) |
mydomain = Rectangle(l0=1.,l1=1.,n0=40, n1=20) |
6 |
# define characteristic function of Gamma^D |
# define characteristic function of Gamma^D |
7 |
x = mydomain.getX() |
x = mydomain.getX() |
8 |
gammaD = x[0].whereZero()+x[1].whereZero() |
gammaD = x[0].whereZero()+x[1].whereZero() |