Vectorial constant fe space

NEWS.md

@@ -15,6 +15,8 @@ and this project adheres to [Semantic Versioning](https://semver.org/spec/v2.0.0
   was using default, it's now using the spetialized implementation for the underlying triangulation type.
   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)
+
 ## [0.17.17] - 2023-02-24
 
 ### Added

src/FESpaces/ConstantFESpaces.jl

@@ -30,7 +30,7 @@ cell_basis=SingleFieldFEBasis(
 cell_dof_basis_array=lazy_map(get_dof_basis,cell_reffe)
 cell_dof_basis=CellDof(cell_dof_basis_array,Triangulation(model),ReferenceDomain())
 
-cell_dof_ids=Fill(Int32[1],num_cells(model))
+cell_dof_ids=Fill(Int32(1):Int32(num_components(field_type)),num_cells(model))
 A=typeof(cell_basis)
 B=typeof(cell_dof_basis)
 C=typeof(cell_dof_ids)

test/FESpacesTests/ConstantFESpaceTests.jl

@@ -1,4 +1,3 @@
-
 module ConstantFESpacesTests
 
 using Gridap
@@ -30,4 +29,13 @@ uh = solve(op)
 @assert sum(∫((uh[1]-u)*(uh[1]-u))dΩ) < 1.0e-14
 abs(sum(∫(uh[2])dΩ)) < 1.0e-12
 
+Λ2=ConstantFESpace(model,field_type=VectorValue{2,Float64})
+Gridap.FESpaces.test_fe_space(Λ2)
+M2=TrialFESpace(Λ2)
+a2(μ,λ) = ∫(λ⋅μ)dΩ
+l2(λ) = ∫(VectorValue(0.0,0.0)⋅λ)dΩ
+op2 = AffineFEOperator(a2,l2,M2,Λ2)
+μ2h = solve(op2)
+@assert sum(∫(μ2h⋅μ2h)dΩ) < 1.0e-12
+
 end # module