Final revision to polytopal drivers

Author: JordiManyer

.gitignore

@@ -10,3 +10,4 @@ tmp/
 .vscode/
 *.code-workspace
 *.DS_Store
+CLAUDE.md

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.19.0"
+version = "0.19.9"
 
 [deps]
 BlockArrays = "8e7c35d0-a365-5155-bbbb-fb81a777f24e"

src/poisson_hdg.jl

@@ -1,3 +1,4 @@
+#
 # In this tutorial, we will implement a Hybridizable Discontinuous Galerkin (HDG) method
 # for solving the Poisson equation. The HDG method is an efficient variant of DG methods
 # that introduces an auxiliary variable m on mesh interfaces to reduce the global system size.
@@ -163,3 +164,15 @@ writevtk(Ω,"output_path/results",cellfields=["uh"=>uh,"qh"=>qh,"eh"=>eh])
 # solving the same problem on a polytopal mesh is available [in the Gridap repository](https://github.com/gridap/Gridap.jl/blob/75efc9a7a7e286c27e7ca3ddef5468e591845484/test/GridapTests/HDGPolytopalTests.jl). 
 # A tutorial for HHO on polytopal meshes is also available. 
 #
+# You can also check out the following reference:
+#
+# @misc{GridapHybrid,
+#   title={A natural language framework for non-conforming hybrid polytopal methods in Gridap.jl}, 
+#   author={Jordi Manyer and Jai Tushar and Santiago Badia},
+#   year={2026},
+#   eprint={2603.00880},
+#   archivePrefix={arXiv},
+#   primaryClass={cs.MS},
+#   url={https://arxiv.org/abs/2603.00880}, 
+# }
+#

src/poisson_hho.jl

@@ -264,11 +264,13 @@ function patch_weakform()
 end
 
 op = MultiField.StaticCondensationOperator(X,patch_assem,patch_weakform())
-ui, ub = solve(op) 
+xh = solve(op) 
 
-eu  = ui - u
-l2u = sqrt(sum( ∫(eu * eu)dΩp))
-h1u = l2u + sqrt(sum( ∫(∇(eu) ⋅ ∇(eu))dΩp))
+uT, uF = xh
+
+eh = R(xh) - u
+l2err = sqrt(sum( ∫(eh * eh)dΩp))
+h1err = l2err + sqrt(sum( ∫(∇(eh) ⋅ ∇(eh))dΩp))
 
 # ## Going further
 #
@@ -280,3 +282,15 @@ h1u = l2u + sqrt(sum( ∫(∇(eu) ⋅ ∇(eu))dΩp))
 #   - [Incompressible Stokes](https://github.com/gridap/Gridap.jl/blob/75efc9a7a7e286c27e7ca3ddef5468e591845484/test/GridapTests/HHOMixedStokesPolytopal.jl) 
 #   - [Linear Elasticity](https://github.com/gridap/Gridap.jl/blob/75efc9a7a7e286c27e7ca3ddef5468e591845484/test/GridapTests/HHOMixedElasticity.jl)
 #
+# You can also check out the following reference:
+#
+# @misc{GridapHybrid,
+#   title={A natural language framework for non-conforming hybrid polytopal methods in Gridap.jl}, 
+#   author={Jordi Manyer and Jai Tushar and Santiago Badia},
+#   year={2026},
+#   eprint={2603.00880},
+#   archivePrefix={arXiv},
+#   primaryClass={cs.MS},
+#   url={https://arxiv.org/abs/2603.00880}, 
+# }
+#