X-Git-Url: http://erislabs.net/gitweb/?a=blobdiff_plain;f=lib%2Fgcd.c;h=5401fd65eec94e9dcb64947ad6d1a3d3724e4583;hb=d1fb15ba12436afdd742eaff1ec46f75ee397bda;hp=539c84d840b5c570519bbca5f46e029f718caf20;hpb=a0a86efb418fd595d9154cf6ca70c78ea777afb9;p=gnulib.git diff --git a/lib/gcd.c b/lib/gcd.c index 539c84d84..5401fd65e 100644 --- a/lib/gcd.c +++ b/lib/gcd.c @@ -1,11 +1,11 @@ /* Arithmetic. - Copyright (C) 2001-2002, 2006 Free Software Foundation, Inc. + Copyright (C) 2001-2002, 2006, 2009-2010 Free Software Foundation, Inc. Written by Bruno Haible , 2001. - This program is free software; you can redistribute it and/or modify + This program is free software: you can redistribute it and/or modify it under the terms of the GNU General Public License as published by - the Free Software Foundation; either version 2, or (at your option) - any later version. + the Free Software Foundation; either version 3 of the License, or + (at your option) any later version. This program is distributed in the hope that it will be useful, but WITHOUT ANY WARRANTY; without even the implied warranty of @@ -13,8 +13,7 @@ GNU General Public License for more details. You should have received a copy of the GNU General Public License - along with this program; if not, write to the Free Software Foundation, - Inc., 51 Franklin Street, Fifth Floor, Boston, MA 02110-1301, USA. */ + along with this program. If not, see . */ #include @@ -48,39 +47,39 @@ GCD (WORD_T a, WORD_T b) if (a & c) { if (b & c) - goto odd_odd; + goto odd_odd; else - goto odd_even; + goto odd_even; } else { if (b & c) - goto even_odd; + goto even_odd; else - abort (); + abort (); } for (;;) { odd_odd: /* a/c and b/c both odd */ if (a == b) - break; + break; if (a > b) - { - a = a - b; - even_odd: /* a/c even, b/c odd */ - do - a = a >> 1; - while ((a & c) == 0); - } + { + a = a - b; + even_odd: /* a/c even, b/c odd */ + do + a = a >> 1; + while ((a & c) == 0); + } else - { - b = b - a; - odd_even: /* a/c odd, b/c even */ - do - b = b >> 1; - while ((b & c) == 0); - } + { + b = b - a; + odd_even: /* a/c odd, b/c even */ + do + b = b >> 1; + while ((b & c) == 0); + } } /* a = b */