计算偏导数时出现分段错误

Question

我正在尝试计算函数 f(x,y)=x+y w.r.t x 和 y 在 (3, 2) 点的偏导数。我编写了代码，但它给我一个名为 Segmentation fault (core dumped) 的错误。我不知道那是什么。谁能指出这有什么问题吗？

#include<iostream>
#include<cmath>
#include<vector>
using namespace std;

class Der{

private:
    double f;  // function value at x
    vector <double> df; // derivative of function at x
    
public:
    Der();
    Der(vector <double>* v);
    Der operator+(Der); // f + g
    
    void print();
    
};

Der :: Der(){}

Der :: Der(vector <double>* v){
    this->f = v->at(0);
    for(int i=1; i < v->size(); i++){
        this -> df[i-1] = v->at(i);
    }
}


Der Der :: operator+(Der g){
    Der h;
    h.f = this->f + g.f;
    for(int i=0; i< df.size(); i++){
        h.df[i] = this->df[i] + g.df[i];
    }
    return h;
}

void Der :: print(){
    cout<<"Function value at given point is : "<<f<<endl;
    for(int i=0; i< df.size(); i++){
    cout<<"Derivative at given point is : "<<df[i]<<endl;
    }
}

int main()
{
    Der f;
    vector <double> t {3,1,0};
    vector <double> k {2,0,1};
    Der x(&t),y(&k);
    
    f = x+y;
    f.print();
}

非常感谢！

Answer 1

在您的构造函数中，您的向量 df 是空的，但您尝试使用

将值放入向量中

this -> df[i-1] = v->at(i);

这将导致未定义的行为，因为您正在访问空向量的边界之外。一种修复方法是用以下内容替换该行：

df.push_back(v->at(i));

这将在每次调用时增加向量的大小。

同样在 Der Der :: operator+(Der g){ 你有：

h.df[i] = this->df[i] + g.df[i];

与构造函数有相同的错误。 h.df 为空。您可以使用以下方法解决该问题：

h.df.push_back(this->df[i] + g.df[i]);

也像 c++ 中提到的 @πìντα ῥεῖ 一样，我们更喜欢通过指针传递 const 引用：

 Der :: Der(const vector <double>& v){

计算偏导数时出现分段错误

Segmentation error while calculating partial derivative

c++

scientific-computing