File Coverage

blib/lib/Math/MVPoly/Integer.pm

Criterion	Covered	Total	%
statement	18	21	85.7
branch			n/a
condition			n/a
subroutine	3	3	100.0
pod	1	1	100.0
total	22	25	88.0

package Math::MVPoly::Integer;

# Copyright (c) 1998 by Brian Guarraci. All rights reserved.
# This program is free software; you can redistribute it and/or modify it
# under the same terms as Perl itself.

use strict;
use Math::MVPoly::Monomial;

sub
gcd
{
	my $u;
	my $v;
	my $q;
	my $t;
	my $g;

	my $a = $_[0];
	my $b = $_[1];

	$u = [1,0,abs($a)];
	$v = [0,1,abs($b)];

	while($$v[2])
	{
		$q = int($$u[2]/$$v[2]);
		$t = [  $$u[0] - $$v[0]*$q,
			$$u[1] - $$v[1]*$q,
			$$u[2] - $$v[2]*$q];
		$u = $v;
		$v = $t;
	}

	$g = $$u[2];

	return $g;
}

1;

__END__