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
_______________________________________________
Mailing list: https://launchpad.net/~ffc
Post to : [email protected]
Unsubscribe : https://launchpad.net/~ffc
More help : https://help.launchpad.net/ListHelp
