23.12.2012 Views

MATLAB Function Reference Volume 1: A - E - Bad Request

MATLAB Function Reference Volume 1: A - E - Bad Request

MATLAB Function Reference Volume 1: A - E - Bad Request

SHOW MORE
SHOW LESS

You also want an ePaper? Increase the reach of your titles

YUMPU automatically turns print PDFs into web optimized ePapers that Google loves.

2dmperm<br />

Purpose Dulmage-Mendelsohn decomposition<br />

Syntax p = dmperm(A)<br />

[p,q,r] = dmperm(A)<br />

[p,q,r,s] = dmperm(A)<br />

dmperm<br />

Description If A is a reducible matrix, the linear system Ax =<br />

b can be solved by permuting<br />

A to a block upper triangular form, with irreducible diagonal blocks, and then<br />

performing block backsubstitution. Only the diagonal blocks of the permuted<br />

matrix need to be factored, saving fill and arithmetic in the blocks above the<br />

diagonal.<br />

See Also sprank<br />

p = dmperm(A) returns a row permutation p so that if A has full column rank,<br />

A(p,:) is square with nonzero diagonal. This is also called a maximum<br />

matching.<br />

[p,q,r] = dmperm(A) where A is a square matrix, finds a row permutation p<br />

and a column permutation q so that A(p,q) is in block upper triangular form.<br />

The third output argument r is an integer vector describing the boundaries of<br />

the blocks: The kth block of A(p,q) has indices r(k):r(k+1)-1.<br />

[p,q,r,s] = dmperm(A), where A is not square, finds permutations p and q<br />

and index vectors r and s so that A(p,q) is block upper triangular. The blocks<br />

have indices (r(i):r(i+1)-1, s(i):s(i+1)-1).<br />

If A is not square, the first block may have more columns and the last block may<br />

have more rows. All other blocks are square. dmperm permutes nonzeros to the<br />

diagonals of square blocks, but does not do this for non-square blocks.<br />

In graph theoretic terms, the diagonal blocks correspond to strong Hall<br />

components of the adjacency graph of A.<br />

2-455

Hooray! Your file is uploaded and ready to be published.

Saved successfully!

Ooh no, something went wrong!