← Back to team overview

ffc team mailing list archive

Re: [HG FFC] Restrictions appears to be working for the Poisson demo.

 



Kristian Oelgaard wrote:
Quoting "Garth N. Wells" <gnw20@xxxxxxxxx>:


Kristian Oelgaard wrote:
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]

Do we have point, edge, face and facet as well?

No, only interval, triangle, tetrahedron, quadrilateral and hexahedron

is predefined in UFL. I assume simplices and restrictions currently only
supports a restriction to one topological dimension lower i.e.,
tetrahedron -> triangle
triangle -> interval


We need a 'facet', which is a triangle in 3D and interval in 2D, so we write dimension-independent forms (at least for 2D and 3D).

Garth

Kristian


Garth

will not have the dof which is inside the cell, it is restricted to the
facets
(which is interval in this case).

Kristian

Kent

v = TestFunction(element)
u = TrialFunction(element)

a = inner(grad(v), grad(u))*dx

They should all work.

evaluate_basis and evaluate_basis_derivatives is not implemented for
restricted
elements. The reason is that I'm not sure about what should happen when I
throw
away basisfunctions.

Kristian

Garth

changeset:   1611:283c1799fe643b297668dfa5347a7de50fee8555
user:        "Kristian Oelgaard <k.b.oelgaard@xxxxxxxxxx>"
date:        Tue Aug 18 17:40:05 2009 +0100
files:       ffc/fem/finiteelement.py ffc/fem/quadratureelement.py
description:
Bug fix to make buildbot happy.


changeset:   1610:2d95d63e8281e4f06cb9a24f1895b6e70563d1c5
user:        "Kristian Oelgaard <k.b.oelgaard@xxxxxxxxxx>"
date:        Tue Aug 18 17:09:17 2009 +0100
files:       ffc/compiler/evaluatebasis.py
ffc/compiler/finiteelement.py
ffc/fem/createelement.py ffc/fem/finiteelement.py
description:
A first implementation of restrictions work for FiniteElements.

----------------------------------------------------------------------
For more details, visit http://www.fenics.org/hg/ffc
_______________________________________________
FFC-dev mailing list
FFC-dev@xxxxxxxxxx
http://www.fenics.org/mailman/listinfo/ffc-dev
_______________________________________________
FFC-dev mailing list
FFC-dev@xxxxxxxxxx
http://www.fenics.org/mailman/listinfo/ffc-dev

_______________________________________________
FFC-dev mailing list
FFC-dev@xxxxxxxxxx
http://www.fenics.org/mailman/listinfo/ffc-dev








Follow ups

References