Eu tenho um algoritmo MCMC específico que gostaria de portar para C / C ++. Grande parte da computação dispendiosa já está em C via Cython, mas quero que todo o amostrador seja escrito em uma linguagem compilada, para que eu possa escrever wrappers para Python / R / Matlab / o que for.
Depois de bisbilhotar, estou inclinado para C ++. Conheço algumas bibliotecas relevantes: Armadillo (http://arma.sourceforge.net/) e Scythe (http://scythe.wustl.edu/). Ambos tentam emular alguns aspectos do R / Matlab para facilitar a curva de aprendizado, que eu gosto muito. Foice quadrada um pouco melhor com o que eu quero fazer, eu acho. Em particular, seu RNG inclui muitas distribuições nas quais o Tatu possui apenas uniforme / normal, o que é inconveniente. O tatu parece estar em desenvolvimento bastante ativo, enquanto o Scythe viu seu último lançamento em 2007.
Então, o que eu estou querendo saber é se alguém tem experiência com essas bibliotecas - ou outras que eu quase certamente perdi - e se sim, se há algo para recomendar uma sobre as outras para um estatístico muito familiarizado com Python / R / Matlab mas menos ainda com idiomas compilados (não completamente ignorantes, mas não exatamente proficientes ...).