C++平方根と再帰

Question

私の平方根プログラムは何らかの理由で、ほとんどの入力に対して何を得るべきであるかというだけの答えです。なぜこのことが分かりませんか？特定の入力だけが間違っています。私はまた、答えが与えられた後、非常に最後にセグメンテーションフォルトを取得し、それがどちらかであるかどうかわからない。あなたがreturn文

double newton(double num, double a) 
{ 


    if ((abs(a*a - num) <= FLT_EPSILON)) 
    { 
     return a; 
    } 
    else 
    { 
     newton(num, (a+num/a)/2); 
     ^^^^^^^^^^^^^^^^^^^^^^^^^^^ 
    } 


} 

を配置するのを忘れたので

#include<iostream> 
#include<cmath> 
#include<cfloat> 
#include<string> 
#include <cstdlib> 
using namespace std; 

//declare sqroot function calls recursive function newton 
double sqroot1(double num); 
double newton(double num, double guess); 
int main(int argc, char **argv) 
{ 

    for(int i = 0 ; i < argc ; i++) 
    { 
     cout<<"sqroot("<<argv[i+1]<<") is "<< sqroot1(atoi(argv[i+1]))  <<endl; 
    } 
} 

double newton(double num, double a) 
{ 
    if ((abs(a*a - num) <= FLT_EPSILON)) 
    { 
     return a; 
    } 
    else 
    { 
     newton(num, (a+num/a)/2); 
    } 
} 

double sqroot1(double num) 
{ 
    double sqrt = newton(num,num/2); 
    return sqrt; 
} 

Answer 1

関数が未定義の動作をしているがなければならない

return newton(num, (a+num/a)/2); 

Answer 2

あなたがargv[argc]にアクセスしようとするので、あなたがクラッシュを取得しています、これはargv配列の最後を過ぎています。

FLT_EPSILON（floatの値）と比較しているため、数字が間違っている可能性があります。は2倍に使用されます。しかし、これらのいずれかと比較することは、この問題の正しいアプローチではありません。 aの値に比例してスケールされた差異を調べる必要があります。

Answer 3

私が最初に固定されなければならない二つの基本的なエラーを参照してください。

1. をあなたの主な機能は、リターンを持っていません。あなたの最後の}の前に、あなたの主な機能の最後にreturn 0;を置くべきです。

2. 存在しないargv[argc]にアクセスしようとしています。その配列の最後の要素はargv[argc - 1]です。ループステートメントを変更して、配列の存在しない部分を読み込もうとしないようにしてください。

Answer 4

1201ProgramAlarmの回答に加えて、なぜ次のように悪いですか？

if ((abs(a*a - num) <= eps)) 

たとえば、epsは1e-15で、DBL_EPSILONより少し大きいとします。もし、aがやや大きければ（5以上のような）、最後の場所（ulp）のユニットがepsより大きく、条件が真であるのはa*a == numの場合のみです。これには到達せず、無限の再帰とクラッシュに終わるリスクがあります。今

のはnumは、初期推定のため

abs(num/2*num/2 - num) <= eps 

だから、例えば、sqroot1は（1E-20）が真である

abs(1e-40/4 - 1e-20) <= 1e-15 

を評価し、あなたがするような非常に小さいとしましょう5e-21を1e-20の平方根とする。これはひどい精度です。