Simplex Algorithm

[ftneek]

The changes you propose to canonicalForm causes test_simplex() to fail on my end. Looks like case 3 fails.
> simplex_test()
["[[-2,-4,1,0,0,-22],[-3,-2,0,1,0,-20],[-4,-5,0,0,1,-40],[6,5,0,0,0,0]] Error: Bad Argument Value"]

Check by reverting to:

Code:

canonicalForm(M, BV):=
BEGIN
LOCAL aij;
FOR I FROM 1 TO SIZE(BV) DO
 aij:=M(I,BV(I));
 IF aij !=0 THEN
  M:=mRow(aij^-1,M,I);
  M:=pivot(M,I,BV(I));
 END;
END;
RETURN M;
END;

> test_simplex()
{1,1,1,1,1,1}
I am also using sum(abs(expected - result)) == 0 as you suggested.