A flexible algorithm for calculating pair interactions on SIMD architectures