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.
- Peter On Mon, Sep 14, 2009 at 12:28 PM, Garth N. Wells <[email protected]> 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 <[email protected]> > > 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 <[email protected]> > > 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 <[email protected]> > > 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 > > [email protected] > > http://www.fenics.org/mailman/listinfo/ffc-dev > > > _______________________________________________ > DOLFIN-dev mailing list > [email protected] > http://www.fenics.org/mailman/listinfo/dolfin-dev >
_______________________________________________ DOLFIN-dev mailing list [email protected] http://www.fenics.org/mailman/listinfo/dolfin-dev
