annotate carfac/carfac_util.cc @ 662:7e18c84ca2b7

Small cleanup of eigen usage in SAI implementation.
author ronw@google.com
date Tue, 16 Jul 2013 19:56:11 +0000
parents 16dfff1de47a
children
rev   line source
alexbrandmeyer@609 1 //
alexbrandmeyer@643 2 // carfac_util.cc
alexbrandmeyer@609 3 // CARFAC Open Source C++ Library
alexbrandmeyer@609 4 //
alexbrandmeyer@643 5 // Created by Alex Brandmeyer on 6/3/13.
alexbrandmeyer@609 6 //
alexbrandmeyer@609 7 // This C++ file is part of an implementation of Lyon's cochlear model:
alexbrandmeyer@609 8 // "Cascade of Asymmetric Resonators with Fast-Acting Compression"
alexbrandmeyer@609 9 // to supplement Lyon's upcoming book "Human and Machine Hearing"
alexbrandmeyer@609 10 //
alexbrandmeyer@609 11 // Licensed under the Apache License, Version 2.0 (the "License");
alexbrandmeyer@609 12 // you may not use this file except in compliance with the License.
alexbrandmeyer@609 13 // You may obtain a copy of the License at
alexbrandmeyer@609 14 //
alexbrandmeyer@609 15 // http://www.apache.org/licenses/LICENSE-2.0
alexbrandmeyer@609 16 //
alexbrandmeyer@609 17 // Unless required by applicable law or agreed to in writing, software
alexbrandmeyer@609 18 // distributed under the License is distributed on an "AS IS" BASIS,
alexbrandmeyer@609 19 // WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied.
alexbrandmeyer@609 20 // See the License for the specific language governing permissions and
alexbrandmeyer@609 21 // limitations under the License.
alexbrandmeyer@609 22
alexbrandmeyer@643 23 #include "carfac_util.h"
alexbrandmeyer@609 24
ronw@644 25 ArrayX CARFACDetect(const ArrayX& x) {
alexbrandmeyer@643 26 ArrayX conductance, z, set;
alexbrandmeyer@609 27 FPType a = 0.175;
alexbrandmeyer@610 28 // This offsets the low-end tail into negative x territory.
alexbrandmeyer@610 29 // The parameter is adjusted for the book, to make the 20% DC response
alexbrandmeyer@610 30 // threshold at 0.1.
alexbrandmeyer@609 31 z = x + a;
alexbrandmeyer@610 32 // Zero is the final answer for many points.
alexbrandmeyer@626 33 conductance = (z < 0).select(0.0, (z*z*z) / (z*z*z + z*z + 0.1));
alexbrandmeyer@609 34 return conductance;
ronw@644 35 }