Use of the log for the computation of the partial derivative of Lauricella. This allows to increase the number of iterations.
The elements of the sum are not any more computed several times.

Fixed: the calculation of the data frame containing the indices for the nested sums was wrong for dimension p>3.