  We prove optimal boundary regularity for bounded positive weak solutions of
fast diffusion equations in smooth bounded domains. This solves a problem
raised by Berryman and Holland in 1980 for these equations in the subcritical
and critical regimes. Our proof of the a priori estimates uses a geometric type
structure of the fast diffusion equations, where an important ingredient is an
evolution equation for a curvature-like quantity.
