code: purgatorio

ref: db1eb844461b07a25ca49117851fa874fd88e065
dir: /libmp/port/mpmod.c/

View raw version
#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);
}