← Back to team overview

ffc team mailing list archive

Re: Forms

 



Martin Sandve Alnæs wrote:
By default, the variables "M", "L", and "a" are checked
and used if they are Form instances. Otherwise, any set
of forms (and elements) can be exported by defining lists
with the instances in:

  V = FiniteElement(...)
  ...
  a = ...*dx
  forms = (a, L, F, J)
  elements = (V, P, TH)

I'll update the manual.


OK.

Still, I would expect an error for the below form which is mixes linear and bilinear terms.

Garth


Martin



On Sun, Jun 7, 2009 at 9:59 AM, Garth N. Wells<gnw20@xxxxxxxxx> wrote:
I couldn't find in the UFL manual whether or not there is something
special about 'a' and 'L' denoting bilinear and linear forms,
respectively, or if this determined by the form compiler?

Either way, it seems that some checks are missing in either UFL or FFC
since

  element = FiniteElement("Lagrange", triangle, 1)
  v = TestFunction(element)
  u = TrialFunction(element)
  f = Function(element)

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

doesn't result in an error. I would expect this to give an error.

Garth
_______________________________________________
FFC-dev mailing list
FFC-dev@xxxxxxxxxx
http://www.fenics.org/mailman/listinfo/ffc-dev



References