archive/share/trash/kach.dem

   1 /*      Demo for Purtilo's implementation of
   2         Hacijan's algorithm....                 */
   3
   4 load("kach")$
   5 load("haipart")$
   6 fpprec:16$
   7 if status(feature,ITS)=true
   8 then translate(hach,check,iterate,fx,fq,fh,getl,criterion,accelerate)$
   9
  10 /* dynamalloc:true$ */
  11
  12 /*      now get a matrix a and a matrix b for the problem
  13         of solving  a.x < b   ...               */
  14
  15 a:bfloat(matrix([1,0],[0,1],[-1,-1]))$
  16
  17 b:bfloat(matrix([-1],[-1],[4.00000001]))$
  18
  19 n:2$
  20
  21 m:3$
  22
  23 /*      the user can find the L as in Hacijan's abstract by
  24         doing ... */
  25
  26 getl(a,b,n,m);
  27
  28 /*      and can find a solution vector x for
  29
  30                 x[1]            <       -1
  31                         x[2]    <       -1
  32                 x[1] +  x[2]    >=      -4
  33
  34         by doing                */
  35
  36 show%all:true$
  37
  38 hach(a,b,n,m,1);
  39
  40 /*      Note that this last example does not pay attention
  41         to the L of the theory.   Doing so, taking a bit longer,
  42         is                      */
  43
  44 show%all:false$
  45
  46 hach(a,b,n,m,9.58);