5 template <
typename Scalar>
7 Matrix< Scalar, Dynamic, Dynamic > &r,
9 Scalar tol = std::sqrt(NumTraits<Scalar>::epsilon()) )
12 typedef DenseIndex Index;
15 Index i, j, k, l, ii, jj;
20 const Index n = r.cols();
21 const Scalar tolr = tol * abs(r(0,0));
22 Matrix< Scalar, Dynamic, 1 > wa(n);
23 eigen_assert(ipvt.size()==n);
27 for (k = 0; k < n; ++k)
28 if (abs(r(k,k)) > tolr) {
30 for (j = 0; j <= k-1; ++j) {
31 temp = r(k,k) * r(j,k);
33 r.col(k).head(j+1) -= r.col(j).head(j+1) * temp;
40 for (k = 0; k <= l; ++k) {
41 for (j = 0; j <= k-1; ++j)
42 r.col(j).head(j+1) += r.col(k).head(j+1) * r(j,k);
43 r.col(k).head(k+1) *= r(k,k);
48 for (j = 0; j < n; ++j) {
51 for (i = 0; i <= j; ++i) {
64 r.topLeftCorner(n,n).template triangularView<StrictlyUpper>() = r.topLeftCorner(n,n).transpose();
Namespace containing all symbols from the Eigen library.
Definition: AdolcForward:45