In a phase field Kim-Kim-Suzuki (KKS) model, the variables are the order parameters, global concentration, and phase concentrations defined by the Allen-Cahn equation, Cahn-Hilliard equation, and the KKS local equilibrium constraints, respectively. I am wondering how the phase concentration variables are solved in FEniCS. Are they solved on the mesh nodes or on the Gaussian quadrature points? Thanks.