Skip to content

Commit

Permalink
quadrature degree fixed to p+1
Browse files Browse the repository at this point in the history
  • Loading branch information
@tlroy
tlroy committed Feb 8, 2024
1 parent 0ee8c76 commit 87374f1
Showing 1 changed file with 12 additions and 6 deletions.
18 changes: 12 additions & 6 deletions echemfem/solver.py
View file
Original file line number Diff line number Diff line change
Expand Up @@ -211,34 +211,40 @@ def _internal_dx(subdomain_id=None, domain=None):
return _dx(
subdomain_id=subdomain_id,
domain=domain,
scheme=quadrature_rule)
scheme=quadrature_rule,
degree=p+1)

if mesh.layers:
# This doesn't include top and bottom surfaces
def _internal_ds(subdomain_id=None, domain=None):
return ds_v(
subdomain_id=subdomain_id,
domain=domain,
scheme=quadrature_rule_face)
scheme=quadrature_rule_face,
degree=p+1)

def _internal_dS(subdomain_id=None, domain=None):
return dS_v(subdomain_id=subdomain_id,
domain=domain,
scheme=quadrature_rule_face) + dS_h(subdomain_id=subdomain_id,
scheme=quadrature_rule_face,
degree=p+1) + dS_h(subdomain_id=subdomain_id,
domain=domain,
scheme=quadrature_rule_face)
scheme=quadrature_rule_face,
degree=p+1)
else:
def _internal_ds(subdomain_id=None, domain=None):
return _ds(
subdomain_id=subdomain_id,
domain=domain,
scheme=quadrature_rule_face)
scheme=quadrature_rule_face,
degree=p+1)

def _internal_dS(subdomain_id=None, domain=None):
return _dS(
subdomain_id=subdomain_id,
domain=domain,
scheme=quadrature_rule_face)
scheme=quadrature_rule_face,
degree=p+1)

self.ds = _internal_ds
self.dx = _internal_dx
Expand Down

0 comments on commit 87374f1

Please sign in to comment.