1 #include "os.h"
2 #include <mp.h>
3 #include "dat.h"
4 
5 // remainder = b mod m
6 //
7 // knuth, vol 2, pp 398-400
8 
9 void
10 mpmod(mpint *b, mpint *m, mpint *remainder)
11 {
12 	mpdiv(b, m, nil, remainder);
13 	if(remainder->sign < 0)
14 		mpadd(m, remainder, remainder);
15 }