dolfin team mailing list archive

[Joachim Haga <jobh@xxxxxxxxx>]

** Attachment added: "0006-Store-reference-to-DirichletBC-input-parameters-to-a.patch"
   https://bugs.launchpad.net/bugs/747289/+attachment/1963779/+files/0006-Store-reference-to-DirichletBC-input-parameters-to-a.patch

-- 
You received this bug notification because you are a member of DOLFIN
Team, which is subscribed to DOLFIN.
https://bugs.launchpad.net/bugs/747289

Title:
  Python: DirichletBC arguments are deleted early

Status in DOLFIN:
  New

Bug description:
  The following code fails with a strange error:
  -----
  from dolfin import *
  class Boundary(SubDomain):
      def inside(self, x, on_boundary):
          return on_boundary
  class BoundaryFunction(Expression):
      def eval(self, values, x):
          values[0] = 1 if near(x[1],1) else 0
  mesh = UnitSquare(2,2)
  V = FunctionSpace(mesh, "CG", 1)
  v,u = TestFunction(V), TrialFunction(V)
  A = assemble(v*u*dx)

  bc = DirichletBC(V, BoundaryFunction(), Boundary())
  bc.apply(A)
  -----
  Replacing the bc = line as follows also fails, with a different strange error:
  -----
  bf = BoundaryFunction()
  bc = DirichletBC(V, BoundaryFunction(), Boundary())
  ----
  Only this works:
  -----
  bf = BoundaryFunction()
  bc = DirichletBC(V, bf, Boundary())
  -----

  I have attached a patch that works around this by saving a reference
  to the arguments in DirichletBC. It's not the proper fix of course,
  but I have no time to find the root cause.