This paper presents an algorithm to accelerate the evaluation of inspiral-merger-ringdown waveform models for gravitational wave data analysis. While the idea can also be applied in the time domain, here we focus on the frequency domain, which is most typically used to reduced computational cost in gravitational wave data analysis. Our work extends the idea of multibanding, which has been developed to accelerate frequency domain waveforms, to include the merger and ringdown and spherical harmonics beyond the dominant quadrupole spherical harmonic. The original method is based on a heuristic algorithm based on the inspiral to de-refine the equi-spaced frequency grid used for data analysis where a coarser grid is sufficient for accurate evaluation of a waveform model. Here we use a different criterion, based on the local interpolation error, which is more flexible and can easily be adapted to general waveforms, if their phenomenology is understood. We discuss our implementation in the LIGO Algorithm Library for the PhenomXHM frequency domain model, and report the acceleration in different parts of the parameter space of compact binary systems.