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