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

rename for clarification

Browse Source
This commit is contained in:
conrad
2022-02-24 15:40:43 +10:00
parent c0f8c37805
commit 5cd8718b85
1 changed files with 6 additions and 6 deletions
Download Patch File Download Diff File Expand all files Collapse all files
include/armadillo_bits/op_inv_rcond_meat.hpp
+6 -6
View File
@@ -102,29 +102,29 @@ op_inv_rcond::apply_direct_spd(Mat<typename T1::elem_type>& out, typename T1::po
eT* colmem = out.memptr();
T max_abs_src_val = T(0);
T max_abs_out_val = T(0);
T max_abs_inv_val = T(0);
for(uword i=0; i<N; ++i)
{
eT& out_ii = colmem[i];
const eT src_val = out_ii;
const eT out_val = eT(1) / src_val;
const eT inv_val = eT(1) / src_val;
if( (src_val == eT(0)) || (access::tmp_real(src_val) <= T(0)) ) { return false; }
out_ii = out_val;
out_ii = inv_val;
const T abs_src_val = std::abs(src_val);
const T abs_out_val = std::abs(out_val);
const T abs_inv_val = std::abs(inv_val);
max_abs_src_val = (abs_src_val > max_abs_src_val) ? abs_src_val : max_abs_src_val;
max_abs_out_val = (abs_out_val > max_abs_out_val) ? abs_out_val : max_abs_out_val;
max_abs_inv_val = (abs_inv_val > max_abs_inv_val) ? abs_inv_val : max_abs_inv_val;
colmem += N;
}
out_rcond = T(1) / (max_abs_src_val * max_abs_out_val);
out_rcond = T(1) / (max_abs_src_val * max_abs_inv_val);
return true;
}