comparison src/CAR.cpp @ 489:20de0b60b694 carfac_cpp

Added the initialization of g0, but parting from the MATLAB way of doing it. (temporary measure)

           (CARFAC::ERB_Hz(pole_freqs[i])/pole_freqs[i]-car_params->min_zeta_));
       r1_coeffs_[i] = 1-zr_coeffs_[i]*min_zeta_mod;
 
       h_coeffs_[i] = c0_coeffs_[i] * f;
 
-      g0_coeffs_[i] = 0; // TODO: matlab design a bit hacky - think for bit
+      //TODO: g0_coeffs_ calculated here for now. Let's talk about this, need
+      // to see the whole picture of what makes sense.
+      // Related: I believe CAR, IHC and AGC classes would be beneficial!
+      float tmp = 1 - 2*r1_coeffs_[i]*a0_coeffs_[i] + r1_coeffs_[i]*r1_coeffs_[i];
+      g0_coeffs_[i] = tmp / ( tmp + h_coeffs_[i]*r1_coeffs_[i]*c0_coeffs_[i]);
   }
 
 }
 
 CAR_coefficients::~CAR_coefficients(){