17  * along with this program; if not, write to the Free Software
42 	if (mpi_copy(&u, a) < 0)  in mpi_invm()
44 	if (mpi_copy(&v, n) < 0)  in mpi_invm()
48 		if (mpi_rshift(u, u, 1) < 0)  in mpi_invm()
50 		if (mpi_rshift(v, v, 1) < 0)  in mpi_invm()
56 	if (!u1)  in mpi_invm()
58 	if (!odd) {  in mpi_invm()
60 		if (!u2)  in mpi_invm()
63 	if (mpi_copy(&u3, u) < 0)  in mpi_invm()
65 	if (mpi_copy(&v1, v) < 0)  in mpi_invm()
67 	if (!odd) {  in mpi_invm()
69 		if (!v2)  in mpi_invm()
71 		if (mpi_sub(v2, u1, u) < 0)  in mpi_invm()
74 	if (mpi_copy(&v3, v) < 0)  in mpi_invm()
76 	if (mpi_test_bit(u, 0)) {	/* u is odd */  in mpi_invm()
78 		if (!t1)  in mpi_invm()
80 		if (!odd) {  in mpi_invm()
82 			if (!t2)  in mpi_invm()
86 		if (mpi_copy(&t3, v) < 0)  in mpi_invm()
92 		if (!t1)  in mpi_invm()
94 		if (!odd) {  in mpi_invm()
96 			if (!t2)  in mpi_invm()
99 		if (mpi_copy(&t3, u) < 0)  in mpi_invm()
104 			if (!odd) {  in mpi_invm()
105 				if (mpi_test_bit(t1, 0) || mpi_test_bit(t2, 0)) {	/* one is odd */  in mpi_invm()
106 					if (mpi_add(t1, t1, v) < 0)  in mpi_invm()
108 					if (mpi_sub(t2, t2, u) < 0)  in mpi_invm()
111 				if (mpi_rshift(t1, t1, 1) < 0)  in mpi_invm()
113 				if (mpi_rshift(t2, t2, 1) < 0)  in mpi_invm()
115 				if (mpi_rshift(t3, t3, 1) < 0)  in mpi_invm()
118 				if (mpi_test_bit(t1, 0))  in mpi_invm()
119 					if (mpi_add(t1, t1, v) < 0)  in mpi_invm()
121 				if (mpi_rshift(t1, t1, 1) < 0)  in mpi_invm()
123 				if (mpi_rshift(t3, t3, 1) < 0)  in mpi_invm()
130 		if (!t3->sign) {  in mpi_invm()
131 			if (mpi_set(u1, t1) < 0)  in mpi_invm()
133 			if (!odd)  in mpi_invm()
134 				if (mpi_set(u2, t2) < 0)  in mpi_invm()
136 			if (mpi_set(u3, t3) < 0)  in mpi_invm()
139 			if (mpi_sub(v1, v, t1) < 0)  in mpi_invm()
143 			if (!odd)  in mpi_invm()
144 				if (mpi_sub(v2, u, t2) < 0)  in mpi_invm()
149 			if (mpi_set(v3, t3) < 0)  in mpi_invm()
153 		if (mpi_sub(t1, u1, v1) < 0)  in mpi_invm()
155 		if (!odd)  in mpi_invm()
156 			if (mpi_sub(t2, u2, v2) < 0)  in mpi_invm()
158 		if (mpi_sub(t3, u3, v3) < 0)  in mpi_invm()
160 		if (t1->sign) {  in mpi_invm()
161 			if (mpi_add(t1, t1, v) < 0)  in mpi_invm()
163 			if (!odd)  in mpi_invm()
164 				if (mpi_sub(t2, t2, u) < 0)  in mpi_invm()
175 	if (!odd) {  in mpi_invm()