함수와 부동 소수점 비교

Question

#include<stdio.h> 
#include<stdlib.h> 
#define abs(a) ((a)>0 ? a: -a) 
#define eps_sqrt 0.00000000000001 
#define it 100 

float sqrt(float x) 
/*The Square Root Function using the Newton's Method*/ 
{ 
    int it_sqrt=0; 
    float a_sqrt = x/2; 
    while ((abs((a_sqrt*a_sqrt)-(x))>=eps_sqrt) && (2.0*a_sqrt != 0) && (it_sqrt<=it)) 
    { 
     a_sqrt = a_sqrt - ((a_sqrt*a_sqrt)-(x)/(2.0*a_sqrt)); 
     it_sqrt++; 
    } 
    return a_sqrt; 
} 

int main() 
{ 
    printf("%.5f\n", sqrt(5)); 
    system ("pause"); 
} 

필자는 Newton의 반복 방법을 사용하여 파이썬에서 제곱근을 찾았으며 완벽하게 작동했습니다. 저는 C 언어에 익숙하지 않아이 기능이 왜 작동하지 않는지 이해할 수 없습니다. 실행할 때마다 "-1. # INF0A"가 반환됩니다. 도움이 될 것입니다.

---

편집 : 나는 0.000001에 EPS를 바꾸는 시도하고 또한 작동하지 않았다.

Answer 1

double mysqrt(double x){ 
    double eps=pow(10,-10); 
    double x0 = 0.0; 
    double x1 = x/2.0; 
    while(fabs(x1 - x0)>eps){ 
     x0 = x1; 
     x1 = x0 + (x - x0*x0)/x0/ 2.0; 
    } 
    return x1; 
} 

매크로 확장
abs((a_sqrt*a_sqrt)-(x))
확장 (((a_sqrt*a_sqrt)-(x))>0 ? (a_sqrt*a_sqrt)-(x): -(a_sqrt*a_sqrt)-(x))
NG : -(a_sqrt*a_sqrt)-(x)

abs((a_sqrt*a_sqrt- x))
확장 (((a_sqrt*a_sqrt- x))>0 ? (a_sqrt*a_sqrt- x): -(a_sqrt*a_sqrt- x))

재
#define abs(a) ((a)>0 ? a: -a)

#define abs(a) ((a)>0 ? a: -(a))

Answer 2

더 큰 엡실론을 사용해보십시오. 아마도 파이썬은 수레 대신에 더블을 사용합니다.

Answer 3

이 줄을 변경 :

   a_sqrt = a_sqrt - ((a_sqrt*a_sqrt)-(x)/(2.0*a_sqrt)); 

나를 위해 작품

   a_sqrt = a_sqrt - ((a_sqrt*a_sqrt - x)/(2.0*a_sqrt)); 

에.

Answer 4

double을 실제로 사용하는 경우가 드문 경우 중 하나입니다. 부동 소수점의 정밀도가 eps_sqrt보다 훨씬 낮은 참고 :

[mic@mic-nb tmp]$ cat tmp2.c 
#include <stdio.h> 
#include <math.h> 

int main() { 
    double a = sqrtl(2.0); 
    printf("%1.20f\n", a - (float) a); 
} 
[mic@mic-nb tmp]$ gcc tmp2.c; ./a.out 
0.00000002420323430563 
vs. your value of: 
0.00000000000001 

그래서 프로그램이, 대부분의 경우, 종료하지 않습니다.