[BACK]Return to integer_linear_inequalities.adb CVS log [TXT][DIR] Up to [local] / OpenXM_contrib / PHC / Ada / Math_Lib / Supports
File: [local] / OpenXM_contrib / PHC / Ada / Math_Lib / Supports / integer_linear_inequalities.adb (download)

Revision 1.1.1.1 (vendor branch), Sun Oct 29 17:45:28 2000 UTC (23 years, 8 months ago) by maekawa
Branch: PHC, MAIN
CVS Tags: v2, maekawa-ipv6, RELEASE_1_2_3, RELEASE_1_2_2_KNOPPIX_b, RELEASE_1_2_2_KNOPPIX, RELEASE_1_2_2, RELEASE_1_2_1, HEAD
Changes since 1.1: +0 -0 lines 

Import the second public release of PHCpack.

OKed by Jan Verschelde.

with Standard_Floating_Numbers;          use Standard_Floating_Numbers;
with Standard_Floating_Vectors;          use Standard_Floating_Vectors;
with Standard_Floating_Matrices;
with Floating_Linear_Inequalities;       use Floating_Linear_Inequalities;

package body Integer_Linear_Inequalities is

  procedure Integer_Complementary_Slackness
                  ( tableau : in out Matrix; feasible : out boolean ) is
  begin
    Integer_Complementary_Slackness(tableau,tableau'last(2)-1,feasible);
  end Integer_Complementary_Slackness;

  procedure Integer_Complementary_Slackness
                  ( tableau : in out Matrix; lastcol : in integer;
                    feasible : out boolean ) is

    tab : Standard_Floating_Matrices.Matrix
                 (tableau'range(1),tableau'first(2)..lastcol);
    rhs,sol : Standard_Floating_Vectors.Vector(tab'range(1));
    tol : constant double_float := 10.0**(-12);  -- double precision
    columns : Standard_Integer_Vectors.Vector(sol'range);

  begin
    for i in tab'range(1) loop
      for j in tab'range(2) loop
        tab(i,j) := double_float(tableau(i,j));
      end loop;
    end loop;
    for i in rhs'range loop
      rhs(i) := double_float(tableau(i,tableau'last(2)));
    end loop;
    Complementary_Slackness(tab,lastcol,rhs,tol,sol,columns,feasible);
  end Integer_Complementary_Slackness;

end Integer_Linear_Inequalities;