I try to answer the question whether spacetime of Schwarzschild vacuum solution is geodesically complete by analyzing its null geodesics. The infinitesimal length element is $$ds^2=(1-\frac{r_s}{r})~c^2 dt^2-(1-\frac{r_s}{r})^{-1} dr^2-r^2 d\Omega^{2}. \tag{1}$$ The solving of corresponding differential equation for $ds=0$ results in $$\pm~ c t=r+r_{s} \log{(r-r_s)}. \tag{2}$$
The "coordinate time" $t$ is coevally affine parameter which value for infalling ($+$) geodesic ranges from $-\infty$ to $+\infty$ for $r\in (+\infty,~ r_{s})$.
As the answers below have proved the statement in bold letters is not true! Coordinate time is a non-affine parameter. However, one can always re-parametrize $t$ to affine parameter $\tau(t)$ (see Appendix in Quasi-geodesics in relativistic gravity) which function is $$\pm~ c \tau=r+2 r_{s} \log{(r-r_s)}-\frac{r_{s}^2}{r-r_s}. \tag{3}$$
For ingoing geodesic the affine parameter $\tau \in (-\infty$,$+\infty$) generates geodesic $r(\tau)\in (+\infty,~ r_{s})$.
Thus, Schwarzschild vacuum solution represents a geodesically complete manifold with null geodesics that never cross the event horizon ($r=r_{s}$) which is not part of that manifold. Do I interpret it correctly (this time)?
Remark
The derived equation (3) is not correct as well because the used equation A.7 from the mentioned reference has a wrong sign (it should be plus instead of minus there). The correct formula for affine parameter is $$\tau\equiv A\int{e^{\nu+\lambda}} dr+B= A r +B \tag{4}.$$
My interpretation based on the at first derived equation (3) was therefore not correct.