tiny change to gammamixture. I tried vectorizing the loop over k, and the performance went down slightly. i think openmp over is possible though