ffc team mailing list archive

[Kristian Oelgaard <k.b.oelgaard@xxxxxxxxxx>]

Quoting Anders Logg <logg@xxxxxxxxx>:

> I've made a first attempt to remove the old stuff. It mostly works but
> there were some problems I couldn't figure out for the quadrature
> representation.
>
> The MixedPoisson demo fails for quadrature with the following error:
>
>   File
>
>
"/usr/lib/python2.5/site-packages/ffc/compiler/codegeneration/quadrature/uflquadraturegenerator_utils.py",
> line 842, in __create_mapping_basis
>     name, non_zeros, zeros, ones = self.name_map[name]
> KeyError: 'FE0_C3'
>
> Any clues what goes wrong?

Yes, after some digging I discovered that it's because someone has changed

if ffc_element.value_mapping(component) == Mapping.AFFINE:
to
if ffc_element.space_mapping(component) == Mapping.AFFINE:

in many places.

However, value_mapping produces the following error now:

  File "/home/oelgaard/fenics/clone/ffc/ffc/fem/mixedelement.py", line 81, in
value_mapping
    i = Index(component)
NameError: global name 'Index' is not defined

This is probably because Index has been removed from FFC? So we should use
space_mapping?

The problem is that for:

BDM = FiniteElement("Brezzi-Douglas-Marini", triangle, q)
DG  = FiniteElement("Discontinuous Lagrange", triangle, q - 1)

mixed_element = BDM + DG

(tau, w) = TestFunctions(mixed_element)
(sigma, u) = TrialFunctions(mixed_element)

a = w*u*dx

The FFC element for the mixed basis function is:
Mixed finite element: [FiniteElement('Brezzi-Douglas-Marini', 'triangle', 1),
FiniteElement('Discontinuous Lagrange', 'triangle', 0)]

But when trying to get the mapping for component 2 (the DG element)

ffc_element.space_mapping(2) =  1

which is not equal to Mapping.AFFINE = 0

Is there a simple way of fixing this using space_mapping, or should we fix the
value_mapping?

Kristian

> --
> Anders
>
>
>
>
> On Mon, Apr 20, 2009 at 10:44:40PM +0200, 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:   1506:fbeeaf084b14549672b1f71a129e958c047a8979
> > tag:         tip
> > user:        Anders Logg <logg@xxxxxxxxx>
> > date:        Mon Apr 20 22:44:36 2009 +0200
> > files:       ffc/__init__.py ffc/compiler/analysis/elementdata.py
> ffc/compiler/analysis/formdata.py
> ffc/compiler/codegeneration/common/finiteelement.py
> ffc/compiler/codegeneration/quadrature/uflquadraturegenerator_utils.py
> ffc/compiler/codegeneration/tensor/tensorgenerator.py
> ffc/compiler/compiler.py ffc/compiler/language/README
> ffc/compiler/language/__init__.py ffc/compiler/language/algebra.py
> ffc/compiler/language/builtins.py ffc/compiler/language/index.py
> ffc/compiler/language/indexcall.py ffc/compiler/language/indexreordering.py
> ffc/compiler/language/integral.py ffc/compiler/language/operators.py
> ffc/compiler/language/reassignment.py ffc/compiler/language/restriction.py
> ffc/compiler/language/signature.py ffc/compiler/language/tokens.py
> ffc/compiler/representation/tensor/facetmap.py
> ffc/compiler/representation/tensor/factorization.py
> ffc/compiler/representation/tensor/geometrytensor.py
> ffc/compiler/representation/tensor/monomialextraction.py
> ffc/compiler/representation/tensor/monomialintegration.py
> ffc/compiler/representation/tensor/referencetensor.py
> ffc/compiler/representation/tensor/tensorreordering.py
> ffc/compiler/representation/tensor/tensorrepresentation.py
> ffc/compiler/representation/tensor/ufltensorrepresentation.py
> ffc/compiler/uflcompiler.py ffc/fem/finiteelement.py ffc/fem/mixedelement.py
> ffc/jit/jit.py ffc/jit/jitobject.py scripts/ffc test/unit/test.py
> > description:
> > Delete a lot of old stuff.
> >
> >
> > changeset:   1505:b6be99e11da21ba89d11651447f28815e591b115
> > user:        Anders Logg <logg@xxxxxxxxx>
> > date:        Mon Apr 20 21:29:22 2009 +0200
> > files:       ChangeLog demo/Constant.form demo/Elasticity.form
> demo/EnergyNorm.form demo/Equation.form demo/FacetIntegrals.form
> demo/FunctionOperators.form demo/Heat.form demo/Mass.form
> demo/MixedMixedElement.form demo/MixedPoisson.form demo/NavierStokes.form
> demo/NeumannProblem.form demo/Optimization.form demo/P5tet.form
> demo/P5tri.form demo/Poisson.form demo/PoissonDG.form demo/PoissonSystem.form
> demo/QuadratureElement.form demo/README demo/Stokes.form demo/SubDomain.form
> demo/SubDomains.form demo/TensorWeightedPoisson.form
> demo/VectorLaplaceGradCurl.form test/regression/reference/Constant.h
> test/regression/reference/Elasticity.h test/regression/reference/EnergyNorm.h
> test/regression/reference/Equation.h
> test/regression/reference/FacetIntegrals.h
> test/regression/reference/FunctionOperators.h
> test/regression/reference/Heat.h test/regression/reference/Mass.h
> test/regression/reference/MixedMixedElement.h
> test/regression/reference/MixedPoisson.h
> test/regression/reference/NavierStokes.h
> test/regression/reference/NeumannProblem.h
> test/regression/reference/Optimization.h test/regression/reference/P5tet.h
> test/regression/reference/P5tri.h test/regression/reference/Poisson.h
> test/regression/reference/PoissonDG.h
> test/regression/reference/PoissonSystem.h
> test/regression/reference/QuadratureElement.h
> test/regression/reference/Stokes.h test/regression/reference/SubDomain.h
> test/regression/reference/SubDomains.h
> test/regression/reference/TensorWeightedPoisson.h
> test/regression/reference/VectorLaplaceGradCurl.h test/regression/test.py
> test/simple_verify_tensors/reference.pickle
> test/simple_verify_tensors/test.py
> > description:
> > Remove .form files and update references
> >
> >
> > changeset:   1504:cc48e5d8e03a5d3c15c20085ae682379081e2740
> > user:        Anders Logg <logg@xxxxxxxxx>
> > date:        Mon Apr 20 19:55:45 2009 +0200
> > files:       scripts/ffc
> > description:
> > Remove support for old style .form files.
> >
> > ----------------------------------------------------------------------
> > 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
>