diff options
Diffstat (limited to 'unsupported/Eigen/src/NonLinearOptimization/covar.h')
| -rw-r--r-- | unsupported/Eigen/src/NonLinearOptimization/covar.h | 10 |
1 files changed, 6 insertions, 4 deletions
diff --git a/unsupported/Eigen/src/NonLinearOptimization/covar.h b/unsupported/Eigen/src/NonLinearOptimization/covar.h index 104898a30..6c77916f5 100644 --- a/unsupported/Eigen/src/NonLinearOptimization/covar.h +++ b/unsupported/Eigen/src/NonLinearOptimization/covar.h @@ -1,9 +1,10 @@ +namespace internal { template <typename Scalar> -void ei_covar( +void covar( Matrix< Scalar, Dynamic, Dynamic > &r, const VectorXi &ipvt, - Scalar tol = ei_sqrt(NumTraits<Scalar>::epsilon()) ) + Scalar tol = sqrt(NumTraits<Scalar>::epsilon()) ) { typedef DenseIndex Index; @@ -14,14 +15,14 @@ void ei_covar( /* Function Body */ const Index n = r.cols(); - const Scalar tolr = tol * ei_abs(r(0,0)); + const Scalar tolr = tol * abs(r(0,0)); Matrix< Scalar, Dynamic, 1 > wa(n); assert(ipvt.size()==n); /* form the inverse of r in the full upper triangle of r. */ l = -1; for (k = 0; k < n; ++k) - if (ei_abs(r(k,k)) > tolr) { + if (abs(r(k,k)) > tolr) { r(k,k) = 1. / r(k,k); for (j = 0; j <= k-1; ++j) { temp = r(k,k) * r(j,k); @@ -61,3 +62,4 @@ void ei_covar( r.diagonal() = wa; } +} // end namespace internal |