I think Stan internally decouples the model from the algorithm. In other words, a model is just an implementation of an interface, providing log probability val

I think Stan internally decouples the model from the algorithm. In other words, a model is just an implementation of an interface, providing log probability values (and gradient) for a given proposal. To wit it supports optimization in addition to sampling.