halfnorm_lddd function

Third derivative tensor of the log-likelihood