← Back to team overview

ffc team mailing list archive

Re: [DOLFIN-dev] [HG FFC] Fixed computation of direction for facet normals

 



Peter Brune wrote:
I think I see what's going on from looking at the generated UFC code from ffc/demo/PoissonDG.ufl. It's only generating and using the normal for one side of the interior jump terms. This is going to cause the "first" cell here to have the proper normal, and the other to have a flipped normal. I can look at what happens in the interior facet integral code and have it generate and use the proper normals for both cells.


Could you register on Launchpad so the bug can be assigned to you.

Garth

- Peter

On Mon, Sep 14, 2009 at 12:28 PM, Garth N. Wells <gnw20@xxxxxxxxx <mailto:gnw20@xxxxxxxxx>> wrote:

    The new facet normal computation seems to give some seriously wrong
    results for the Poisson DG demo. Take a looks at the DOLFIN demo

        demo/pde/dg/poisson/python

    Peter: Could you please take a look?

    Garth

    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:   1651:09868489fc408b1b0024bbd43e901112076d1e46
     > tag:         tip
     > user:        Peter Brune <brune@xxxxxxxxxxxxxxx
    <mailto:brune@xxxxxxxxxxxxxxx>>
     > date:        Mon Sep 14 11:08:25 2009 -0500
     > files:       ffc/compiler/codesnippets.py
     > description:
     > Fixed computation of direction for facet normals
     >
     >
     > changeset:   1650:aa7511f9cd9c4784418e76ad7b027b8bc1d30b39
     > user:        Anders Logg <logg@xxxxxxxxx <mailto:logg@xxxxxxxxx>>
     > date:        Mon Sep 14 17:20:00 2009 +0200
     > files:       demo/Biharmonic.ufl ffc/compiler/codesnippets.py
    ffc/compiler/removeunused.py
    test/regression/reference/quadrature/Biharmonic.h
    test/regression/reference/quadrature/Constant.h
    test/regression/reference/quadrature/Elasticity.h
    test/regression/reference/quadrature/ElementRestriction.h
    test/regression/reference/quadrature/EnergyNorm.h
    test/regression/reference/quadrature/Equation.h
    test/regression/reference/quadrature/FacetIntegrals.h
    test/regression/reference/quadrature/FunctionOperators.h
    test/regression/reference/quadrature/Heat.h
    test/regression/reference/quadrature/Mass.h
    test/regression/reference/quadrature/MetaData.h
    test/regression/reference/quadrature/MixedMixedElement.h
    test/regression/reference/quadrature/MixedPoisson.h
    test/regression/reference/quadrature/NavierStokes.h
    test/regression/reference/quadrature/NeumannProblem.h
    test/regression/reference/quadrature/Normals.h
    test/regression/reference/quadrature/Optimization.h
    test/regression/reference/quadra
    ture/P5tet.h test/regression/reference/quadrature/P5tri.h
    test/regression/reference/quadrature/Poisson.h
    test/regression/reference/quadrature/PoissonDG.h
    test/regression/reference/quadrature/PoissonSystem.h
    test/regression/reference/quadrature/QuadratureElement.h
    test/regression/reference/quadrature/ReactionDiffusion.h
    test/regression/reference/quadrature/StabilisedStokes.h
    test/regression/reference/quadrature/Stokes.h
    test/regression/reference/quadrature/SubDomain.h
    test/regression/reference/quadrature/SubDomains.h
    test/regression/reference/quadrature/TensorWeightedPoisson.h
    test/regression/reference/quadrature/VectorLaplaceGradCurl.h
    test/regression/reference/tensor/Constant.h
    test/regression/reference/tensor/Elasticity.h
    test/regression/reference/tensor/ElementRestriction.h
    test/regression/reference/tensor/EnergyNorm.h
    test/regression/reference/tensor/Equation.h
    test/regression/reference/tensor/FacetIntegrals.h
    test/regression/reference/tensor/Heat.h test/regression/refere
    nce/tensor/Mass.h test/regression/reference/tensor/MetaData.h
    test/regression/reference/tensor/MixedMixedElement.h
    test/regression/reference/tensor/MixedPoisson.h
    test/regression/reference/tensor/NavierStokes.h
    test/regression/reference/tensor/NeumannProblem.h
    test/regression/reference/tensor/Normals.h
    test/regression/reference/tensor/Optimization.h
    test/regression/reference/tensor/P5tet.h
    test/regression/reference/tensor/P5tri.h
    test/regression/reference/tensor/Poisson.h
    test/regression/reference/tensor/PoissonSystem.h
    test/regression/reference/tensor/ReactionDiffusion.h
    test/regression/reference/tensor/StabilisedStokes.h
    test/regression/reference/tensor/Stokes.h
    test/regression/reference/tensor/SubDomain.h
    test/regression/reference/tensor/SubDomains.h
    test/regression/reference/tensor/VectorLaplaceGradCurl.h
    test/regression/test.py test/regression/update-references.sh
     > description:
     > Merge Peter's patch, add suppot for bool in remove_unused and
    update references.
     >
     >
     > changeset:   1649:2ef63eb5e7e7d0d24e81dee86eecb983ff044389
     > user:        Peter Brune <brune@xxxxxxxxxxxxxxx
    <mailto:brune@xxxxxxxxxxxxxxx>>
     > date:        Fri Sep 11 12:21:19 2009 -0500
     > files:       demo/Normals.ufl
     > description:
     > Added a demo for the normals
     >
     >
    ----------------------------------------------------------------------
     > For more details, visit http://www.fenics.org/hg/ffc
     > _______________________________________________
     > FFC-dev mailing list
     > FFC-dev@xxxxxxxxxx <mailto:FFC-dev@xxxxxxxxxx>
     > http://www.fenics.org/mailman/listinfo/ffc-dev


    _______________________________________________
    DOLFIN-dev mailing list
    DOLFIN-dev@xxxxxxxxxx <mailto:DOLFIN-dev@xxxxxxxxxx>
    http://www.fenics.org/mailman/listinfo/dolfin-dev





References