6template <
typename Scalar>
17 typedef DenseIndex
Index;
24 JacobiRotation<Scalar> givens;
35 s.topLeftCorner(n,n).template triangularView<StrictlyLower>() = s.topLeftCorner(n,n).transpose();
38 for (j = 0; j < n; ++j) {
52 for (k = j; k < n; ++k) {
55 givens.makeGivens(-s(k,k), sdiag[k]);
59 s(k,k) = givens.c() * s(k,k) + givens.s() * sdiag[k];
60 temp = givens.c() * wa[k] + givens.s() * qtbpj;
61 qtbpj = -givens.s() * wa[k] + givens.c() * qtbpj;
65 for (i = k+1; i<n; ++i) {
66 temp = givens.c() * s(i,k) + givens.s() * sdiag[i];
67 sdiag[i] = -givens.s() * s(i,k) + givens.c() * sdiag[i];
76 for(nsing=0; nsing<n && sdiag[nsing]!=0; nsing++) {}
78 wa.tail(n-nsing).setZero();
79 s.topLeftCorner(nsing, nsing).transpose().template triangularView<Upper>().solveInPlace(wa.head(nsing));
86 for (j = 0; j < n; ++j) x[ipvt[j]] = wa[j];
Derived & setZero(Index size)
Namespace containing all symbols from the Eigen library.
EIGEN_DEFAULT_DENSE_INDEX_TYPE Index