ufl team mailing list archive

[Anders Logg <logg@xxxxxxxxx>]

Yes, most of all. :-)

-- 
Anders


On Mon, Jun 08, 2009 at 10:25:30AM +0200, Martin Sandve Alnæs wrote:
> Facet normals should always be restricted as well, right?
> 
> Martin
> 
> 
> 
> On Fri, Jun 5, 2009 at 8:54 PM, Garth N. Wells<gnw20@xxxxxxxxx> wrote:
> >
> >
> > Martin Sandve Alnæs wrote:
> >>
> >> I have zero knowledge of discontinuous Galerkin methods,
> >> so there are no checks at all. What should we check for?
> >>
> >> Do you mean (implementation must be different):
> >>
> >> if integral is on interior facets:
> >>    for f in basis function and functions in integrand:
> >>        assert that f is restricted
> >> ?
> >>
> >
> > I would be nice to have the above. It would have saved me a few hours this
> > afternoon . . .
> >
> > Garth
> >
> >
> >> We should probably propagate ('+') and ('-') up
> >> the tree to reach the leaves (terminals), right?
> >>
> >> Martin
> >>
> >>
> >>
> >> On Thu, May 14, 2009 at 7:52 AM, Garth N. Wells <gnw20@xxxxxxxxx> wrote:
> >>>
> >>> FFC used to throw an error if functions were integrated on interior
> >>> facets but not restricted to the '+' or '-' side, but UFL doesn't. If I
> >>> do
> >>>
> >>>    element = FiniteElement("Discontinuous Lagrange", triangle, 1)
> >>>    v = TestFunction(element)
> >>>    u = TrialFunction(element)
> >>>
> >>>    # Integrate on interior facet
> >>>    a = v*u*dS
> >>>
> >>> what result should I expect? I had anticipated an error.
> >>>
> >>> Garth
> >>> _______________________________________________
> >>> UFL-dev mailing list
> >>> UFL-dev@xxxxxxxxxx
> >>> http://fenics.org/mailman/listinfo/ufl-dev
> >>>
> >
> >
> >
> _______________________________________________
> UFL-dev mailing list
> UFL-dev@xxxxxxxxxx
> http://fenics.org/mailman/listinfo/ufl-dev

Attachment: signature.asc
Description: Digital signature