norm_gg function

Second derivative matrix of the expected per-observation log-likelihood