private SimplexAlgorithm GetResult(SimplexTable simplexTable, int startRowForTransform)
{
SimplexTable simplexT = (SimplexTable)simplexTable.Clone();
CanonicalTransformation.TransformForSimplex(simplexT, startRowForTransform);
Console.WriteLine(simplexT.ToString());
SimplexAlgorithm smpAlg = new SimplexAlgorithm(simplexT, true);
smpAlg.GetResultForSimplex();
if (smpAlg.Result == null)
{
return(smpAlg);
}
int nOfDouble = IsDouble(smpAlg.Result);
if (nOfDouble == -1)
{
return(smpAlg);
}
Console.WriteLine("double result({0})={1}", nOfDouble, smpAlg.Result[nOfDouble]);
//Branch 1
Console.WriteLine("start branch 1");
SimplexTable simplexT1 = (SimplexTable)simplexT.Clone();
AddLimitation(simplexT1, nOfDouble, smpAlg.Result[nOfDouble], "<=");
SimplexAlgorithm smpAlg1;
smpAlg1 = GetResult(simplexT1, simplexT1.nRows - 1);
//Brackh 2
Console.WriteLine("start branch 2");
SimplexTable simplexT2 = (SimplexTable)simplexT.Clone();
AddLimitation(simplexT2, nOfDouble, smpAlg.Result[nOfDouble], ">=");
SimplexAlgorithm smpAlg2;
smpAlg2 = GetResult(simplexT2, simplexT2.nRows - 1);
if (smpAlg1.Result == null)
{
return(smpAlg2);
}
if (smpAlg2.Result == null)
{
return(smpAlg1);
}
if (smpAlg1.FunctionValue > smpAlg2.FunctionValue)
{
return(smpAlg1);
}
else
{
return(smpAlg2);
}
}
public SimplexAlgorithm(SimplexTable simplexTable, bool taskForMax)
{
this.simplexTable = (SimplexTable)simplexTable.Clone();
d = new Fraction[simplexTable.nColumns];
basis = new int[simplexTable.nRows];
TaskForMax = taskForMax;
}
