$V_i \subset \mathbb{R}^n$ and $V_0 \supset V_1 \supset ... \supset V_i \supset ...$,

$A_i, B_i: V_i \rightarrow V_i$ be square non-symmetric positive definite matrices,

$Q_i:V_{i-1}\rightarrow V_i$ (respectively $Q_i^ T:V_i\rightarrow V_{i-1}$) be rectangular matrices representing $L^2$-projections, $I_i$ be the identity in $V_i$ and

$A_i=Q_i A_{i-1} Q_i^T$.

Is there a way to represent the following matrix series as a vector-of-matrices times bloc-matrix times vector-of-matrices multiplication of the form $\overline{v}^ T\overline{M} \overline{v}$? Does this matrix series sound familiar to anyone?

\begin{align} (I_0-B_0 A_0)&(I_0-Q_1^T A_1^{-1} Q_1 A_0)(I_0-B_0 A_0) + \\ +(I_0-B_0 A_0)Q_1^T(I_1-B_1 A_1)&(I_1-Q_2^T A_2^{-1} Q_2 A_1)(I_1-B_1 A_1)A_1^{-1}Q_1 A_0(I_0-B_0 A_0)+ \\ +(I_0-B_0 A_0)Q_1^T(I_1-B_1 A_1)Q_2^T(I_2-B_2 A_2)&(I_2-Q_3^T A_3^{-1} Q_3 A_2)(I_2-B_2 A_2)A_2^{-1}Q_2 A_1(I_1-B_1 A_1) A_1^{-1} Q_1 A_0(I_0-B_0 A_0) \\ &... \end{align}

Should there be several representations, I am particularly interested in $\overline{M}$ to be as full as possible.


  • $\begingroup$ Well for sure there are trivial representations where the "vector of matrices" is $\begin{bmatrix}I\\I\\\vdots\\I\end{bmatrix}$ and $M$ is diagonal... $\endgroup$ Nov 16, 2016 at 19:12
  • $\begingroup$ Yes, of course, unfortunately that leaves the complexity too concentrated in the diagonal to use theory of infinite matrices of operators... $\endgroup$
    – Astor
    Nov 16, 2016 at 22:06

1 Answer 1


(cannot comment)

It resembles something which can arise in multigrid methods. If we assume that you are solving $Ax = b$,$V_i$ correspond to grid functions for different meshes, and $B_i$ are smoothers. But usually it is written as a recursive procedure and one stops at some $i=I$. You can try to look for something useful in classical multigrid books I guess.

  • $\begingroup$ Thanks a lot!! You actually hit the nail on the head, this expression comes from a multigrid context, unfortunately because $A_i$ and $B_i$ are non-symmetric, I am outside "classical" multigrid methods though :-( I was thinking on how to take this up to matrices of matrices, such that I can apply generalized matrices theory. $\endgroup$
    – Astor
    Nov 16, 2016 at 18:26

Your Answer

By clicking “Post Your Answer”, you agree to our terms of service, privacy policy and cookie policy

Not the answer you're looking for? Browse other questions tagged or ask your own question.