18 Jan
2016
18 Jan
'16
2:28 p.m.
On 18/01/16 14:07, Angwenyi David wrote:
f = Function(V).interpolate(Expression("x[0]*x[1]*x[2]*t"))
This looks wrong. Expressions by default have access to the coordinates (x[0..2]), but if you want to provide other values, you need to do so like so: expr = Expression("x[0]*x[1]*x[2]*t", t=value_for_t) f = Function(V).interpolate(expr) Although it looks like you don't use the function "f" later at all. note that if you want to change the value used for t in this expression, you will then also need to re-interpolate it: while t < T: expr.t = new_value_for_t f.interpolate(expr) Cheers, Lawrence