Message181793
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. |
|
Date |
User |
Action |
Args |
2013-02-10 12:19:23 | mark.dickinson | set | recipients:
+ mark.dickinson, rhettinger, serhiy.storchaka |
2013-02-10 12:19:23 | mark.dickinson | set | messageid: <1360498763.41.0.309503113821.issue17149@psf.upfronthosting.co.za> |
2013-02-10 12:19:23 | mark.dickinson | link | issue17149 messages |
2013-02-10 12:19:23 | mark.dickinson | create | |
|