We apply nonsmooth analysis to a well-known optical inverse problem, phase retrieval. The phase retrieval problem arises in many different modalities of electromagnetic imaging and has been studied in the optics literature for over forty years. The state of the art for this problem in two dimensions involves iterated projections for solving a nonconvex feasibility problem. Despite widespread use of these algorithms, current mathematical theory cannot explain their success. At the heart of projection algorithms is a nonconvex, nonsmooth optimization problem. We obtain some insight into these algorithms by applying techniques from nonsmooth analysis. In particular, we show that the weak closure of the set of directions toward the projection generate the subdifferential of the corresponding squared set distance function. Following a pattern of proof described in [F. H. Clarke, Yu. S. Ledyaev, R. J. Stern, and P. R. Wolenski, Nonsmooth Analysis and Control Theory, Springer-Verlag, New York, 1998], this result is generalized to provide conditions under which the subdifferential of an integral function equals the integral of the subdifferential.