python - 勒让德多项式的根c++

我正在编写一个程序，使用c++查找n阶勒让德多项式的根。我的代码附在下面:

double* legRoots(int n)
{
 double myRoots[n];
 double x, dx, Pi = atan2(1,1)*4;
 int iters = 0;
 double tolerance = 1e-20;
 double error = 10*tolerance;
 int maxIterations = 1000;

 for(int i = 1; i<=n; i++)
 {
  x = cos(Pi*(i-.25)/(n+.5));
  do
  {
   dx -= legDir(n,x)/legDif(n,x);
   x += dx;
   iters += 1;
   error = abs(dx);
  } while (error>tolerance && iters<maxIterations);
  myRoots[i-1] = x;
 }
 return myRoots;
}

假设我确实有起作用的勒让德多项式和勒让德多项式导数生成函数，但我认为这样做会使代码文本难以理解。从返回数组计算值的意义上讲，此函数起作用，但是它们疯狂地关闭了，输出以下内容:

3.95253e-323
6.94492e-310
6.95268e-310
6.42285e-323
4.94066e-323
2.07355e-317

我在Python中编写的等效函数给出了以下内容:

[-0.90617985 -0.54064082  0.          0.54064082  0.90617985]

我希望另一组眼睛可以帮助我了解C++代码中的问题是什么导致值急剧下降。我在Python代码中所做的工作与在C++中所做的没有任何不同，因此，感谢任何人对此提供的任何帮助。作为参考，我主要尝试模仿Rosetta代码中关于高斯求积:http://rosettacode.org/wiki/Numerical_integration/Gauss-Legendre_Quadrature的方法。

最佳答案

您正在将地址返回到堆栈中的临时变量

{
    double myRoots[n];
    ...
    return myRoots; // Not a safe thing to do
}

我建议将您的函数定义更改为

void legRoots(int n, double *myRoots)

省略return语句，并在调用函数之前定义myroots

double myRoots[10];
legRoots(10, myRoots);

选项2是使用new或malloc动态分配myRoots。