I'm trying to use periodic boundary conditions within a Raviart-Thomas finite element space in Fenics (dolfin 1.2.0) in a Ubuntu 12.04 (amd64) machine (with python 2.7). If other FE space is used, everything works fine, but if a Raviart-Thomas is considered, python crashes (due to libdolfin errors). In what follows, it is attached a test code to reproduce the error:
from dolfin import * # Periodic boundary condition class class PeriodicWall(SubDomain): def inside(self, x, on_boundary): return bool(on_boundary and near(x, 0.5)) def map(self, x, y): y = x y = x + 1.0 # Mesh and periodic boundary conditions mesh = RectangleMesh(-0.5,-0.5,0.5,0.5,5,5,"left") periodicwall = PeriodicWall() # Raviart-Thomas space with periodic boundaries V = FunctionSpace(mesh, "RT", 1, constrained_domain=periodicwall)
Any kind of help would be helpful: are there any errors in the test code? or is it a known bug in Dolfin?