In "Lecture Notes in Quantum Chemistry II" from the European Summer School in Quantum Chemistry, the authors note that when discussing the LYP functional and other GGA functionals, that
Thus one can immediately see that an evaluation of $v_{xc}$ demands an evaluation fo the second derivatives of basis functions. (pp102).
Indeed, when (for a functional like B88) I start with the functional itself and compute the functional derivative (either by hand or with a CAS), I arrive at a (messy) expression that involves second derivatives with respect to space. However, libraries like libxc only need scalar functions of the gradient to be passed, in particular $\nabla \rho \cdot \nabla \rho$ , in order to compute $v_{xc}$. Now libxc contains the expressions for its functionals in Maple and then performs symbolic differentiation on them and then uses the maple-to-C translator; thus the actual C code in libxc is uninterpretable by a human.
Is there a known expression for $v_{xc}$ for B88 or another GGA-like functional that might give me a hint as to how to simplify the functional derivative to only depend on these parameters?