Auxil function

An auxiliary function used in calculating the negative log-likelehood and its gradient