Method compute_derivative

tools/Wires/Parameters/VertexOffsetParameter.cpp:42–65

}

MatrixFr VertexOffsetParameter::compute_derivative() const {
    const VectorF bbox_min = m_wire_network->get_bbox_min();
    const VectorF bbox_max = m_wire_network->get_bbox_max();
    const VectorF center = 0.5 * (bbox_min + bbox_max);
    const VectorF half_bbox_size = bbox_max - center;

    const size_t dim = m_wire_network->get_dim();
    const size_t num_vertices = m_wire_network->get_num_vertices();
    const size_t roi_size = m_roi.size();
    const MatrixFr& vertices = m_wire_network->get_vertices();
    assert(m_axis < dim);

    MatrixFr derivative = MatrixFr::Zero(num_vertices, dim);

    for (size_t i=0; i<roi_size; i++) {
        size_t v_idx = m_roi[i];
        assert(v_idx < num_vertices);
        const VectorF& v = vertices.row(v_idx);
        Float sign = v[m_axis] - center[m_axis] > 0.0 ? 1.0 : -1.0;
        derivative(m_roi[i], m_axis) = sign * half_bbox_size[m_axis];
    }

    return derivative;
}