The result of divide(a,b) is the unique element q such that a = q.b + r, with 0 <= r < b. Here b is assumed to be > 0, the sign of a may be arbitrary. For instance, divide(-1,2) should be -1, so that -1 = -1.2 + 1. On my machine, -1/2 is 0 (which is the other value accepted by the C standard.) [Note that the correct symmetry to apply to |a|, one that maps classes with the same quotient to each other, is not 
but 
, where the latter value can be conveniently written as |~a| in C or C++. Amazingly Fokko's incorrect original expresion |-(-a/b -1)| never did any harm. MvL]
Definition at line 48 of file intutils_def.h.
Referenced by atlas::smithnormal::columnReduce(), atlas::matrix::columnReduce(), atlas::smithnormal::rowReduce(), and atlas::matrix::rowReduce().
