Segfault when using grid function in ufl boundary integrals
@gridFunction(mesh, name="u_gamma", order=1)
def u_gamma(xg):
return q_interp(xg[0])
eps = 1e-8
bc_top = DirichletBC(space, Constant(310.), x[1] > - eps)
k = Constant(1)
dt = Constant(0.01, name='dt')
## WORKS
A = u * v * ufl.dx + dt * k * ufl.dot(ufl.grad(u), ufl.grad(v)) * ufl.dx
b = uold * v * ufl.dx - u_gamma * v * ufl.dx
scheme = solutionScheme([A == b, bc_top], solver='gmres', parameters={})
scheme.solve(target=unew)
print("Solved first step")
## DOES NOT WORK !!! Using u_gamma on the boundary
A = u * v * ufl.dx + dt * k * ufl.dot(ufl.grad(u), ufl.grad(v)) * ufl.dx
b = uold * v * ufl.dx - u_gamma * v * ufl.ds(3)
scheme = solutionScheme([A == b, bc_top], solver='gmres', parameters={})
scheme.solve(target=unew)
```[bndbug.py](/uploads/48edb48b3c8d8641c88557084ccddd44/bndbug.py)
Edited by Robert K