Function TEST

test/unit/MathUtilsTest.cpp:43–92

{

TEST(MathUtilsTest, barycentric)
{
    using namespace math;

    // Triangle (1, 1), (5, 4), (3, 6) -- z's 55, 25, 100
    double x1 = 1;
    double y1 = 1;
    double z1 = 55;
    double x2 = 5;
    double y2 = 4;
    double z2 = 25;
    double x3 = 3;
    double y3 = 6;
    double z3 = 100;

    struct point
    {
        double x;
        double y;
    };

    std::vector<point> points { { 3, 3 }, {4, 2}, {5, 4}, {4, 6}, {1, 5},
        {4, 5}, {2, 7.0 / 4}, {3, 4}, {3, 6} };

    double inf = std::numeric_limits<double>::infinity();
    std::vector<double> results { 48.57142857142, inf, 25, inf, inf,
        62.5, 47.5, 65.7142857142, 100 };

    double z;
    for (size_t i = 0; i < points.size(); ++i)
    {
        const point& p = points[i];
        z = barycentricInterpolation(x1, y1, z1, x2, y2, z2, x3, y3, z3, p.x, p.y);
        if (std::isinf(results[i]))
            EXPECT_TRUE(std::isinf(z));
        else
            EXPECT_NEAR(z, results[i], .0000000001);
    }

    // Re-order triangle points (x2 is before x1 in input).  Results should be the same.
    for (size_t i = 0; i < points.size(); ++i)
    {
        const point& p = points[i];
        z = barycentricInterpolation(x2, y2, z2, x1, y1, z1, x3, y3, z3, p.x, p.y);
        if (std::isinf(results[i]))
            EXPECT_TRUE(std::isinf(z));
        else
            EXPECT_NEAR(z, results[i], .0000000001);
    }
}

TEST(MathUtilsTest, bary_issue_4694)
{