Return to integer_linear_inequalities.adb CVS log

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

Annotation of OpenXM_contrib/PHC/Ada/Math_Lib/Supports/integer_linear_inequalities.adb, Revision 1.1

1.1     ! maekawa     1: with Standard_Floating_Numbers;          use Standard_Floating_Numbers;
        !             2: with Standard_Floating_Vectors;          use Standard_Floating_Vectors;
        !             3: with Standard_Floating_Matrices;
        !             4: with Floating_Linear_Inequalities;       use Floating_Linear_Inequalities;
        !             5:
        !             6: package body Integer_Linear_Inequalities is
        !             7:
        !             8:   procedure Integer_Complementary_Slackness
        !             9:                   ( tableau : in out Matrix; feasible : out boolean ) is
        !            10:   begin
        !            11:     Integer_Complementary_Slackness(tableau,tableau'last(2)-1,feasible);
        !            12:   end Integer_Complementary_Slackness;
        !            13:
        !            14:   procedure Integer_Complementary_Slackness
        !            15:                   ( tableau : in out Matrix; lastcol : in integer;
        !            16:                     feasible : out boolean ) is
        !            17:
        !            18:     tab : Standard_Floating_Matrices.Matrix
        !            19:                  (tableau'range(1),tableau'first(2)..lastcol);
        !            20:     rhs,sol : Standard_Floating_Vectors.Vector(tab'range(1));
        !            21:     tol : constant double_float := 10.0**(-12);  -- double precision
        !            22:     columns : Standard_Integer_Vectors.Vector(sol'range);
        !            23:
        !            24:   begin
        !            25:     for i in tab'range(1) loop
        !            26:       for j in tab'range(2) loop
        !            27:         tab(i,j) := double_float(tableau(i,j));
        !            28:       end loop;
        !            29:     end loop;
        !            30:     for i in rhs'range loop
        !            31:       rhs(i) := double_float(tableau(i,tableau'last(2)));
        !            32:     end loop;
        !            33:     Complementary_Slackness(tab,lastcol,rhs,tol,sol,columns,feasible);
        !            34:   end Integer_Complementary_Slackness;
        !            35:
        !            36: end Integer_Linear_Inequalities;