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

fix

Browse Source
This commit is contained in:
conrad
2021-03-22 10:55:49 +10:00
parent 0309803163
commit 9e66f44bc2
1 changed files with 8 additions and 8 deletions
Download Patch File Download Diff File Expand all files Collapse all files
include/armadillo_bits/op_inv_meat.hpp
+8 -8
View File
@@ -64,13 +64,13 @@ op_inv::apply_direct(Mat<typename T1::elem_type>& out, const Base<typename T1::e
arma_debug_check( (out.is_square() == false), "inv(): given matrix must be square sized" );
bool status = false;
if((out.n_rows <= 4) && is_cx<eT>::no)
{
const Mat<eT> tmp = out;
Mat<eT> tmp(out.n_rows, out.n_rows);
status = op_inv::apply_tiny_noalias(out, tmp);
const bool status = op_inv::apply_tiny_noalias(tmp, out);
if(status) { arrayops::copy(out.memptr(), tmp.memptr(), tmp.n_elem); return true; }
}
else
if(out.is_diagmat())
@@ -320,14 +320,14 @@ op_inv_sympd::apply_direct(Mat<typename T1::elem_type>& out, const Base<typename
if((out.n_rows <= 4) && is_cx<eT>::no)
{
const Mat<eT> tmp = out;
Mat<eT> tmp(out.n_rows, out.n_rows);
const bool status = op_inv::apply_tiny_noalias(out, tmp);
const bool status = op_inv::apply_tiny_noalias(tmp, out);
if(status) { return true; }
if(status) { arrayops::copy(out.memptr(), tmp.memptr(), tmp.n_elem); return true; }
}
return auxlib::inv_sympd(out, expr.get_ref());
return auxlib::inv_sympd(out);
}