dolfin team mailing list archive

[Pietro Maximoff <question95954@xxxxxxxxxxxxxxxxxxxxx>]

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

Hello

I have a ufl file like so (for the Heat Equation with time-stepping via the theta-method):

element = FiniteElement("Lagrange", "triangle", 2)

v         = TestFunction(element)
u         = TrialFunction(element)
kappa = Coefficient(element)
u_old  = Coefficient(element)
f          = Coefficient(element)

dt = ....

a = u*v*dx + theta*dt*kappa*inner(grad(u), grad(v))*dx
L = ((u_old + dt*f)*v - (1-theta)*dt*kappa*inner(grad(u), grad(v)))*dx

The problem is that I can't assign kappa to L.kappa in the linear form, i.e.,

L.kappa = kappa

This is fine though for the Bilinear form.

On examining the generated header file, it turns out that the LinearForm class has no member kappa which is exactly the error message given. 

What am I doing wrong?

Best wishes 

Pietro


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