ufl team mailing list archive
-
ufl team
-
Mailing list archive
-
Message #01320
Re: Restrictions at facets
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
>>>
>
>
>
Follow ups
References