ffc team mailing list archive

[Kristian Oelgaard <k.b.oelgaard@xxxxxxxxx>]

On 25 March 2010 16:44, Mehdi Nikbakht <m.nikbakht@xxxxxxxxxx> 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     : ffc@xxxxxxxxxxxxxxxxxxx
> Unsubscribe : https://launchpad.net/~ffc
> More help   : https://help.launchpad.net/ListHelp

Attachment: signature.asc
Description: OpenPGP digital signature