Return to dft.h CVS log

Up to [local] / OpenXM_contrib2 / asir2000 / fft

Annotation of OpenXM_contrib2/asir2000/fft/dft.h, Revision 1.1

1.1     ! noro        1: /* $OpenXM: OpenXM/src/asir99/fft/dft.h,v 1.1.1.1 1999/11/10 08:12:27 noro Exp $ */
        !             2: #if LONG_IS_32BIT
        !             3: #if defined(VISUAL)
        !             4: typedef unsigned _int64 UL;
        !             5: #else
        !             6: typedef unsigned long long UL;
        !             7: #endif
        !             8: #else
        !             9: typedef unsigned long UL;
        !            10: #endif
        !            11:
        !            12: struct PrimesS {
        !            13:        int     prime, primroot;
        !            14:        int  bitwid, d;
        !            15: };
        !            16:
        !            17: #if USE_FLOAT
        !            18: #define NPrimes 536
        !            19: #else
        !            20: #define NPrimes 13681
        !            21: #endif
        !            22:
        !            23: #if USE_FLOAT
        !            24: typedef float ModNum;
        !            25:
        !            26: #define AxBplusCmodP(ans,anstype,a,b,c,p,pinv) \
        !            27:   { double _w_ = ((double)(a))*((double)(b))+((double)(c)), \
        !            28:       _dp_ = (double)(p); \
        !            29:       if ( _w_ >= _dp_ ) _w_ -= (_dp_*(double)((int)(_w_*(pinv)))); \
        !            30:       ans = (anstype) (_w_>=_dp_?_w_-_dp_: _w_); }
        !            31:
        !            32: #define AxBmodP(ans,anstype,a,b,p,pinv) \
        !            33:   { double _w_ = ((double)(a))*((double)(b)), _dp_ = (double)(p); \
        !            34:       _dp_ = (double)(p); \
        !            35:       if ( _w_ >= _dp_ ) _w_ -= (_dp_*(double)((int)(_w_*(pinv)))); \
        !            36:       ans = (anstype) (_w_>=_dp_?_w_-_dp_: _w_); }
        !            37:
        !            38: #define AxBmodPnostrchk(ans,anstype,a,b,p,pinv) \
        !            39:   { double _w_ = ((double)(a))*((double)(b)), _dp_ = (double)(p); \
        !            40:       if ( _w_ >= _dp_ ) _w_ -= (_dp_*(double)((int)(_w_*(pinv)))); \
        !            41:       ans = (anstype) _w_; }
        !            42:
        !            43: #define AxBplusCmodPnostrchk(ans,anstype,a,b,c,p,pinv) \
        !            44:   { double _w_ = ((double)(a))*((double)(b))+((double)(c)), \
        !            45:       _dp_ = (double)(p); \
        !            46:       if ( _w_ >= _dp_ ) _w_ -= (_dp_*(double)((int)(_w_*(pinv)))); \
        !            47:       ans = (anstype) _w_; }
        !            48: #else
        !            49:
        !            50: #include "inline.h"
        !            51:
        !            52: typedef unsigned int ModNum;
        !            53:
        !            54: #define AxBplusCmodP(ans,anstype,a,b,c,p,pinv) \
        !            55:   { DMAR(a,b,c,p,ans) }
        !            56:
        !            57: #define AxBmodP(ans,anstype,a,b,p,pinv) \
        !            58:   { DMAR(a,b,0,p,ans) }
        !            59:
        !            60: #define AxBmodPnostrchk(ans,anstype,a,b,p,pinv) \
        !            61:   { DMAR(a,b,0,p,ans) }
        !            62:
        !            63: #define AxBplusCmodPnostrchk(ans,anstype,a,b,c,p,pinv) \
        !            64:   { DMAR(a,b,c,p,ans) }
        !            65: #endif
        !            66:
        !            67: #define AstrictmodP(a,p) ((a)>=(p)?((a)-(p)):(a))
        !            68:
        !            69: #define AplusBmodP(a,b,p) (((p)-(b))>(a)?((a)+(b)): ((a)-((p)-(b))))
        !            70: #define A_BmodP(a,b,p) ((a)>=(b)?((a)-(b)): (((p)-(b))+(a)))
        !            71:
        !            72: void C_DFT_FORE(ModNum *,int,int,int,ModNum *,
        !            73: #ifdef POWA_STRIDE
        !            74:        int
        !            75: #else
        !            76:        ModNum *,int,ModNum,double,ModNum *
        !            77: #endif
        !            78: );
        !            79:
        !            80: void C_DFT_BACK(ModNum *,int,int,int,ModNum *,
        !            81: #ifdef POWA_STRIDE
        !            82:                 int
        !            83: #else
        !            84:                 ModNum *,int,int,int,ModNum,ModNum,double,ModNum *
        !            85: #endif
        !            86: );
        !            87:
        !            88: void C_PREP_ALPHA(ModNum,int,int,int,ModNum *,ModNum,double);