Return to standard_integer_norms.adb CVS log

Up to [local] / OpenXM_contrib / PHC / Ada / Math_Lib / Matrices

Annotation of OpenXM_contrib/PHC/Ada/Math_Lib/Matrices/standard_integer_norms.adb, Revision 1.1

1.1     ! maekawa     1: with Standard_Common_Divisors;           use Standard_Common_Divisors;
        !             2:
        !             3: package body Standard_Integer_Norms is
        !             4:
        !             5:   function gcd ( v : Vector ) return integer is
        !             6:
        !             7:     tmp : Vector(v'range);
        !             8:     res : integer;
        !             9:
        !            10:   begin
        !            11:     for i in v'range loop
        !            12:       if v(i) < 0
        !            13:        then tmp(i) := -v(i);
        !            14:        else tmp(i) :=  v(i);
        !            15:       end if;
        !            16:     end loop;
        !            17:     res := tmp(tmp'first);
        !            18:     for i in (tmp'first+1)..tmp'last loop
        !            19:       res := gcd(res,tmp(i));
        !            20:       exit when (res = 1);
        !            21:     end loop;
        !            22:     return res;
        !            23:   end gcd;
        !            24:
        !            25:   procedure Normalize ( v : in out Vector ) is
        !            26:
        !            27:     g : constant integer := gcd(v);
        !            28:
        !            29:   begin
        !            30:     if (g /= 0) and then (g /= 1)
        !            31:      then for i in v'range loop
        !            32:             v(i) := v(i)/g;
        !            33:           end loop;
        !            34:     end if;
        !            35:   end Normalize;
        !            36:
        !            37: end Standard_Integer_Norms;