PyMC computes the element wise log likelihood (per observation) at the end because this is generally needed for model comparison.
This quantity is not computed at all during NUTS sampling, where a combined scalar log likelihood + log prior (and gradients) is computed for all the parameters at once