blob: 89d017c6cedbdff23115e83b5ac21fa2e036191a (
plain)
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
|
/*****
* mod.h
* Andy Hammerlindl 2005/03/16
*
* Definition of implementation-independent mod function.
*****/
#ifndef MOD_H
#define MOD_H
#include <cmath>
#include "common.h"
using std::fmod;
inline Int Mod(Int x, Int y) {return x % y;}
inline double Mod(double x, double y) {return fmod(x,y);}
template<typename T>
inline T portableMod(T x,T y)
{
// Implementation-independent definition of mod; ensure that result has
// same sign as divisor
T val=Mod(x,y);
if((y > 0 && val < 0) || (y < 0 && val > 0)) val += y;
return val;
}
inline Int imod(Int x, Int y)
{
return portableMod<Int>(x,y);
}
inline Int imod(Int i, size_t n) {
i %= (Int) n;
if(i < 0) i += (Int) n;
return i;
}
#endif
|