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. Is a similar scaling needed for quadrature? If so, we can have a common scaling. It looks like it's needed for quadrature. Check the results for the last tabulate_tensor in the output for Normals.ufl. -- Anders
signature.asc
Description: Digital signature
_______________________________________________ Mailing list: https://launchpad.net/~ffc Post to : [email protected] Unsubscribe : https://launchpad.net/~ffc More help : https://help.launchpad.net/ListHelp
