ufl team mailing list archive
-
ufl team
-
Mailing list archive
-
Message #00198
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
-
Re: Dirichlet boundary conditions
From: Martin Sandve Alnæs, 2008-06-19
-
Re: Dirichlet boundary conditions
From: Kent-Andre Mardal, 2008-06-19
-
Re: Dirichlet boundary conditions
From: Martin Sandve Alnæs, 2008-06-19
-
Re: Dirichlet boundary conditions
From: Kent-Andre Mardal, 2008-06-19
-
Re: Dirichlet boundary conditions
From: Martin Sandve Alnæs, 2008-06-19
-
Re: Dirichlet boundary conditions
From: kent-and, 2008-06-20
-
Re: Dirichlet boundary conditions
From: Anders Logg, 2008-06-20
-
Re: Dirichlet boundary conditions
From: kent-and, 2008-06-23
-
Re: Dirichlet boundary conditions
From: Martin Sandve Alnæs, 2008-06-23
-
Re: Dirichlet boundary conditions
From: kent-and, 2008-06-23