Mercurial > hg > aimc
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 } |