Improve error messages and add bound violation tests

src/analysis/Polynomial1D.hpp

@@ -202,14 +202,14 @@ class [[nodiscard]] Polynomial1D
   PUGS_INLINE
   double& coefficient(const size_t i)
   {
-    Assert(i < m_coefficients.size());
+    Assert(i < m_coefficients.size(), "invalid coefficient number");
     return m_coefficients[i];
   }
 
   PUGS_INLINE
   const double& coefficient(const size_t i) const
   {
-    Assert(i < m_coefficients.size());
+    Assert(i < m_coefficients.size(), "invalid coefficient number");
     return m_coefficients[i];
   }
 

tests/test_Polynomial1D.cpp

@@ -254,6 +254,15 @@ TEST_CASE("Polynomial1D", "[analysis]")
   }
 
 #ifndef NDEBUG
-  SECTION("checking for bounds validation") {}
+  SECTION("checking for bounds validation")
+  {
+    Polynomial1D p({-1, 0, 0, 1});
+    REQUIRE_THROWS_AS(p.coefficient(4), AssertError);
+
+    const Polynomial1D q({-1, 0, 0, 1});
+    REQUIRE_THROWS_AS(q.coefficient(4), AssertError);
+
+    REQUIRE_THROWS_AS(Polynomial1D{std::vector<double>{}}, AssertError);
+  }
 #endif   // NDEBUG
 }