Computing with the COM-Poisson distribution