Function compute_min_angles

tools/Wires/Attributes/WireVertexMinAngleAttribute.cpp:46–67

    template<typename T>
    void compute_min_angles(const WireNetwork& network, MatrixFr& values,
            T& angle) {
        const size_t num_vertices = network.get_num_vertices();

        // By default, min angle is the max angle possible: 180 degrees.
        values = MatrixFr::Ones(num_vertices, 1) * M_PI;
        const MatrixFr& vertices = network.get_vertices();
        for (size_t i=0; i<num_vertices; i++) {
            const VectorF& v = vertices.row(i);
            VectorI neighbors = network.get_vertex_neighbors(i);

            const size_t valance = neighbors.size();
            for (size_t j=0; j<valance; j++) {
                VectorF ej = vertices.row(neighbors[j]) - v.transpose();
                for (size_t k=j+1; k<valance; k++) {
                    VectorF ek = vertices.row(neighbors[k]) - v.transpose();
                    values.coeffRef(i,0) =
                        std::min(values.coeffRef(i,0), angle(ej, ek));
                }
            }
        }
    }
}

using namespace WireVertexMinAngleAttributeHelper;