Merge pull request #918 from gridap/patch-in-modal-C0-bases-for-1.9

Update ModalC0Bases.jl for Julia 1.9

Author: amartinhuertas

NEWS.md

@@ -24,6 +24,7 @@ and this project adheres to [Semantic Versioning](https://semver.org/spec/v2.0.0
   Since PR [#884](https://github.com/gridap/Gridap.jl/pull/884).
 
 - Fixed `cell_dof_ids` for the case of vectorial `ConstantFESpace`. Since PR [#888](https://github.com/gridap/Gridap.jl/pull/888)
+- Fixed generation of Modal C0 bases for Julia 1.9. Since PR [#918](https://github.com/gridap/Gridap.jl/pull/918).
 
 ## [0.17.17] - 2023-02-24
 

src/Polynomials/ModalC0Bases.jl

@@ -266,8 +266,8 @@ function _sort_by_nfaces!(terms::Vector{CartesianIndex{D}},orders) where D
   # Generate indices of n-faces and order s.t.
   # (1) dimension-increasing (2) lexicographic
   bin_rang_nfaces = tfill(0:1,Val{D}())
-  bin_ids_nfaces = collect(Iterators.product(bin_rang_nfaces...))
-  sum_bin_ids_nfaces = [sum(bin_ids_nfaces[i]) for i in eachindex(bin_ids_nfaces)]
+  bin_ids_nfaces = vec(collect(Iterators.product(bin_rang_nfaces...)))
+  sum_bin_ids_nfaces = sum.(bin_ids_nfaces)
   bin_ids_nfaces = permute!(bin_ids_nfaces,sortperm(sum_bin_ids_nfaces))
 
   # Generate LIs of basis funs s.t. order by n-faces
@@ -529,4 +529,4 @@ function _hessian_nd_mc0!(
 
   end
 
-end
+end