Found additional juice left... ncol(P) and colSums(P) were inside the
gradient computation... removed a few lines, combined, debugged. Also
cut n.integration to 500 from 1000 (halves execution time, appears to
be sufficient for very long tailed distributions)