Abstract: We consider trace functions generated by the functional calculus of a real function and obtain a complete characterization of the class of geodesically convex trace functions of this form. We then apply these results to the theory of multivariate operator means. Of particular interest is the calculation of the barycenter of relative entropy for a cluster of cubic states.