ufl team mailing list archive

[Kristian Oelgaard <k.b.oelgaard@xxxxxxxxxx>]

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)

Also, what is the strategy for handling the form

a = (inner(grad(v), grad(u)) + v*u)*dx ?

will the order of the form be determined by the highest order, or will it be
possible to split the integrand according to order? The latter could get
complicated for forms like:

a = f*(inner(grad(v), grad(u)) + v*u))*dx

I do notice the # TODO in formdata.py, I'm just curious.

Kristian