3f310e7 is a commit that has the specialized reverse mode functions for the explicit laplace distributions. for instance the file here that has reverse mode specializations for neg_binomial_2. Right now we use embedded reverse and higher order autodiff to calculate these from the likelihood function. If we need to speedup laplace one of the first places we should look is adding these back.