Kayan Nokta Modulo Çalışması

Question

Trigonometri için aralık azaltma işlemini uygulamaya çalışıyorum. Ancak bunun yerine, gelen veriler üzerinde sadece bir modulo pi/2 işlemi gerçekleştirmenin daha iyi olabileceğini düşünüyorum. 32-bit IEEE 754 kayan nokta için hangi algoritmaların var olduğunu ve bu işlem için ne kadar verimli olduğunu merak ediyordum?

Bunu montajda uygulamak zorundayım, bu yüzden fmod, bölme, çarpma, vb. Sadece bir talimatla mevcut değildir. İşlemcim 16 bit kelimeler kullanıyor ve 32 bit kayan nokta ekleme, çıkarma, çarpma, bölme, karekök, kosinüs ve sinüs uyguladık. Sadece kosinüs ve sinüs değerlerini girmek için aralık küçültmeye (modül) ihtiyacım var.

Answer 1

Ben standart kütüphane en fmod() çoğu durumda en iyi seçim olacağını düşünüyorum. İşte birkaç basit algoritmanın bir tartışması için bir link. Benim makinede

, fmod() kullanımları optimize satır içi montaj kodu (/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 

Yani aslında hesaplanması için özel bir işlemci talimatı (fprem) kullanır.

Answer 2

Belki buradaki noktayı kaçırıyorum, ama sadece fmod kullanarak bir şey var mı?

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

Answer 3

istediğiniz algoritma, 0 ve n bazı katsayıları arasındaki kayan nokta value sınırlamak için: Örneğin 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 için

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

e = 0.42331082513074800