function relaxG()

global v
global LpU
global D
global omega;

    v = omega*(D\(LpU*v)) + (1 - omega)*v;
%    s1(i) = norm(v2,inf);
