Any reason holds Julia from having a proper Pochhammer function?


#1

Overflow? Inaccuracy?


#2

Nobody has needed it and so nobody has implemented it yet? This would be useful for the Combinatorics.jl package. Why don’t you try making a pull request?


#3

I see. Thanks for the encouragement. But I’m just a beginner and dare not write any thing for Julia or any libraries for serious use. :sweat_smile:


#4

It would be a pretty good introduction to how go about contributing, actually!


#5

Good point. I’ll start to work on it though slowly and issue a pull request when ready. Irresistible to make my own contribution to Julia…


#6

You can find an implementation in FastTransforms.jl:


#7

Thanks! Is this something that we can add to Base? I think people would be reluctant to load a whole package just for using a (simple) function.


#8

All the special functions in Base have been moved to SpecialFunctions.jl, so it could be moved there


#9

SpecialFunctions is not the right place for an integer only function, though of course you could extend it via the Gamma function.


#10

It’s easy to add a fallback that calls lgamma (Need to use lgamma Instead of gamma to avoid overflow)