logmultinomial implements log(Combinatorics.multinomial)
It's likely I haven't met proper guidelines for style/documentation, so please verify those when reviewing the pull request.
Is there any other library (numpy, mathematica, etcetera) that implements this function?
Is there any other library (numpy, mathematica, etcetera) that implements this function?
As far as I know, no. Numpy definitely doesn't have it, but numpy doesn't even have multinomial, which I was very glad to find Julia's Combinatorics does have. Multinomials are useful for various discrete probability calculations, and when those calculations are done under a log transform, logmultinomial would be handy. I implemented it because I needed it for game-theory research code where I'm computing expected utilities in large symmetric games, and I think there's a good chance it could be useful to others.