diff --git a/include/armadillo_bits/op_inv_meat.hpp b/include/armadillo_bits/op_inv_meat.hpp index e8af684c..5edd6644 100644 --- a/include/armadillo_bits/op_inv_meat.hpp +++ b/include/armadillo_bits/op_inv_meat.hpp @@ -318,6 +318,9 @@ op_inv_sympd::apply_direct(Mat& out, const Base::yes) { arma_debug_warn_level(1, "inv_sympd(): given matrix is not hermitian"); } } + // TODO: for complex matrices, warn if imaginary components of diagonal are greater than eps? + // TODO: however, this may be inconsistent with the check in Base::is_sympd(), which uses data derived threshold + if((out.n_rows <= 4) && is_cx::no) { Mat tmp(out.n_rows, out.n_rows, arma_nozeros_indicator());