log(0) warning about arspy

[MFreidank]

Hey @bsxfan,
Thanks for reporting this.
Currently, the only alternative I see would be to cap inputs to log at a small value to avoid the error, but that is guaranteed to affect sample statistics so I'm not sure it's a better solution.

I will double-check that the formula computed is correct (it's been a while since I looked at it) and if yes, I'll probably just add a few lines of documentation and suggest your workaround for such cases.

Update: on a separate note: do you know how the matlab implementation handles this? I don't have access to Matlab right now, but ideally we should handle this case in the same way.

[bsxfan]

@MFreidank

• Sorry, my MATLAB license stayed behind with my previous employer. But as far as I remember from MATLAB, the default behaviour would be to quietly return -inf for log(0). So if the same thing happened in the MATLAB implementation, nobody would have noticed.
• I can mention that the problem seems rare. As you can see in my test code, I had to do 1000 invocations to make it happen. If I do just 10, chances are it does not occur.