clear all
% FEM-Berechnungsmodell:
[xy km ep kr bk] = finmod_m ('Rendbrem.dat') ;
% Systemsteifigkeitsmatrix und -belastungsvektor:
[succ Ksys fq]   = syswbc_m (xy , km , ep , kr , bk) ;

Ks = SymBand2Sparse (Ksys) ;       % "Sparse matrix"
Kf = SymBand2Quad   (Ksys) ;       % "Voll besetzte Matrix"

t1 = cputime ;
x  = Ks \ fq ;
disp (['Zeit für Standard-Solver (Sparse): ' num2str(cputime-t1)]) ;

t1 = cputime ;
R  = chol (Ks) ;                   % Cholesky-Zerlegung
disp (['Zeit für Function chol (Sparse):   ' num2str(cputime-t1)]) ;
t1 = cputime   ;
x  = R\(R'\fq) ;                   % Vorwaert-Rueckwaertseinsetzen
disp (['Zeit für VorwRueckwEins (Sparse):  ' num2str(cputime-t1)]) ;

t1 = cputime ;
x  = Kf \ fq ;
disp (['Zeit für Standard-Solver (Full):   ' num2str(cputime-t1)]) ;

t1 = cputime ;
R  = chol (Kf) ;                   % Cholesky-Zerlegung
disp (['Zeit für Function chol (Full):     ' num2str(cputime-t1)]) ;
t1 = cputime   ;
x  = R\(R'\fq) ;                   % Vorwaert-Rueckwaertseinsetzen
disp (['Zeit für VorwRueckwEins (Full):    ' num2str(cputime-t1)]) ;

t1 = cputime ;
opts.SYM = true ; opts.POSDEF = true ;
x  = linsolve (Kf , fq , opts) ;
disp (['Zeit für linsolve (opts):          ' num2str(cputime-t1)]) ;