On 29 January 2010 20:21, Anders Logg <l...@simula.no> wrote:
On Fri, Jan 29, 2010 at 04:13:13PM +0100, Kristian Oelgaard wrote:I think it is almost there, but we need to rethink the way we manipulate the degrees of finite elements and forms, I think we have discussed this before but I don't recall if we ever reached a conclusion. Try running the QuadratureElement.ufl demo.I don't see exactly what the problem is. Is it the automatic degree selection? I've added some more debug output and the following is now printed for the QuadratureElement demo: Found quadrature element(s) with the following degree(s): [3] Selecting quadrature degree based on quadrature element: 3 quadrature_degree: auto --> 3 That looks correct to me. Shouldn't this result in the same points being chosen for both the QuadratureElement and the quadrature loop?
I see now that it is only a problem for tensor representation. For the linear form the first term v*u*dx makes it a form of degree 4, and tensor representation will use that degree for the entire form if I remember correctly? Kristian
-- Anders -----BEGIN PGP SIGNATURE----- Version: GnuPG v1.4.9 (GNU/Linux) iEYEARECAAYFAktjNTMACgkQTuwUCDsYZdHwxACcCTzBr0hzp1nXowRvQKnIsrmg /zAAn36IpBEEvbxC2lVEh1y0TP+F7eDZ =/PfD -----END PGP SIGNATURE-----
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
