improve computation of outer product gradient for univariate kalman filters allow to use hessian from optimizer == 5, without computing it with the hessian.m function