Method lookup

tools/Wires/Inflator/GeometryCorrectionTable.cpp:196–213

}

Vector2F GeometryCorrectionTable::lookup(Float half_width, Float half_height) {
    const size_t num_entries = m_measurements.size();
    assert(m_design_sizes.size() == num_entries);
    Vector2F target(half_width*2, half_height*2);

    std::vector<Float> dist(num_entries);
    std::vector<size_t> indices(num_entries);
    for (size_t i=0; i<num_entries; i++) {
        const Vector2F& measure = m_measurements[i];
        dist[i] = (target - measure).squaredNorm();
        indices[i] = i;
    }

    auto index_compare = [=](size_t i, size_t j) { return dist[i] < dist[j]; };
    std::sort(indices.begin(), indices.end(), index_compare);

    return interpolate(target, indices[0], indices[1], indices[2]) * 0.5;
}

Vector2F GeometryCorrectionTable::interpolate(const Vector2F& target,
        size_t idx1, size_t idx2, size_t idx3) {