Skip to content

Mimetic spectral elements #5

New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

Sign up for GitHub

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Merged
chrisrichardson merged 3 commits into from
Sep 10, 2019
Merged

Mimetic spectral elements #5

Show file tree
Hide file tree
Changes from all commits
Commits
Show all changes
3 commits
Select commit Hold shift + click to select a range
7f8ab0a
Added mse primal and dual elements
celdred Jun 27, 2019
b7208e8
Merge branch 'master' into mse
chrisrichardson Sep 10, 2019
2d5b853
Merge branch 'master' into mse
chrisrichardson Sep 10, 2019
File filter

Filter by extension

Filter by extension
Conversations
Failed to load comments.
Loading
Jump to
Jump to file
Failed to load files.
Loading
Diff view
Diff view
22 changes: 22 additions & 0 deletions test/test_elements.py
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -199,3 +199,25 @@ def test_radau():

element = FiniteElement("Radau", cell, degree)
assert element == eval(repr(element))


def test_mse():
for degree in (2, 3, 4, 5):
element = FiniteElement('EGL', interval, degree)
assert element == eval(repr(element))

element = FiniteElement('EGL-Edge', interval, degree - 1)
assert element == eval(repr(element))

element = FiniteElement('EGL-Edge L2', interval, degree - 1)
assert element == eval(repr(element))

for degree in (1, 2, 3, 4, 5):
element = FiniteElement('GLL', interval, degree)
assert element == eval(repr(element))

element = FiniteElement('GLL-Edge', interval, degree - 1)
assert element == eval(repr(element))

element = FiniteElement('GLL-Edge L2', interval, degree - 1)
assert element == eval(repr(element))
14 changes: 14 additions & 0 deletions ufl/finiteelement/elementlist.py
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -258,6 +258,20 @@ def show_elements():
register_alias("Q- L2", lambda family, dim, order,
degree: feec_element_l2(family, dim, order, degree))

# mimetic spectral elements - primal and dual complexs
register_element("Extended-Gauss-Legendre", "EGL", 0, H1, "identity", (2, None),
("interval",))
register_element("Extended-Gauss-Legendre Edge", "EGL-Edge", 0, L2, "identity", (1, None),
("interval",))
register_element("Extended-Gauss-Legendre Edge L2", "EGL-Edge L2", 0, L2, "L2 Piola", (1, None),
("interval",))
register_element("Gauss-Lobatto-Legendre Edge", "GLL-Edge", 0, L2, "identity", (0, None),
("interval",))
register_element("Gauss-Lobatto-Legendre Edge L2", "GLL-Edge L2", 0, L2, "L2 Piola", (0, None),
("interval",))


# NOTE- the edge elements for primal mimetic spectral elements are accessed by using variant='mse' in the appropriate places

def feec_element(family, n, r, k):
"""Finite element exterior calculus notation
Expand Down