Mercurial > octave-nkf
diff liboctave/UMFPACK/UMFPACK/MATLAB/umfpack.m @ 5164:57077d0ddc8e
[project @ 2005-02-25 19:55:24 by jwe]
author | jwe |
---|---|
date | Fri, 25 Feb 2005 19:55:28 +0000 |
parents | |
children |
line wrap: on
line diff
--- /dev/null Thu Jan 01 00:00:00 1970 +0000 +++ b/liboctave/UMFPACK/UMFPACK/MATLAB/umfpack.m Fri Feb 25 19:55:28 2005 +0000 @@ -0,0 +1,44 @@ +function [out1, out2, out3, out4, out5] = umfpack (in1, in2, in3, in4, in5) +% UMFPACK v4.4 is a MATLAB mexFunction for solving sparse linear systems. +% +% UMFPACK v4.4: | MATLAB approximate equivalent: +% --------------------------------------------------------------------- +% x = umfpack (A, '\', b) ; | x = A \ b +% | +% x = umfpack (b, '/', A) ; | x = b / A +% | +% [L,U,P,Q] = umfpack (A) ; | [m,n] = size (A) ; +% | I = speye (n) ; +% | Q = I (:, colamd (A)) ; +% | [L,U,P] = lu (A*Q) ; +% | +% [L,U,P,Q,R] = umfpack (A) ; | [m,n] = size (A) ; +% | I = speye (n) ; +% | Q = I (:, colamd (A)) ; +% | r = full (sum (abs (A), 2)) ; +% | r (find (r == 0)) = 1 ; +% | R = spdiags (r, 0, m, m) ; +% | [L,U,P] = lu ((R\A)*Q) ; +% | +% [P,Q,F,C] = umfpack (A, 'symbolic') | [m,n] = size (A) ; +% | I = speye (n) ; +% | Q = I (:, colamd (A)) ; +% | [count,h,parent,post] = ... +% | symbfact (A*Q, 'col') ; +% +% A must be sparse. It can be complex, singular, and/or rectangular. A must be +% square for '/' or '\'. b must be a full real or complex vector. For +% [L,U,P,Q,R] = umfpack (A), the factorization is L*U = P*(R\A)*Q. If A has a +% mostly symmetric nonzero pattern, then replace "colamd" with "amd" in the +% MATLAB-equivalent column in the table above. Type umfpack_details for more +% information. +% +% See also: lu_normest, colamd, amd. +% To use UMFPACK for an arbitrary b, see umfpack_solve. + +% UMFPACK Version 4.4, Copyright (c) 2005 by Timothy A. Davis. +% All Rights Reserved. Type umfpack_details for License. + +help umfpack +error ('umfpack mexFunction not found! Use umfpack_make to compile umfpack.') ; +