The function factorizes the symmetric (or Hermitian)
matrix as:
is the outer factor
is the middle term
is the set of eigenvectors of P
is the vector
, where
,
, and
are the number of positive, zero, and negative eigenvalues respectively.