--- rpl/lapack/lapack/dlaed1.f 2010/12/21 13:53:28 1.7 +++ rpl/lapack/lapack/dlaed1.f 2011/07/22 07:38:06 1.8 @@ -26,9 +26,9 @@ * of a full symmetric matrix (which was reduced to tridiagonal form) * are desired. * -* T = Q(in) ( D(in) + RHO * Z*Z' ) Q'(in) = Q(out) * D(out) * Q'(out) +* T = Q(in) ( D(in) + RHO * Z*Z**T ) Q**T(in) = Q(out) * D(out) * Q**T(out) * -* where Z = Q'u, u is a vector of length N with ones in the +* where Z = Q**T*u, u is a vector of length N with ones in the * CUTPNT and CUTPNT + 1 th elements and zeros elsewhere. * * The eigenvectors of the original matrix are stored in Q, and the