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

proposed fix

Browse Source
This commit is contained in:
conrad
2025-04-23 01:17:26 +10:00
parent bd912040a5
commit eb8ea83e11
1 changed files with 7 additions and 3 deletions
Download Patch File Download Diff File Expand all files Collapse all files
include/armadillo_bits/op_expmat_meat.hpp
+7 -3
View File
@@ -124,11 +124,15 @@ op_expmat::apply_direct(Mat<typename T1::elem_type>& out, const Base<typename T1
if(arma_isfinite(norm_val) == false) { return false; }
const double log2_val = (norm_val > T(0)) ? double(eop_aux::log2(norm_val)) : double(0);
// const double log2_val = (norm_val > T(0)) ? double(eop_aux::log2(norm_val)) : double(0);
//
// int exponent = int(0); std::frexp(log2_val, &exponent);
//
// const uword s = uword( (std::max)(int(0), exponent + int(1)) );
int exponent = int(0); std::frexp(log2_val, &exponent);
int exponent = int(0); std::frexp(norm_val, &exponent);
const uword s = uword( (std::max)(int(0), exponent + int(1)) );
const uword s = uword(exponent);
A /= eT(eop_aux::pow(double(2), double(s)));