109   MatrixType matA(matrix);  in compute()
124       m_sigma[k] = matA.col(k).end(m-k).norm();  in compute()
127         if (matA(k,k) < 0.0)  in compute()
129         matA.col(k).end(m-k) /= m_sigma[k];  in compute()
130         matA(k,k) += 1.0;  in compute()
140 …Scalar t = matA.col(k).end(m-k).eigen2_dot(matA.col(j).end(m-k)); // FIXME dot product or cwise pr…  in compute()
141         t = -t/matA(k,k);  in compute()
142         matA.col(j).end(m-k) += t * matA.col(k).end(m-k);  in compute()
147       e[j] = matA(k,j);  in compute()
152       m_matU.col(k).end(m-k) = matA.col(k).end(m-k);  in compute()
170         work.end(m-k-1) = matA.corner(BottomRight,m-k-1,n-k-1) * e.end(n-k-1);  in compute()
172           matA.col(j).end(m-k-1) += (-e[j]/e[k+1]) * work.end(m-k-1);  in compute()
185     m_sigma[nct] = matA(nct,nct);  in compute()
189     e[nrt] = matA(nrt,p-1);  in compute()