Function domainReductionTest

test/bsplinetestingutilities.cpp:100–141

}

bool domainReductionTest(BSpline &bs, const BSpline &bs_orig)
{
    if (bs.getNumVariables() != 2 || bs_orig.getNumVariables() != 2)
        return false;

    // Check for error
    if (!compareBSplines(bs, bs_orig))
        return false;

    auto lb = bs.getDomainLowerBound();
    auto ub = bs.getDomainUpperBound();

    bool flag = false;
    unsigned int index = 0;
    for (; index < lb.size(); index++)
    {
        if (ub.at(index)-lb.at(index) > 1e-1)
        {
            flag = true;
            break;
        }
    }

    if (flag)
    {
        auto split = (ub.at(index) + lb.at(index))/2;

        auto lb2 = lb;
        auto ub2 = ub; ub2.at(index) = split;
        BSpline bs2(bs);
        bs2.reduceSupport(lb2, ub2);

        auto lb3 = lb; lb3.at(index) = split;
        auto ub3 = ub;
        BSpline bs3(bs);
        bs3.reduceSupport(lb3, ub3);

        return (domainReductionTest(bs2, bs_orig) && domainReductionTest(bs3, bs_orig));
    }

    return true;
}

bool runRecursiveDomainReductionTest()
{