Merge branch 'Ferrite-FEM:master' into L2-pyr-wed

Author: AbdAlazezAhmed

.github/workflows/Check.yml

@@ -10,7 +10,7 @@ jobs:
   pre-commit:
     runs-on: ubuntu-latest
     steps:
-      - uses: actions/checkout@v4
+      - uses: actions/checkout@v6
       - uses: pre-commit/action@2c7b3805fd2a0fd8c1884dcaebf91fc102a13ecd # v3.0.1
         env:
           # Skip runic-pre-commit since we use runic-action below instead
@@ -20,7 +20,7 @@ jobs:
     name: "Runic"
     runs-on: ubuntu-latest
     steps:
-      - uses: actions/checkout@v4
+      - uses: actions/checkout@v6
       - uses: julia-actions/setup-julia@v2
         with:
           version: '1.11'
@@ -33,7 +33,7 @@ jobs:
     runs-on: ubuntu-latest
     name: "ExplicitImports.jl"
     steps:
-      - uses: actions/checkout@v4
+      - uses: actions/checkout@v6
       # - uses: julia-actions/setup-julia@v2
       #   with:
       #     version: '1'

.github/workflows/Documentation.yml

@@ -6,14 +6,14 @@ on:
     tags: ['*']
 jobs:
   docs:
-    name: Julia 1.11 - ubuntu-latest
+    name: Julia 1.12 - ubuntu-latest
     runs-on: ubuntu-latest
     timeout-minutes: 60
     steps:
-      - uses: actions/checkout@v4
+      - uses: actions/checkout@v6
       - uses: julia-actions/setup-julia@v2
         with:
-          version: '1.11'
+          version: '1.12'
       - uses: julia-actions/cache@v2
       - name: Install dependencies
         run: julia --project=docs -e 'using Pkg; Pkg.instantiate(); Pkg.precompile()'

.github/workflows/Downstream.yml

@@ -29,14 +29,14 @@ jobs:
           - {user: 'Ferrite-FEM', repo: 'FerriteInterfaceElements.jl'}
           # - {user: 'Ferrite-FEM', repo: 'FerriteViz.jl'} # Requires release
     steps:
-      - uses: actions/checkout@v4
+      - uses: actions/checkout@v6
       - uses: julia-actions/setup-julia@v2
         with:
           version: ${{ matrix.julia-version }}
       - uses: julia-actions/cache@v2
       # - uses: julia-actions/julia-buildpkg@v1
       - name: Clone Downstream
-        uses: actions/checkout@v4
+        uses: actions/checkout@v6
         with:
           repository: ${{ matrix.package.user }}/${{ matrix.package.repo }}
           path: downstream

.github/workflows/Test.yml

@@ -9,18 +9,18 @@ jobs:
     name: Julia ${{ matrix.julia-version }} - ${{ matrix.os }}
     runs-on: ${{ matrix.os }}
     timeout-minutes: 40
-    continue-on-error: ${{ matrix.julia-version == 'nightly' }}
+    continue-on-error: ${{ matrix.julia-version == 'pre' }}
     strategy:
       matrix:
-        julia-version: ['1.10', '1', 'nightly']
+        julia-version: ['lts', '1', 'pre']
         os: ['ubuntu-latest']
         include:
           - os: windows-latest
             julia-version: '1'
           - os: macOS-latest
             julia-version: '1'
     steps:
-      - uses: actions/checkout@v4
+      - uses: actions/checkout@v6
       - uses: julia-actions/setup-julia@v2
         with:
           version: ${{ matrix.julia-version }}
@@ -34,7 +34,7 @@ jobs:
           Pkg.add("JET")
           Pkg.develop(path = pwd())
           include(joinpath(pwd(), "test/jet.jl"))
-        if: ${{ matrix.os == 'ubuntu-latest' && matrix.julia-version != 'nightly' }}
+        if: ${{ matrix.os == 'ubuntu-latest' && matrix.julia-version != 'pre' }}
       - uses: julia-actions/julia-processcoverage@v1
         with:
           directories: 'src,ext'

CHANGELOG.md

@@ -7,11 +7,41 @@ and this project adheres to [Semantic Versioning](https://semver.org/spec/v2.0.0
 
 ## [Next] - xxxx-xx-xx
 
+## [v1.3.0] - 2026-01-09
+
 ### Added
-- Support for directly assembling to `SparseMatrixCSR` (from `SparseMatricesCSR.jl`). ([#864])
+ - New function `Ferrite.write_facetset` for exporting facet sets to VTK. ([#1132], [#1259])
+
+### Documentation
+ - Fix errors in incompressible elasticity tutorial. ([#1251])
+ - Add utility for downloading meshes (and other assets) for tutorials. ([#1037])
+
+## [v1.2.0] - 2025-11-14
+
+### Added
+ - Support for directly assembling to `SparseMatrixCSR` (from `SparseMatricesCSR.jl`). ([#864])
+ - Enhance `generate_grid` to support outputting line meshes embedded in two and three
+   spatial dimentions. ([#1122], [#1214])
+
+### Fixes
+ - Fix L2 projection of tensor fields on discontinuous interpolations. ([#1197], [#1198])
+ - Fix VTK output of tensor data in embedded domains. ([#1212], [#1213])
+ - Fix VTK nodal output for discontinuous output. ([#1194])
+ - Fix `facetskeleton` for mixed grids. ([#1223])
+ - Fix normal vector computation for boundary integration in embedded domains. ([#1226],
+   [#1235])
+ - Relax vector input type from `Vector` to `AbstractVector` in `evaluate_at_points!`.
+   ([#1183])
+- The result of grid coloring is now not dependent on iteration order of dictionaries and sets. This may change the coloring obtained from `create_coloring`
+  compared to previous Ferrite versions.
 
 ### Documentation updates
-- Extended assembly docs with information on how to support direct assembly into new matrix types. ([#864])
+ - Extended assembly docs with information on how to support direct assembly into new matrix
+   types. ([#864])
+ - Add a list of reserach papers where Ferrite was used for numbeir simulations. Please add
+   your paper! ([#1221])
+ - Add section on global and local DoF numbering. ([#1089])
+ - Fix some typos and grammar ([#1210], [#1224], [#1228])
 
 ## [v1.1.0] - 2025-05-01
 
@@ -873,6 +903,8 @@ poking into Ferrite internals:
 [v0.3.13]: https://github.com/Ferrite-FEM/Ferrite.jl/releases/tag/v0.3.13
 [v0.3.14]: https://github.com/Ferrite-FEM/Ferrite.jl/releases/tag/v0.3.14
 [v1.0.0]: https://github.com/Ferrite-FEM/Ferrite.jl/releases/tag/v1.0.0
+[v1.1.0]: https://github.com/Ferrite-FEM/Ferrite.jl/releases/tag/v1.1.0
+[v1.2.0]: https://github.com/Ferrite-FEM/Ferrite.jl/releases/tag/v1.2.0
 [#352]: https://github.com/Ferrite-FEM/Ferrite.jl/issues/352
 [#363]: https://github.com/Ferrite-FEM/Ferrite.jl/issues/363
 [#378]: https://github.com/Ferrite-FEM/Ferrite.jl/issues/378
@@ -1045,6 +1077,8 @@ poking into Ferrite internals:
 [#779]: https://github.com/Ferrite-FEM/Ferrite.jl/issues/779
 [#835]: https://github.com/Ferrite-FEM/Ferrite.jl/issues/835
 [#855]: https://github.com/Ferrite-FEM/Ferrite.jl/issues/855
+[#864]: https://github.com/Ferrite-FEM/Ferrite.jl/issues/864
+[#867]: https://github.com/Ferrite-FEM/Ferrite.jl/issues/867
 [#880]: https://github.com/Ferrite-FEM/Ferrite.jl/issues/880
 [#888]: https://github.com/Ferrite-FEM/Ferrite.jl/issues/888
 [#914]: https://github.com/Ferrite-FEM/Ferrite.jl/issues/914
@@ -1054,7 +1088,31 @@ poking into Ferrite internals:
 [#949]: https://github.com/Ferrite-FEM/Ferrite.jl/issues/949
 [#953]: https://github.com/Ferrite-FEM/Ferrite.jl/issues/953
 [#974]: https://github.com/Ferrite-FEM/Ferrite.jl/issues/974
+[#1044]: https://github.com/Ferrite-FEM/Ferrite.jl/issues/1044
+[#1045]: https://github.com/Ferrite-FEM/Ferrite.jl/issues/1045
 [#1058]: https://github.com/Ferrite-FEM/Ferrite.jl/issues/1058
 [#1059]: https://github.com/Ferrite-FEM/Ferrite.jl/issues/1059
+[#1063]: https://github.com/Ferrite-FEM/Ferrite.jl/issues/1063
 [#1083]: https://github.com/Ferrite-FEM/Ferrite.jl/issues/1083
+[#1089]: https://github.com/Ferrite-FEM/Ferrite.jl/issues/1089
 [#1096]: https://github.com/Ferrite-FEM/Ferrite.jl/issues/1096
+[#1122]: https://github.com/Ferrite-FEM/Ferrite.jl/issues/1122
+[#1146]: https://github.com/Ferrite-FEM/Ferrite.jl/issues/1146
+[#1151]: https://github.com/Ferrite-FEM/Ferrite.jl/issues/1151
+[#1162]: https://github.com/Ferrite-FEM/Ferrite.jl/issues/1162
+[#1176]: https://github.com/Ferrite-FEM/Ferrite.jl/issues/1176
+[#1178]: https://github.com/Ferrite-FEM/Ferrite.jl/issues/1178
+[#1183]: https://github.com/Ferrite-FEM/Ferrite.jl/issues/1183
+[#1194]: https://github.com/Ferrite-FEM/Ferrite.jl/issues/1194
+[#1197]: https://github.com/Ferrite-FEM/Ferrite.jl/issues/1197
+[#1198]: https://github.com/Ferrite-FEM/Ferrite.jl/issues/1198
+[#1210]: https://github.com/Ferrite-FEM/Ferrite.jl/issues/1210
+[#1212]: https://github.com/Ferrite-FEM/Ferrite.jl/issues/1212
+[#1213]: https://github.com/Ferrite-FEM/Ferrite.jl/issues/1213
+[#1214]: https://github.com/Ferrite-FEM/Ferrite.jl/issues/1214
+[#1221]: https://github.com/Ferrite-FEM/Ferrite.jl/issues/1221
+[#1223]: https://github.com/Ferrite-FEM/Ferrite.jl/issues/1223
+[#1224]: https://github.com/Ferrite-FEM/Ferrite.jl/issues/1224
+[#1226]: https://github.com/Ferrite-FEM/Ferrite.jl/issues/1226
+[#1228]: https://github.com/Ferrite-FEM/Ferrite.jl/issues/1228
+[#1235]: https://github.com/Ferrite-FEM/Ferrite.jl/issues/1235

Project.toml

@@ -1,6 +1,6 @@
 name = "Ferrite"
 uuid = "c061ca5d-56c9-439f-9c0e-210fe06d3992"
-version = "1.1.0"
+version = "1.3.0"
 
 [deps]
 EnumX = "4e289a0a-7415-4d19-859d-a7e5c4648b56"

benchmark/benchmark-cellvalues.jl

@@ -0,0 +1,67 @@
+using Ferrite, BenchmarkTools
+
+function element_routine(cv::CellValues)
+    s = 0.0
+    for q_point in 1:getnquadpoints(cv)
+        dΩ = getdetJdV(cv, q_point)
+        for i in 1:getnbasefunctions(cv)
+            for j in 1:getnbasefunctions(cv)
+                s += norm(shape_value(cv, q_point, i)) * norm(shape_gradient(cv, q_point, j)) * dΩ
+            end
+        end
+    end
+    return s
+end
+
+function element_routine(cv::MultiFieldCellValues)
+    s = 0.0
+    for q_point in 1:getnquadpoints(cv)
+        dΩ = getdetJdV(cv, q_point)
+        for i in 1:getnbasefunctions(cv.s1)
+            for j in 1:getnbasefunctions(cv.s1)
+                s += norm(shape_value(cv.s1, q_point, i)) * norm(shape_gradient(cv.s1, q_point, j)) * dΩ
+            end
+        end
+    end
+    return s
+end
+
+
+CT = QuadraticHexahedron
+ip_geo = geometric_interpolation(CT)
+x = Ferrite.reference_coordinates(ip_geo)
+dim = Ferrite.getrefdim(ip_geo)
+RefShape = getrefshape(ip_geo)
+qr = QuadratureRule{RefShape}(2)
+ip1 = Lagrange{RefShape, 1}()
+ip2 = Lagrange{RefShape, 2}()
+
+update_gradients = true
+update_hessians = false
+
+cv_s1 = CellValues(qr, ip1, ip_geo; update_gradients, update_hessians)
+cv_v1 = CellValues(qr, ip1^dim, ip_geo; update_gradients, update_hessians)
+cv_s2 = CellValues(qr, ip2, ip_geo; update_gradients, update_hessians)
+cv_v2 = CellValues(qr, ip2^dim, ip_geo; update_gradients, update_hessians)
+
+cmv_s1 = MultiFieldCellValues(qr, (s1 = ip1,), ip_geo; update_gradients, update_hessians)
+cmv_s1_square = MultiFieldCellValues(qr, (s1 = ip1, s2 = ip1), ip_geo; update_gradients, update_hessians)
+cmv_s1_v1 = MultiFieldCellValues(qr, (s1 = ip1, v1 = ip1^dim), ip_geo; update_gradients, update_hessians)
+cmv_s2_v2 = MultiFieldCellValues(qr, (s2 = ip2, v2 = ip2^dim), ip_geo; update_gradients, update_hessians)
+cmv_s1_v1_s2_v2 = MultiFieldCellValues(qr, (s1 = ip1, v1 = ip1^dim, s2 = ip2, v2 = ip2^dim), ip_geo; update_gradients, update_hessians)
+
+for key in (
+        :cv_s1, :cv_v1, :cv_s2, :cv_v2,
+        :cmv_s1, :cmv_s1_square, :cmv_s1_v1, :cmv_s2_v2, :cmv_s1_v1_s2_v2,
+    )
+    cv = getproperty(Main, key)
+    print(key, ": ")
+    @btime reinit!($cv, $x)
+end
+
+# Check that the indexing inside the element routine doesn't affect the performance
+for key in (:cv_s1, :cmv_s1_v1_s2_v2)
+    cv = getproperty(Main, key)
+    reinit!(cv, x)
+    @btime element_routine($cv)
+end