کد مقاله | کد نشریه | سال انتشار | مقاله انگلیسی | نسخه تمام متن |
---|---|---|---|---|
4598624 | 1631093 | 2016 | 36 صفحه PDF | دانلود رایگان |
The paper presents a new algorithm to compute the LU-factorization of a matrix represented in a quasiseparable or semiseparable form (i.e., using generators). It obtains the quasiseparable representations of the factors L and U of an N×NN×N block matrix via O(N)O(N) arithmetic operations on the block entries. The algorithm uses recursions based exclusively on unitary transformations which provide numerical stability even in singular cases. The method of the paper is based on the theory developed in [1] and provides an alternative to the approach proposed in [7] for strongly regular matrices. The algorithm presented here works also for some matrices with possibly singular principle submatrices. The results of numerical tests show that also for strongly regular matrices the new algorithm is comparable with the previous methods.
Journal: Linear Algebra and its Applications - Volume 502, 1 August 2016, Pages 5–40