Mark Dickinson added the comment:
I suspect that this is simply an error in the original code: the docstring
says that mu should be in the range [0, 2*pi), so reducing mu modulo 2*pi makes
little sense. I guess the lines at the end of the method were intended to be
written:
if u3 >= 0.5:
theta = (mu + _acos(f)) % TWOPI
else:
theta = (mu - _acos(f)) % TWOPI
instead of:
if u3 >= 0.5:
theta = (mu % TWOPI) + _acos(f)
else:
theta = (mu % TWOPI) - _acos(f)
That would then give consistent results, at least.
----------
_______________________________________
Python tracker <[email protected]>
<http://bugs.python.org/issue17149>
_______________________________________
_______________________________________________
Python-bugs-list mailing list
Unsubscribe:
http://mail.python.org/mailman/options/python-bugs-list/archive%40mail-archive.com
