pow ( )【C语言库函数源代码】

【C语言库函数源代码】

【本程序在Dev C++ 4.9.9.2 下编译通过】

double my_pow(double x, double y)

{

   register double ret, value;

   double r = 1.0;

   long p = (long) y;

 

   if (x == 0.0 && y > 0.0)

      return 0.0;

   if (y == (double) p)

   {

     

      if (p == 0)

        return 1.0;

      if (p < 0)

      {

        p = -p;

        x = 1.0 / x;

      }

      while (1)

      {

        if (p & 1)

           r *= x;

        p >>= 1;

        if (p == 0)

           return r;

        x *= x;

      }

   }

   __asm__(

      "fmul  %%st(1);"

      "fst   %%st(1);"

      "frndint;\n\t"

      "fxch;\n\t"

      "fsub %%st(1);\n\t"

      "f2xm1;\n\t"

      : "=t" (ret), "=u" (value)

      :  "0" (log2 (x)), "1" (y)

   );

   ret += 1.0;

   __asm__(

      "fscale"

      : "=t" (ret)

      : "0" (ret), "u" (value)

   );

   return ret;

}

//这是一个求x的y次幂的函数，因为要求是浮点数类型的，所以代码复杂了许多。

//如果仅仅是求整数的x的y次幂，那么就简单了许多。

int main()

{

   printf("%f\n",my_pow(0.0,0.0));

   printf("%f\n",my_pow(1024.0,0.0));

   printf("%f\n",my_pow(2.0,10.0));

   printf("%f\n",pow(2.5,12.6));

   printf("%f\n",my_pow(2.5,12.6));

   system("pause");

   return 0;

}