DG PeriodicBoundary
Asked by
Miro
Hi, I tried to make DG work in periodic domain following the code in demo/undocument
The program crashes when constructing the function space. Can someone point out, what I am doing
wrong? I am running trunk version of dolfin. Thanks.
from dolfin import *
mesh = UnitSquareMesh(4,4)
class PeriodicBoundar
def inside(self, x, on_boundary):
return bool(x[0] < DOLFIN_EPS and x[0] > -DOLFIN_EPS and on_boundary)
def map(self, x, y):
y[0] = x[0] - 1.0
y[1] = x[1]
pbc = PeriodicBoundary()
V = FunctionSpace(
Question information
- Language:
- English Edit question
- Status:
- Answered
- For:
- DOLFIN Edit question
- Assignee:
- No assignee Edit question
- Last query:
- Last reply:
Can you help with this problem?
Provide an answer of your own, or ask Miro for more information if necessary.
To post a message you must log in.