ffc team mailing list archive

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

Quoting Marie Rognes <meg@xxxxxxxxxxx>:

> Kristian Oelgaard wrote:
> > Quoting Fredrik Bengzon <fredrik.bengzon@xxxxxxxxxxx>:
> >
> >   
> >> Hello
> >> I'm wondering if there is some simple way (i.e., other than component 
> >> wise) of writing the term
> >>
> >> integral ( [ outer(n, v) ] : < grad u > ) dS
> >>     
> >
> > If I do the following:
> >
> > element = FiniteElement("BDM", "triangle", 1)
> >
> > v = TrialFunction(element)
> > u = TestFunction(element)
> > n = FacetNormal("triangle")
> >
> > print "jump", jump(outer(n,v))
> > print "avg", avg(grad(u))
> >
> > a = dot(jump(outer(n,v)), avg(grad(u)))*dS
> >
> > print "a", a
> >
> > at least the output (the print statements) looks OK as far as I can tell,
> but
> > have a look. The above form will however produce an error in simplify.py.
> Maybe
> > you can have a look and see what is going wrong?
> >   
> 
> A proper bug in simplify. Should be fixed now.
> 

Exellent

Kristian
 
> 
> -- 
> Marie E. Rognes
> Ph.D Fellow, 
> Centre of Mathematics for Applications, 
> 
> University of Oslo
> 
> http://folk.uio.no/meg
>