分段错误（核心已转储）- 向量

Question

我正在研究可以使用基本段计算曲线长度的算法。所以如果我有曲线的 x 和 y 坐标向量，我需要计算这个基本段的数量。我为此发明了自己的递归算法。有一个代码：

#include <iostream>
#include <vector>
#include <cmath>

using namespace std;
//***********************************************
inline double Pitagoras(double xp, double yp, double xf, double yf)
{
    return sqrt((xp - xf)*(xp - xf) + (yp - yf)*(yp - yf));
}
//***************************************************
inline double calculateX(double xp, double yp, double a, double b, const double eps)
{
   double delta;
   double x1, x2;

     delta = (-2.0 * xp + 2.0 *a*b - 2.0 *a*yp)*(-2.0 * xp + 2.0 *a*b - 2.0 *a*yp)
        - 4.0* (1.0 + a*a)*(xp*xp + yp*yp + b*b - 2.0*b*yp - eps*eps);

     x1 = (-(-2.0 * xp + 2.0 *a*b - 2.0 *a*yp) - sqrt(delta))/(2.0 * (1.0 + a*a));
     x2 = (-(-2.0 * xp + 2.0 *a*b - 2.0 *a*yp) + sqrt(delta))/(2.0 * (1.0 + a*a));

   if(x1 >= xp)
     return x1;
     else
       return x2;
}
//***************************************************
inline double calculateY(double x, double a, double b)
{
  return a*x + b;
}
//***********************************************
unsigned long algorithmKolmogorow(double xp, double yp, double xf,
double yf, const double eps, vector<double> &vectorX, vector<double> &vectorY);
//***********************************************
int main()
{
    vector<double> vctrY; //vector of value of function
    vector<double> vctrX; //vector of x
    double xP,yP,xF,yF; //coordinates of two points on the curve
    const double Eps = 0.0001; //length of elementary line

      for(double x=1.0; x<=5 ;x +=0.001)
        {
          vctrX.push_back(x);
          vctrY.push_back(x*x); //f(x) = x^2
        }

    xP = vctrX[0];
    yP = vctrY[0];
    xF = vctrX[1];
    yF = vctrY[1]; //set beginning value

    cout<<algorithmKolmogorow(xP, yP, xF, yF, Eps, vctrX, vctrY)*Eps;

    return 0;
}
//***************************************************
unsigned long algorithmKolmogorow(double xp, double yp, double xf,
 double yf, const double eps, vector<double> &vectorX, vector<double> &vectorY)
{
    static unsigned long N; //licznik
    static unsigned long i = 1;
    double d;
    double a,b;

      d = Pitagoras(xp, yp, xf, yf);

         if(d >= eps){
            a = (yf - yp)/(xf - xp);
            b = yp - a*xp;
            xp = calculateX(xp, yp, a, b, eps);
            yp = calculateY(xp, a, b);
            N++;
         }

           else{
             i++;
             xf = vectorX[i];
             yf = vectorY[i];
             //cout<<i<<"\t"<<vectorX[i]<<"\t"<<vectorY[i]<<endl;
           }

             if(i < vectorX.size())
              N =  algorithmKolmogorow(xp, yp, xf, yf, eps, vectorX, vectorY);

return N;
}

如您所见，我正在为抛物线函数设置 x,y 坐标。当 Eps 很大时，例如 Eps = 0.001，一切正常。如果我设置较小的值，如 Eps = 0.0001，那么我会在主题和程序停止运行中出现错误。我完全不知道为什么。

我可以添加您需要的任何新信息（关于我的编译器、IDE、OS 等）。

拉法尔

Answer 1

         i++;
         xf = vectorX[i];
         yf = vectorY[i];
         //cout<<i<<"\t"<<vectorX[i]<<"\t"<<vectorY[i]<<endl;
       }

         if(i < vectorX.size())

这将超出向量的边界 1。也许您需要 if(i < vectorX.size()-1) 或者（注意 AndyG 的评论）也许您需要 i++ 低于 [= 的这两种用法13=].

Answer 2

由于递归，您可能运行出栈 space。

解决这个问题的一个快速技巧是在 shell 中使用 ulimit 增加堆栈大小，例如

(ulimit -s unlimited; ./my_program)

然而，真正的解决方法是删除递归。 algorithmKolmogorow 看起来它只做尾递归，它总是可以转换为循环：

unsigned long algorithmKolmogorow(double xp, double yp, double xf,
 double yf, const double eps, vector<double> &vectorX, vector<double> &vectorY)
{
    static unsigned long N; //licznik
    static unsigned long i = 1;
    double d;
    double a,b;

    while(true) {
      d = Pitagoras(xp, yp, xf, yf);

         if(d >= eps){
            a = (yf - yp)/(xf - xp);
            b = yp - a*xp;
            xp = calculateX(xp, yp, a, b, eps);
            yp = calculateY(xp, a, b);
            N++;
         }

           else{
             i++;
             if(i >= vectorX.size())
               return N;
             xf = vectorX[i];
             yf = vectorY[i];
             //cout<<i<<"\t"<<vectorX[i]<<"\t"<<vectorY[i]<<endl;
           }

    }

}

还解决了访问一个元素超出其边界的向量的问题。

还有一些问题：

N 应该初始化
algorithmKolmogorow 每次程序调用只能使用一次。从 N 和 i 中删除 static 来解决这个问题。

Answer 3

您 Eps 的值太小了。它会导致算法进展太慢，这意味着您运行在计算完成之前就出栈了 space。我将 Eps 增加到 0.005，代码运行没问题。虽然不确定它是否会产生正确的答案；你必须检查一下。

此外，您需要更改这些行：

xf = vectorX[i];
yf = vectorY[i];

为此：

xf = vectorX[i-1];
yf = vectorY[i-1];

否则，您将访问最后一次迭代中的最后一个向量元素。您可能还想对 Eps 进行公差研究，以了解它必须有多小才能获得准确的计算。

分段错误（核心已转储）- 向量

Segmentation Fault (core dumped) - vectors

c++

coredump

segmentation-fault

stdvector