test_Polynomial.cpp

#include <catch2/catch_test_macros.hpp>

#include <Kokkos_Core.hpp>

#include <utils/PugsAssert.hpp>
#include <utils/Types.hpp>

#include <algebra/TinyMatrix.hpp>
#include <analysis/Polynomial.hpp>

// Instantiate to ensure full coverage is performed
template class Polynomial<0>;
template class Polynomial<1>;
template class Polynomial<2>;
template class Polynomial<3>;
template class Polynomial<4>;
template class Polynomial<5>;

// clazy:excludeall=non-pod-global-static

TEST_CASE("Polynomial", "[analysis]")
{
  SECTION("construction")
  {
    REQUIRE_NOTHROW(Polynomial<2>{TinyVector<3>{2, 3, 4}});
  }
  SECTION("degree")
  {
    Polynomial<2> P(TinyVector<3>{2, 3, 4});
    REQUIRE(P.degree() == 2);
  }
  SECTION("equality")
  {
    Polynomial<2> P(TinyVector<3>{2, 3, 4});
    Polynomial<2> Q(TinyVector<3>{2, 3, 4});
    Polynomial<2> S(TinyVector<3>{2, 3, 5});

    REQUIRE(P == Q);
    REQUIRE(P != S);
  }
  SECTION("addition")
  {
    Polynomial<2> P(TinyVector<3>{2, 3, 4});
    Polynomial<2> Q(TinyVector<3>{-1, -3, 2});
    Polynomial<2> S(TinyVector<3>{1, 0, 6});
    Polynomial<3> T(TinyVector<4>{0, 3, 1, -2});
    Polynomial<3> U(TinyVector<4>{2, 6, 5, -2});
    REQUIRE(S == (P + Q));
    REQUIRE((T + P) == U);
  }
  SECTION("opposed")
  {
    Polynomial<2> P(TinyVector<3>{2, 3, 4});
    Polynomial<2> Q = -P;
    REQUIRE(Q == Polynomial<2>(TinyVector<3>{-2, -3, -4}));
  }
  SECTION("difference")
  {
    Polynomial<2> P(TinyVector<3>{2, 3, 4});
    Polynomial<2> Q(TinyVector<3>{3, 4, 5});
    Polynomial<2> D(TinyVector<3>{-1, -1, -1});
    REQUIRE(D == (P - Q));
    Polynomial<3> R(TinyVector<4>{2, 3, 4, 1});
    REQUIRE(D == (P - Q));
    REQUIRE((P - R) == Polynomial<3>(TinyVector<4>{0, 0, 0, -1}));
    R -= P;
    REQUIRE(R == Polynomial<3>(TinyVector<4>{0, 0, 0, 1}));
  }
  SECTION("product_by_scalar")
  {