ffc team mailing list archive

[Mehdi Nikbakht <m.nikbakht@xxxxxxxxxx>]

Hello,

It seems that something is going wrong in handling Exponential function.
If I compile ,

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

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

L = exp(f)*v*f*dx

the code for computing element tensor is as following,

   A[j] += FE0[ip][j]*std::exp(F0)*std::exp(F0)*W4[ip]*det;

which is obviously wrong. It should be something like,

   A[j] += FE0[ip][j]*std::exp(F0)*F0*W4[ip]*det;

Mehdi