Function TEST

test/unit/filters/RadiusAssignFilterTest.cpp:12–65

{

TEST(RadiusAssignFilterTest, basic_usage)
{
    Options ro;
    ro.add("filename", Support::datapath("filters/radiusassign/grid_20x10x10.txt"));
    StageFactory factory;
    Stage& r = *(factory.createStage("readers.text"));
    r.setOptions(ro);

    std::vector<std::string> is3ds = {"false", "false", "true", "true"};
    std::vector<double> radVals = {1, 3.5, 1, 3.5};
    std::vector<unsigned int> nbUpdatedPoints = {0, 0, 0, 0};
    std::vector<unsigned int> nbExpectedUpdatedPoints = {10, 370, 1, 179};

    for (size_t ii = 0; ii < radVals.size(); ii++)
    {
        Options ao;
        ao.add("value", "Classification=1");
        ao.add("where", "(X == 5) && (Y == 5) && (Z == 5)");

        Stage& a = *(factory.createStage("filters.assign"));
        a.setInput(r);
        a.setOptions(ao);

        Options fo;
        fo.add("is3d", is3ds[ii]);
        fo.add("reference_domain", "Classification[1:1]");
        fo.add("update_expression", "Classification = 2");
        fo.add("radius", radVals[ii]);

        Stage& f = *(factory.createStage("filters.radiusassign"));
        f.setInput(a);
        f.setOptions(fo);

        PointTable table;
        f.prepare(table);
        PointViewSet viewSet = f.execute(table);
        PointViewPtr view = *viewSet.begin();

        EXPECT_EQ(1u, viewSet.size());

        PointRef point(*view, 0);

        for (PointId id = 0; id < view->size(); ++id)
        {
            point.setPointId(id);
            if (view->getFieldAs<unsigned int>(Dimension::Id::Classification, id) == 2)
                nbUpdatedPoints[ii] += 1;
        }
        // Check that some points are updated
        EXPECT_TRUE(nbUpdatedPoints[ii] == nbExpectedUpdatedPoints[ii])
            << "Found: '" << nbUpdatedPoints[ii] << "'" << std::endl
            << "expected: '" << nbExpectedUpdatedPoints[ii] <<"'" << std::endl;
    }
}