ffc team mailing list archive

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

I don't know if it is a bug in the new quadrature code, it depends on how you look at it. :)
The generated code in tabulate_tensor is exactly identical (apart from formatting of float values and some comments).
However, this means that I no longer reset the values of A before going crazy with the += operator. This might explain the different results. I guess I should add back the reset A code since we decided to move away form allowing different representations on the same subdomain.

Kristian

On 29 January 2010 21:42, Anders Logg <logg@xxxxxxxxx> wrote:
> I've tracked down the error in the DOLFIN fem unit test to the
> following simple test:
>
>  from dolfin import *
>
>  mesh = UnitSquare(1, 1)
>  V = FunctionSpace(mesh, "DG", 1)
>
>  v = TestFunction(V)
>  u = TrialFunction(V)
>  h = CellSize(mesh)
>
>  a = (1/h)*v*u*ds
>
>  A = assemble(a)
>  info(A, True)
>
> With the old FFC, the result is
>
> row 0: (0, 0.235702)  (1, 0.117851)  (2, 0)
> row 1: (0, 0.117851)  (1, 0.707107)  (2, 0.235702)
> row 2: (0, 0)  (1, 0.235702)  (2, 0.471405)
> row 3: (3, 0.707107)  (4, 0.353553)  (5, 0)
> row 4: (3, 0.353553)  (4, 1.64992)  (5, 0.471405)
> row 5: (3, 0)  (4, 0.471405)  (5, 0.942809)
>
> With the new FFC, the result is instead
>
> row 0: (0, 0.235702)  (1, 0.117851)  (2, 0)
> row 1: (0, 0.117851)  (1, 0.471405)  (2, 0.117851)
> row 2: (0, 0)  (1, 0.117851)  (2, 0.235702)
> row 3: (3, 0.235702)  (4, 0.117851)  (5, 0)
> row 4: (3, 0.117851)  (4, 0.471405)  (5, 0.117851)
> row 5: (3, 0)  (4, 0.117851)  (5, 0.235702)
>
> --
> Anders
>
> -----BEGIN PGP SIGNATURE-----
> Version: GnuPG v1.4.9 (GNU/Linux)
>
> iEYEARECAAYFAktjSCkACgkQTuwUCDsYZdFBmwCdFX69dDe+T61eY1XePeKmGVaz
> wMoAn0nyAR4BSNj6MyplYCOsNq88ZTV/
> =xJMD
> -----END PGP SIGNATURE-----
>
> _______________________________________________
> 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