alexbrandmeyer@609: //
alexbrandmeyer@643: //  carfac_util.cc
alexbrandmeyer@609: //  CARFAC Open Source C++ Library
alexbrandmeyer@609: //
alexbrandmeyer@643: //  Created by Alex Brandmeyer on 6/3/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@643: #include "carfac_util.h"
alexbrandmeyer@609: 
ronw@644: ArrayX CARFACDetect(const ArrayX& x) {
alexbrandmeyer@643:   ArrayX conductance, z, set;
alexbrandmeyer@609:   FPType a = 0.175;
alexbrandmeyer@610:   // This offsets the low-end tail into negative x territory.
alexbrandmeyer@610:   // The parameter is adjusted for the book, to make the 20% DC response
alexbrandmeyer@610:   // threshold at 0.1.
alexbrandmeyer@609:   z  = x + a;
alexbrandmeyer@610:   // Zero is the final answer for many points.
alexbrandmeyer@626:   conductance = (z < 0).select(0.0, (z*z*z) / (z*z*z + z*z + 0.1));
alexbrandmeyer@609:   return conductance;
ronw@644: }