2 hnrmod - compute mod h * 2^n +r
26 This builtin in faster than the standard mod in that is makes use
27 of shifts and additions when h == 1. When h > 1, a division by h
31 > print hnrmod(2^177-1, 1, 177, -1), hnrmod(10^40, 17, 51, 1)
40 void zhnrmod(ZVALUE v, ZVALUE h, ZVALUE zn, ZVALUE zr, ZVALUE *res)