← Back to team overview

ufl team mailing list archive

Re: Dirichlet boundary conditions

 

On Mon, Jun 23, 2008 at 01:22:13PM +0200, kent-and@xxxxxxxxx wrote:

> > Maybe we can do:
> >
> > V = FiniteElement("Lagrange", "triangle", 1)
> > Vg = Restricted(V, g, i) # subspace of V with u=g on dS_i
> > u = TrialFunction(Vg)
> > f = Function(V)
> >
> >
> 
> I like the abstraction/name trace better than restriction, but no hard
> feelings.
> It could also be 'Constraints'.
> 
> Kent

How about this:

A FunctionSpace is created from a FiniteElement and a set of
Constraints, so we may do things like

  element = FiniteElement("Lagrange", "triangle", 1)
  constraints = (Trace(u, s0) == g, Trace(u, s1) == h)

  V = FunctionSpace(element, constraints)

We can let __eq__ in Trace return a Constraint object.

See here also:

  http://www.geuz.org/getdp/doc/slides/getdp-10.html

-- 
Anders

Attachment: signature.asc
Description: Digital signature


Follow ups

References