On Thu, Nov 26, 2009 at 03:43:51PM +0100, Kristian Oelgaard wrote: > > > On Wed, Nov 25, 2009 at 10:56 PM, Anders Logg <[email protected]> wrote: > > I think I've tracked down the bug but don't know how it should be > handled. We need input from Kristian on this. > > I've added some debug printing in FFC that should make things clear > when trying to compile the following simple form with FFC: > > P1 = FiniteElement("Lagrange", triangle, 1) > Q1 = FiniteElement("Quadrature", triangle, 2) > > v = TestFunction(P1) > g = Function(Q1) > > L = v*g*ds > > What happens is that FFC creates quadrature elements in compiler.py by > calling create_element. The constructor of QuadratureElement then > calls make_quadrature in FIAT just to get the number of quadrature > points. In this call, it uses the cell, not the facet although it > should be integrated on the facet. > > Then later make_quadrature is also called from quadraturerepresentation.py > (strange that this is done twhice), but that time with the facet > argument. > > The results is that in the first call to make_quadrature, the number > of points is 4, but in the second case only 2. > > Kristian, could you take a quick look? > > > Yes, I'll take a look. > > Kristian
Great! -- Anders
signature.asc
Description: Digital signature
_______________________________________________ Mailing list: https://launchpad.net/~dolfin Post to : [email protected] Unsubscribe : https://launchpad.net/~dolfin More help : https://help.launchpad.net/ListHelp
