Splitting mixed functions with deep copy fails for periodic bcs

Below code reproduces crash

from dolfin import *
class PeriodicBoundary(SubDomain):
    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]

mesh = UnitSquareMesh(32, 32)
V = FunctionSpace(mesh, "CG", 1, constrained_domain=PeriodicBoundary())
W = V*V
u = Function(W)

u1 = u.split(deepcopy=True)[0]

Comments (3)