ecearth3/sources/lpjg/tests/math_test.cpp

///////////////////////////////////////////////////////////////////////////////////////
/// \file math_test.cpp
/// \brief Unit tests functionality in guessmath.h
///
/// \author Joe Siltberg
/// $Date: 2013-10-10 10:20:33 +0200 (Thu, 10 Oct 2013) $
///
///////////////////////////////////////////////////////////////////////////////////////

#include "config.h"
#include "catch.hpp"

#include "guessmath.h"

TEST_CASE("Historic/add", "Some basic tests of adding values to a Historic") {
	Historic<double, 3> history;

	REQUIRE(history.size() == 0);

	history.add(1);

	REQUIRE(history.size() == 1);
	REQUIRE(history.sum() == Approx(1));
	REQUIRE(history[0] == 1);
	REQUIRE(history.mean() == Approx(1));

	history.add(2);

	REQUIRE(history.size() == 2);
	REQUIRE(history.sum() == Approx(3));
	REQUIRE(history[0] == 1);
	REQUIRE(history[1] == 2);
	REQUIRE(history.mean() == Approx(1.5));

	history.add(3);

	REQUIRE(history.size() == 3);
	REQUIRE(history.sum() == Approx(6));
	REQUIRE(history[0] == 1);
	REQUIRE(history[1] == 2);
	REQUIRE(history[2] == 3);
	REQUIRE(history.mean() == Approx(2));

	history.add(4);

	REQUIRE(history.size() == 3);
	REQUIRE(history.sum() == Approx(9));
	REQUIRE(history[0] == 2);
	REQUIRE(history[1] == 3);
	REQUIRE(history[2] == 4);
	REQUIRE(history.mean() == Approx(3));
}

TEST_CASE("variation_coefficient", "Tests of variation coefficient") {

	double single_value[] = { 7 };

	REQUIRE(variation_coefficient(single_value, 1) == Approx(-1));

	double two_values[] = { 7 , 7};

	REQUIRE(variation_coefficient(two_values, 2) == Approx(0));

	double values[] = { 2, 4, 4, 4, 5, 5, 7, 9 };

	REQUIRE(variation_coefficient(values, 8) == Approx(0.427618));
}
TEST_CASE("Historic/periodic", "Some basic tests of the period functions in Historic") {
	Historic<double, 3> history;

	REQUIRE(history.size() == 0);

	history.add(1);
	history.add(2);
	history.add(3);

	REQUIRE(history.periodicsum(2) == Approx(5));
	REQUIRE(history.periodicmean(2) == Approx(2.5));
	REQUIRE(history.min() == Approx(1));
	REQUIRE(history.max() == Approx(3));
	history.add(4);
	REQUIRE(history.periodicsum(2) == Approx(7));
	REQUIRE(history.periodicmean(2) == Approx(3.5));
	REQUIRE(history.lastadd() == Approx(4));
	REQUIRE(history.min() == Approx(2));
	REQUIRE(history.max() == Approx(4));
}