Describir: Inexact primal-dual active set iteration for optimal distribution control of stationary heat or cold source