phynexis-deps/armadillo-code
1
0
Fork 0
Code Issues Pull Requests Packages Projects Releases Wiki Activity

use the same setting for singular value threshold as pinv()

Browse Source
This commit is contained in:
conrad
2022-03-01 14:19:42 +10:00
parent 934e24012c
commit ec1892e29d
1 changed files with 4 additions and 2 deletions
Download Patch File Download Diff File Expand all files Collapse all files
include/armadillo_bits/auxlib_meat.hpp
+4 -2
View File
@@ -5078,7 +5078,8 @@ auxlib::solve_approx_svd(Mat<typename T1::pod_type>& out, Mat<typename T1::pod_t
blas_int nrhs = blas_int(B.n_cols);
blas_int lda = blas_int(A.n_rows);
blas_int ldb = blas_int(tmp.n_rows);
eT rcond = eT(-1); // -1 means "use machine precision"
//eT rcond = eT(-1); // -1 means "use machine precision"
eT rcond = (std::max)(A.n_rows, A.n_cols) * std::numeric_limits<eT>::epsilon();
blas_int rank = blas_int(0);
blas_int info = blas_int(0);
@@ -5199,7 +5200,8 @@ auxlib::solve_approx_svd(Mat< std::complex<typename T1::pod_type> >& out, Mat< s
blas_int nrhs = blas_int(B.n_cols);
blas_int lda = blas_int(A.n_rows);
blas_int ldb = blas_int(tmp.n_rows);
T rcond = T(-1); // -1 means "use machine precision"
//T rcond = T(-1); // -1 means "use machine precision"
T rcond = (std::max)(A.n_rows, A.n_cols) * std::numeric_limits<T>::epsilon();
blas_int rank = blas_int(0);
blas_int info = blas_int(0);