gumbel_p1_lddd function

Third derivative tensor of the normalized log-likelihood