9512fe850e98fdd448c638ca63fdd92a8a510255ahl/* Copyright (c) 1984, 1986, 1987, 1988, 1989 AT&T */
9512fe850e98fdd448c638ca63fdd92a8a510255ahl/* All Rights Reserved */
9512fe850e98fdd448c638ca63fdd92a8a510255ahl * Copyright (c) 1980 Regents of the University of California.
9512fe850e98fdd448c638ca63fdd92a8a510255ahl * All rights reserved. The Berkeley software License Agreement
9512fe850e98fdd448c638ca63fdd92a8a510255ahl * specifies the terms and conditions for redistribution.
9512fe850e98fdd448c638ca63fdd92a8a510255ahl/* Portions Copyright(c) 1988, Sun Microsystems Inc. */
9512fe850e98fdd448c638ca63fdd92a8a510255ahl/* All Rights Reserved */
9512fe850e98fdd448c638ca63fdd92a8a510255ahl * Copyright (c) 1997, by Sun Microsystems, Inc.
9512fe850e98fdd448c638ca63fdd92a8a510255ahl * All rights reserved.
9512fe850e98fdd448c638ca63fdd92a8a510255ahl/* LINTLIBRARY */
9512fe850e98fdd448c638ca63fdd92a8a510255ahl MINT x, y, z, w;
9512fe850e98fdd448c638ca63fdd92a8a510255ahl while (y.len != 0) {
9512fe850e98fdd448c638ca63fdd92a8a510255ahl mp_mdiv(&x, &y, &w, &z);
9512fe850e98fdd448c638ca63fdd92a8a510255ahl * Minimize calls to allocators. Don't use pointers for local
9512fe850e98fdd448c638ca63fdd92a8a510255ahl * variables, for the one "initialized" multiple precision
9512fe850e98fdd448c638ca63fdd92a8a510255ahl * variable, do it just once.
23b5c241225a8ade2b6b9f06ebb891ee459e3b02tomee /* invariant: x0*u1 + x1*u2 = u3 */
23b5c241225a8ade2b6b9f06ebb891ee459e3b02tomee /* invariant: x0*v1 + x2*v2 = v3 */
23b5c241225a8ade2b6b9f06ebb891ee459e3b02tomee /* invariant: x(n+1) = x(n-1) % x(n) */
9512fe850e98fdd448c638ca63fdd92a8a510255ahl /* now x0*u1 + x1*u2 == 1, therefore, (u2*x1) % x0 == 1 */
_mp_xfree(&q);
_mp_xfree(&r);
_mp_xfree(&t);