File: mpinvert.c

package info (click to toggle)
drawterm 20091003-1
  • links: PTS
  • area: main
  • in suites: squeeze
  • size: 2,920 kB
  • ctags: 5,505
  • sloc: ansic: 55,149; makefile: 566; asm: 20
file content (21 lines) | stat: -rw-r--r-- 394 bytes parent folder | download | duplicates (4) 
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
 
#include "os.h"
#include <mp.h>

#define iseven(a)	(((a)->p[0] & 1) == 0)

// use extended gcd to find the multiplicative inverse
// res = b**-1 mod m
void
mpinvert(mpint *b, mpint *m, mpint *res)
{
	mpint *dc1, *dc2;	// don't care

	dc1 = mpnew(0);
	dc2 = mpnew(0);
	mpextendedgcd(b, m, dc1, res, dc2);
	if(mpcmp(dc1, mpone) != 0)
		abort();
	mpmod(res, m, res);
	mpfree(dc1);
	mpfree(dc2);
}