def _compute_proj(A):
 (U, _, _) = linalg.svd(A, full_matrices=False)
 return (np.identity(A.shape[0]) - np.dot(U, U.T.conjugate()))

def _compute_proj ( A ) : ( U , _ , _ ) = linalg . svd ( A , full_matrices = False ) return ( np . identity ( A . shape [ 0 ] ) - np . dot ( U , U . T . conjugate ( ) ) )
compute the orthogonal projection operation for a manifold vector a .