procedure ssum(var x1,x2,x3: scalar);
begin { ssum }
  underflow:=false;
  r_sum(down,x1.l,x2.l,x3.l);
  r_sum(up,x1.u,x2.u,x3.u);
  if underflow then snormalize(x3);
end { ssum };

