s_h
double y, z, r, s, t, w, s_h, s_l, x_h, x_l, zz[3], ee[2], z_h,
s_h = (double) ((float) s);
s_l = -v * (s_h * r - (w - s_h * t));
z_l = (s_h * s_h - z_h);
zz[2] = z_l + s_l * (s + s_h);
v = s_h * ee[0];
long double y, z, r, s, t, w, s_h, s_l, x_h, x_l, zz[3], ee[2], z_h,
s_h = s; HALF(s_h);
s_l = -v * (s_h * r - (w - s_h * t));
z_l = (s_h * s_h - z_h);
zz[2] = z_l + s_l * (s + s_h);
v = s_h * ee[0];