diff --git a/include/armadillo_bits/auxlib_meat.hpp b/include/armadillo_bits/auxlib_meat.hpp index 8a70da4c..72bd2e75 100644 --- a/include/armadillo_bits/auxlib_meat.hpp +++ b/include/armadillo_bits/auxlib_meat.hpp @@ -28,6 +28,10 @@ auxlib::inv(Mat& A) { arma_debug_sigprint(); + // NOTE: given a matrix with NaN values, getrf() and getri() do not necessarily fail, + // NOTE: and can produce matrices with NaN values. + // NOTE: we're not checking for non-finite values to avoid breaking existing user code. + if(A.is_empty()) { return true; } #if defined(ARMA_USE_LAPACK)