diff --git a/hog/forms.py b/hog/forms.py
index 064e0a91847d455accf95a8cf711a73440406149..8ae046393a7a8ad2b0330b12c960bb6eccfa7c1b 100644
--- a/hog/forms.py
+++ b/hog/forms.py
@@ -770,7 +770,10 @@ The resulting matrix must be multiplied with a vector of ones to be used as the
 def divdiv(
     trial: TrialSpace,
     test: TestSpace,
+    component_trial: int,
+    component_test: int,
     geometry: ElementGeometry,
+    quad: Quadrature,
     symbolizer: Symbolizer,
     blending: GeometryMap = IdentityMap(),
 ) -> Form:
@@ -783,6 +786,8 @@ divdiv operator which is used as a stabilization term that is taken from
 
 for the P1-P0 stabilized Stokes discretization.
 
+Component trial: {component_trial}
+Component test:  {component_test}
 Geometry map:    {blending}
 
 Weak formulation
@@ -806,6 +811,7 @@ Weak formulation
         docstring=docstring,
     )
 
+
 def k_divdiv(
     trial: TrialSpace,
     test: TestSpace,
@@ -815,13 +821,7 @@ def k_divdiv(
     blending: GeometryMap = IdentityMap(),
 ) -> Form:
     docstring = f"""
-divdiv operator which is used as a stabilization term that is taken from
-
-    Blank, L. (2014).
-    On Divergence-Free Finite Element Methods for the Stokes Equations (Freie Universität Berlin).
-    p. 84, eq. (6.2)
-
-for the P1-P0 stabilized Stokes discretization.
+divdiv operator which is the compressible part of full Stokes operator
 
 Geometry map:    {blending}