Thanks to the basis-function as in Zaglmayr’s thesis, we can use standard smoothers (point Gauss-Seidel, block Gauss-Seidel, IC,…) for the high-order system (plus something good for the lowest-order sub-system).
See also this talk, around page 18.
Joachim