comp.lang.idl-pvwave archive: archive » Some MAKE

Home » Public Forums » archive » Some MAKE_DLL questions

Show: Today's Messages :: Show Polls :: Message Navigator
E-mail to friend

Re: Some MAKE_DLL questions [message #92148 is a reply to message #92143]

Tue, 20 October 2015 08:13

natha
Messages: 482
Registered: October 2007

Senior Member

OK, so it seems to work...
I took the example code of Michael's book and I wrote the following OpenMP version:

#include <omp.h>
#include <stdio.h>
#include <stdlib.h>

float callex_total(float arr[], int *n) {

int nthreads, i, tid;
float total = 0.;

#pragma omp parallel shared(arr) private(i,tid)
{
tid = omp_get_thread_num();
if (tid == 0)
{
nthreads = omp_get_num_threads();
printf("Number of threads = %d\n", nthreads);
}

printf("Thread %d starting...\n",tid);

#pragma omp for
for (i = 0; i < *n; i++)
{
total += arr[i];
printf("Thread %d: total += arr[%d]\n",tid,i);
}

} /* end of parallel section */

return(total);
}

Then, I had to play a lot with the compiler but at the end I've found the way:

cfile='callex_total.c'
cfile_noext=file_basename(cfile, '.c')

srcdir = file_dirname(file_expand_path(cfile))

cc='gcc -fopenmp -lgomp -fPIC %c -c -o %o'
ld='gcc -fopenmp -shared -o %L %O %X'

make_dll, cfile_noext, 'IDL_Load', input_directory=srcdir, output_directory=srcdir, /verbose, /show_all_output, cc=cc, ld=ld ;;extra_cflag='-fopenmp -c'

result = call_external(cfile_noext+'.so', cfile_noext, findgen(10), 10, /f_value, /auto_glue)
print, result

And here is the output:
Thread 6 starting...
Thread 4 starting...
Thread 4: total += arr[8]
Thread 4: total += arr[9]
Thread 2 starting...
Thread 2: total += arr[4]
Thread 2: total += arr[5]
Thread 7 starting...
Thread 5 starting...
Thread 1 starting...
Thread 1: total += arr[2]
Thread 1: total += arr[3]
Thread 3 starting...
Thread 3: total += arr[6]
Thread 3: total += arr[7]
Number of threads = 8
Thread 0 starting...
Thread 0: total += arr[0]
Thread 0: total += arr[1]
45.0000

I am very glad to know that we can link OpenMP to IDL.

Report message to a moderator

[Message index]

		Some MAKE_DLL questions By: natha on Mon, 19 October 2015 06:30
		Re: Some MAKE_DLL questions By: Russell[1] on Mon, 19 October 2015 08:04
		Re: Some MAKE_DLL questions By: natha on Mon, 19 October 2015 08:44
		Re: Some MAKE_DLL questions By: Michael Galloy on Mon, 19 October 2015 09:58
		Re: Some MAKE_DLL questions By: natha on Mon, 19 October 2015 10:36
		Re: Some MAKE_DLL questions By: Jim Pendleton on Mon, 19 October 2015 18:40
		Re: Some MAKE_DLL questions By: natha on Tue, 20 October 2015 08:13

Previous Topic:	Read ENVI file
Next Topic:	"% Unable to allocate memory" when using histogram.

-=] Back to Top [=-

[ Syndicate this forum (XML) ] [

] [

]

Current Time: Fri Oct 10 06:11:57 PDT 2025

Total time taken to generate the page: 1.11922 seconds