Fix issue #203

Author: JordiManyer

Project.toml

@@ -1,7 +1,7 @@
 name = "Tutorials"
 uuid = "ab87341f-2653-5384-8dc6-691ea82e91b3"
-authors = ["Santiago Badia <santiago.badia@monash.edu>", "Francesc Verdugo <fverdugo@cimne.upc.edu>"]
 version = "0.20.0"
+authors = ["Santiago Badia <santiago.badia@monash.edu>", "Francesc Verdugo <fverdugo@cimne.upc.edu>"]
 
 [deps]
 BlockArrays = "8e7c35d0-a365-5155-bbbb-fb81a777f24e"
@@ -40,12 +40,13 @@ BlockArrays = "1"
 DataStructures = "0.19"
 Gridap = "0.20"
 GridapDistributed = "0.4"
-GridapEmbedded = "0.9.9"
+GridapEmbedded = "0.9"
 GridapGmsh = "0.7"
 GridapP4est = "0.3"
 GridapPETSc = "0.5"
 GridapSolvers = "0.7"
 MPI = "0.20"
+PartitionedArrays = "0.3"
 SpecialFunctions = "1"
 julia = "1.10"
 

src/poisson_hho.jl

@@ -191,7 +191,7 @@ R  = reconstruction_operator(ptopo,order,Y,Ωp,Γp,dΩp,dΓp)
 #   - The stabilization term `s`
 #   - The rhs term `l`
 
-hTinv = CellField(1 ./ collect(get_array(∫(1)dΩp)), Ωp)
+hTinv = CellField(1 ./ collect(get_array(∫(1)dΩp) .^ (1/2)), Ωp)
 
 function a(u,v)
   Ru_Ω, Ru_Γ = R(u)