SCAM 기상모델의 성능향상을 위한 LAPACK BLAS 라이브러리의 활용

Performance Improvements of SCAM Climate Model using LAPACK BLAS Library

Abstract

슈퍼 컴퓨팅 기술 및 하드웨어 기술의 발달로 수치 연산 방식 또한 고도화되고 있다. 그에 따라 이전 대비 향상된 기상 예측 또한 가능해진다. 본 논문에서는 SCAM(Single-Columns Atmospheric Model, CESM(Community Earth System Model)을 간소화 한 버전)에 포함되어 있으며 대기 연산을 수행하는 적운 모수화 코드, Unicon(A Unified Convection Scheme)의 성능을 향상하기 위하여 소스 코드 내의 선형대수 수치적 연산 부분에 고밀도 선형대수 연산을 위한 라이브러리인 LAPACK(Linear Algebra PACKage) BLAS(Basic Linear Algebra Subprograms)의 level1 함수를 적용할 것을 제안한다. 이를 분석하기 위하여 SCAM의 전체적인 실행 구조도를 제시하고 해당 실행환경에서 테스트를 진행하였다. 기존 소스 코드 대비 SCOPY 함수는 0.4053%, DSCAL 함수는 0.7812%, DDOT 함수는 0.0469%의 성능 향상을 이끌어 내었으며 이를 모두 적용한 결과 기존 소스 코드 대비 0.8537%의 성능 향상을 보였다. 이는 본 논문에서 제안한 고밀도 선형대수 연산을 위한 라이브러리인 LAPACK BLAS 적용 방법이 동일한 CPU 환경에서 추가적인 하드웨어의 개입 없이 성능을 향상시킬 수 있음을 의미한다.

With the development of supercomputing technology and hardware technology, numerical computation methods are also being advanced. Accordingly, improved weather prediction becomes possible. In this paper, we propose to apply the LAPACK(Linear Algebra PACKage) BLAS(Basic Linear Algebra Subprograms) library to the linear algebraic numerical computation part within the source code to improve the performance of the cumulative parametric code, Unicon(A Unified Convection Scheme), which is included in SCAM(Single-Columns Atmospheric Model, simplified version of CESM(Community Earth System Model)) and performs standby operations. In order to analyze this, an overall execution structure diagram of SCAM was presented and a test was conducted in the relevant execution environment. Compared to the existing source code, the SCOPY function achieved 0.4053% performance improvement, the DSCAL function 0.7812%, and the DDOT function 0.0469%, and all of them showed a 0.8537% performance improvement. This means that the LAPACK BLAS application method, a library for high-density linear algebra operations proposed in this paper, can improve performance without additional hardware intervention in the same CPU environment.

1. 서론

슈퍼 컴퓨팅 기술의 발전에 따라 각국의 기후 예측모델 또한 고도화되고 있다. 한국 기상청(KMA, Korea Meteorological Administration) 또한 기후모델(Climate Mode)을 도입하여 기상 예보를 수행하고 있다. 한국 기상청은 2004년 미국 국립대기 과학연구소(NCAR, National Center for Atmospheric Research)와 기술 협력을 체결하였고 2014년에는 전지구 기상모델(GCM, Global Climate Model)인 GloSea5(Global Seasonal Forecast System)를 도입하여 슈퍼컴퓨터 5호기에 이식하였다. 또한 해당 모델을 GloSea6로 버전 업그레이드를 진행한 바 있다.

CESM(Community Earth System Model)은 주로 NSF(National Science Foundation)에서 지원하는 대기, 해양, 얼음, 지표면, 탄소 순환 및 기타 구성요소로 구성된 지구의 과거, 현재 및 미래 기후 상태에 대한 최첨단 컴퓨터 시뮬레이션을 제공하는 지구 기후모델이다. 이는 그림 1과 같이 지구 모델을 일정한 수평, 수직의 격자 크기로 나눈 후 각 격자에 대한 대기의 상태와 운동에 대한 방정식을 계산하여 얻은 결과를 통해 시뮬레이션 한다[1].

그림 1. CESM, SCAM의 격자 모델링 형태

Fig. 1. Grid modeling of CESM, SCAM

SCAM(Single-Columns Atmospheric Model)[1]은 특정 기상 연구를 위해 CESM을 간소화 한 버전이다. 전 지구가 아닌 단일 격자(Single-Columns)에 대하여 CESM이 수행할 수 있는 4가지 요소 중 오직 대기(Atmospheric) 모델에만 집중하여 모델을 수행한다. SCAM은 2088번의 step (1step=1,200초)을 수행할 수 있으며 이를 Day(일)로 환산할 시 최대 29일까지 모델을 수행할 수 있다[2].

LAPACK(Linear Algebra PACKage)은 선형대 수학의 수치적 계산을 수행하는 라이브러리로, 일차 연립방정식, 행렬 분해, 선형 최소제곱법 등을 수행하는 함수를 제공한다. LAPACK은 1992년에 Fortran77로 작성되었으며 2008년에 Fortran90언어로 변경되었다[3]. BLAS(Basic Linear Algebra Subroutines)는 기본 선형 대수 하위 프로그램에서 벡터 덧셈, 스칼라 곱셈, 내적, 선형 조합 및 행렬 곱셈과 같은 일반적인 선형 대수 연산을 수행하는 선형 대수 라이브러리를 위한 표준 저 수준 루틴이다[4]. 그림 2와 같이 LAPACK은 BLAS 라이브러리의 호출을 통해 가능한 한 많은 계산이 수행되도록 작성되었다[5]. 따라서 고성능 컴퓨터에서 LAPACK BLAS를 사용할 경우 성능 효율성 증가를 기대할 수 있다[6].

그림 2. LAPACK의 개요

Fig. 2. Overview of LAPACK

SCAM 모델에서 적운 모수화 연산을 담당하는 코드는 Unicon(A Unified Convection Scheme)[7]이며 해당 코드에 LAPACK BLAS 라이브러리를 적용함으로써 Unicon 코드를 포함한 SCAM 모델 전체의 성능 향상을 끌어내어 관련 기술을 축적하고 이를 한국형 기상 모델(KIM)에 적용하여 한국의 계절 예측 시스템의 고도화 방안을 연구하고자 한다.

2. 관련 연구

Merchant, et al는 BLAS 및 LAPACK이 여러 고성능 컴퓨팅(HPC) 애플리케이션의 기본 구성 요소를 형성하므로 성능 결정과 BLAS/LAPACK 루틴의 성능 튜닝을 위한 FPU(Floating Point Unit) 마이크로 아키텍처의 영향을 조사하였다[8]. Kim, et a l.는 행렬의 대규모 그룹에 적용된 BLAS/LAPACK 루틴에서 SIMD 벡터 길이에 따라 행렬을 블록으로 인터리브하는 새로운 소형 데이터 레이아웃을 제안하고 이를 계층적 병렬 애플리케이션 내에서 사용할 수 있는 새로운 인터페이스와 결합하여 아키텍처의 속도를 향상시켰다[9]. Phillips는 공유 메모리 멀티 프로세서를 기반으로 각 처리 요소가 스칼라 특성만 갖는 기본 선형 대수 하위 프로그램(BL AS)의 성능을 고려하여 동일한 환경에서 실행되었을 때 특정 LAPACK 루틴의 성능에 대한 블록 크기의 영향을 조사하였다[10]. Bjørstad는 BLAS 루틴이 컴파일러 제한에 의해 고성능이 억제된다는 것을 조사하였으며 그에 따라 최적의 데이터 이동을 위해 단일 블록에 적용되는 명시적 LU 인수 분해의 필요성을 제기하고 블록 크기의 확장을 통해 LAPACK 루틴의 블록 구조를 수정하지 않고도 고성능에 달성할 수 있음을 보였다[11].

3. SCAM 실행 구조

SCAM 모델을 수행한 방법은 그림 3과 같이 Local PC Client에서 Gate로 SSH 원격 접속을 이용하여 실제로 모델을 수행하는 Cluster Server에 연산을 요청하는 형태로 진행하였다. 해당 방식으로 사용자가 SCAM 모델을 수행시키게 되면 Gate가 USER 계정의 SCAM 소스 코드를 Cluster Server로 전달하게 된다. 전달받은 소스 코드는 Batch System(pbs)에 SCAM Running Job으로 전달되어 Job scheduler(torque)를 거쳐 실행된다. 모델을 실행할 때마다 컴파일은 새로 진행되며 모델 수행이 종료되면 최종 출력값인 nc 확장자 파일 4개와 컴파일 로그 및 러닝 로그가 생성되고 이를 rynch를 이용하여 User가 전달받을 수 있다[12].

그림 3. SCAM 실행 구조도

Fig. 3. Running Process of SCAM System

4. 성능 향상을 위한 LAPACK BLAS 적용 방법

4.1 LAPACK BLAS

LAPACK BLAS는 고밀도 선형대수 연산을 위해 만들어진 고성능 라이브러리이다. 그중 Level1 함수는 벡터 – 벡터의 일반화된 연산 y = αx + y (여기서 x, y는 벡터, α는 스칼라를 의미함)을 수행한다. 따라서 본 논문에서는 대기 모델의 연산에 사용되는 방정식 풀이를 LAPACK BLAS의 Level1 함수를 통해 연산하도록 하여 성능을 향상시킬 것을 제시한다.

BLAS function의 명명 규칙은 그림 4와 같이 첫 번째 글자는 Number type을 의미하며 각 함수의 기능에 따라 명명된다[13]. 본 논문에서는 Unicon 소스코드 연산에 사용되는 변수들의 Data type을 고려하여 SCOPY, DSCAL, DDOT의 함수를 적용하였으며 3가지의 함수가 모두 적용되어 연산을 수행하도록 하였다.

그림 4. BLAS routines 명명 규칙

Fig. 4. BLAS routines: naming conventions

4.2 SCOPY 함수 적용

SCOPY 함수는 벡터x를 벡터y로 unrolled loops를 사용하여 복사하는 함수이다. SCOPY 함수에 사용되는 각 파라미터가 의미하는 바는 아래 표 1과 같다[14]. Unicon.F90에서 행렬을 복사하는 수식이 loop 안에 존재할 경우 수동으로 loop 문을 분리하고 SCOPY 함수를 call 하여 연산하도록 하였다. 그림 5는 SCOPY가 적용된 Unicon.F90 소스코드의 일부이다.

표 1. Parameters of SCOPY

Table 1. Parameters of SCOPY

그림 5. SCOPY 함수 적용

Fig. 5. Apply of SCOPY function

4.3 DSCAL 함수 적용

DSCAL 함수는 벡터 x를 스칼라 α만큼의 곱연산을 unrolled loops를 사용하여 결괏값을 연산하고 벡터 x를 반환하는 함수이다. DSCAL에 사용되는 파라미터가 의미하는 바는 표 2와 같다[15]. 해당 함수 또한 마찬가지로 unrolled loops를 사용하여 연산하므로 적용하려는 연산식이 loop문 안에 존재할 경우 수동으로 loop문에서 분리한 후 function call하여야 한다. 그림 6은 DSCAL이 적용된 Unicon.F90 소스코드의 일부이다.

표 2. Parameters of DSCAL

Table 2. Parameters of DSCAL​​​​​​​

그림 6. DSCAL 함수 적용

Fig. 6. Apply of DSCAL function​​​​​​​

4.4 DDOT 함수 적용

DDOT 함수는 벡터x 및 벡터y의 내적, 즉 x ∙ y = x₁y₁ + x₂y₂ + … + x_ny_n 를 연산하는 함수이다. DDOT에 사용되는 파라미터가 의미하는 바는 표 3과 같다[16]. 앞선 두 함수의 경우 parameter에 out값 또한 선언되어 있었으나 DDOT 함수는 return할 값의 타입이 parameter에 정해져있지 않아 그림 7과 같이 fortran에서 제공하는 External이라는 사용자 정의 함수를 사용하여 수동으로 정의한 후 적용하였다. 그림 8은 DDOT이 적용된 Unicon.F90 소스코드의 일부이다.

그림 7. External DDOT 선언

Fig. 7. External DDOT declaration

표 3. Parameters of DDOT

Table 3. Parameters of DDOT​​​​​​​

그림 8. DDOT 함수 적용

Fig. 8. Apply of DDOT function

5. 실험 결과 및 성능 평가

성능 분석을 위하여 3장에서 제시한 SCAM 실행방법에 따라 모델 수행을 30회 실행하고 그에 대한 평균 시간을 기준으로 하였다. SCAM의 step은 2088, 즉 29일로 수행하였다. 또한 1회당 수행시간은 컴파일시간(compile time)과 모델 수행시간(running time)을 합한 것으로 한다. 함수를 개별적으로 적용한 경우와 모두 적용한 경우의 평균값을 각각 도출하였다. 실험에 수행된 시스템 스펙은 표4와 같다.

표 4. 실험에 수행된 시스템 스펙

Table 4. System specs used for SCAM​​​​​​​

5.1 SCOPY 성능평가

적운 모수화 코드의 기존 연산 대신 SCOPY 함수를 적용하여 연산한 결과는 그림 9와 같다. 기존 소스코드는 평균 5.1069초 , SCOPY가 적용된 코드는 5.0862초의 성능을 보이며 이는 기존 대비 0.4053% 성능 향상을 보인다.

그림 9. SCOPY 성능평가

Fig. 9. Evaluation of SCOPY​​​​​​​

5.2 DSCAL 성능평가

적운 모수화 코드의 기존 연산 대신 DSCAL 함수를 적용하여 연산한 결과는 그림 10과 같다. 기존 소스코드는 평균 5.1069초 , DSCAL이 적용된 코드는 5.067초의 성능을 보이며 이는 기존 대비 0.7812% 성능향상을 보인다.

그림 10. DSCAL 성능평가

Fig. 10. Evaluation of DSCAL​​​​​​​

5.3 DDOT 성능평가

적운 모수화 코드의 기존 연산 대신 DDOT함수를 적용하여 연산한 결과는 그림 11과 같다. 기존 소스코드는 평균 5.1069초 , DDOT이 적용된 코드는 5.1045초의 성능을 보이며 이는 기존 대비 0.0469% 성능향상을 보인다.

그림 11. DDOT 성능평가

Fig. 11. Evaluation of DDOT​​​​​​​

5.4 SCOPY+DSCAL+DDOT 성능평가

적운 모수화 코드의 기존 연산 대신 앞서 제시한 세 함수를 모두 적용하여 연산한 결과는 그림 12와 같다. 기존 소스코드는 평균 5.1069초 , 세 함수가 모두 적용된 코드는 5.0633초의 성능을 보이며 이는 기존 대비 0.8537% 성능이 향상되었음을 보인다.

그림 12. DSCAL+SCOPY+DDOT 성능평가

Fig. 12. Evaluation of DSCAL+SCOPY+DDOT​​​​​​​

6. 결론

본 논문에서는 SCAM에 포함되어 대기 연산을 담당하고 있는 적운 모수화 코드 Unicon의 성능 향상을 위하여 LAPACK BLAS 라이브러리 중 level1의 function 세 가지를 적용하는 방안을 제시하였다. 그리고 이를 구축된 환경에서 SCAM step 2088 (즉, 29일)을 기준으로 실행 및 분석한 결과 SCOPY 함수 적용은 0.4053% 성능 향상을, DSCAL 함수 적용은 0.7812%의 성능 향상을, DDOT 함수 적용은 0.0469%의 성능 향상을 보여주었으며 모든 함수를 종합하여 적용한 결과 0.8537%의 성능 향상을 보였다. 이는 동일한 CPU 환경에서 추가적인 하드웨어의 개입 없이 고밀도 선형대수 연산을 위한 고성능 라이브러리인 BLAS를 사용하여 성능을 향상시킬 수 있음을 의미한다. 또한, 본 논문에서 제안한 LAPACK BLAS 라이브러리 적용 방법이 SCAM의 성능 향상을 이끌어내었음을 보여준다.