Ignore:
Timestamp:
12/10/06 19:48:05 (18 years ago)
Author:
r2d
Message:
  • variable name correction
File:
1 edited

Legend:

Unmodified
Added
Removed
  • libmpc/branches/r2d/libmpcenc/quant.c

    r156 r159  
    163163    float  signal = 1.e-30f;
    164164    float  error = 1.e-30f;
    165     float  sig, err;
     165    float  err, sig;
    166166
    167167    // Summation of the absolute power and the quadratic error
    168168    for ( k = 0; k < 36; k++ ) {
    169         err = input[k] * NoiseInjectionCompensation1D [res];
    170                 sig = mpc_nearbyintf(err * fac) * invfac - err;
    171 
    172         error += sig * sig;
    173         signal += err * err;
     169        sig = input[k] * NoiseInjectionCompensation1D [res];
     170                err = mpc_nearbyintf(sig * fac) * invfac - sig;
     171
     172        error += err * err;
     173        signal += sig * sig;
    174174    }
    175175
     
    185185    float  fac    = A [res];
    186186    float  invfac = C [res];
    187         float  signal, error, ret, sig, err;
     187        float  signal, error, ret, err, sig;
    188188
    189189    // Summation of the absolute power and the quadratic error
     
    191191        signal = error = 1.e-30f;
    192192    for ( ; k < 12; k++ ) {
    193         err = input[k] * NoiseInjectionCompensation1D [res];
    194                 sig = mpc_nearbyintf(err * fac) * invfac - err;
    195 
    196         error += sig * sig;
    197         signal += err * err;
    198     }
    199     sig = signal > error  ?  error / (SNRcomp * signal)  :  error / signal;
    200     ret = sig;
     193        sig = input[k] * NoiseInjectionCompensation1D [res];
     194                err = mpc_nearbyintf(sig * fac) * invfac - sig;
     195
     196        error += err * err;
     197        signal += sig * sig;
     198    }
     199    err = signal > error  ?  error / (SNRcomp * signal)  :  error / signal;
     200    ret = err;
    201201    signal = error = 1.e-30f;
    202202    for ( ; k < 24; k++ ) {
    203         err = input[k] * NoiseInjectionCompensation1D [res];
    204                 sig = mpc_nearbyintf(err * fac) * invfac - err;
    205 
    206         error += sig * sig;
    207         signal += err * err;
    208     }
    209     sig = signal > error  ?  error / (SNRcomp * signal)  :  error / signal;
    210         ret = maxf(ret, sig);
     203        sig = input[k] * NoiseInjectionCompensation1D [res];
     204                err = mpc_nearbyintf(sig * fac) * invfac - sig;
     205
     206        error += err * err;
     207        signal += sig * sig;
     208    }
     209    err = signal > error  ?  error / (SNRcomp * signal)  :  error / signal;
     210        ret = maxf(ret, err);
    211211
    212212    signal = error = 1.e-30f;
    213213    for ( ; k < 36; k++ ) {
    214         err = input[k] * NoiseInjectionCompensation1D [res];
    215                 sig = mpc_nearbyintf(err * fac) * invfac - err;
    216 
    217         error += sig * sig;
    218         signal += err * err;
    219     }
    220     sig = signal > error  ?  error / (SNRcomp * signal)  :  error / signal;
    221         ret = maxf(ret, sig);
     214        sig = input[k] * NoiseInjectionCompensation1D [res];
     215                err = mpc_nearbyintf(sig * fac) * invfac - sig;
     216
     217        error += err * err;
     218        signal += sig * sig;
     219    }
     220    err = signal > error  ?  error / (SNRcomp * signal)  :  error / signal;
     221        ret = maxf(ret, err);
    222222
    223223    return ret;
Note: See TracChangeset for help on using the changeset viewer.