alexbrandmeyer@609: //
alexbrandmeyer@609: //  car_params.cc
alexbrandmeyer@609: //  CARFAC Open Source C++ Library
alexbrandmeyer@609: //
alexbrandmeyer@609: //  Created by Alex Brandmeyer on 5/10/13.
alexbrandmeyer@609: //
alexbrandmeyer@609: // This C++ file is part of an implementation of Lyon's cochlear model:
alexbrandmeyer@609: // "Cascade of Asymmetric Resonators with Fast-Acting Compression"
alexbrandmeyer@609: // to supplement Lyon's upcoming book "Human and Machine Hearing"
alexbrandmeyer@609: //
alexbrandmeyer@609: // Licensed under the Apache License, Version 2.0 (the "License");
alexbrandmeyer@609: // you may not use this file except in compliance with the License.
alexbrandmeyer@609: // You may obtain a copy of the License at
alexbrandmeyer@609: //
alexbrandmeyer@609: //     http://www.apache.org/licenses/LICENSE-2.0
alexbrandmeyer@609: //
alexbrandmeyer@609: // Unless required by applicable law or agreed to in writing, software
alexbrandmeyer@609: // distributed under the License is distributed on an "AS IS" BASIS,
alexbrandmeyer@609: // WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied.
alexbrandmeyer@609: // See the License for the specific language governing permissions and
alexbrandmeyer@609: // limitations under the License.
alexbrandmeyer@609: 
alexbrandmeyer@609: #include "car_params.h"
alexbrandmeyer@609: 
alexbrandmeyer@609: CARParams::CARParams(){
alexbrandmeyer@609:   velocity_scale_ = 0.1; //for the velocity nonlinearity
alexbrandmeyer@609:   v_offset_ = 0.04; //offset gives quadratic part
alexbrandmeyer@609:   min_zeta_ = 0.1; //minimum damping factor in mid-freq channels
alexbrandmeyer@609:   max_zeta_ = 0.35; //maximum damping factor in mid-freq channels
alexbrandmeyer@609:   first_pole_theta_ = 0.85 * PI;
alexbrandmeyer@609:   zero_ratio_ = 1.4142; //how far zero is above pole
alexbrandmeyer@609:   high_f_damping_compression_ = 0.5; //0 to 1 to compress theta
alexbrandmeyer@609:   erb_per_step_ = 0.5; //assume G&M's ERB formula
alexbrandmeyer@609:   min_pole_hz_ = 30;
alexbrandmeyer@609:   erb_break_freq_ = 165.3; //Greenwood map's break frequency
alexbrandmeyer@609:   erb_q_ = 1000/(24.7*4.37);
alexbrandmeyer@609: };
alexbrandmeyer@609: 
alexbrandmeyer@609: //This method has been created for debugging purposes and depends on <iostream>.
alexbrandmeyer@609: //Could possibly be removed in the final version to reduce dependencies.
alexbrandmeyer@609: void CARParams::OutputParams(){
alexbrandmeyer@609:   std::cout << "CARParams Values" << std::endl;
alexbrandmeyer@609:   std::cout << "****************" << std::endl;
alexbrandmeyer@609:   std::cout << "velocity_scale_ = " << velocity_scale_ << std::endl;
alexbrandmeyer@609:   std::cout << "v_offset_ = " << v_offset_ << std::endl;
alexbrandmeyer@609:   std::cout << "min_zeta_ = " << min_zeta_ << std::endl;
alexbrandmeyer@609:   std::cout << "max_zeta_ = " << max_zeta_ << std::endl;
alexbrandmeyer@609:   std::cout << "first_pole_theta_ = " << first_pole_theta_ << std::endl;
alexbrandmeyer@609:   std::cout << "zero_ratio_ = " << zero_ratio_ << std::endl;
alexbrandmeyer@609:   std::cout << "high_f_damping_compression_ = " << high_f_damping_compression_
alexbrandmeyer@609:             << std::endl;
alexbrandmeyer@609:   std::cout << "erb_per_step_ = " << erb_per_step_ << std::endl;
alexbrandmeyer@609:   std::cout << "min_pole_hz_ = " << min_pole_hz_ << std::endl;
alexbrandmeyer@609:   std::cout << "erb_break_freq_ = " << erb_break_freq_ << std::endl;
alexbrandmeyer@609:   std::cout << "erb_q_ = " << erb_q_ << std::endl << std::endl;
alexbrandmeyer@609: }
alexbrandmeyer@609: 
alexbrandmeyer@609: void CARParams::SetParams(FPType vs, FPType voff, FPType min_z, FPType max_z,
alexbrandmeyer@609:                           FPType fpt, FPType zr, FPType hfdc, FPType erbps,
alexbrandmeyer@609:                           FPType mph, FPType erbbf, FPType erbq) {  
alexbrandmeyer@609:   velocity_scale_ = vs;
alexbrandmeyer@609:   v_offset_ = voff;
alexbrandmeyer@609:   min_zeta_ = min_z;
alexbrandmeyer@609:   max_zeta_ = max_z;
alexbrandmeyer@609:   first_pole_theta_ = fpt;
alexbrandmeyer@609:   zero_ratio_ = zr;
alexbrandmeyer@609:   high_f_damping_compression_ = hfdc;
alexbrandmeyer@609:   erb_per_step_ = erbps;
alexbrandmeyer@609:   min_pole_hz_ = mph;
alexbrandmeyer@609:   erb_break_freq_ = erbbf;
alexbrandmeyer@609:   erb_q_ = erbq;
alexbrandmeyer@609: };