public FullSolution Run()
{
Solver solver = new Solver();
Random rand = new Random();
if (A is null)
{
if (_generator is null)
{
throw new System.InvalidOperationException("Generating values parameters in GeneratorTaskCondition must be assigned to Problem");
}
(A, B) = _generator.GenerateAB(N, rand);
L = _generator.GenerateL(R, rand);
Alpha = _generator.GenerateAlpha(R);
Cs = new double[R][];
for (var i = 0; i < R; i++)
{
Cs[i] = _generator.GenerateMatrix(N, rand);
}
}
(double[][] xs, double[] fsForXs) = solver.GetSolutions(Cs, A, B);
(double[] solution, _) = solver.SolveSeveral(Cs, A, B, L, Alpha, fsForXs);
double[] newAlpha = new double[R];
Alpha.CopyTo(newAlpha, 0);
double[] optimalX = solver.UpdateX(Cs, A, B, L, ref newAlpha, fsForXs, solution);
double[] fsForX = solver.CalculateFs(Cs, optimalX);
double[] deltas = solver.CalculateDeltas(fsForX, fsForXs);
double[] ys = solver.CalculateYs(deltas, L);
double functionValue = solver.CalculateOptimalFunc(ys, newAlpha);
double[] distances = solver.CalculateDistances(xs, optimalX);
Solution solutionWithoutChange = new Solution(optimalX, functionValue, Alpha, newAlpha, Cs, fsForX, xs, fsForXs, deltas, ys, distances, B.Length);
double[][] newCs = new double[R][];
PercentFinder.CopyMultidimensional(Cs, ref newCs);
ParametersForDefined parameters = new ParametersForDefined(optimalX, newCs, A, B, L, Alpha, _cChangeParameters);
ParametersForDefined = parameters;
double persentOfChange = PercentFinder.FindPercentOfChange(parameters, solver, rand);
//double persentOfChange = PercentFinder.SearchMeanPercent(PercentFinder.FindPercentOfChange, parameters, 0.01, solver, rand);
double[] newX = parameters.NewX;
double[] newNewAlpha = new double[R];
newAlpha.CopyTo(newNewAlpha, 0);
(double[][] newXs, double[] newFsForXs) = solver.GetSolutions(newCs, A, B);
newX = solver.UpdateX(newCs, A, B, L, ref newNewAlpha, newFsForXs, newX);
double[] newfsForX = solver.CalculateFs(newCs, newX);
double[] newDeltas = solver.CalculateDeltas(newfsForX, newFsForXs);
double[] newYs = solver.CalculateYs(newDeltas, L);
double newFunctionValue = solver.CalculateOptimalFunc(newYs, newNewAlpha);
double[] newDistances = solver.CalculateDistances(newXs, newX);
Solution solutionWithChange = new Solution(newX, newFunctionValue, newAlpha, newNewAlpha, newCs, newfsForX, newXs, newFsForXs, newDeltas, newYs, newDistances, B.Length);
return(new FullSolution(A, B, L, solutionWithoutChange, persentOfChange, solutionWithChange));
}
public override SearchParameters Copy()
{
var newCs = new double[_csBase.Length][];
PercentFinder.CopyMultidimensional(_csBase, ref newCs);
var newParams = new ParametersForDefined(OldX, newCs, A, B, L, Alpha, CChangeParameters);
return(newParams);
}
public ParametersForDefined(double[] oldX, double[][] cs, double[] a, double[] b, double[] l, double[] alpha, CChangeParameters parameters) : base(parameters)
{
A = a;
B = b;
L = l;
Alpha = alpha;
OldX = oldX;
_csBase = new double[cs.Length][];
_cs = new double[cs.Length][];
PercentFinder.CopyMultidimensional(cs, ref _csBase);
PercentFinder.CopyMultidimensional(cs, ref _cs);
}
public override void Clear()
{
PercentFinder.CopyMultidimensional(_csBase, ref _cs);
}
