ffc team mailing list archive

[Anders Logg <logg@xxxxxxxxx>]

On Thu, Jan 28, 2010 at 10:16:44AM +0100, Kristian Oelgaard wrote:
>
>
> On 27 January 2010 23:39, Anders Logg <logg@xxxxxxxxx> wrote:
> >Some results from tabulate_tensor seemed to be off by a factor 2 (in 2D).
> >I tracked this down to a scaling issue where we need to scale the
> >weights on a facet such that the sum of the weights is 1.
> >
> >I have added this in monomialintegration.py where I multiply by 2 in
> >2D and by 6 in 3D.
>
> You don't need any scaling, but if we did shouldn't it be 8 in 3D? We also had a scaling factor of 0.125 in the old FFC.
> This is why some results are 'still' wrong for quadrature I think.

You're correct. I hadn't pulled the FIAT update.

Looks much better now! :-)

I think the scaling should have been a factor 6 which is the inverse
area of the reference tetrahedron but I guess it's not important.

--
Anders

Attachment: signature.asc
Description: Digital signature