double a, b;
  double diff = a -b;

if (abs(diff)/(abs(a) + abs(b)) < EPSILON)
      my_diff = 0;
  else
      my_diff = sqrt(diff)

double diff = a-b;
if (diff <= 0)
    my_diff = 0;
else
    my_diff = sqrt(diff);