Skip to content

Commit

Permalink
Merge pull request #179 from firedrakeproject/fix-facet-avg-round-two
Browse files Browse the repository at this point in the history
* fix-facet-avg-round-two:
  Possible more right?
  • Loading branch information
wence- committed Dec 5, 2018
2 parents e4cadde + 8904c3e commit ed84cbc
Showing 1 changed file with 3 additions and 2 deletions.
5 changes: 3 additions & 2 deletions tsfc/fem.py
Original file line number Diff line number Diff line change
Expand Up @@ -285,11 +285,12 @@ def facet_avg(self, o):
integrand, = o.ufl_operands
domain = o.ufl_domain()
measure = ufl.Measure(self.context.integral_type, domain=domain)
integrand, degree, argument_multiindices = entity_avg(integrand / CellVolume(domain), measure, self.context.argument_multiindices)
integrand, degree, argument_multiindices = entity_avg(integrand / FacetArea(domain), measure, self.context.argument_multiindices)

config = {name: getattr(self.context, name)
for name in ["ufl_cell", "precision", "index_cache",
"integration_dim", "entity_ids"]}
"integration_dim", "entity_ids",
"integral_type"]}
config.update(quadrature_degree=degree, interface=self.context,
argument_multiindices=argument_multiindices)
expr, = compile_ufl(integrand, point_sum=True, **config)
Expand Down

0 comments on commit ed84cbc

Please sign in to comment.