Cauchy matrix
From Wikipedia, the free encyclopedia
In mathematics, the Cauchy matrix is an m×n matrix A, whose elements are given by
where xi and yj are elements of a field , and where (xi) and (yj) are injective sequences (they do not contain repeated elements; elements are distinct).
Contents |
[edit] Properties
- When m=n, the determinant, known as a Cauchy determinant, is given explicitly by
- As a consequence of the injectivity of (xi) and (yj), all square Cauchy matrices are invertible.
- Every submatrix of a Cauchy matrix is itself a Cauchy matrix.
[edit] Examples
The Hilbert matrix is a special case of the Cauchy matrix, where
[edit] Generalization
A matrix C is called Cauchy-like if it is of the form
Defining X = diag(xi), Y = diag(yi), one sees that both Cauchy and Cauchy-like matrices satisfy the displacement equation
- XC + CY = rsT
(with for the Cauchy one). Hence Cauchy-like matrices have a common displacement structure, which can be exploited while working with the matrix. For example, there are known algorithms in literature for
- approximate Cauchy matrix-vector multiplication with O(nlogn) ops (e.g. the fast multipole method),
- (pivoted) LU factorization with O(n2) ops (GKO algorithm), and thus linear system solving,
- approximated or unstable algorithms for linear system solving in O(nlog2n).
Here n denotes the size of the matrix (one usually deals with square matrices, though all algorithms can be easily generalized to rectangular matrices).
[edit] References
- A. Gerasoulis, A fast algorithm for the multiplication of generalized Hilbert matrices with vectors, Mathematics of Computation, 1988; vol. 50, no. 181, pp. 179-188.
- I. Gohberg, T. Kailath, V. Olshevsky, Fast Gaussian elimination with partial pivoting for matrices with displacement structure. Mathematics of Computation, 1995; vol. 64, no. 212, pp. 1557-1576.
- P.G. Martinsson, M. Tygert, V. Rokhlin, An O(Nlog2N) algorithm for the inversion of general Toeplitz matrices, Computers & Mathematics with Applications, 2005; 50, pp. 741-752.