ref: 072c45d4633d5f03e3b79f9c60c3655a6d4a1149
dir: /sys/src/libmp/port/mpmod.c/
#include "os.h" #include <mp.h> #include "dat.h" // remainder = b mod m // // knuth, vol 2, pp 398-400 void mpmod(mpint *b, mpint *m, mpint *remainder) { mpdiv(b, m, nil, remainder); if(remainder->sign < 0) mpadd(m, remainder, remainder); }