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<dreal>(
        ccwAngle(Vector2(1, 0), Vector2(0, 1))) - 90 < 0.001);
    REQUIRE(radiansToDegrees<dreal>(
        ccwAngle(Vector2(0, 1), Vector2(-1, 0))) - 90 < 0.001);
    REQUIRE(radiansToDegrees<dreal>(
        ccwAngle(Vector2(-1, 0), Vector2(0, -1))) - 90 < 0.001);
    REQUIRE(radiansToDegrees<dreal>(
        ccwAngle(Vector2(0, -1), Vector2(1, 0))) - 90 < 0.001);

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

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

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

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

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

        dreal t0 = 0;
        dreal t1 = 0;

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

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

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

    REQUIRE(failings <= 10);
}