On 25 March 2010 16:44, Mehdi Nikbakht <[email protected]> wrote:
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;
Obviously, I pushed a fix. Kristian
Mehdi _______________________________________________ Mailing list: https://launchpad.net/~ffc Post to : [email protected] Unsubscribe : https://launchpad.net/~ffc More help : https://help.launchpad.net/ListHelp
signature.asc
Description: OpenPGP digital signature
_______________________________________________ Mailing list: https://launchpad.net/~ffc Post to : [email protected] Unsubscribe : https://launchpad.net/~ffc More help : https://help.launchpad.net/ListHelp
