% Comparaison : vitesse de resolution matlab/moncode
%
function [Tmat Tper] = SpeedChol(Dim,n0,nbIt)
% nbIt=nb de resolutions
% Dim: dimension du pb
% N= nb de DL
Tmat=zeros(nbIt,2);
Tper=zeros(nbIt,2);

for iter = 1:nbIt
  n=n0*iter;
  N=n^Dim;
% construction du 2nd mbre aleatoire:
  f=rand(N,1);
  switch Dim
    case 1
        A=Laplacien1D(n);
    case 2
        A=Laplacien2D(n);
    case 3
        A=Laplacien3D(n);
    otherwise
        A=zeros(n,n);
        msg='mauvaise valeur de Dim';
  end
  Tper(iter,1)=N;
  Tmat(iter,1)=N;
  t=cputime;
  SolChol(A,f);
  Tper(iter,2)=cputime-t;
  
  t=cputime;
  A\f;
  Tmat(iter,2)=cputime-t;
end