I have some questions related to Schauder theory for PDEs, as in, given that $F\in C^{0,\alpha}(\Omega)$ and $A\in C^{0,\alpha}(\Omega)$, the solution $u$ to the pde \begin{equation} \begin{cases} -\nabla\cdot(A(x)\nabla u)=-\nabla\cdot F&\text{ in }\Omega\\ u=0 &\mbox{ on } \partial \Omega \end{cases} \end{equation} is in $C^{1,\alpha}(\omega)$ for $\omega\Subset\Omega$ where $\Omega\subset\mathbb{R}^n$. My questions are:
The standard perturbation method proves $C^{0,\alpha}$ estimates for the gradient of u using the Campanato seminorm. In showing $u\in C^{1,\alpha}$, would one not first need to prove that $u$ is $C^\alpha$ before being able to talk about the gradient, or does this follow from holder regularity of the gradient? I remember seeing counterexamples for this sort of thing but I can't recall them.
The solution u is compared to the solution of equation with constant coefficients and boundary condition $u$ on a ball $B_R$ with radius $R$. This involves using the known decay estimates on these functions that look like \begin{equation} \int_{B_r}|u(y)|^2\,dy \leq C (r/R)^n \int_{B_R}|u(y)|^2\,dy. \end{equation} My question is why is this estimate true up to $B_R$? The proofs that I saw stop at a set smaller than the set on which the PDE is posed. Would someone know about this?