On Mon, Jun 09, 2014 at 11:30:09AM +0200, Jan Blechta wrote:
> On Mon, 9 Jun 2014 11:10:12 +0200
> Anders Logg <[email protected]> wrote:
>
> > For vector elements, the jump() operator in UFL is defined as follows:
> >
> >   dot(v('+'), n('+')) + dot(v('-'), n('-'))
> >
> > I'd like to argue that it should instead be implemented like so:
> >
> >   outer(v('+'), n('+')) + outer(v('-'), n('-'))
>
> This inconsistency has been already encountered by users
> http://fenicsproject.org/qa/359/discontinuous-galerkin-jump-operators

Interesting! I hadn't noticed.

Are there any objections to changing this definition in UFL?

--
Anders
_______________________________________________
fenics mailing list
[email protected]
http://fenicsproject.org/mailman/listinfo/fenics

Reply via email to