new_exc[i] = exc[i] + comb_gain * (
fact += step;  new_exc[i] = exc[i] + comb_gain * fact * (
void comb_filter(float *exc,          /*decoded excitation*/float *new_exc,      /*enhanced excitation*/float *ak,           /*LPC filter coefs*/int p,               /*LPC order*/int nsf,             /*sub-frame size*/int pitch,           /*pitch period*/float *pitch_gain,   /*pitch gain (3-tap)*/float  comb_gain     /*gain of comb filter*/){   int i;   float exc_energy=0, new_exc_energy=0;   float gain;   /*Compute excitation energy prior to enhancement*/   for (i=0;i<nsf;i++)      exc_energy+=exc[i]*exc[i];   /*Apply pitch comb-filter (filter out noise between pitch harmonics)*/   for (i=0;i<nsf;i++)   {      new_exc[i] = exc[i] + comb_gain * (                                         pitch_gain[0]*exc[i-pitch+1] +                                         pitch_gain[1]*exc[i-pitch] +                                         pitch_gain[2]*exc[i-pitch-1]                                         );   }      /*Gain after enhancement*/   for (i=0;i<nsf;i++)      new_exc_energy+=new_exc[i]*new_exc[i];   /*Compute scaling factor and normalize energy*/   gain = sqrt(exc_energy)/sqrt(.1+new_exc_energy);   for (i=0;i<nsf;i++)      new_exc[i] *= gain;}