Regularity of solutions to the Dirichlet problem for fast diffusion equations

We prove global H"older gradient estimates for bounded positive weak solutions of fast diffusion equations in smooth bounded domains with the homogeneous Dirichlet boundary condition, which then lead us to establish their optimal global regularity. This solves a problem raised by Berryman and Holland in 1980.