为何sqrt不能正常工作？

winnie

环境：mplab6.61+picc8.35pl1
现象：在c程序中使用double变量，做了初始化, + - * \ 都正常，但是sqrt会出现一个莫名其妙的结果
，在build的option中尝试过选择32位或24位格式的浮点数，但都不行
pow()一样不能工作
比sqrt还要糟，debug时跑了半天，跳出来一个对话框说：no error message
本来还想自己编一个sqrt的，可是现在觉得肯定要从编译器方面找原因，否则发现一个不对且不就要自己写一个？

winnie

CCS 的庫函數 math.h 裡面的源程序

////////////////////////////////////////////////////////////////////////////
// float sqrt(float x)
////////////////////////////////////////////////////////////////////////////
// Description : returns the square root of x
// Date : N/A
//
float sqrt(float x)
{
   float y, res;
   BYTE *p;

   #ifdef _ERRNO
   if(x < 0)
   {
      errno=EDOM;
   }
   #endif

   if( x<=0.0)
      return(0.0);

   y=x;
   p=&y;
   (*p)=(BYTE)((((int16)(*p)) + 127) >> 1);

   do {
      res=y;
      y+=(x/y);
      (*p)--;
   } while(res != y);

   return(res);
}




小心，CCS 和 Hi-Tech 的浮點格式好像不一樣喔，可能沒法在 PICC 上使。

winnie

ANSI-C 都可以啦

float sqrt ( float x )
{
   float y, z;
   unsigned char *p;

   if(x<=0.0)  return(0.0);

   y = x;
// p = &y; (*p)=(int8)((((int16)(*p)) + 127) >> 1);   // CCS 直接把指數除以二，減少運算次數

   while (1)
   {
      z = y;
      y =(y + x/y) / 2.0;

      if (y == z) break;
   }
   return ( y );
}