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 <l...@simula.no> 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@lists.launchpad.net 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 : ffc@lists.launchpad.net Unsubscribe : https://launchpad.net/~ffc More help : https://help.launchpad.net/ListHelp
