Stan 's developers in particular use it for hierarchical models, but it general anything with highly correlated parameters works better with HMC than MCMC, IIRC.
Michael Betancourt (à stan dev working on the HMC parts) has a pair of YouTube videos which go into details.
That said, I switched to pymc3 so that I could compute logp via opencl more easily and if there were better ways to do this, I'm happy to see them.
How's that? The algorithms have the same goal - find the posterior distribution. The time to get to that distribution is what is important and what is compared, provided that both algorithms get proper results. How they do it underneath is irrelevant for the user who waits.
That's like saying that comparing a horse cart and an automobile is comparing apples and oranges.
That being said, there are other things where Stan might fare better. User familiarity, or maturity, or personal taste...
Michael Betancourt (à stan dev working on the HMC parts) has a pair of YouTube videos which go into details.
That said, I switched to pymc3 so that I could compute logp via opencl more easily and if there were better ways to do this, I'm happy to see them.