c++ - 运行简单的 .oct 文件时 Octave 崩溃，我的代码有什么问题？

Question

我在编写一个更长、更复杂的程序(FEM 求解器)时发现了这个问题，但为了解决它，我开始删除代码，直到我有了这个以同样方式崩溃的简单示例。

下面的 .oct 文件 大约每 4 次崩溃 1 次，但它是随机的，有时它可以运行 50 次而不会崩溃。当它崩溃时，它会在没有警告的情况下关闭 Octave，因此我无法调试它。 (我拥有更多变量的较大代码大约有 80% 的时间会崩溃)。

我已经在适用于 Windows 10 x64 的 Octave 5.1.0 和 5.1.1(非官方) 以及至少一台其他计算机上尝试过此操作。我在 CLI 或 GUI 模式下运行 Octave 也没关系。

我希望我犯了一些简单的错误......但我似乎找不到它。

.m 脚本来检查它何时崩溃:

for i = 1:50
  a = test7(5)
end 

编译成.oct的test7.cc文件:

#include <octave/oct.h>

DEFUN_DLD (test7, args, nargout,
           "Test code that crashes...")
{
  if (args.length () != 1) {
    print_usage ();    
    return octave_value();
  }

  const int N = args(0).int_value();        
  const int savesTot = 500000;
  int iter = 0; 

//Declaration method 1 (also crashes)
  //OCTAVE_LOCAL_BUFFER (double, Q, N); 
//Declaration method 2 (also crashes)
  //double *Q = new double(N);
//Declaration method 3 (crashes)
  NDArray Q_o (dim_vector (N,1));
  double *Q = Q_o.fortran_vec();

  //Q_o(0) = 1.1;   // Calling NDArray class (also crashes)
  Q[0] = 1.1;       // Calling pointer to raw data directly (faster)

  while (iter < savesTot) {

    for (int i = 0; i<N; i++) {    
      if (i == 0) {        
        //Q_o(i+1)      = 1.0;  
          Q[i+1]        = 1.0;        
      } else if ( i == (N-1) ) { 
        //Q_o(i+1)      = 2.0;
          Q[i+1]        = 2.0;  
      } else {        
        //Q_o(i+1)      = 3.0;  
          Q[i+1]        = 3.0;         
      }
    }  

    iter++;
    OCTAVE_QUIT;        // Check if user has pressed Ctrl+C
  }

  octave_value_list retval;    // Define return list (can contain many variables)
  retval(0) = iter; 

  if (! error_state) { return retval; }

  return octave_value_list ();
}

test7.cc 使用以下命令在 Windows 10 x64 上的 Octave 5.1.1 或 5.1.0 中编译为 test7.oct:

mkoctfile test7.cc

Answer 1

你的问题是你写的越界了:

    for (int i = 0; i<N; i++) {    
       //...
       Q[i+1]        = 3.0;
       //...
    }

你可能想循环 i=1; i<N; i++和索引Q[i] , 或循环 i=0; i<N-1; i++ .