[논문 리뷰] smt: a Matlab structured matrices toolbox

Michela Redivo‐Zaglia, G. Rodríguez|arXiv (Cornell University)|2009. 03. 30.

Matrix Theory and Algorithms참고 문헌 13인용 수 4

한 줄 요약

이 논문은 객체 지향 프로그래밍을 통해 순환 행렬과 토플리츠 행렬에 대한 최적화된 저장 방식과 빠른 산술 연산을 제공하는 MATLAB 도구상자 smt를 소개한다. 이는 투명성과 확장성을 유지하면서도 효율적인 계산을 가능하게 하며, 구조적 행렬 연산, 테스트 행렬, 순환 행렬 조절자, 토플리츠 선형 시스템을 위한 두 가지 빠른 해법기를 지원한다. 이는 공개적으로 이용 가능한, 생산용으로 사용 가능한 구조적 행렬 소프트웨어 분야의 핵심적 공백을 메운다.

ABSTRACT

We introduce the smt toolbox for Matlab. It implements optimized storage and fast arithmetics for circulant and Toeplitz matrices, and is intended to be transparent to the user and easily extensible. It also provides a set of test matrices, computation of circulant preconditioners, and two fast algorithms for Toeplitz linear systems.

연구 동기 및 목표

• 순환 및 토플리츠 행렬과 같은 구조적 행렬을 위한 공개적으로 이용 가능한 최적화된 소프트웨어의 부족을 해결하기 위해.
• 스parser 행렬을 초월하여 MATLAB의 기본 기능을 구조적 행렬 유형에 대한 1급 지원으로 확장하기 위해.
• 스퍼스 행렬을 다룰 때와 유사하게 투명하고 사용자 우호적이며 확장 가능한 프레임워크를 제공하기 위해.
• 연구자들과 엔지니어들이 구조적 행렬 알고리즘을 다시 구현하지 않고도 빠르게 사용할 수 있도록 하기 위해.
• 이미지 복원, 신호 처리, 대규모 선형 시스템 해법과 같은 성능이 중요한 응용 분야를 지원하기 위해.

제안 방법

• 도구상자는 MATLAB의 객체 지향 프로그래밍 시스템을 사용하여 새로운 데이터 유형을 정의한다: 순환 행렬을 위한 smcirc와 토플리츠 행렬을 위한 smtoep.
• 기본 산술 연산자와 함수를 오버로드하여 기존 MATLAB 워크플로우와의 투명한 통합을 보장한다.
• 메모리 사용량을 줄이고 행렬 연산을 가속화하기 위해 최적화된 저장 형식을 구현한다.
• 표준 테스트 행렬을 구조적 형식으로 생성하기 위한 갤러 기능(smtgallery)을 포함한다.
• 고성능을 위해 MEX 파일을 통해 인터페이스된 두 가지 빠른 토플리츠 시스템 해법기(toms729(Fortran 기반) 및 tlls(C 기반))를 통합한다.
• 사용자가 표시 및 동작 설정을 커스터마이징할 수 있도록 구성 기능(smtconfig)을 제공하며, 일관성을 유지하기 위해 경고 메시지를 사용한다.

실험 결과

연구 질문

• RQ1 어떻게 하면 MATLAB에서 스퍼스 행렬 연산만큼 사용자 투명성과 효율성을 갖춘 구조적 행렬 연산을 실현할 수 있는가?
• RQ2 MATLAB 환경 내에서 순환 및 토플리츠 행렬에 대한 최적화된 저장 및 산술 연산을 효과적으로 구현하는 가장 좋은 방법은 무엇인가?
• RQ3 순환 및 토플리츠를 초월한 새로운 유형의 구조적 행렬을 지원할 수 있는 모듈식이고 확장 가능한 도구상자 설계가 가능한가?
• RQ4 기존의 빠른 구조적 행렬 알고리즘을 고수준의, 생산용 소프트웨어 프레임워크로 노출시킬 수 있는가?
• RQ5 이러한 도구상자의 정확성, 이식성, 유지보수성을 보장하기 위해 필요한 핵심 소프트웨어 공학 실천 방법은 무엇인가?

주요 결과

• smt 도구상자는 MATLAB의 객체 지향 시스템을 활용하여 순환 및 토플리츠 행렬에 대한 최적화된 저장 및 빠른 산술 연산을 성공적으로 구현하였다.
• 도구상자는 행렬 생성, 산술 연산, 선형 시스템 해법을 포함한 전반적인 연산 범위를 지원하며, 수작업 최적화 코드와 비교해도 성능이 유사하다.
• 두 가지 컴파일된 MEX 기반 해법기(toms729 및 tlls)의 통합으로 대규모 토플리츠 시스템의 고성능 해법이 가능해졌다.
• smtgallery를 통해 표준 테스트 행렬의 포괄적인 세트를 제공하며, MATLAB의 gallery 기능을 모델링하여 실수 및 복소수 행렬 모두에 대한 문서화 및 지원을 제공한다.
• MATLAB 버전 6.5부터 호환되며, isfloat 호환성을 위한 패치가 포함되어 있으며, 여러 플랫폼용 사전 컴파일된 바이너리 파일도 포함되어 있다.
• 도구상자는 오픈소스이며 공개적으로 http://bugs.unica.it/smt/ 에서 이용 가능하며, 향후 추가적인 구조적 행렬 클래스 지원을 위한 계획이 수립되어 있다.