function [Aq , Bq , D] = SkalSymmGlsSparse (A , B)

[m ,n ] = size (A) ;
[mb,nb] = size (B) ;

if (m ~= mb)  
    error ('Zeilenanzahl von A und B muessen uebereinstimmen') 
end

D  = sparse(diag(1./sqrt(diag(A)))) ;
Aq = D*A*D ;
Bq = D*B   ;
Aq = (Aq+Aq')*0.5 ;