return (t >= 0)?t:slowpow(x,y,z);
return (t >= 0)?t:__slowpow(x,y,z);
static double power1(double x, double y) {  double z,a,aa,error, t,a1,a2,y1,y2;  z = log2(x,&aa,&error);  t = y*134217729.0;  y1 = t - (t-y);  y2 = y - y1;  t = z*134217729.0;  a1 = t - (t-z);  a2 = z - a1;  a = y*z;  aa = ((y1*a1-a)+y1*a2+y2*a1)+y2*a2+aa*y;  a1 = a+aa;  a2 = (a-a1)+aa;  error = error*ABS(y);  t = __exp1(a1,a2,1.9e16*error);  return (t >= 0)?t:slowpow(x,y,z);}