ufl team mailing list archive

[Kent Andre <kent-and@xxxxxxxxx>]

On fr., 2009-02-13 at 13:32 +0100, Martin Sandve Alnæs wrote:
> On Fri, Feb 13, 2009 at 1:20 PM, Kristian Oelgaard
> <k.b.oelgaard@xxxxxxxxxx> wrote:
> > Quoting Martin Sandve Alnæs <martinal@xxxxxxxxx>:
> >
> >> On Fri, Feb 13, 2009 at 11:53 AM, Kristian Oelgaard
> >> <k.b.oelgaard@xxxxxxxxxx> wrote:
> >> >
> >> >
> >> > Hello,
> >> >
> >> > Any thoughts on how to implement a better rule for determining the
> >> polynomial
> >> > order of a form?
> >> >
> >> > as it is now the forms:
> >> >
> >> > a = v*u*dx
> >> > a = f*v*u*dx
> >> > a = f*f*f*.....*v*u*dx
> >> > a = inner(grad(v), grad(u))*dx
> >> >
> >> > all result in the same order (2 for linear basis functions)
> >>
> >> What do we want? Add the degrees of all functions and basis functions
> >> that are multiplied? Should be possible to do.

In general I don't think this is what you 'want' to do. 
It is common to use the same quadrature rules for 
a = inner(grad(v), grad(u))*dx
and 
a = f*inner(grad(v), grad(u))*dx

even though f is a function. 

It would however be nice to choose quadrature rule by setting some kind
of meta data. 

Kent