10
$\begingroup$

In my calculations I need to use something which is "between" a matrix and its inverse. That is, I invert only some dimensions. I am interested if it has an established name.

That is, a matrix (here 2x2 real, but it is more general) $$ \begin{bmatrix} u' \\ v' \end{bmatrix} = M \begin{bmatrix} u \\ v \end{bmatrix} $$ defines a hyperplane in coordinates $(u,v,u',v')$. Its inverse (if exists) can be defined as a linear operator such that $$ \begin{bmatrix} u \\ v \end{bmatrix} = M^{-1} \begin{bmatrix} u' \\ v' \end{bmatrix}. $$

I am interested in inverting only some coordinates, e.g. $$ \begin{bmatrix} u \\ v' \end{bmatrix} = M^{(-1,1)} \begin{bmatrix} u' \\ v \end{bmatrix}. $$

I know it is a relatively simple thing related to the implicit function theorem, with simple formulas. Yet, I use it a lot and I need to call it somehow. So:

  • does it have its own name?
  • if not, is "partial inverse" fine? (not colliding with other names, not (too) confusing, etc)

If you are curious, I use it in physics (optics) to relate a scattering matrix (relating input to output) to a transfer matrix (relating left/right of an interface).

$\endgroup$
1

3 Answers 3

12
$\begingroup$

It is a principal pivot transform, also known as sweep operator or gyration. You can check the linked review paper.

$\endgroup$
6
  • $\begingroup$ So, as I see, it has many names. Thanks for this review! $\endgroup$ Commented Nov 3, 2014 at 9:06
  • $\begingroup$ @PiotrMigdal Yes, it's one of those operations that pop up unexpectedly in lots of different fields. It is interesting to find out that it has an application in optics, too! $\endgroup$ Commented Nov 3, 2014 at 13:20
  • 1
    $\begingroup$ It is used in optics, just with no name (as often it is used only once). My paper bases on jumping back and forth (e.g. unitary matrix -> partial inverse -> diagonalization -> partial inverse of $C$, $D$ and $C^{-1}$). In any case, all other names than partial inverse seem to me a bit contrived (i.e. there is no way I would have guessed what do they mean). $\endgroup$ Commented Nov 3, 2014 at 13:31
  • 1
    $\begingroup$ @PiotrMigdal Another one of its many names is exchange operator. Maybe you find that more enlightening? $\endgroup$ Commented Nov 3, 2014 at 14:02
  • 3
    $\begingroup$ beware of the other en.wikipedia.org/wiki/Exchange_operator $\endgroup$ Commented Nov 3, 2014 at 18:21
8
$\begingroup$

Partial Inversion For Linear Systems And Partial Closure Of Independence Graphs

$\endgroup$
2
  • $\begingroup$ Just curious, why inversion not inverse? $\endgroup$ Commented Nov 3, 2014 at 8:59
  • 3
    $\begingroup$ the operation is called a partial inversion, the outcome a partial inverse $\endgroup$ Commented Nov 3, 2014 at 9:00
2
$\begingroup$

Partial inversion is a valid name, and yeah, it is still a hot topic in S-matrix theory and Pseudo-Unitary Quantum Mechanics.

The partial inversion transformation $\hat{Y}_{ik}$ over a general rectangular matrix $Z$, with entries $Z_{\rho\sigma}$, is basically the application of the entry-wise transformations below.

\begin{align} Z_{ik} &\to Z_{ik}^{-1} ; \\ Z_{i\sigma} &\to -Z_{ik}^{-1} Z_{i\sigma}, \quad \text{for all} \quad \sigma \neq k; \\ Z_{\rho k} &\to Z_{\rho k}Z_{ik}^{-1}, \quad \text{for all} \quad \rho \neq i; \\ Z_{\rho \sigma} &\to Z_{\rho \sigma} - Z_{\rho k} Z_{ik}^{-1} Z_{i\sigma}, \quad \text{for all} \quad \sigma \neq k, \rho \neq i. \end{align}

It has a diagrammatic representation in this paper similar to permutation paths.

If you use it as a matrix transformation that is simultaneous to the permutation of vector entries between the input and output of a $\mathbf{a}=M\mathbf{b}$ relation, then the equation system defined by this relation is conserved.

According to this paper, it can be used to:

  1. reshape a linear equation system;
  2. build a matrix representation for a renormalized diagram of a scattering process;
  3. change pseudo-unitarity of a matrix, grow a quantum circuit, then come back to the original unitary group;
  4. reversibly "invert" singular matrices (you can apply the regular matrix inversion via partial inversion for singular matrices if you consider special conventions for zero-handling, but it may lack basic properties of inverses);
  5. define exact symbolic inverses for arbitrary dimensions of square matrices;
  6. define a flexibilization of matrix inverse for arbitrary rectangular matrices;
  7. define matrix inversion as a non-linear matrix operator that uses both (i) a commutator operator and (ii) a special commutator that resembles a Yang-Baxter equation;
  8. define generalized inverses and their relations to generalized transpositions along (or around) any diagonal other than the main diagonal.

Also, according to the graphs in soaklot's answer to this other question, it can:

  1. invert fully symbolic matrices on sympy (python package) for dimensions greater than (4,4) faster than sympy's default inversion algorithm.
$\endgroup$
1
  • 1
    $\begingroup$ Thanks for your answer! It is good to see even more applications. A remark: if I am not mistaken, inversion via PPTs is equivalent to Gauss--Jordan elimination; it is fine to use this algorithm with symbolic matrices, but it has worse backward stability properties if implemented in floating-point arithmetic; readers beware. $\endgroup$ Commented Aug 27 at 8:18

You must log in to answer this question.

Not the answer you're looking for? Browse other questions tagged .