Allow discontinuous elements (#293)

* add discontinuous elements

* ++version

* mypy

Author: mscroggs

.github/workflows/test-packages.yml

@@ -18,7 +18,7 @@ jobs:
           python-version: ${{ matrix.python-version }}
       - uses: actions/checkout@v4
         with:
-          ref: v2024.10.0
+          ref: v2025.3.0
         name: Get latest Symfem version
       - run: rm -rf symfem VERSION
         name: Remove downloaded symfem
@@ -51,7 +51,7 @@ jobs:
           activate-environment: symfem
       - uses: actions/checkout@v4
         with:
-          ref: v2024.10.0
+          ref: v2025.3.0
         name: Get latest Symfem version
       - run: rm -rf symfem VERSION
         name: Remove downloaded symfem
@@ -83,15 +83,15 @@ jobs:
           python-version: ${{ matrix.python-version }}
       - uses: actions/checkout@v4
         with:
-          ref: v2024.10.0
+          ref: v2025.3.0
           path: symfem-src
       - name: Move tests and delete Symfem source
         run: |
           mv symfem-src/test .
           mv symfem-src/README.md .
           rm -r symfem-src
       - run: |
-          wget -O symfem.tar.gz https://pypi.io/packages/source/s/symfem/symfem-2024.10.0.tar.gz
+          wget -O symfem.tar.gz https://pypi.io/packages/source/s/symfem/symfem-2025.3.0.tar.gz
           tar -xvzf symfem.tar.gz
         name: Download and unpack latest version of Symfem
       - name: Install LaTeΧ
@@ -101,7 +101,7 @@ jobs:
       - run: pip install pytest
         name: Install pytest
       - run: |
-          cd symfem-2024.10.0
+          cd symfem-2025.3.0
           pip install .[optional]
         name: Install requirements
       - run: |

CHANGELOG_SINCE_LAST_VERSION.md

@@ -3,3 +3,4 @@
 - Correct Guzman-Neilan second kind element
 - Add Guzman-Neilan first kind element
 - Correct Guzman-Neilan bubbles
+- Add discontinuous elements

CITATION.cff

@@ -5,8 +5,8 @@ authors:
     given-names: Matthew W.
     orcid: 0000-0002-4658-2443
 title: Symfem
-version: 2024.10.0
-date-released: 2024-10-05
+version: 2025.3.0
+date-released: 2025-03-03
 license: MIT
 url: https://github.com/mscroggs/symfem
 preferred-citation:

LICENSE

@@ -1,6 +1,6 @@
 MIT License
 
-Copyright (c) 2020-2024 Matthew Scroggs
+Copyright (c) 2020-2025 Matthew Scroggs
 
 Permission is hereby granted, free of charge, to any person obtaining a copy
 of this software and associated documentation files (the "Software"), to deal

VERSION

@@ -1 +1 @@
-2024.10.0
+2025.3.0

codemeta.json

@@ -1 +1 @@
-{"@context": "https://raw.githubusercontent.com/codemeta/codemeta/master/codemeta.jsonld", "@type": "Code", "author": [{"@id": "0000-0002-4658-2443", "@type": "Person", "email": "[email protected]", "name": "Matthew Scroggs", "affiliation": "Department of Engineering, University of Cambridge"}], "identifier": "", "codeRepository": "https://github.com/mscroggs/symfem", "datePublished": "2021-01-23", "dateModified": "2024-10-05", "dateCreated": "2021-01-23", "description": "A symbolic finite element definition library", "keywords": "Python, finite element method, numerical analysis", "license": "MIT", "title": "Symfem", "version": "2024.10.0"}
+{"@context": "https://raw.githubusercontent.com/codemeta/codemeta/master/codemeta.jsonld", "@type": "Code", "author": [{"@id": "0000-0002-4658-2443", "@type": "Person", "email": "[email protected]", "name": "Matthew Scroggs", "affiliation": "Department of Engineering, University of Cambridge"}], "identifier": "", "codeRepository": "https://github.com/mscroggs/symfem", "datePublished": "2021-01-23", "dateModified": "2025-03-03", "dateCreated": "2021-01-23", "description": "A symbolic finite element definition library", "keywords": "Python, finite element method, numerical analysis", "license": "MIT", "title": "Symfem", "version": "2025.3.0"}

pyproject.toml

@@ -4,7 +4,7 @@ requires = ["flit_core >=3.8.0,<4"]
 
 [project]
 name = "symfem"
-version = "2024.10.0"
+version = "2025.3.0"
 description = "a symbolic finite element definition library"
 readme = "README.md"
 requires-python = ">=3.8.0"

symfem/create.py

@@ -6,6 +6,7 @@
 
 import symfem.references as _references
 from symfem.finite_element import FiniteElement as _FiniteElement
+from symfem.finite_element import CiarletElement as _CiarletElement
 from symfem.geometry import SetOfPointsInput as _SetOfPointsInput
 from symfem.geometry import parse_set_of_points_input as _parse_set_of_points_input
 
@@ -196,6 +197,16 @@ def create_element(
             raise ValueError(f"Order {order} {element_type} element cannot be created.")
         return element_class(reference, order, **kwargs)
 
+    if element_type.startswith("discontinuous "):
+        e = create_element(cell_type, element_type[14:], order, vertices, **kwargs)
+        if isinstance(e, _CiarletElement):
+            for d in e.dofs:
+                d.entity = (reference.tdim, 0)
+            e.continuity = "L2"
+            e.names = ["discontinuous " + i for i in e.names]
+
+            return e
+
     raise ValueError(f"Unsupported element type: {element_type}")
 
 

symfem/finite_element.py

@@ -585,6 +585,7 @@ def name(self) -> str:
     last_updated = version
     cache = True
     value_type = "unknown"
+    continuity: typing.Optional[str] = None
     _max_continuity_test_order = 4
 
 

symfem/version.py

@@ -1,3 +1,3 @@
 """Version number."""
 
-version = "2024.10.0"
+version = "2025.3.0"