Merge pull request #110 from gridap/polytopal

Gridap 0.19

Author: JordiManyer

NEWS.md

@@ -4,6 +4,12 @@ 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.9.7] - 2025-6-11
+
+### Added
+
+- Added support for Gridap v0.19. Since PR[#110](https://github.com/gridap/GridapEmbedded.jl/pull/110).
+
 ## [0.9.6] - 2025-04-19
 
 ### Added

Project.toml

@@ -1,23 +1,7 @@
-authors = ["Francesc Verdugo <[email protected]>", "Eric Neiva <[email protected]>", "Pere Antoni Martorell <[email protected]>", "Santiago Badia <[email protected]>"]
 name = "GridapEmbedded"
 uuid = "8838a6a3-0006-4405-b874-385995508d5d"
-version = "0.9.6"
-
-[compat]
-AbstractTrees = "0.3.3, 0.4"
-Algoim = "0.2.2"
-Combinatorics = "1"
-CxxWrap = "0.16"
-FillArrays = "0.10, 0.11, 0.12, 0.13, 1"
-FiniteDiff = "2.27.0"
-ForwardDiff = "0.10.38, 1"
-Graphs = "1.12.0"
-Gridap = "0.18.12"
-GridapDistributed = "0.3, 0.4"
-MPI = "0.20"
-MiniQhull = "0.1.0, 0.2, 0.3, 0.4"
-PartitionedArrays = "0.3.4"
-julia = "1.3"
+authors = ["Francesc Verdugo <[email protected]>", "Eric Neiva <[email protected]>", "Pere Antoni Martorell <[email protected]>", "Santiago Badia <[email protected]>"]
+version = "0.9.7"
 
 [deps]
 AbstractTrees = "1520ce14-60c1-5f80-bbc7-55ef81b5835c"
@@ -37,6 +21,22 @@ Random = "9a3f8284-a2c9-5f02-9a11-845980a1fd5c"
 Test = "8dfed614-e22c-5e08-85e1-65c5234f0b40"
 algoimWrapper_jll = "3c43aa7b-5398-51f3-8d75-8f051e6faa4d"
 
+[compat]
+AbstractTrees = "0.3.3, 0.4"
+Algoim = "0.2.2"
+Combinatorics = "1"
+CxxWrap = "0.16"
+FillArrays = "0.10, 0.11, 0.12, 0.13, 1"
+FiniteDiff = "2.27.0"
+ForwardDiff = "0.10.38, 1"
+Graphs = "1.12.0"
+Gridap = "0.18.12, 0.19"
+GridapDistributed = "0.3, 0.4"
+MPI = "0.20"
+MiniQhull = "0.1.0, 0.2, 0.3, 0.4"
+PartitionedArrays = "0.3.4"
+julia = "1.3"
+
 [extras]
 FiniteDiff = "6a86dc24-6348-571c-b903-95158fe2bd41"
 MPIPreferences = "3da0fdf6-3ccc-4f1b-acd9-58baa6c99267"

src/Interfaces/EmbeddedFacetDiscretizations.jl

@@ -215,12 +215,15 @@ end
 
 function _restrict_boundary_triangulation(model,facets,bgfacet_to_mask)
   facet_to_bgfacet = facets.glue.face_to_bgface
-  facet_to_mask = lazy_map(Reindex(bgfacet_to_mask),facet_to_bgfacet)
+  
   n_bgfacets = length(bgfacet_to_mask)
-  bgfacet_to_mask2 = fill(false,n_bgfacets)
-  bgfacet_to_mask2[facet_to_bgfacet] .= facet_to_mask
+  new_bgfacet_to_mask = fill(false,n_bgfacets)
+  new_bgfacet_to_mask[facet_to_bgfacet] .= view(bgfacet_to_mask,facet_to_bgfacet)
 
-  BoundaryTriangulation(model,bgfacet_to_mask2,facets.glue.bgface_to_lcell)
+  new_bgfacet_to_lcell = fill(Int8(1),n_bgfacets)
+  new_bgfacet_to_lcell[facet_to_bgfacet] .= facets.glue.face_to_lcell
+
+  BoundaryTriangulation(model,new_bgfacet_to_mask,new_bgfacet_to_lcell)
 end
 
 function compute_bgfacet_to_inoutcut(cut::EmbeddedFacetDiscretization,geo::CSG.Geometry)