Difference between revisions of "Remncomm.m"

From Spinach Documentation Wiki
Jump to: navigation, search
(Normalise sequential blank lines)
m (Rename Arguments section heading to Parameters)
Line 6: Line 6:
 
     C=remncomm(A,EvB)
 
     C=remncomm(A,EvB)
  
==Arguments==
+
==Parameters==
  
 
     A    -  a square matrix
 
     A    -  a square matrix

Revision as of 18:51, 5 June 2026

Removes from the Hermitian operator A the part that does not commute with the Hermitian operator B.

Syntax

    C=remncomm(A,EvB)

Parameters

   A     -  a square matrix

   EvB   -  a square matrix containing eigenvectors 
            of B in columns

Outputs

   C     -  a square matrix

Notes

When the matrix B is diagonal, the part of A that commutes with it is the diagonal part, so just use diag(diag(A)).

See also

Numerical infrastructure

Version 2.5, authors: Ilya Kuprov