 *          END DO * *        Apply PERMUTATIONS P and P**T: *        P * inv(U**T) * inv(D) * inv(U) * P**T. *        Interchange rows and columns I and IPIV(I) in reverse order *        from the formation order of IPIV vector for 