ceres-solver - Ceres-Solver 因“终止：残差和雅可比评估失败”而失败。当非静态变量用于定义雅可比 for 循环时

Question

我有一个非线性系统，我正在使用 Ceres 来解决。这是一个具有稀疏块结构的稀疏系统。由于我也在处理图像数据，因此我的代码基于“去噪.cc”示例。

我遇到的问题是我的代码因“终止：残差和雅可比评估失败。”而失败。我可以通过在 Evaluate 中硬编码变量 'num_weights' 来解决这个问题。

当我在一个或多个像素上调用此函数时，问题仍然存在。对于每个像素，我的权重都不同。

任何关于为什么会这样的见解都会有所帮助。

谢谢！

Cost::Cost(const std::vector<double> &weights) : _weights(weights)
{
    set_num_residuals(1);
    mutable_parameter_block_sizes()->push_back(1); //has more parameters than weights

    for (int i = 0; i < _weights.size(); ++i)
        mutable_parameter_block_sizes()->push_back(1);
}

bool Cost::Evaluate(double const* const* parameters,
                             double *residuals,
                             double **jacobians) const
{
    int num_weights = (int)_weights.size();

    float d0 = parameters[0][0];
    residuals[0] = d0;

    for (int i = 0; i < num_weights; ++i)
    {
        residuals[0] += parameters[i+1][0];
    }

    if (jacobians != NULL)
    {
        for (int i = 0; i < num_weights+1; ++i)
        {
            if (jacobians[i] != NULL)
            {
                jacobians[i][0] = 0;
            }
        }
    }

    return true;
}

Answer 1

刚刚想通了。

基本上，我的向量“std::vector weights”超出了范围。ceres 成本函数不保留向量的所有权，因此“num_weights”评估为 0。