In response to the questioner's comment: one standard viewpoint is that distributions (or tempered distributions, similarly) are "weak-dual-topology" limits of test functions (or of Schwartz functions), where the weak-dual topology (also called "weak-star topology) is given by semi-norms $\mu_\varphi(u)=|u(\varphi)|$ for test functions (or Schwartz) $\varphi$, at first for other test functions (or Schwartz $\varphi$) and then on the (quasi-) completion, by continuity. Thus, in this viewpoint, the integral representation of $\delta$ is not a literal integral, but is an extension-by-continuity of literal integrals giving approximations to $\delta$. I think this viewpoint, of having integrals not being literal, but extensions-by-continuity, may be most consonant with a common physics-y viewpoint... and does (to my mind) convey a more natural sense of "generalized functions/distributions", not as "functionals", but as an extended notion of function.
Of course, the second viewpoint, perhaps more standard in mathematics per se, is that distributions are continuous linear functionals, and that the integrals written about don't make sense... again, at least not as literal integrals, but perhaps more... haughtily said? That is, the same expressions are obtained, but if/because one is insisting on literal integrals, implicitly disallowing writing integral signs for extensions-by-continuity, some new symbols must be made up. (As you can tell, I am not deeply moved by this viewpoint.)
Either way, the underlying mathematics and mechanisms are the same. That is, e.g., for each Schwartz function $\varphi$, the integrals against the truncated Fourier inversion expression for $\delta$ go to $\varphi(0)$ as the cut-off goes to infinity.
The non-uniformity issue/feature is that there is (provably) no estimate for this limit-approach uniform in $\varphi$. It depends, as it turns out, on $\varphi$.
Is this to the point?
EDIT: I forgot to make explicit the point that, whether or not the integrals are pointwise/numerically valued, literal-or-not, etc., the operator/functional indicated by the integral (as extended by continuity, say) has certain expected properties, so that the usual manipulations which would be correct for numerical/literal integrals are also correct for the non-literal (extended-by-continuity) maps written as integrals... because they extend those integrals, and have the same properties... apart, perhaps, from non-critical numerical/pointwise evaluability.
This is why many physics-y heuristics succeed, after all: what is written is not what is literally meant (from the most exclusionary viewpoint), and it is "still ok" because things extend by continuity, with the same correct transformation properties of whatever-it-is operator denoted by "integral" (and extended from a literal integral).
Yes, on good days, we have our cake and eat it, too. This is (to my mind) the marvelous-ness of mathematics, both internally and in applications.
True, on a bad day bad things happen which will be inexplicable to the (naive?) optimist. Yet, crazily enough, even those bad things can often be repurposed to interesting ends. That's the second miracle in this business, I think.
Edit-edit: following up on the edit to the original question, yes, in that example, for $f$ represented by its Taylor-Maclaurin series in a fixed interval $[-\varepsilon_o,\varepsilon_o]$, as $\varepsilon\to 0$ we can estimate the rate of approach to $f(0)$ uniformly.
But if we look at the aggregate of the classes of such functions over all $\varepsilon_o>0$, this uniform-ness is lost.
Also, although I'd not advocate making this one's first worry, there is an old theorem of E. Borel that says that an arbitrary sequence of real (or complex) numbers can be the Taylor-Maclaurin coefficients at $0$ of a smooth function (obviously, then, not necessarily represented by the infinite series, but the finite-series-with-error-term still works...). This shows, strikingly, that there's no way to hope for true uniformity.
Also, as a sort of paraphrase of the Borel resuls: in the classification of distributions supported at $0$ (essentially proven by Taylor-Maclaurin expansions), we find that they are only the finite linear combinations of $\delta$ and its derivatives. Thus, there is no sequence $c_n$ such that $\sum_{n\ge 1} c_n\cdot \delta^{(n)}$ converges as a distribution.
Best Answer
In terms of the modulus function $M_\nu (z)$ and phase function $\theta _\nu (z)$, we can write $$ H_\nu ^{(1)} (z) = M_\nu (z)\mathrm{e}^{\mathrm{i}\theta _\nu (z)}, \quad J_\nu (z) = M_\nu (z)\cos \theta _\nu (z). $$ It is known that \begin{align*} M_{\nu}^{2}\left(z\right)\sim \;&\frac{2}{\pi z}\left(1 +\frac{1}{2}\frac{\mu-1}{(2z)^{2}}+\frac{1\cdot 3}{2\cdot 4}\frac{(\mu-1)(\mu-9)}{(2z)^{4}}\right.\\ & \qquad\qquad\qquad\quad\;\;\,\left.+\frac{1 \cdot 3\cdot 5}{2\cdot 4\cdot 6}\frac{(\mu-1)(\mu-9)(\mu-25)}{(2z)^{6}}+\ldots \right) \end{align*} and \begin{align*} \theta_{\nu}\left(z\right)\sim z-\left(\frac{1}{2}\nu+\frac{1}{4}\right)\pi & + \frac{\mu-1}{2(4z)}+\frac{(\mu-1)(\mu-25)}{6(4z)^{3}}+\frac{(\mu-1)(\mu^{2}-11 4\mu+1073)}{5(4z)^{5}}\\& +\frac{(\mu-1)(5\mu^{3}-1535\mu^{2}+54703\mu-3 75733)}{ 14(4z)^{7}}+\ldots, \end{align*} as $z\to \infty$ in the sector $|\arg z|\leq\pi-\delta<\pi$, with $\mu =4\nu^2$. It is seen that in order for the terms to decay for large $z$, we need only $\nu =o(z)$. See Section 5.11.4 in Y. L. Luke's book The Special Functions and their Approximations, Vol. 1 and Section 13.75 in G. N. Watson's book A Treatise on the Theory of Bessel Functions.
You may also try using the known uniform asymptotic expansions.