- Library Home /
- Search Collections /
- Open Collections /
- Browse Collections /
- UBC Theses and Dissertations /
- Ildl factorizations for sparse skew-symmetric matrices
Open Collections
UBC Theses and Dissertations
UBC Theses and Dissertations
Ildl factorizations for sparse skew-symmetric matrices He, Shiwen
Abstract
Our goal is to solve a sparse skew-symmetric linear system efficiently. We propose a slight modification to the Bunch LDL\T factorization with partial pivoting strategy for skew-symmetric matrices, which saves approximately one third of the overall number of swaps of rows and columns on average. We also develop a rook pivoting strategy for this LDL\T factorization in Crout order. We derive an incomplete LDL\T factorization based on the full LDL\T factorization, with both partial pivoting strategy and rook pivoting strategy. The incomplete factorization can be used as a preconditioner for Krylov subspace solvers. Various column versions of ILUTP factorizations are investigated. We show that the approach taken saves approximately half of the work and storage compared to standard ILU factorizations that do not exploit the structure of the matrix, and the rook pivoting strategy is often superior to the partial pivoting strategy, in particular when the matrix is ill-conditioned and the overall number of elements in the rows of the factors is restricted.
Item Metadata
Title |
Ildl factorizations for sparse skew-symmetric matrices
|
Creator | |
Publisher |
University of British Columbia
|
Date Issued |
2012
|
Description |
Our goal is to solve a sparse skew-symmetric linear system efficiently. We propose a slight modification to the Bunch LDL\T factorization with partial pivoting strategy for skew-symmetric matrices, which saves approximately one third of the overall number of swaps of rows and columns on average. We also develop a rook pivoting strategy for this LDL\T factorization in Crout order. We derive an incomplete LDL\T factorization based on the full LDL\T factorization, with both partial pivoting strategy and rook pivoting strategy. The incomplete factorization can be used as a preconditioner for Krylov subspace solvers. Various column versions of ILUTP factorizations are investigated. We show that the approach taken saves approximately half of the work and storage compared to standard ILU factorizations that do not exploit the structure of the matrix, and the rook pivoting strategy is often superior to the partial pivoting strategy, in particular when the matrix is ill-conditioned and the overall number of elements in the rows of the factors is restricted.
|
Genre | |
Type | |
Language |
eng
|
Date Available |
2012-04-23
|
Provider |
Vancouver : University of British Columbia Library
|
Rights |
Attribution-NonCommercial-NoDerivatives 4.0 International
|
DOI |
10.14288/1.0103449
|
URI | |
Degree | |
Program | |
Affiliation | |
Degree Grantor |
University of British Columbia
|
Graduation Date |
2012-05
|
Campus | |
Scholarly Level |
Graduate
|
Rights URI | |
Aggregated Source Repository |
DSpace
|
Item Media
Item Citations and Data
Rights
Attribution-NonCommercial-NoDerivatives 4.0 International