I noticed that #253 omitted the normalization step when refactoring the sample() method of the Dirichlet distribution. The sum local variable is not used after it is computed. For large dimensions of alpha, the magnitude of sampled vectors is noticeably shorter than those generated by statrs@0.17.1 or scipy.
I noticed that #253 omitted the normalization step when refactoring the
sample()method of the Dirichlet distribution. Thesumlocal variable is not used after it is computed. For large dimensions of alpha, the magnitude of sampled vectors is noticeably shorter than those generated by statrs@0.17.1 or scipy.