  We present an algorithm for holographic shaping of partially coherent light,
bridging the gap between traditional coherent and geometric optical approaches.
The description of partially coherent light relies on a mode expansion
formalism, with possibly thousands of individual modes, and the inverse-design
optimization algorithm is based on gradient descent with the help of
algorithmic differentiation. We demonstrate numerically and experimentally that
an optical system consisting of two phase patterns can potentially achieve any
intensity profile transformation with good accuracy.
