Method compute

tools/Wires/Parameters/VertexOffsetParameterDerivative.cpp:17–60

using namespace VertexOffsetParameterDerivativeHelper;

MatrixFr VertexOffsetParameterDerivative::compute() {
    WireNetwork::Ptr wire_network = m_parameter->get_wire_network();

    const size_t dim = m_mesh->get_dim();
    const size_t num_mesh_vertices = m_mesh->get_num_vertices();
    const size_t num_mesh_faces = m_mesh->get_num_faces();
    const size_t num_wire_vertices = wire_network->get_num_vertices();

    assert(dim == wire_network->get_dim());

    VectorI roi = m_parameter->get_roi();
    BoolVector in_roi(num_wire_vertices, false);
    const size_t roi_size = roi.size();
    for (size_t i=0; i<roi_size; i++) {
        in_roi[roi[i]] = true;
    }

    VectorF weights = VectorF::Zero(num_mesh_vertices);
    MatrixFr derivative_v = MatrixFr::Zero(num_mesh_vertices, dim);
    MatrixFr wire_derivative = m_parameter->compute_derivative();

    for (size_t i=0; i< num_mesh_faces; i++) {
        int source = m_face_source[i];
        if (source < 0) {
            // Source is edge
            size_t edge_idx = -source - 1;
            compute_derivative_on_edge(
                    edge_idx, i, in_roi, wire_derivative,
                    derivative_v, weights);
        } else if (source > 0) {
            // Source is vertex
            size_t vertex_idx = source - 1;
            compute_derivative_on_vertex(
                    vertex_idx, i, in_roi, wire_derivative,
                    derivative_v, weights);
        }
    }

    for (size_t i=0; i<num_mesh_vertices; i++) {
        if (weights[i] > 0)
            derivative_v.row(i) /= weights[i];
    }
    return derivative_v;
}

void VertexOffsetParameterDerivative::compute_derivative_on_edge(
        size_t wire_edge_index,