diff --git a/include/armadillo_bits/op_log_det_meat.hpp b/include/armadillo_bits/op_log_det_meat.hpp index 69c3250b..dd1a6d25 100644 --- a/include/armadillo_bits/op_log_det_meat.hpp +++ b/include/armadillo_bits/op_log_det_meat.hpp @@ -204,6 +204,8 @@ op_log_det_sympd::apply_direct(typename T1::pod_type& out_val, const Base::yes) { arma_debug_warn_level(1, "log_det_sympd(): given matrix is not hermitian"); } } + // TODO: for complex matrices, ensure imaginary components on the diagonal are zero (within tolerance?) + return auxlib::log_det_sympd(out_val, A); }