Polytopal methods

.github/workflows/ci.yml

@@ -10,7 +10,7 @@ jobs:
       fail-fast: false
       matrix:
         version:
-          - '1.9'
+          - '1.10'
         os:
           - ubuntu-latest
         arch:
@@ -45,7 +45,7 @@ jobs:
       fail-fast: false
       matrix:
         version:
-          - '1.9'
+          - '1.10'
         os:
           - ubuntu-latest
         arch:
@@ -85,7 +85,7 @@ jobs:
       - uses: actions/checkout@v4
       - uses: julia-actions/setup-julia@v2
         with:
-          version: '1.8'
+          version: '1.10'
       - run: |
           julia --project=docs -e '
             using Pkg

NEWS.md

@@ -5,6 +5,13 @@ All notable changes to this project will be documented in this file.
 The format is based on [Keep a Changelog](https://keepachangelog.com/en/1.0.0/),
 and this project adheres to [Semantic Versioning](https://semver.org/spec/v2.0.0.html).
 
+## [0.4.8] 2025-06-11
+
+### Added
+
+- Added support for Gridap v0.19, with distributed counterparts for the new feaures introduced. This includes support for polytopal meshes, polytopal methods and patch assembly. Since PR[#175](https://github.com/gridap/GridapDistributed.jl/pull/175).
+- Added `MacroDiscreteModel`, which gives a global numbering and classification of the interfaces between processors. Since PR[#175](https://github.com/gridap/GridapDistributed.jl/pull/175).
+
 ## [0.4.7] 2025-03-04
 
 ### Added

Project.toml

@@ -1,7 +1,7 @@
 name = "GridapDistributed"
 uuid = "f9701e48-63b3-45aa-9a63-9bc6c271f355"
 authors = ["S. Badia <santiago.badia@monash.edu>", "A. F. Martin <alberto.f.martin@anu.edu.au>", "F. Verdugo <f.verdugo.rojano@vu.nl>"]
-version = "0.4.7"
+version = "0.4.8"
 
 [deps]
 BlockArrays = "8e7c35d0-a365-5155-bbbb-fb81a777f24e"
@@ -18,8 +18,8 @@ WriteVTK = "64499a7a-5c06-52f2-abe2-ccb03c286192"
 [compat]
 BlockArrays = "1"
 FillArrays = "1"
-ForwardDiff = "0.10"
-Gridap = "0.18.7"
+ForwardDiff = "0.10, 1"
+Gridap = "0.19.1"
 LinearAlgebra = "1"
 MPI = "0.16, 0.17, 0.18, 0.19, 0.20"
 PartitionedArrays = "0.3.3"

docs/src/Geometry.md

@@ -3,4 +3,13 @@
 ```@autodocs
 Modules = [GridapDistributed]
 Pages   = ["Geometry.jl"]
-```
+```
+
+## MacroDiscreteModels
+
+These are functionalities to select and globally number the interfaces between processors:
+
+```@autodocs
+Modules = [GridapDistributed]
+Pages   = ["MacroDiscreteModels.jl"]
+```

src/Adaptivity.jl

@@ -1155,7 +1155,7 @@ function refine_cell_gids(
   cmodel::DistributedDiscreteModel{Dc},
   fmodels::AbstractArray{<:DiscreteModel{Dc}}
 ) where Dc
-  cgids = get_cell_gids(cmodel)
+  cgids = partition(get_cell_gids(cmodel))
   f_own_to_local = map(cgids,fmodels) do cgids,fmodel
     glue = get_adaptivity_glue(fmodel)
     f2c_map = glue.n2o_faces_map[Dc+1]

src/Algebra.jl

@@ -117,6 +117,15 @@ function to_parray_of_arrays(a::AbstractArray{<:DebugArray})
   end
 end
 
+function to_parray_of_arrays(a::AbstractArray{<:Vector})
+  indices = linear_indices(first(a))
+  map(indices) do i
+    map(a) do aj
+      aj[i]
+    end
+  end
+end
+
 # This type is required because MPIArray from PArrays 
 # cannot be instantiated with a NULL communicator
 struct MPIVoidVector{T} <: AbstractVector{T}

src/Autodiff.jl

@@ -209,7 +209,7 @@ function distributed_autodiff_array_gradient(a, i_to_x, j_to_i::AbstractArray{<:
     # where each entry is a 2-block BlockVector with the first block being the
     # contribution of the plus side and the second, the one of the minus side
     is_single_field = eltype(eltype(j_to_result_plus)) <: Number
-    k = is_single_field ? BlockMap(2,[1,2]) : Fields.BlockBroadcasting(BlockMap(2,[1,2]))
+    k = is_single_field ? BlockMap(2,[1,2]) : Arrays.BlockBroadcasting(BlockMap(2,[1,2]))
     lazy_map(k,j_to_result_plus,j_to_result_minus)
   end
 
@@ -244,7 +244,7 @@ function distributed_autodiff_array_jacobian(a, i_to_x, j_to_i::AbstractArray{<:
       [(CartesianIndex(1,), CartesianIndex(1, 2)), (CartesianIndex(2,), CartesianIndex(2, 2))]  # Minus -> Second column
     ]
     is_single_field = eltype(eltype(j_to_result_plus)) <: AbstractArray
-    k = is_single_field ? Fields.MergeBlockMap((2,2),I) : Fields.BlockBroadcasting(Fields.MergeBlockMap((2,2),I))
+    k = is_single_field ? Arrays.MergeBlockMap((2,2),I) : Arrays.BlockBroadcasting(Arrays.MergeBlockMap((2,2),I))
     lazy_map(k,j_to_result_plus,j_to_result_minus)
   end
 

src/Geometry.jl

@@ -419,7 +419,7 @@ function Geometry.DiscreteModel(
   gids = PRange(partition)
 
   models = map(lcell_to_cell) do lcell_to_cell
-    DiscreteModelPortion(model,lcell_to_cell)
+    Geometry.restrict(model,lcell_to_cell)
   end
 
   GenericDistributedDiscreteModel(models,gids)

src/GridapDistributed.jl

@@ -65,4 +65,8 @@ include("Adaptivity.jl")
 
 include("Autodiff.jl")
 
+include("MacroDiscreteModels.jl")
+
+include("PatchAssemblers.jl")
+
 end # module