Computing the aggregate loss distribution based on numerical inversion of the compound empirical characteristic function of frequency and severity