dolfin team mailing list archive

Thread
Date

Re: [Question #95954]: Absent Coefficient in Linear Form

To: dolfin@xxxxxxxxxxxxxxxxxxx
From: Garth Wells <question95954@xxxxxxxxxxxxxxxxxxxxx>
Date: Sat, 02 Jan 2010 18:30:26 -0000
Reply-to: question95954@xxxxxxxxxxxxxxxxxxxxx
Sender: bounces@xxxxxxxxxxxxx

Question #95954 on DOLFIN changed:
https://answers.launchpad.net/dolfin/+question/95954

    Status: Open => Answered

Garth Wells proposed the following answer:

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

-- 
You received this question notification because you are a member of
DOLFIN Team, which is an answer contact for DOLFIN.