diff --git a/src/algebra/TinyMatrix.hpp b/src/algebra/TinyMatrix.hpp index 38d8c28393bb059b636ce8885b0cb68da5ad79a1..9ca3746c7bf5acb6cfa78ca9dc3de0da224d88e7 100644 --- a/src/algebra/TinyMatrix.hpp +++ b/src/algebra/TinyMatrix.hpp @@ -20,8 +20,7 @@ class TinyMatrix static_assert((N > 0), "TinyMatrix size must be strictly positive"); PUGS_FORCEINLINE - constexpr size_t - _index(const size_t& i, const size_t& j) const noexcept + constexpr size_t _index(const size_t& i, const size_t& j) const noexcept // LCOV_EXCL_LINE (due to forced inline) { return i * N + j; } @@ -343,9 +342,7 @@ det(const TinyMatrix<N, T>& A) determinent *= -1; } } - const size_t I = index[i]; - if (M(I, i) == 0) - return 0; + const size_t I = index[i]; const T inv_Mii = 1. / M(I, i); for (size_t k = i + 1; k < N; ++k) { const size_t K = index[k];