윈디하나의 솔라나라: OpenMP

• 분류: 솔라리스
• 문서정보: OpenMP의 소개 및 솔라리스에서 OpenMP를 사용하는 방법에 대해 설명
• 업데이트: 2009-03-27
• 이 문서는 윈디하나의 솔라나라, http://www.solanara.net/에서 최근에 갱신된 문서를 찾을 수 있다.

목차

• OpenMP 소개
• OpenMP 예제
  • OpenMP 지원 확인
  • OpenMP 환경변수
  • OpenMP를 이용한 행렬 계산
  • OpenMP를 이용한 QuickSort

OpenMP 소개

• CPU의 코어수가 늘어남에 따라 쉽고 빠른 쓰레드 프로그래밍이 각광받고 있는 이때, OpenMP가 나왔다. 현재 표준이 되었기 때문에 사용하는데 부담은 없을듯 하다.
• 솔라리스 8 이상, Sun Studio 10 이상에서 작동한다. (테스트한 시스템은 솔라리스 10, Sun Studio 12이다) GCC 3.4이상에서도 작동하고, 타 플랫폼(예를 들면 윈도우즈 + VC환경)에서도 작동한다.
• 이 문서에서는 CPU의 코어가 2개인 경우를 가정한다. (HyperThread가 지원되어서 2개로 나오는 경우는 안된다. 코어가 2개 이상 있어야 한다)
• OpenMP 공식 홈페이지: http://openmp.org/
• 참고
  • Sun Studio 12 OpenMP API User's Guide
  • OpenMP Tutorial
  • OpenMP Application Program Interface Version 3.0

OpenMP 예제

OpenMP 지원 확인

/*
	Print OpenMP status.
	WindyHana's Solanara: OpenMP http://www.solanara.net/solanara/openmp
	cc -O3 -xopenmp=parallel -o openmpstatus openmp_status.c
*/
#include <stdio.h>
#include <omp.h>

int main(int argc, char *argv[]) {
	int mt = 1, iam = 0, np = 1, omp = 0;

	#pragma omp parallel default(shared) private(mt, iam, np, omp)
	{
		#if defined (_OPENMP)
			mt = omp_get_max_threads();
			np = omp_get_num_threads();
			iam = omp_get_thread_num();
			omp = 1;
		#endif
		if (omp) {
			printf("MaxThreads: %d, UsedThreads: %d, ThreadNum: %d\n", mt, np, iam);
		} else {
			printf("No OpenMP used.\n");
		}
	}
	return 0;
}

root@wl ~ # export OMP_NUM_THREADS=2
root@wl ~ # cc -O3 -xopenmp=parallel -o openmpstatus openmp_status.c
root@wl ~ # ./openmpstatus
MaxThreads: 2, UsedThreads: 2, ThreadNum: 0
MaxThreads: 2, UsedThreads: 2, ThreadNum: 1
root@wl ~ # 

OpenMP 환경변수

• OMP_NUM_THREADS: OpenMP에서 사용할 쓰래드의 개수. 기본값은 1이며, 보통 코어의 개수를 지정해준다.
• OMP_DYNAMIC: FALSE로 설정하면 시스템의 코어개수에 상관 없이 OMP_NUM_THREADS에 설정된 개수를 사용한다.

root@wl ~ # export OMP_NUM_THREADS=4
root@wl ~ # export OMP_DYNAMIC=FALSE
root@wl ~ # ./matrix

• SUNW_MP_PROCBIND: 각각의 쓰레드를 어느 가상 프로세서에 바인딩 할 것인지 결정한다. 번호는 가상 프로세서의 순서이다. [psrinfo | grep "on-line"]을 이용했을때 나오는 순서를 적어준다. (프로세서 아이디를 말하는게 아니다)
  

  root@wl ~ $ /usr/sbin/psrinfo | grep "on-line"
  0       on-line   since 10/07/09 21:00:00
  1       on-line   since 10/07/09 21:00:00
  10      on-line   since 10/07/09 21:00:00
  11      on-line   since 10/07/09 21:00:00
  20      on-line   since 10/07/09 21:00:00
  21      on-line   since 10/07/09 21:00:00
  root@wl ~ $ export SUNW_MP_PROCBIND=false
  root@wl ~ $ export SUNW_MP_PROCBIND=0
  root@wl ~ $ export SUNW_MP_PROCBIND="0 2 3 4 5"
  root@wl ~ $ export SUNW_MP_PROCBIND="0-6"
  

• STACKSIZE: 각각의 쓰레드의 스택 크기를 지정한다. 단위는 킬로바이트이다.
  예) 스택사이즈를 16MiB로 지정: export STACKSIZE=16384

OpenMP를 이용한 행렬 계산

/*
	Calc simpel matrix
	WindyHana's Solanara: OpenMP http://www.solanara.net/solanara/openmp
	cc -O3 -mt -lmtmalloc -xopenmp=parallel -o openmpmatrix openmp_matrix.c
*/
#include <stdio.h>
#include <omp.h>
#include <mtmalloc.h>
const int nrows =      10;
const int ncols = 1000000;

int main(int argc, char *argv[]) {
	int* array[nrows];

	printf("Initializing...\n");
	#pragma omp parallel for
		for (int i = 0; i < nrows; i++) {
			array[i] = (int*) malloc(sizeof(int) * ncols);
		}

	#pragma omp parallel for
		for (int i = 0; i < nrows; i++) {
			printf("Calc %d column by ThreadNum %d\n", i, omp_get_thread_num());
			for (int j = 0; j < ncols; j++) {
				for (int k = 0; k < nrows; k++) {
					array[i][j] = array[i][k] * array[k][j];
				}
			}
		}
	return 0
}

root@wl ~ # export OMP_NUM_THREADS=2
root@wl ~ # cc -O3 -mt -lmtmalloc -xopenmp=parallel -o openmpmatrix openmp_matrix.c
root@wl ~ # ./openmpmatrix
Initializing...
Calc 0 column by ThreadNum 0
Calc 5 column by ThreadNum 1
Calc 6 column by ThreadNum 1
Calc 1 column by ThreadNum 0
Calc 2 column by ThreadNum 0
Calc 7 column by ThreadNum 1
Calc 3 column by ThreadNum 0
Calc 8 column by ThreadNum 1
Calc 4 column by ThreadNum 0
Calc 9 column by ThreadNum 1
root@wl ~ #

OpenMP를 이용한 QuickSort

/*
	QuickSort for Solaris, OpenMP
	WindyHana's Solanara: OpenMP http://www.solanara.net/solanara/openmp
	cc -O3 -xopenmp=parallel -o openmpquick openmp_quick.c
*/
#include <stdio.h>
#include <stdlib.h> 
#include <string.h>
#include <omp.h>
#include <sys/time.h>

#define MAX_ELEMENTS 500000
#define MEMALLOCSIZE MAX_ELEMENTS * sizeof(int)

void quicksort(unsigned int * data, int lo, int hi);
void quicksort_openmp(unsigned int * data, int lo, int hi);
void printdata(unsigned int * data, int size);
unsigned long long gettimestamp();

int main(int argc, char *argv[]){
	unsigned int * data;
	unsigned int * data2;
	unsigned long long e1, e2, t1, t2;
	unsigned int seed = (unsigned int)gettimestamp();
	printf("Seed is %d\n", seed);
	srand(seed);
	
	printf("Alloc Memory: %d * 2 bytes\n", MEMALLOCSIZE);
	data = (unsigned int *) malloc(MEMALLOCSIZE);
	data2 = (unsigned int *) malloc(MEMALLOCSIZE);

	printf("Initializing Data [%d]\n", MAX_ELEMENTS);
	for (int i = 0; i < MAX_ELEMENTS; i++) { 
		int r1 = rand() & 0x0000FFFF;
		int r2 = rand() & 0x0000FFFF;
		unsigned int t = (r1 << 16) | r2;
		data[i] = t;
	}
	memcpy((void *)data2, (void *)data, MEMALLOCSIZE);

    #if defined (_OPENMP)
		printf("Max OpenMP Threads: %d\n", omp_get_max_threads());
    #endif

	t1 = gettimestamp();
	quicksort(data, 0, MAX_ELEMENTS - 1);
	t2 = gettimestamp();
	printf("QuickSort\n");
	printdata(data, MAX_ELEMENTS);
	e1 = t2 - t1;

	t1 = gettimestamp();
	quicksort_openmp(data2, 0, MAX_ELEMENTS - 1);
	t2 = gettimestamp();
	printf("QuickSort(OpenMP)\n");
	printdata(data2, MAX_ELEMENTS);
	e2 = t2 - t1;

	printf("       Elapsed: %lld\n", e1);
	printf("OpenMP Elapsed: %lld\n", e2);
	free(data);
	free(data2);
	return 0;
}

void quicksort_openmp(unsigned int * data, int lo, int hi) {
	int i = lo, j = hi, temp;
	int x = data[(lo + hi) / 2];
	do { 
		while (data[i] < x) i++; 
		while (data[j] > x) j--;
		if (i <= j) {
			temp = data[i]; data[i] = data[j]; data[j] = temp; // SWAP
			i++; j--;
		}
	} while (i <= j);
	#pragma omp parallel sections
	{
		#pragma omp section
		if (lo < j) quicksort(data, lo, j);
		#pragma omp section
		if (i < hi) quicksort(data, i, hi);
	}
}

// Just remove openmp progma from quicksort_openmp()
void quicksort(unsigned int * data, int lo, int hi) {
	int i = lo, j = hi, temp;
	int x = data[(lo + hi) / 2];
	do { 
		while (data[i] < x) i++; 
		while (data[j] > x) j--;
		if (i <= j) {
			temp = data[i]; data[i] = data[j]; data[j] = temp; // SWAP
			i++; j--;
		}
	} while (i <= j);
	if (lo < j) quicksort(data, lo, j);
	if (i < hi) quicksort(data, i, hi);
}

// print integer array (for debug)
void printdata(unsigned int * data, int size) {
	if (size > 100) {
		for (int i = 0; i < 2; i ++) {
			printf("data[%d] = %d\n", i, data[i]);
		}
		printf("...\n");
		for (int i = size - 3; i < size; i ++) {
			printf("data[%d] = %d\n", i, data[i]);
		}
	} else {
		for (int i = 0; i < size; i ++) {
			printf("data[%d] = %d\n", i, data[i]);
		}
	}
}

// get current timestamp (micro seconds)
unsigned long long gettimestamp() {
	struct timeval t;
	gettimeofday(&t, NULL);
	return t.tv_sec * 1000000 + t.tv_usec;
}

root@wl ~ # export OMP_NUM_THREADS=2
root@wl ~ # cc -O3 -xopenmp=parallel -o openmpquick openmp_quick.c
root@wl ~ # ./openmpquick
Seed is 1389706669
Alloc Memory: 2000000 * 2 bytes
Initializing Data [500000]
Max OpenMP Threads: 2
QuickSort
data[0] = 636
data[1] = 2446
...
data[499997] = 2147448088
data[499998] = 2147448460
data[499999] = 2147449545
QuickSort(OpenMP)
data[0] = 636
data[1] = 2446
...
data[499997] = 2147448088
data[499998] = 2147448460
data[499999] = 2147449545
       Elapsed: 84537
OpenMP Elapsed: 77274
root@wl ~ #

Copyright © 2004-2010 Jo HoSeok. All rights reserved.