  Splitting functions are universal functions describing the collinear dynamics
of gauge theories, and as such are crucial ingredients for a wide variety of
calculations in perturbative QCD. We present analytic results for the triple
collinear splitting functions in QCD with a single massive parton. We derive
the splitting functions using two distinct methods; first by expanding the
squared matrix elements in the collinear limit, and secondly by using
soft-collinear effective theory with massive quarks. We find agreement between
these two approaches, providing a strong check of our results. Additionally, we
also check all iterated and soft limits of our results, finding agreement with
predictions from factorization. Our results provide an important ingredient for
higher order perturbative calculations involving massive partons, and for the
description of the collinear dynamics of heavy flavor jets.
