tdecomp function

Root of singularity-corrected eigenvalue decomposition