The LR Cholesky algorithm for symmetric hierarchical matrices

P Benner, T Mach - Linear Algebra and its Applications, 2013 - Elsevier
Linear Algebra and its Applications, 2013Elsevier
We investigate the application of the LR Cholesky algorithm to symmetric hierarchical
matrices, symmetric simple structured hierarchical matrices and symmetric hierarchically
semiseparable (HSS) matrices. The data-sparsity of these matrices make the otherwise
expensive LR Cholesky algorithm applicable, as long as the data-sparsity is preserved. We
will see in an example that the ranks of the low rank blocks grow and the data-sparsity gets
lost. We will explain this behavior by applying a theorem on the structure preservation of …
We investigate the application of the LR Cholesky algorithm to symmetric hierarchical matrices, symmetric simple structured hierarchical matrices and symmetric hierarchically semiseparable (HSS) matrices. The data-sparsity of these matrices make the otherwise expensive LR Cholesky algorithm applicable, as long as the data-sparsity is preserved. We will see in an example that the ranks of the low rank blocks grow and the data-sparsity gets lost. We will explain this behavior by applying a theorem on the structure preservation of diagonal plus semiseparable matrices under LR Cholesky transformations. Therefore we have to give a new more constructive proof for the theorem. We will show that the structure of H-matrices is almost preserved and so the LR Cholesky algorithm is of almost quadratic complexity for H-matrices.
Elsevier
以上显示的是最相近的搜索结果。 查看全部搜索结果