1function x = cs_ltsolve (L,b) %#ok 2%CS_LTSOLVE solve a sparse upper triangular system L'*x=b. 3% x = cs_ltsolve(L,b) computes x = L'\b, L must be lower triangular with a 4% zero-free diagonal. b must be a full vector. 5% 6% Example: 7% Prob = ssget ('HB/bcsstk01') ; L = cs_chol (Prob.A) ; n = size (L,1) ; 8% b = rand (n,1) ; x = cs_ltsolve (L,b) ; norm (L'*x-b) 9% 10% See also CS_LSOLVE, CS_USOLVE, CS_UTSOLVE, MLDIVIDE. 11 12% Copyright 2006-2012, Timothy A. Davis, http://www.suitesparse.com 13 14error ('cs_ltsolve mexFunction not found') ; 15