Octave аварийно завершает работу при запуске простого файла .oct, что не так в моем коде? - PullRequest
Новая
       30

Octave аварийно завершает работу при запуске простого файла .oct, что не так в моем коде?

0 голосов
/ 15 октября 2019

Я обнаружил проблему во время написания более длинной и более сложной программы (FEM Solver), но для ее устранения я начал удалять код, пока у меня не появился этот простой пример, который вылетает таким же образом.

СледующееФайл .oct дает сбой примерно 1 в 4 раза, но это случайный случай, и иногда он может работать 50 раз без сбоя . Когда он падает, он закрывает Octave без предупреждений, поэтому я не смог его отладить. (Большой код, который у меня есть с большим количеством переменных, вылетает примерно в 80% случаев).

Я пробовал это в Octave 5.1.0 и 5.1.1 (неофициально) для Windows 10 x64 и хотя бы на одном другом компьютере. Также не имеет значения, запускаю ли я Octave в режиме CLI или GUI.

Я надеюсь, что совершил простую ошибку ... но я не могу ее найти.

.m скрипт, чтобы проверить, когда он падает: 

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

test7.cc файл, скомпилированный в .oct: 

#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 скомпилирован в test7.oct внутри Octave 5.1.1 или 5.1.0 в Windows 10 x64 с помощью следующей команды: 

mkoctfile test7.cc

1 Ответ

1 голос
/ 15 октября 2019

Ваша проблема в том, что вы пишете вне границ: 

    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++.

...