浮動小数点モジュロ演算

Question

私は三角法の範囲削減操作を実装しようとしています。しかし、代わりに、着信データに対してモジュロπ/ 2演算を実行するほうが良いかもしれないと思います。私はどのアルゴリズムが存在し、32ビットIEEE 754浮動小数点でこの演算に効率的であるのか疑問に思っていましたか？

これはアセンブリで実装する必要があります。したがって、fmod、除算、乗算などは1つの命令では利用できません。私のプロセッサは16ビットワードを使用し、32ビット浮動小数点加算、減算、乗算、除算、平方根、余弦、サインを実装しています。私は余弦と正弦に値を入力するための範囲の縮小（モジュラス）が必要です。

Answer 1

私は、標準ライブラリのfmod()がほとんどの場合に最良の選択と考えます。いくつかの簡単なアルゴリズムについては、linkをご覧ください。私のマシン上で

、fmod()用途最適化されたインラインアセンブリコード（/usr/include/bits/mathinline.h）：

#if defined __FAST_MATH__ && !__GNUC_PREREQ (3, 5) 
__inline_mathcodeNP2 (fmod, __x, __y, \ 
    register long double __value;       \ 
    __asm __volatile__         \ 
    ("1: fprem\n\t"       \ 
    "fnstsw %%ax\n\t"        \ 
    "sahf\n\t"         \ 
    "jp 1b"        \ 
    : "=t" (__value) : "0" (__x), "u" (__y) : "ax", "cc");   \ 
    return __value) 
#endif 

は、だから、実際に計算する専用のCPU命令（FPREM）を使用しています。

Answer 2

多分私はここでポイントを逃しているが、 あなたは単にfmodを使用して何かに対して何かを持っていますか？

double theta = 10.4; 
const double HALF_PI = 2 * atan(1); 
double result = fmod(theta, HALF_PI); 

Answer 3

したいアルゴリズム、0とnいくつかの係数の間の浮動小数点valueを制限する例：pi mod e（3.14159265358979 MOD 2.718281828459045）

3.14159265358979/2.718281828459045 
    = 1.1557273497909217179 

Trunc(1.1557273497909217179) 
    = 1 

1.1557273497909217179 - 1 
    = 0.1557273497909217179 

0.1557273497909217179 * e 
    = 0.1557273497909217179 * 2.718281828459045 
    = 0.42331082513074800 

PI modのための

Double fmod(Double value, Double modulus) 
{ 
    return value - Trunc(value/modulus)*modulus; 
} 

e = 0.42331082513074800