The authors consider the wave equation as a boundary integral equation with a retarded potential. They discretize the equation making use of the convolution quadrature in time and of the Galerkin boundary element method in space variables. As this method implies densely populated matrices, and consequently at least quadratic complexity, they use the panel clustering approximation in order to reduce the complexity. According to this approximation, they carry out upper bounds for both storage requirements and computational complexity.