MATLAB Answers

Mex OpenMP C code

23 views (last 30 days)
Yangyang Xu
Yangyang Xu on 22 Jan 2019
Edited: sjhstone on 19 Aug 2020
Hi,
I tried to use OpenMP in Matlab 2018a. My code is very simple, called testOpenMP.c.
I successfully compiled the code by mex -v testOpenMP.c CFLAGS="$CFLAGS -fopenmp" LDFLAGS="$LDFLAGS -fopenmp"
However, we I run the compile mex extension by just typing testOpenMP in Matlab's command window. It always crashed.
Do you know how to correctly compile such an openmp c code and run it?
Thanks,
#include "mex.h"
#include "math.h"
#include "stdio.h"
#include <omp.h>
void mexFunction(
int nargout,
mxArray *pargout [ ],
int nargin,
const mxArray *pargin [ ]
) {
int thread_count = 32;
#pragma omp parallel for default(none)
for(int i = 0; i < 32; i++)
{
printf("Hello from %d\n",i);
}
}
  5 Comments
Sushant Mahajan
Sushant Mahajan on 4 May 2019
Did you find any solution? I am facing the exact same issue with a FORTRAN subroutine using MEX configured with a gfortran-6 compiler.
It runs fine when it is compiled without the -fopenmp flag, but MATLAB shuts down when the subroutine is called after using the -fopenmp flag during compilation. This happens even if I do not use any OMP_PARALLEL (equivalent to C's #pragma omp parallel) construct in the code.

Sign in to comment.

Answers (2)


sjhstone
sjhstone on 19 Aug 2020
Edited: sjhstone on 19 Aug 2020
Both printf and mexPrintf are not thread-safe and should be avoided in multi-thread OpenMP MEX C codes.
You can get more information at https://walkingrandomly.com/?p=1795

Tags

Products

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!