OpenMP в Mex: StackOverflow ошибка

Question

я получил следующую часть кода, получая мне ошибку переполнения стека

#pragma omp parallel shared(Mo1, Mo2, sum_normalized_p_gn, Data, Mean_Out,Covar_Out,Prior_Out, det) private(i) num_threads(number_threads) 
{ 

     //every thread has a new copy 
     double* normalized_p_gn = (double*)malloc(NMIX*sizeof(double)); 
     #pragma omp critical 
     { 
      int id = omp_get_thread_num(); 
      int threads = omp_get_num_threads(); 
      mexEvalString("drawnow"); 
     } 

     #pragma omp for 
     //some parallel process..... 

} 

переменные, объявленные в совместно создаются таНос. и они потребляют с большим объемом памяти

есть 2 вопроса относительно вышеуказанного кода. 1) почему это привело бы к возникновению ошибки переполнения стека (т. Е. Ошибки сегментации), прежде чем он перейдет в параллельный цикл? он отлично работает, когда он работает в последовательном режиме .... 2) Я прав на динамическое выделение памяти для каждого потока, такого как «normalized_p_gn» выше?

С уважением Эдвин

Answer 1

Одна из возможностей, которую мы не можем исключить, ваш фрагмент кода не выявил каких-либо цифр, является то, что вы просто пытаетесь выделить слишком много памяти при запуске этого параллельно. Если вы можете подтвердить, что это не ваша проблема, прокомментируйте или отредактируйте свой вопрос, и я еще раз посмотрю.

Answer 2

Вместо malloc используйте файлы mxMalloc в файлах mex (см. here). Не забудьте сделать mxFree, когда закончите с памятью.