TY - JOUR T1 - Fast Evaluation of the Room Transfer Function Using Multipole Expansion JF - Audio, Speech, and Language Processing, IEEE Transactions on Y1 - 2007 A1 - Duraiswami, Ramani A1 - Zotkin,Dmitry N A1 - Gumerov, Nail A. KW - acoustics;computational KW - algorithm;frequency-domain KW - Alien-Berkley KW - analysis;reverberation; KW - complexity;frequency-domain KW - computations;image KW - expansion;receivers;reverberant KW - field;architectural KW - fields;reverberation;room KW - function;acoustic KW - method;image KW - potential;multipole KW - sound KW - source KW - sources;monopole KW - transfer AB - Reverberation in rooms is often simulated with the image method due to Allen and Berkley (1979). This method has an asymptotic complexity that is cubic in terms of the simulated reverberation length. When employed in the frequency domain, it is relatively computationally expensive if there are many receivers in the room or if the source or receiver positions are changing with time. The computational complexity of the image method is due to the repeated summation of the fields generated by a large number of image sources. In this paper, a fast method to perform such summations is presented. The method is based on multipole expansion of the monopole source potential. For offline computation of the room transfer function for N image sources and M receiver points, use of the Allen-Berkley algorithm requires O(NM) operations, whereas use of the proposed method requires only O(N+M) operations, resulting in significantly faster computation of reverberant sound fields. The proposed method also has a considerable speed advantage in situations where the room transfer function must be rapidly updated online in response to source/receiver location changes. Simulation results are presented, and algorithm accuracy, speed, and implementation details are discussed. For problems that require frequency-domain computations, the algorithm is found to generate sound fields identical to the ones obtained with the frequency-domain version of the Allen-Berkley algorithm at a fraction of computational cost VL - 15 SN - 1558-7916 CP - 2 M3 - 10.1109/TASL.2006.876753 ER -