I am reading the paper of [Convergence of the Yamabe flow for arbitrary initial energy][1]


  [1]: https://projecteuclid.org/journals/journal-of-differential-geometry/volume-69/issue-2/Convergence-of-the-Yamabe-flow-for-arbitrary-initial-energy/10.4310/jdg/1121449107.full

I am stuck by one part of the paper. Suppose $u_\infty>0$ is a smooth function on $(M, g_0)$ and 
$$L_0=\frac{4(n-1)}{n-2}\Delta_{g_0}-R_{g_0}$$
is the conformal laplacian, where $R_{g_0}$ is the scalar curvature. $u_\infty$ satisfies
$$L_0u_\infty+r_\infty u_\infty^{\frac{n+2}{n-2}}=0$$ 
consider the linear operator
$$\psi\mapsto T\psi=u_\infty^{-\frac{4}{n-2}}L_0\psi$$
which is symmetric with respect to the inner product
$$(\psi_1,\psi_2)_\infty=\int_M u_\infty^{\frac{4}{n-2}}\psi_1\psi_2d\mu_0$$
on $L^2(M)$. By the spectral theorem, eigenvalue and eigenvector of $T$, suppose they are $\{-\lambda_i,i\in\mathbb{N}\}$ and $\{\psi_i, i\in\mathbb{N}\}$, satisfies
$$L_0\psi_i+\lambda_i u_\infty^{\frac{4}{n-2}}\psi_i=0$$
$$\int_M u_\infty^{\frac{4}{n-2}}\psi_i\psi_jd\mu_0=\begin{cases}1\quad i=j\\0\quad i\neq j\end{cases}$$
$$\lambda_i\to\infty \text{ as } i\to \infty$$
Let $A=\{i|\lambda_i\leq \frac{n+2}{n-2}r_\infty\}$. He constructs the projection operator $\Pi$ which is
$$\Pi f=\sum_{i\not\in A}\left(\int_M \psi_ifd\mu_0\right)u_\infty^{\frac{4}{n-2}}\psi_i=f-\sum_{i\in A}\left(\int_M \psi_ifd\mu_0\right)u_\infty^{\frac{4}{n-2}}\psi_i$$
He said that by using the **implicit function theorem** one can get the following significance: for every vecrtor $z\in \mathbb{R}^A$ sufficiently small, there exists a smooth function $\bar{u}_z$ such that
$$\int_M u^{\frac{4}{n-2}}(\bar{u}_z-u_\infty)\psi_id\mu_0=z_i \quad \forall\,i\in A$$
$$\Pi\left(L_0\bar{u}_z+r_\infty\bar{u}_z^\frac{n+2}{n-2}\right)=0$$
futhermore the map $z\mapsto \bar{u}_z$ is real analytic.

Here comes my questions, how can we apply and justify the assumption of implicit function theorem? I can not get the motivation of the projection operator either.