The documentation of boomeramg mention that it's possible to solve also matrices arising from the discretization of system of PDEs. But there is no more information on it. What should I do to make use of it in PETSc? Thomas