dolfin team mailing list archive

["Garth N. Wells" <gnw20@xxxxxxxxx>]

Pietro Maximoff wrote:
> Question #95954 on DOLFIN changed:
> https://answers.launchpad.net/dolfin/+question/95954
> 
>     Status: Answered => Open
> 
> Pietro Maximoff is still having a problem:
> I intend to use different theta values but in this particular instance,
> theta = 1.
> 
> You're right,  L should read:
> 
> L = ((u_old + dt*f)*v - (1-theta)*dt*kappa*inner(grad(u_old),
> grad(v)))*dx
> 
> Even then, kappa is absent from the LinearForm class and hence, can't be
> set by
> 
> L.kappa = kappa
> 
> I wonder why this is the case.
> 

(1-theta) = 0 if theta=1 -> (1-theta)*dt*kappa*inner(grad(u_old) = 0 ->
UFL throws it away.

Garth