Ulf@477: #include "CAR.h"
Ulf@484: #include "CARFAC_common.h"
Ulf@487: #include "CARFAC.h"
Ulf@486: #include <cmath>
Ulf@477: 
Ulf@486: CAR_coefficients::CAR_coefficients(CAR_parameters* car_params, float fs,
Ulf@484:                                    FloatArray pole_freqs){
Ulf@486: 
Ulf@486:   float f = pow(car_params->zero_ratio_, 2) + 1;
Ulf@486: 
Ulf@488:   // TODO: dirty FloatArray initialisation. Redo this later
Ulf@487:   r1_coeffs_= pole_freqs;
Ulf@487:   a0_coeffs_= pole_freqs;
Ulf@487:   c0_coeffs_= pole_freqs;
Ulf@487:   h_coeffs_= pole_freqs;
Ulf@487:   g0_coeffs_= pole_freqs;
Ulf@487:   zr_coeffs_= pole_freqs;
Ulf@487: 
Ulf@486:   FloatArray theta = pole_freqs;
Ulf@486: 
Ulf@487:   float min_zeta_mod;
Ulf@487:   float x;
Ulf@487:   float ff = car_params->high_f_damping_compression_;
Ulf@487: 
Ulf@488:   // TODO: temp loop until we get eigen in  (or we just leave it like this)
Ulf@488:   for(float i=0; i<pole_freqs.size(); i++){
Ulf@486:       theta[i] *= (2*kPi*fs); // scalar mult.
Ulf@487:       c0_coeffs_[i] = sin(theta[i]);
Ulf@487:       a0_coeffs_[i] = cos(theta[i]);
Ulf@487: 
Ulf@487:       x = theta[i]/kPi;
Ulf@487:       zr_coeffs_[i] = kPi * (x - ff * x*x*x);
Ulf@487: 
Ulf@487:       min_zeta_mod = (car_params->min_zeta_ + 0.25 *
Ulf@487:           (CARFAC::ERB_Hz(pole_freqs[i])/pole_freqs[i]-car_params->min_zeta_));
Ulf@487:       r1_coeffs_[i] = 1-zr_coeffs_[i]*min_zeta_mod;
Ulf@487: 
Ulf@487:       h_coeffs_[i] = c0_coeffs_[i] * f;
Ulf@488: 
Ulf@488:       g0_coeffs_[i] = 0; // TODO: matlab design a bit hacky - think for bit
Ulf@486:   }
Ulf@486: 
Ulf@477: }
Ulf@477: 
Ulf@484: CAR_coefficients::~CAR_coefficients(){
Ulf@483:   // TODO Auto-generated destructor stub
Ulf@477: }