[Date Prev][Date Next][Thread Prev][Thread Next]
- Subject: RE: Modulo infinity
- From: Dan Tull <dtull@...>
- Date: Mon, 31 Mar 2014 21:46:14 +0000
> > Just check this code in your favorite C compiler:
> > http://codepad.org/ebn79KuY
> That code does not calculate N%inf. It calculates fmod(x,Inf).
> But fmod is not %...
Out of curiosity (apologies if I'm missing something obvious), why isn't % the same as fmod?
#define luai_nummod(a,b) ((a) - floor((a)/(b))*(b))
In some cursory checking, they return results that are almost always equal (or different only at the 12th decimal place or so) aside from this discrepancy around infinity as the second operand. Meanwhile, pow (on the very next line in luaconf.h) just delegates to the C level pow function.
Is there some compatibility or performance issue with fmod? (Note that I'm looking at 5.1, so this predates 5.3's integer handling)
I suppose regardless of the reason, anybody that cares about this case can tweak luaconf.h to swap in fmod for luai_nummod's definition easily enough.
Just curious -- DT