I know that Paul Malliavin was heavily influenced by ideas from differential geometry while developing his calculus on Wiener space. But what are the concrete analogies between both areas of mathematics? What has this to do with Hörmander's theorem (if so)?

I can't speak for Paul Malliavin's influences, but I do know a bit about Hormander's theorem (by no means, an expert), and it is naturally suited to differentiable manifolds involving largely the idea of pullbacks of vectors. Malliavin calculus was apparently initiated to give a probabilistic proof of Hormander's theorem. Following Rogers and Williams, Hormander's theorem concerns stochastic differential equations of the form $$ \partial X = \sum_q U_q(X) \partial B^q + W(X)\partial t. $$ Here, X is a stochastic process taking values in R^{n}, U_{q} and W are smooth vector fields, B^{q} are Brownian motions and ∂ represents the Stratonovich integral. As Stratonovich integration satisfies the standard change of variables formula, this SDE makes sense on an arbitrary differentiable manifold. Next, according to the statement of Hormander's theorem, let [.,.] be the usual Lie Bracket for vector fields, and let A_{0},A_{1},... be the sequence of Lie algebras defined as follows. $$ \begin{align} & A_0={\rm Lie}(U_1,U_2,...),\\ & A_k={\rm Lie}([U,W]\colon U\in A_{k1}) \end{align} $$ Then Hormander's theorem states that if ∪_{n}A_{n} spans the tangent space at each point of R^{n}, then X has smooth transition densities. Hormander's theorem is naturally a statement concerning diffusions on differentiable manifolds, as everything I said above makes perfect sense, and is true, if R^{n} is replaced by any differentiable manifold. The idea behind the Malliavin proof is to consider differentiating with respect to perturbations of the Brownian motions B^{q}. The point is, that if B^{q} has a small bump applied at time t, then this creates a small bump in the solution X proportional to U_{q} at this time, which will then be propagated along the solution. In fact, solutions to the original SDE with all different points give rise to stochastically moving frames, and bumps in the solution X are transported along with these frames in a similar as vector fields give rise to transport of vectors along these fields. The solution for X is smooth with respect to smooth bumps in the Brownian motions, which can be shown by converting bumps in the Brownian motion into changes in the probability measure, using Girsanov transforms. So, according to Malliavin calculus, you can always differentiate the solution with respect to the Brownian motions. The idea behind the proof of Hormanders theorem, is to invert this process of varying the Brownian motion→bumps in the final position of X. To do this, it is necessary to invert the process of transporting along the moving frames. That is, it must be a 11 map on the tangent spaces. Then, by a stochastic "pullback" on the manifold (or R^{n}), you can interpret differentiating the solution with respect to its position at any time in terms of differentiating with respect to the Brownian motions. So, this method of proving Hormander's theorem requires you to be able to differentiate with respect to the Brownian motion (i.e., Malliavin calculus) and the rest is all differential geometry (i.e., pullbacks of tangent vectors). 


Some of Malliavin's ideas are wellexplained in his own book: Stochastic Analysis. I think the main geometric idea behind his proof of Hormander's theorem is the idea of submersion. More precisely, if we consider a stochastic differential equation in Stratonovitch form $dX_t=V_0(X_t) dt +\sum_{i=1}^n V_i(X_t) \circ dW^i_t $ where $W$ is a Wiener process and the $V_i$'s vector fields on $R^n$ (or a manifold), we can see the solution $X_t$ as function of the path $\{ W_s , s \le t\}$. We thus have a map $\Pi$ from the Wiener space $\mathbb{W}$, which is the path space of $W$ and $\mathbb{R}^n$. A key insight of Malliavin is now to see the Wiener space $\mathbb{W}$ as an infinitedimensional manifold with tangent space the CameronMartin space (space of functions with derivatives in $L^2$) and to realize that Hormander's conditions make the map $\Pi$ a submersion. As a consequence the Wiener measure is mapped into a smooth measure, which means that if Hormander's conditions are satisfied then the distribution of $X_t$ has a smooth density. It implies the hypoellipticity of the generator of the Markov process $(X_t)_{t \ge 0}$ which is $L=V_0+\frac{1}{2} \sum_{i=1}^n V_i^2$. Of course, this has to be made rigorous and some care has to be given because we have to use Sobolev type derivatives instead of Frechet derivatives. This idea that there is a "natural" manifold structure on path spaces is quite fruitful and besides Malliavin's proof of Hormander's theorem explains and leads to several interesting properties of the operator $V_0+\frac{1}{2} \sum_{i=1}^n V_i^2$. 

