test_math_functions.cpp

// Description: Testing for math functions
// DocumentationOf: math_functions.h

#include "test/test_init.h"

#include "pastel/sys/random.h"
#include "pastel/sys/math_functions.h"

TEST_CASE("CcwAngle (CcwAngle)")
{
    REQUIRE(radiansToDegrees<real>(
        ccwAngle(Vector2(1, 0), Vector2(0, 1))) - 90 < 0.001);
    REQUIRE(radiansToDegrees<real>(
        ccwAngle(Vector2(0, 1), Vector2(-1, 0))) - 90 < 0.001);
    REQUIRE(radiansToDegrees<real>(
        ccwAngle(Vector2(-1, 0), Vector2(0, -1))) - 90 < 0.001);
    REQUIRE(radiansToDegrees<real>(
        ccwAngle(Vector2(0, -1), Vector2(1, 0))) - 90 < 0.001);

    REQUIRE(radiansToDegrees<real>(
        ccwAngle(Vector2(1, 0), Vector2(1, 1))) - 45 < 0.001);
    REQUIRE(radiansToDegrees<real>(
        ccwAngle(Vector2(1, 0), Vector2(-1, 1))) - 135 < 0.001);
    REQUIRE(radiansToDegrees<real>(
        ccwAngle(Vector2(1, 0), Vector2(-1, -1))) - 225 < 0.001);
    REQUIRE(radiansToDegrees<real>(
        ccwAngle(Vector2(1, 0), Vector2(1, -1))) - 315 < 0.001);
}

TEST_CASE("Quadratic (Quadratic)")
{
    integer failings = 0;

    for (integer i = 0;i < 10000;++i)
    {
        real aRoot = random<real>() * 100 - 50;
        real bRoot = random<real>() * 100 - 50;

        if (bRoot < aRoot)
        {
            std::swap(aRoot, bRoot);
        }

        const real aCoeff = random<real>() * 100 + 2;
        const real bCoeff = (-aRoot - bRoot) * aCoeff;
        const real cCoeff = (aRoot * bRoot) * aCoeff;

        real t0 = 0;
        real t1 = 0;

        quadratic(aCoeff, bCoeff, cCoeff, t0, t1);

        real aDelta = abs(aRoot - t0);
        real bDelta = abs(bRoot - t1);

        if (aDelta >= 0.1 || bDelta >= 0.1)
        {
            ++failings;
        }
    }

    REQUIRE(failings <= 10);
}