¿Por qué OpenMP en un archivo mex solo produce 1 hilo?

Soy nuevo en OpenMP. Tengo el siguiente código que compila bien usando Matlab mex configurado con MSVS2010. La computadora tiene 8 procesadores disponibles (que revisé también usando matlabpool).¿Por qué OpenMP en un archivo mex solo produce 1 hilo?

#include "mex.h" 
#include <omp.h> 

typedef unsigned char uchar; 
typedef unsigned int uint; 
//Takes a uint8 input array and uint32 index array and preallocated uint8 array the same 
//size as the first one and copies the data over using the indexed mapping 
void mexFunction(int nlhs, mxArray *plhs[], int nrhs, const mxArray*prhs[]) 
{ 
    uint N = mxGetN(prhs[0]); 
    mexPrintf("n=%i\n", N); mexEvalString("drawnow"); 
    uchar *input = (uchar*)mxGetData(prhs[0]); 
    uint *index = (uint*)mxGetData(prhs[1]); 
    uchar *output = (uchar*)mxGetData(prhs[2]); 

    uint nThreads, tid; 
#pragma omp parallel private(tid) shared(input, index, output, N, nThreads) num_threads(8) 
    { 
     tid = omp_get_thread_num(); 

     if (tid==0) { 
      nThreads = omp_get_num_threads(); 

     } 

     for (int i=tid*N/nThreads;i<tid*N/nThreads+N/nThreads;i++){ 
      output[i]=input[index[i]]; 
     } 
    } 
    mexPrintf("nThreads = %i\n",nThreads);mexEvalString("drawnow"); 
}

La salida que se ve es

n=600000000 
nThreads = 1

¿Por qué es sólo un hilo está creando a pesar que solicita 8?

Fuente

2011-11-07 twerdster

suspiro. Típico, pasar horas intentando y fallando y luego encontrar la respuesta 5 minutos después de publicar en SO.

El archivo debe ser mexed con soporte OpenMP

mex mexIndexedCopy.cpp COMPFLAGS="/openmp $COMPFLAGS"

Fuente

2011-11-07 01:16:04 twerdster

Conozco esa sensación, hermano. – CptSupermrkt

¿Cuál es la opción equivalente en Linux con gcc como compilador? '-fopenmp'? – linello

@linello Sí. De hecho, solo perdí un par de horas porque no estaba pasando '-fopenmp' correctamente. Debe pasarlo tanto al compilador como al enlazador. 'mex CXXFLAGS =" \ $ CXXFLAGS -fopenmp "LDFLAGS =" \ $ LDFLAGS -fopenmp "[otras opciones] ' para C++. (Para C, use 'CFLAGS' en lugar de' CXXFLAGS'; para C y C++, use ambos). –

¿Por qué OpenMP en un archivo mex solo produce 1 hilo?

Respuesta

Cuestiones relacionadas