HOME

TheInfoList



OR:

In mathematics, the weakly chained diagonally dominant matrices are a family of
nonsingular matrices In linear algebra, an -by- square matrix is called invertible (also nonsingular or nondegenerate), if there exists an -by- square matrix such that :\mathbf = \mathbf = \mathbf_n \ where denotes the -by- identity matrix and the multiplicati ...
that include the strictly
diagonally dominant matrices In mathematics, a square matrix is said to be diagonally dominant if, for every row of the matrix, the magnitude of the diagonal entry in a row is larger than or equal to the sum of the magnitudes of all the other (non-diagonal) entries in that row ...
.


Definition


Preliminaries

We say row i of a complex matrix A = (a_) is strictly diagonally dominant (SDD) if , a_, >\textstyle, a_, . We say A is SDD if all of its rows are SDD. Weakly diagonally dominant (WDD) is defined with \geq instead. The
directed graph In mathematics, and more specifically in graph theory, a directed graph (or digraph) is a graph that is made up of a set of vertices connected by directed edges, often called arcs. Definition In formal terms, a directed graph is an ordered pa ...
associated with an m \times m complex matrix A = (a_) is given by the vertices \ and edges defined as follows: there exists an edge from i \rightarrow j if and only if a_ \neq 0.


Definition

A complex square matrix A is said to be weakly chained diagonally dominant (WCDD) if * A is WDD and * for each row i_1 that is ''not'' SDD, there exists a
walk Walking (also known as ambulation) is one of the main gaits of terrestrial locomotion among legged animals. Walking is typically slower than running and other gaits. Walking is defined by an 'inverted pendulum' gait in which the body vaults ov ...
i_1 \rightarrow i_2 \rightarrow \cdots \rightarrow i_k in the directed graph of A ending at an SDD row i_k.


Example

The m \times m matrix :\begin1\\ -1 & 1\\ & -1 & 1\\ & & \ddots & \ddots\\ & & & -1 & 1 \end is WCDD.


Properties


Nonsingularity

A WCDD matrix is nonsingular. ''Proof'': Let A=(a_) be a WCDD matrix. Suppose there exists a nonzero x in the null space of A. Without loss of generality, let i_1 be such that , x_, =1\geq, x_j, for all j. Since A is WCDD, we may pick a walk i_1\rightarrow i_2\rightarrow\cdots\rightarrow i_k ending at an SDD row i_k. Taking moduli on both sides of :-a_x_ = \sum_ a_x_j and applying the triangle inequality yields :\left, a_\\leq\sum_\left, a_\\left, x_j\\leq\sum_\left, a_\, and hence row i_1 is not SDD. Moreover, since A is WDD, the above chain of inequalities holds with equality so that , x_, =1 whenever a_\neq0. Therefore, , x_, =1. Repeating this argument with i_2, i_3, etc., we find that i_k is not SDD, a contradiction. \square Recalling that an
irreducible In philosophy, systems theory, science, and art, emergence occurs when an entity is observed to have properties its parts do not have on their own, properties or behaviors that emerge only when the parts interact in a wider whole. Emergence ...
matrix is one whose associated directed graph is
strongly connected In the mathematical theory of directed graphs, a graph is said to be strongly connected if every vertex is reachable from every other vertex. The strongly connected components of an arbitrary directed graph form a partition into subgraphs that ...
, a trivial corollary of the above is that an ''irreducibly diagonally dominant matrix'' (i.e., an irreducible WDD matrix with at least one SDD row) is nonsingular.


Relationship with nonsingular M-matrices

The following are equivalent: * A is a nonsingular WDD
M-matrix In mathematics, especially linear algebra, an ''M''-matrix is a ''Z''-matrix with eigenvalues whose real parts are nonnegative. The set of non-singular ''M''-matrices are a subset of the class of ''P''-matrices, and also of the class of inverse-p ...
. * A is a nonsingular WDD
L-matrix In mathematics, the class of L-matrices are those matrices whose off-diagonal entries are less than or equal to zero and whose diagonal entries are positive; that is, an L-matrix ''L'' satisfies :L=(\ell_);\quad \ell_ > 0; \quad \ell_\leq 0, \quad ...
; * A is a WCDD
L-matrix In mathematics, the class of L-matrices are those matrices whose off-diagonal entries are less than or equal to zero and whose diagonal entries are positive; that is, an L-matrix ''L'' satisfies :L=(\ell_);\quad \ell_ > 0; \quad \ell_\leq 0, \quad ...
; In fact, WCDD L-matrices were studied (by
James H. Bramble James Henry Bramble (December 1, 1930 – July 20, 2021) was an American mathematician known for his fundamental contributions in the development of the finite element methods, including the Bramble–Hilbert lemma,J. H. Bramble and S. R. Hilbert. ...
and B. E. Hubbard) as early as 1964 in a journal article in which they appear under the alternate name of ''matrices of positive type''. Moreover, if A is an n\times n WCDD L-matrix, we can bound its inverse as follows: :\left\Vert A^\right\Vert _\leq\sum_\left _\prod_^(1-u_)\right   where   u_=\frac\sum_^\left, a_\. Note that u_n is always zero and that the right-hand side of the bound above is \infty whenever one or more of the constants u_i is one. Tighter bounds for the inverse of a WCDD L-matrix are known.


Applications

Due to their relationship with M-matrices (see above), WCDD matrices appear often in practical applications. An example is given below.


Monotone numerical schemes

WCDD L-matrices arise naturally from monotone approximation schemes for
partial differential equations In mathematics, a partial differential equation (PDE) is an equation which imposes relations between the various partial derivatives of a multivariable function. The function is often thought of as an "unknown" to be solved for, similarly to ...
. For example, consider the one-dimensional Poisson problem :u^(x) + g(x)= 0   for   x \in (0,1) with
Dirichlet boundary conditions In the mathematical study of differential equations, the Dirichlet (or first-type) boundary condition is a type of boundary condition, named after Peter Gustav Lejeune Dirichlet (1805–1859). When imposed on an ordinary or a partial differential ...
u(0)=u(1)=0. Letting \ be a numerical grid (for some positive h that divides unity), a monotone finite difference scheme for the Poisson problem takes the form of :-\fracA\vec + \vec = 0   where  
vec Vec may mean: Mathematics: * vec(''A''), the vectorization of a matrix ''A''. * Vec denotes the category of vector spaces over the reals. Other: * Venetian language (Vèneto), language code. * Vecuronium, a muscle relaxant. * vec, a sentient mora ...
j = g(jh) and :A = \begin2 & -1\\ -1 & 2 & -1\\ & -1 & 2 & -1\\ & & \ddots & \ddots & \ddots\\ & & & -1 & 2 & -1\\ & & & & -1 & 2 \end. Note that A is a WCDD L-matrix.


References

{{reflist Matrices