function B = ddomaltf(A)
%  B = ddomaltf(A) returns a diag. dom. matrix
%      by modifying the diagonal of A
d = diag(A) ;
a = abs(d) ;
f= sum(abs(A), 2) - a ;
i = find(f >= a) ;
A(i,i) = A(i,i) + diag(f(i)) ;
B = A;
