diff --git a/src/PatchBasedSmoothers/seq/PatchTriangulations.jl b/src/PatchBasedSmoothers/seq/PatchTriangulations.jl
index b16abfd3..fead9f41 100644
--- a/src/PatchBasedSmoothers/seq/PatchTriangulations.jl
+++ b/src/PatchBasedSmoothers/seq/PatchTriangulations.jl
@@ -44,12 +44,12 @@ function Gridap.Geometry.Triangulation(PD::PatchDecomposition)
   return PatchTriangulation(trian,PD,patch_cells,nothing,nothing)
 end
 
-function Gridap.Geometry.BoundaryTriangulation(PD::PatchDecomposition{Dc}) where Dc
-  Df       = Dc -1 
+function Gridap.Geometry.BoundaryTriangulation(PD::PatchDecomposition{Dc};tags="boundary") where Dc
+  Df       = Dc-1 
   model    = PD.model
   labeling = get_face_labeling(model)
 
-  is_boundary = get_face_mask(labeling,["boundary"],Df)
+  is_boundary = get_face_mask(labeling,tags,Df)
   patch_faces = get_patch_faces(PD,Df,is_boundary)
   pfaces_to_pcells = get_pfaces_to_pcells(PD,Df,patch_faces)
 
diff --git a/test/LinearSolvers/GMGTests.jl b/test/LinearSolvers/GMGTests.jl
index ce78dc42..bff9ab83 100644
--- a/test/LinearSolvers/GMGTests.jl
+++ b/test/LinearSolvers/GMGTests.jl
@@ -64,6 +64,7 @@ function gmg_driver(t,parts,mh,spaces,qdegree,smoothers,biform,liform,u)
                         mode=:preconditioner)
   
   solver = CGSolver(gmg;maxiter=20,atol=1e-14,rtol=1.e-6,verbose=i_am_main(parts))
+  #solver = GMRESSolver(5;Pr=gmg,maxiter=20,atol=1e-14,rtol=1.e-6,verbose=i_am_main(parts))
   ns = numerical_setup(symbolic_setup(solver,A),A)
   toc!(t,"GMG setup")