Efficient multipole expansion: Choice of order and density partitioning techniques