ref: 059c85dd75f4790625ea3ba19f1c76cf2f3e6b4a
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);
}