DOI QR코드

DOI QR Code

Two dimensional variable-length vector storage format for efficient storage of sparse matrix in the finite element method

유한요소법에서 희소행렬의 효율적인 저장을 위한 2차원 가변길이 벡터 저장구조

  • Boo, Hee-Hyung (School of Electrical Engineering and Computer Science, Kyungpook National University) ;
  • Kim, Sung-Ho (School of Computer Science and Engineering, Kyungpook National University)
  • 부희형 (경북대학교 전자전기컴퓨터학부) ;
  • 김승호 (경북대학교 컴퓨터학부)
  • Received : 2012.02.23
  • Accepted : 2012.09.07
  • Published : 2012.09.30

Abstract

In this paper, we propose the two dimensional variable-length vector storage format which can be used for efficient storage of sparse matrix in the FEM (finite element method). The proposed storage format is the method storing only actual needed non-zero values of each row on upper triangular matrix with the total rows N, by using two dimensional variable-length vector instead of $N{\times}N$ large sparse matrix of entire equation of finite elements. This method only needs storage spaces of the number of minimum 1 to maximum 5 in 2D grid structure and the number of minimum 1 to maximum 14 in 3D grid structure of analysis target. The number doesn't excess two times although involving index number. From the experimental result, we can find out that the proposed storage format can reduce the memory space more effectively, as the total number of nodes increases, than the existing skyline storage format storing maximum column height.

본 논문에서는 유한요소법에서 희소행렬의 효율적인 저장을 위한 2차원 가변길이 벡터 저장구조를 제안한다. 제안한 저장구조는 유한요소 전체 방정식의 거대희소행렬 $N{\times}N$ 대신, 전체 행의 개수 N의 상삼각행렬에서 0이 아닌 실제 필요한 값들만 2차원 가변길이 벡터를 이용하여 저장하는 방법이다. 이 방법을 이용하면, 해석대상의 2차원 격자구조에서는 각 절점당 최소 1개에서 최대 5개까지의 저장 공간이 필요하게 되고, 3차원 격자구조에서는 각 절점당 최소 1개에서 최대 14개까지의 저장 공간이 필요하게 된다. 인덱스를 포함해도 2배 이상을 넘지 않는다. 본 논문의 실험 결과에 의해, 제안한 저장구조는 총 절점 개수가 많아질수록 기존의 최대칼럼 높이를 저장하는 스카이 라인 저장구조보다 메모리 공간을 효과적으로 줄일 수 있는 구조임을 알 수 있었다.

Keywords

References

  1. G. R. LIU, S. S. QUEK, "THE FINITE ELEMENT METHOD: A practical course," Butterworth- Heinemann, 2003.
  2. R.E. Lewis, J.P. Ward, "The finite element method: principles and applications," Addison-Wesley, 1991.
  3. Takeo Taniguchi, Kohji Fujiwara, "Parallel Skyline Method using Two Dimensional Array," Memoirs of the Faculty of Engineering, Okayama University, VoI.24, No.2, pp.99-1I2, March 1990.
  4. Felippa, C.A., "Solution of Linear Equations with Skyline-stored Symmetric Matrix," Computers and Structures, Vol. 5, pp. 13-29, 1975. https://doi.org/10.1016/0045-7949(75)90016-4
  5. Eun-Jin Im, "An Efficient Computation of Matrix Triple Products," Journal of the Korea Society of Computer and Information, Vol.11, No.3, pp. 141-149, July 2006.
  6. Thomas J. Rudolphi, "Finite Element Method," McGraw-Hill Korea, 2010.
  7. NVIDIA CUDA (Compute Unified Device Architecture ): Programming Guide, Version 2.1, December 2008.
  8. Junghwan Kim, Jinsoo Kim, "Implementation of Efficient Power Method on CUDA GPU ," Journal of the Korea Society of Computer and Information, Vol.16, No.2, pp. 9-16, Feb. 2011. https://doi.org/10.9708/jksci.2011.16.2.009
  9. Zhihui Zhang, Qinghai Miao, Ying Wang, "CUDABased Jacobi's Iterative Method," Computer Science-Technology and Applications, IFCSTA '09. International Forum, IEEE Computer Society, Vol. 1, pp. 259-262, 2009.