Output gap estimation

It is about the initialization of the Kalman filter recursion, see e.g. [Estimation with the diffuse_filter option)