Pretty much everywhere I look it is stated that the full two point Green function (let's say for the Klein-Gordon field) is a geometric series in the one particle irreducible diagrams, ie. in momentum space,
$$G(k) = G_0(k)+G_0(k)\Sigma(k)G_0(k)+G_0(k)\Sigma(k)G_0(k)\Sigma(k)G_0(k) + \dots $$ $$ = G_0(k)\big(1+\Sigma(k)G_0(k)+\Sigma(k)G_0(k)\Sigma(k)G_0(k)+\dots\big)$$
and that the sum of this is, using the geometric series formula, $$G(k)=\frac{G_0(k)}{1-\Sigma(k)G_0(k)}=\frac{1}{k^2+m^2-\Sigma(k)-i\epsilon}$$
(for instance on the first page on the top right hand corner of here, or the bottom/top of pages 56/57 here).
However, what is the justification for summing this geometric series in that way? It never seems to be justified and it doesn't seem that $|\Sigma(k)G_0(k)|<1\,.$ Even if I do renormalize the mass $m$ so that $m_R^2 := m^2-\Sigma(k)$ is finite, if the geometric series summation isn't justified the infinities won't cancel out and everything will still diverge. Is there a resummation going on implicitly? Is this step completely nonperturbative?
$\mathbf{Edit\;1}$: Basically as I see it, the situation is this: I ask what the amplitude for the propagator is (in momentum space), and you say $$\frac{G_0(k)}{1-\Sigma(k)G_0(k)}=\frac{1}{k^2+m_R^2-i\epsilon}$$ where $m_R^2=m^2-\Sigma(k)$ is the (finite) mass. Then I notice that you can calculate it approximately using the perturbation series given by the theory, which is $$G_0(k)\big(1+\Sigma(k)G_0(k)+\Sigma(k)G_0(k)\Sigma(k)G_0(k)+\dots\big)\,.$$ However I notice that $\Sigma(k)G_0(k)>1$ and I know that for $x>1$ the perturbation series given by $$1+x+x^2+\cdots$$ is not a good approximation to $\frac{1}{1-x}\,,$ at any order of the perturbation series. So it would seem that perturbation theory fails because it doesn't give a good approximation, at any order, for the propagator. So then I ask what is the justification for this whole thing in the first place? Are they basically just seeing that the perturbation series for the propagator is the same as the perturbation series for $\frac{1}{1-x}\,,$ but realizing that $|x|>1$ and then assuming that the actual propagator must be $\frac{1}{1-x}\,?$ Because this would seem to be a completely nonperturbative step.
$\mathbf{Edit\;2}$: I am going to do the actual computation which is confusing me, and if someone can point out where it is wrong (if it is) that would be a great help. By the way I did a Wick rotation previously, but I am not going to do that this time:
We have that $$\frac{G_0(k)}{1-\Sigma(k)G_0(K)}=\frac{1}{k^2-m_R^2}$$
where $m_R$ is the physical (and finite) mass, and where $G_0(k)=\frac{1}{k^2-m^2}$ where $m$ is the cut-off dependent bare mass. Rearranging, I get $$\Sigma(k)G_0(k)=1-G_0(k)(k^2-m_R^2)=1-\frac{k^2-m_R^2}{k^2-m^2}\,.$$ Let's look at the region where $m_R^2<k^2<m^2\,.$ $m^2\to\infty$ as the cut-off is removed, so this is a very large region. Now it looks to me like that in the above expression for $\Sigma(k)G_0(k)\,,$ that the right hand side is ALWAYS greater than one, and even possibly near infinity for certain values of $k$ (or maybe it doesn't obtain values near infinity because $k$ has to be restricted to values below the cut-off, but this isn't really important). This would seem to make the whole geometric summation suspicious. Is there a mistake?
If I have a basic misunderstanding of how this works then I'd like to know, this is the most basic case of renormalization in QFT but I don't understand it.