Evaluation of aggregation integrals on nested grids refined toward an arbitrary point. Approximation of density distribution through high order polynomials. Explicit algorithms consisting of manageable amount of code. Evaluation of aggregation source and sink terms in almost linear complexity. Numerical tests leading to (heuristic) strategies for the choice of hp refinement.