Ignore:
Timestamp:
09/28/06 18:15:21 (16 years ago)
Author:
r2d
Message:

added a PsyModel struct, too much members in it. some will be removed later

File:
1 edited

Legend:

Unmodified
Added
Removed
  • mppenc/branches/r2d/libmpcpsy/fft4g.c

    r59 r60  
    1818 */
    1919
    20 #include "mppenc.h"
     20#include "libmpcpsy.h"
    2121
    2222/* F U N C T I O N S */
     
    5757    float  xi;
    5858
    59     ENTER(30);
    6059    if ( n > 4) {
    6160        bitrv2  ( n, ip + 2, a );
     
    6968    a[0] += a[1];
    7069    a[1]  = xi;
    71     LEAVE(30);
    7270    return;
    7371}
     
    8482    double  delta;
    8583
    86     ENTER(31);
    8784    ip[0] = nw;
    8885    ip[1] = 1;
     
    106103        }
    107104    }
    108     LEAVE(31);
    109105    return;
    110106}
     
    118114    double  delta;
    119115
    120     ENTER(32);
    121116    ip[1] = nc;
    122117    if ( nc > 1 ) {
     
    130125        }
    131126    }
    132     LEAVE(32);
    133127    return;
    134128}
     
    142136    float  xr, xi, yr, yi;
    143137
    144     ENTER(33);
    145138    ip[0] = 0;
    146139    l     = n;
     
    235228        }
    236229    }
    237     LEAVE(33);
    238230    return;
    239231}
     
    246238    float  x0r, x0i, x1r, x1i, x2r, x2i, x3r, x3i;
    247239
    248     ENTER(34);
    249240    l = 2;
    250241    if ( n > 8 ) {
     
    291282        } while ( j += 2, j < l );
    292283    }
    293     LEAVE(34);
    294284    return;
    295285}
     
    303293    float  x0r, x0i, x1r, x1i, x2r, x2i, x3r, x3i;
    304294
    305     ENTER(35);
    306295    x0r   = a[ 0] + a[ 2];
    307296    x0i   = a[ 1] + a[ 3];
     
    401390        a[j + 15] = wk3r * x0i + wk3i * x0r;
    402391    } while ( j += 16, j < n );
    403     LEAVE(35);
    404     return;
    405 }
    406 
    407 extern void Cdecl cftmdl_3DNow_1 ( const int n, const int l, float* a, float* w );
    408 extern void Cdecl cftmdl_3DNow_2 ( const int n, const int l, float* a, float* w );
     392    return;
     393}
     394
     395// extern void Cdecl cftmdl_3DNow_1 ( const int n, const int l, float* a, float* w );
     396// extern void Cdecl cftmdl_3DNow_2 ( const int n, const int l, float* a, float* w );
    409397
    410398
     
    416404    float  x0r, x0i, x1r, x1i, x2r, x2i, x3r, x3i;
    417405
    418     ENTER(36);
    419406    m = l << 2;
    420407
     
    467454        a[j3 + 1] = wk1r * (x0i + x0r);
    468455    }
    469     LEAVE(36);
    470 
    471     ENTER(39);
     456
    472457    k1 = 0;
    473458    m2 = 2 * m;
     
    542527        } while ( j += 2, j < l+k+m );
    543528    }
    544     LEAVE(39);
    545     return;
    546 }
    547 
    548 
    549 static void
    550 cftmdl_3DNow ( const int n, const int l, float* a, float* w )
    551 {
    552     int    j, j1, j2, j3, k, k1, m, m2;
    553     float  wk1r, wk1i, wk2r, wk2i, wk3r, wk3i;
    554     float  x0r, x0i, x1r, x1i, x2r, x2i, x3r, x3i;
    555 
    556     ENTER(36);
    557     cftmdl_3DNow_1 (n,l,a,w);
    558     LEAVE(36);
    559 
    560     ENTER(39);
    561     m  = l << 2;
    562     k1 = 0;
    563     m2 = 2 * m;
    564     for ( k = m2; k < n; k += m2 ) {
    565         k1  += 2;
    566         wk2r = w[k1];
    567         wk2i = w[k1 + 1];
    568         wk1r = w[2*k1];
    569         wk1i = w[2*k1 + 1];
    570         wk3r = wk1r - 2 * wk2i * wk1i;
    571         wk3i = 2 * wk2i * wk1r - wk1i;
    572         j    = k;
    573         do {
    574             j1        = j  + l;
    575             j2        = j1 + l;
    576             j3        = j2 + l;
    577             x0r       = a[j]      + a[j1];
    578             x0i       = a[j + 1]  + a[j1 + 1];
    579             x1r       = a[j]      - a[j1];
    580             x1i       = a[j + 1]  - a[j1 + 1];
    581             x2r       = a[j2]     + a[j3];
    582             x2i       = a[j2 + 1] + a[j3 + 1];
    583             x3r       = a[j2]     - a[j3];
    584             x3i       = a[j2 + 1] - a[j3 + 1];
    585             a[j]      = x0r + x2r;
    586             a[j + 1]  = x0i + x2i;
    587             x0r      -= x2r;
    588             x0i      -= x2i;
    589             a[j2]     = wk2r * x0r - wk2i * x0i;
    590             a[j2 + 1] = wk2r * x0i + wk2i * x0r;
    591             x0r       = x1r - x3i;
    592             x0i       = x1i + x3r;
    593             a[j1]     = wk1r * x0r - wk1i * x0i;
    594             a[j1 + 1] = wk1r * x0i + wk1i * x0r;
    595             x0r       = x1r + x3i;
    596             x0i       = x1i - x3r;
    597             a[j3]     = wk3r * x0r - wk3i * x0i;
    598             a[j3 + 1] = wk3r * x0i + wk3i * x0r;
    599         } while ( j += 2, j < l + k );
    600 
    601         wk1r = w[2*k1 + 2];
    602         wk1i = w[2*k1 + 3];
    603         wk3r = wk1r - 2 * wk2r * wk1i;
    604         wk3i = 2 * wk2r * wk1r - wk1i;
    605         j    = k + m;
    606         do {
    607             j1        = j + l;
    608             j2        = j1 + l;
    609             j3        = j2 + l;
    610             x0r       = a[j]      + a[j1];
    611             x0i       = a[j + 1]  + a[j1 + 1];
    612             x1r       = a[j]      - a[j1];
    613             x1i       = a[j + 1]  - a[j1 + 1];
    614             x2r       = a[j2]     + a[j3];
    615             x2i       = a[j2 + 1] + a[j3 + 1];
    616             x3r       = a[j2]     - a[j3];
    617             x3i       = a[j2 + 1] - a[j3 + 1];
    618             a[j]      = x0r + x2r;
    619             a[j + 1]  = x0i + x2i;
    620             x0r      -= x2r;
    621             x0i      -= x2i;
    622             a[j2]     = -wk2i * x0r - wk2r * x0i;
    623             a[j2 + 1] = -wk2i * x0i + wk2r * x0r;
    624             x0r       = x1r - x3i;
    625             x0i       = x1i + x3r;
    626             a[j1]     = wk1r * x0r - wk1i * x0i;
    627             a[j1 + 1] = wk1r * x0i + wk1i * x0r;
    628             x0r       = x1r + x3i;
    629             x0i       = x1i - x3r;
    630             a[j3]     = wk3r * x0r - wk3i * x0i;
    631             a[j3 + 1] = wk3r * x0i + wk3i * x0r;
    632         } while ( j += 2, j < l+k+m );
    633     }
    634     LEAVE(39);
    635     return;
    636 }
     529    return;
     530}
     531
     532
     533// static void
     534// cftmdl_3DNow ( const int n, const int l, float* a, float* w )
     535// {
     536//     int    j, j1, j2, j3, k, k1, m, m2;
     537//     float  wk1r, wk1i, wk2r, wk2i, wk3r, wk3i;
     538//     float  x0r, x0i, x1r, x1i, x2r, x2i, x3r, x3i;
     539//
     540//     cftmdl_3DNow_1 (n,l,a,w);
     541//
     542//     m  = l << 2;
     543//     k1 = 0;
     544//     m2 = 2 * m;
     545//     for ( k = m2; k < n; k += m2 ) {
     546//         k1  += 2;
     547//         wk2r = w[k1];
     548//         wk2i = w[k1 + 1];
     549//         wk1r = w[2*k1];
     550//         wk1i = w[2*k1 + 1];
     551//         wk3r = wk1r - 2 * wk2i * wk1i;
     552//         wk3i = 2 * wk2i * wk1r - wk1i;
     553//         j    = k;
     554//         do {
     555//             j1        = j  + l;
     556//             j2        = j1 + l;
     557//             j3        = j2 + l;
     558//             x0r       = a[j]      + a[j1];
     559//             x0i       = a[j + 1]  + a[j1 + 1];
     560//             x1r       = a[j]      - a[j1];
     561//             x1i       = a[j + 1]  - a[j1 + 1];
     562//             x2r       = a[j2]     + a[j3];
     563//             x2i       = a[j2 + 1] + a[j3 + 1];
     564//             x3r       = a[j2]     - a[j3];
     565//             x3i       = a[j2 + 1] - a[j3 + 1];
     566//             a[j]      = x0r + x2r;
     567//             a[j + 1]  = x0i + x2i;
     568//             x0r      -= x2r;
     569//             x0i      -= x2i;
     570//             a[j2]     = wk2r * x0r - wk2i * x0i;
     571//             a[j2 + 1] = wk2r * x0i + wk2i * x0r;
     572//             x0r       = x1r - x3i;
     573//             x0i       = x1i + x3r;
     574//             a[j1]     = wk1r * x0r - wk1i * x0i;
     575//             a[j1 + 1] = wk1r * x0i + wk1i * x0r;
     576//             x0r       = x1r + x3i;
     577//             x0i       = x1i - x3r;
     578//             a[j3]     = wk3r * x0r - wk3i * x0i;
     579//             a[j3 + 1] = wk3r * x0i + wk3i * x0r;
     580//         } while ( j += 2, j < l + k );
     581//
     582//         wk1r = w[2*k1 + 2];
     583//         wk1i = w[2*k1 + 3];
     584//         wk3r = wk1r - 2 * wk2r * wk1i;
     585//         wk3i = 2 * wk2r * wk1r - wk1i;
     586//         j    = k + m;
     587//         do {
     588//             j1        = j + l;
     589//             j2        = j1 + l;
     590//             j3        = j2 + l;
     591//             x0r       = a[j]      + a[j1];
     592//             x0i       = a[j + 1]  + a[j1 + 1];
     593//             x1r       = a[j]      - a[j1];
     594//             x1i       = a[j + 1]  - a[j1 + 1];
     595//             x2r       = a[j2]     + a[j3];
     596//             x2i       = a[j2 + 1] + a[j3 + 1];
     597//             x3r       = a[j2]     - a[j3];
     598//             x3i       = a[j2 + 1] - a[j3 + 1];
     599//             a[j]      = x0r + x2r;
     600//             a[j + 1]  = x0i + x2i;
     601//             x0r      -= x2r;
     602//             x0i      -= x2i;
     603//             a[j2]     = -wk2i * x0r - wk2r * x0i;
     604//             a[j2 + 1] = -wk2i * x0i + wk2r * x0r;
     605//             x0r       = x1r - x3i;
     606//             x0i       = x1i + x3r;
     607//             a[j1]     = wk1r * x0r - wk1i * x0i;
     608//             a[j1 + 1] = wk1r * x0i + wk1i * x0r;
     609//             x0r       = x1r + x3i;
     610//             x0i       = x1i - x3r;
     611//             a[j3]     = wk3r * x0r - wk3i * x0i;
     612//             a[j3 + 1] = wk3r * x0i + wk3i * x0r;
     613//         } while ( j += 2, j < l+k+m );
     614//     }
     615//     return;
     616// }
    637617
    638618
     
    643623    float  wkr, wki, xr, xi, yr, yi;
    644624
    645     ENTER(37);
    646625    m  = n >> 1;
    647626    ks = 2 * nc / m;
     
    664643        kk       += ks;
    665644    } while ( j += 2, j < m );
    666     LEAVE(37);
    667645    return;
    668646}
Note: See TracChangeset for help on using the changeset viewer.