Using scipy special functions

[dokken]

You should create a UserExpression to do this computation:

How to difine UserExpression that takes the solution u

Consider the following: import dolfin as dl N = 100 X = 1.0 mesh = dl.IntervalMesh(N, 0, X) degree = 1 # linear basis functions Vh = dl.FunctionSpace(mesh, "Lagrange", degree) u = dl.interpolate(dl.Expression("x[0]", degree = 1), Vh) alpha = 2.0 class user_expr(dl.UserExpression): def set_alpha(self, u, alpha): # in the user_expr class self.u = u self.alpha = alpha def eval(self, value, x): # in the user_expr class "set value[0] to value at point x" va…