gausser_impl_3_test.cpp

#include "gausser_impl_3.h"

#include "gausser.h"

// Boost.Test does not play well with -Weffc++


#include <boost/test/unit_test.hpp>

using namespace ribi;

BOOST_AUTO_TEST_CASE(ribi_gausser_impl_3)
{
  const double sd{1.0};
  const gausser_impl_3 g(sd);
  // Must be in range of gausser, which is [-4.0,4.0]
  for (double x = -3.9; x < 3.9; x += 0.001)
  {
    const double expected{gauss(x, sd)};
    const double measured{g(x)};
    BOOST_CHECK_CLOSE(measured, expected, 0.000001);
  }
  // Extremes should also work
  BOOST_CHECK_NO_THROW(g(-5.0));
  BOOST_CHECK_NO_THROW(g(5.0));
}