On Sat, Jan 30, 2010 at 08:32:50PM +0100, Kristian Oelgaard wrote: > > > On 30 January 2010 00:24, Anders Logg <l...@simula.no> wrote: > >On Fri, Jan 29, 2010 at 10:25:52PM +0100, Kristian Oelgaard wrote: > >> > >>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 > > > >I think that might be the problem. You need to set all entries to zero > >in case you happen not to assign to all of them in the code. DOLFIN > >does not reset the element tensor when it assembles so you need to do > >that as part of the generated code. > > Did you check this again? I have pushed a changeset that should fix this > issue. > > Kristian
Yes, it works perfect. The DOLFIN unit tests all pass again. -- Anders
signature.asc
Description: 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