1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17
|
function x = cs_cholsol (A,b,order) %#ok
%CS_CHOLSOL solve A*x=b using a sparse Cholesky factorization.
% x = cs_cholsol(A,b) computes x = A\b, where A sparse symmetric positive
% definite, and b is a full vector. A 3rd input parameter allows the
% ordering to be modified: 0: natural, 1:amd(A), 2: amd(S'*S) where S=A except
% with no dense rows, 3:amd(A'*A). The default ordering option is 1.
%
% Example:
% Prob = ssget ('HB/bcsstk01') ; A = Prob.A ; b = rand (size (A,1),1) ;
% x = cs_cholsol (A,b) ; norm (A*x-b)
%
% See also CS_CHOL, CS_AMD, CS_LUSOL, CS_QRSOL, MLDIVIDE.
% CXSparse, Copyright (c) 2006-2022, Timothy A. Davis. All Rights Reserved.
% SPDX-License-Identifier: LGPL-2.1+
error ('cs_cholsol mexFunction not found') ;
|