Fast numerical solution of nonlinear nonlocal cochlear models