1function x = cs_lsolve (L,b) %#ok 2%CS_LSOLVE solve a sparse lower triangular system L*x=b. 3% x = cs_lsolve(L,b) computes x = L\b, L must be lower triangular with a 4% zero-free diagonal. b must be a column vector. x is full if b is full. 5% If b is sparse, x is sparse but the nonzero pattern of x is NOT sorted (it 6% is returned in topological order). 7% 8% Example: 9% Prob = ssget ('HB/bcsstk01') ; L = cs_chol (Prob.A) ; n = size (L,1) ; 10% b = rand (n,1) ; x = cs_lsolve (L,b) ; norm (L*x-b) 11% 12% See also CS_LTSOLVE, CS_USOLVE, CS_UTSOLVE, MLDIVIDE. 13 14% Copyright 2006-2012, Timothy A. Davis, http://www.suitesparse.com 15 16error ('cs_lsolve mexFunction not found') ; 17