# HG changeset patch
# User Ulf.Hammarqvist@gmail.com
# Date 1333791699 0
# Node ID c23860d91135bd21aa55ff0f346c28a6db98b547
# Parent  ff12d0432d9ca9dd6da9a130b0ba8c69ed3c95e1


diff -r ff12d0432d9c -r c23860d91135 branches/carfac_cpp/src/CAR.cpp
--- a/branches/carfac_cpp/src/CAR.cpp	Sat Apr 07 09:35:10 2012 +0000
+++ b/branches/carfac_cpp/src/CAR.cpp	Sat Apr 07 09:41:39 2012 +0000
@@ -8,7 +8,7 @@
 
   float f = pow(car_params->zero_ratio_, 2) + 1;
 
-  // dirty FloatArray initialisation. Redo this later
+  // TODO: dirty FloatArray initialisation. Redo this later
   r1_coeffs_= pole_freqs;
   a0_coeffs_= pole_freqs;
   c0_coeffs_= pole_freqs;
@@ -22,8 +22,8 @@
   float x;
   float ff = car_params->high_f_damping_compression_;
 
-  // temp loop until we get eigen in
-  for(float i=0; i<theta.size(); i++){
+  // TODO: temp loop until we get eigen in  (or we just leave it like this)
+  for(float i=0; i<pole_freqs.size(); i++){
       theta[i] *= (2*kPi*fs); // scalar mult.
       c0_coeffs_[i] = sin(theta[i]);
       a0_coeffs_[i] = cos(theta[i]);
@@ -33,14 +33,13 @@
 
       min_zeta_mod = (car_params->min_zeta_ + 0.25 *
           (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
   }
 
-  // g0_coeffs_[i] = ... TODO: restructure this!
-
 }
 
 CAR_coefficients::~CAR_coefficients(){