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