loglikelihood_hessian_diag_bicm function

Computes the loglikelihood hessian for the bicm function