Method process

kernels/TileKernel.cpp:198–284

We calculate the origin specially in order to avoid a "first point" check for every point iteration, seeing as we might have BILLIONS of points to process.

// check for every point iteration, seeing as we might have BILLIONS
// of points to process.
void TileKernel::process(const Readers& readers)
{
    using std::placeholders::_1;
    using std::placeholders::_2;
    using std::placeholders::_3;
    SplitterFilter::PointAdder adder =
        std::bind(&TileKernel::adder, this, _1, _2, _3);

    bool haveOrigin(false);
    StageWrapper::ready(m_splitter, m_table);
    for (auto&& rp : readers)
    {
        Streamable& r = *(rp.second);
        std::vector<bool> skips(m_table.capacity());
        PointId idx(0);
        PointRef point(m_table, idx);

        StreamableWrapper::ready(r, m_table);
        if (m_repro)
            StreamableWrapper::spatialReferenceChanged(*m_repro,
                r.getSpatialReference());

        // Read first point.
        bool finished(false);
        finished = !StreamableWrapper::processOne(r, point);
        if (!haveOrigin && !finished)
        {
            if (std::isnan(m_xOrigin))
                m_xOrigin = point.getFieldAs<double>(Dimension::Id::X);
            if (std::isnan(m_yOrigin))
                m_yOrigin = point.getFieldAs<double>(Dimension::Id::Y);
            m_splitter.setOrigin(m_xOrigin, m_yOrigin);
            haveOrigin = true;
        }

        idx++;

        while (!finished)
        {
            // Read subsequent points.
            while (true)
            {
                point.setPointId(idx);
                finished = !StreamableWrapper::processOne(r, point);
                if (finished || idx == m_table.capacity() - 1)
                    break;
                idx++;
            }

            // "end" points just beyond the last valid point ID
            PointId end = finished ? idx : m_table.capacity();

            SpatialReference srs = r.getSpatialReference();
            if (!srs.empty())
                m_table.setSpatialReference(srs);
            // Reproject if necessary.
            if (m_repro)
            {