% TP3
% ---
% Ex 8: Methode de Cholesky
%---------------------------

function u = SolChol(A,f)
% resolution de A.u = f
% A supposee SDP

% Calcul de L triang inf telle que A = L.L'
% ATTENTION: chol(A) -> L' => A=chol(A)' * chol(A)
LT = chol(A);
% resolution descente-remontee:
y = LT'\f;
u = LT\y;