最速化するにはどうしたらいいでしょうか？

質問

#include <stdio.h>
#include <math.h>
#include <time.h>
double abssucc(double x){
return x / (1 - pow(2,-53));
}
double abspred(double x){
return x * (1 - pow(2,-53));
}
double succ (double x){
if(abs(x) < pow(2,-1022)) return x + pow(2,-1074);
if (x >= 0) return abssucc(x);
else return abspred(x);
}
double pred(double x){
if(abs(x) < pow(2,-1022)) return x - pow(2,-1074);
if (x >= 0) return abspred(x);
else return abssucc(x);
}

int main(){
clock_t start, end;

double r= 0;
int i=0;
start = clock();

for(i;i<1000000;i++){

r= succ(1);

printf("%.50f\n",r);

}end = clock();
printf("ループ100万回の時間: %f秒\n", (double)(end - start) / CLOCKS_PER_SEC);

}

このプログラムをpow関数使わないでもっと速く動かしたいのですが、やり方がわかりません。

教えて下さい！お願いします。

回答No.1

kmee · ベストアンサー

使う値が決まっているのだから、先に計算して変数に入れておく。

> #include <time.h>
double pow_2_53 ;
double pow_2_1022 ;
double pow_2_1074 ;
/* とでもしておいて */
> double abssucc(double x){
return x / (1 - pow_2_53);
/* と使う */

> int main(){
pow_2_53 = pow(2,-53);
pow_2_1022 = pow(2,-1022);
pow_2_1074 = pow(2,-1074);
/*とmainの最初の方で代入 */

....って、2の-1074って、doubleの範囲外にならないですか?
doubleは64ビットでIEEE754準拠だとアンダーフローになると思うのですが。

お礼

ありがとうございます！これで試してみます。

回答No.3

kmee

あっと、見落してた。質問には直接関係ないけど
> abs(x)
このプログラムがC言語なら、int abs(int) です。
doubleに使うならfabsです。

お礼

そうなんですか。無知ですみません＞＜直してみます！

回答No.2

Tacosan

nextafter とか nexttoward は?

お礼

そうですね！それも一つの方法だと思うので、試してみます。