ffc team mailing list archive
-
ffc team
-
Mailing list archive
-
Message #02808
Re: [HG FFC] Restrictions appears to be working for the Poisson demo.
Quoting kent-and@xxxxxxxxx:
> > Quoting kent-and@xxxxxxxxx:
> >
> >> > Quoting "Garth N. Wells" <gnw20@xxxxxxxxx>:
> >> >
> >> >>
> >> >>
> >> >> FFC wrote:
> >> >> > One or more new changesets pushed to the primary ffc repository.
> >> >> > A short summary of the last three changesets is included below.
> >> >> >
> >> >> > changeset: 1612:de211a9d88bce134c4abc2b128283ad35415a529
> >> >> > tag: tip
> >> >> > user: "Kristian Oelgaard <k.b.oelgaard@xxxxxxxxxx>"
> >> >> > date: Tue Aug 18 17:57:36 2009 +0100
> >> >> > files: ffc/fem/finiteelement.py
> >> >> > description:
> >> >> > Restrictions appears to be working for the Poisson demo.
> >> >> >
> >> >>
> >> >> Nice. What does the UFL syntax look like?
> >> >
> >> > scalar = FiniteElement("CG", triangle, 3)
> >> > vector = VectorElement("CG", triangle, 3)
> >> >
> >> > #element = scalar
> >> > #element = scalar["interval"]
> >> > #element = scalar[interval] + scalar
> >> > element = vector[interval]
> >>
> >> What does this mean ?
> >
> > It's related to this thread
> >
> > http://www.fenics.org/pipermail/ffc-dev/2009-July/002726.html
> >
> > the 'scalar' element will have 10 dofs, 3 at the vertices 6 on the edges
> > and 1
> > inside the cell.
> >
> > element = scalar[interval]
> >
> > will not have the dof which is inside the cell, it is restricted to the
> > facets
> > (which is interval in this case).
> >
> > Kristian
> >
>
> Nice!!!
>
> Will the 'interval' dofs include the vertex dofs ?
Yes, in the current implementation. We can always change this later once we
agree on the interpretation.
Kristian
> Kent
>
References