public ProjectedCostFunction(CostFunction costFunction, Vector parametersValues, List <bool> parametersFreedoms)
{
numberOfFreeParameters_ = 0;
fixedParameters_ = parametersValues;
actualParameters_ = parametersValues;
parametersFreedoms_ = parametersFreedoms;
costFunction_ = costFunction;
if (!(fixedParameters_.Count == parametersFreedoms_.Count))
{
throw new ApplicationException("fixedParameters_.Count!=parametersFreedoms_.Count");
}
for (int i = 0; i < parametersFreedoms_.Count; i++)
{
if (!parametersFreedoms_[i])
{
numberOfFreeParameters_++;
}
}
if (!(numberOfFreeParameters_ > 0))
{
throw new ApplicationException("numberOfFreeParameters==0");
}
}
public ProjectedCostFunction(CostFunction costFunction, Vector parametersValues, List<bool> parametersFreedoms)
{
numberOfFreeParameters_ = 0;
fixedParameters_ = parametersValues;
actualParameters_ = parametersValues;
parametersFreedoms_ = parametersFreedoms;
costFunction_ = costFunction;
if (!(fixedParameters_.Count==parametersFreedoms_.Count))
throw new ApplicationException("fixedParameters_.Count!=parametersFreedoms_.Count");
for (int i =0; i<parametersFreedoms_.Count; i++)
if(!parametersFreedoms_[i])
numberOfFreeParameters_++;
if (!(numberOfFreeParameters_>0))
throw new ApplicationException("numberOfFreeParameters==0");
}
public ProjectedCostFunction(CostFunction costFunction, Vector parametersValues, List <bool> parametersFreedoms)
{
numberOfFreeParameters_ = 0;
fixedParameters_ = parametersValues;
actualParameters_ = parametersValues;
parametersFreedoms_ = parametersFreedoms;
costFunction_ = costFunction;
Utils.QL_REQUIRE(fixedParameters_.Count == parametersFreedoms_.Count, () =>
"fixedParameters_.Count!=parametersFreedoms_.Count");
for (int i = 0; i < parametersFreedoms_.Count; i++)
{
if (!parametersFreedoms_[i])
{
numberOfFreeParameters_++;
}
}
Utils.QL_REQUIRE(numberOfFreeParameters_ > 0, () => "numberOfFreeParameters==0");
}
//! default constructor
//public Problem(CostFunction costFunction, Constraint constraint, Vector initialValue = Array())
public Problem(CostFunction costFunction, Constraint constraint, Vector initialValue)
{
costFunction_ = costFunction;
constraint_ = constraint;
currentValue_ = (Vector)initialValue.Clone();
}
//! default constructor
//public Problem(CostFunction costFunction, Constraint constraint, Vector initialValue = Array())
public Problem(CostFunction costFunction, Constraint constraint, Vector initialValue) {
costFunction_ = costFunction;
constraint_ = constraint;
currentValue_ = (Vector)initialValue.Clone();
}
protected void crossover(List <Candidate> oldPopulation,
List <Candidate> population,
List <Candidate> mutantPopulation,
List <Candidate> mirrorPopulation,
CostFunction costFunction)
{
if (configuration().crossoverIsAdaptive)
{
adaptCrossover();
}
Vector mutationProbabilities = getMutationProbabilities(population);
List <Vector> crossoverMask = new InitializedList <Vector>(population.Count);
crossoverMask.ForEach((ii, vv) => crossoverMask[ii] = new Vector(population.First().values.size(), 1.0));
List <Vector> invCrossoverMask = new InitializedList <Vector>(population.Count);
invCrossoverMask.ForEach((ii, vv) => invCrossoverMask[ii] = new Vector(population.First().values.size(), 1.0));
getCrossoverMask(crossoverMask, invCrossoverMask, mutationProbabilities);
// crossover of the old and mutant population
for (int popIter = 0; popIter < population.Count; popIter++)
{
population[popIter].values = Vector.DirectMultiply(oldPopulation[popIter].values, invCrossoverMask[popIter])
+ Vector.DirectMultiply(mutantPopulation[popIter].values,
crossoverMask[popIter]);
// immediately apply bounds if specified
if (configuration().applyBounds)
{
for (int memIter = 0; memIter < population[popIter].values.size(); memIter++)
{
if (population[popIter].values[memIter] > upperBound_[memIter])
{
population[popIter].values[memIter] = upperBound_[memIter]
+ rng_.nextReal()
* (mirrorPopulation[popIter].values[memIter]
- upperBound_[memIter]);
}
if (population[popIter].values[memIter] < lowerBound_[memIter])
{
population[popIter].values[memIter] = lowerBound_[memIter]
+ rng_.nextReal()
* (mirrorPopulation[popIter].values[memIter]
- lowerBound_[memIter]);
}
}
}
// evaluate objective function as soon as possible to avoid unnecessary loops
try
{
population[popIter].cost = costFunction.value(population[popIter].values);
if (Double.IsNaN(population[popIter].cost))
{
population[popIter].cost = Double.MaxValue;
}
}
catch
{
population[popIter].cost = Double.MaxValue;
}
}
}
protected void calculateNextGeneration(List <Candidate> population,
CostFunction costFunction)
{
List <Candidate> mirrorPopulation = null;
List <Candidate> oldPopulation = (List <Candidate>)population.Clone();
switch (configuration().strategy)
{
case Strategy.Rand1Standard:
{
population.Shuffle();
List <Candidate> shuffledPop1 = (List <Candidate>)population.Clone();
population.Shuffle();
List <Candidate> shuffledPop2 = (List <Candidate>)population.Clone();
population.Shuffle();
mirrorPopulation = (List <Candidate>)shuffledPop1.Clone();
for (int popIter = 0; popIter < population.Count; popIter++)
{
population[popIter].values = population[popIter].values
+ configuration().stepsizeWeight
*(shuffledPop1[popIter].values - shuffledPop2[popIter].values);
}
}
break;
case Strategy.BestMemberWithJitter:
{
population.Shuffle();
List <Candidate> shuffledPop1 = (List <Candidate>)population.Clone();
population.Shuffle();
Vector jitter = new Vector(population[0].values.size(), 0.0);
for (int popIter = 0; popIter < population.Count; popIter++)
{
for (int jitterIter = 0; jitterIter < jitter.Count; jitterIter++)
{
jitter[jitterIter] = rng_.nextReal();
}
population[popIter].values = bestMemberEver_.values
+ Vector.DirectMultiply(
shuffledPop1[popIter].values - population[popIter].values
, 0.0001 * jitter + configuration().stepsizeWeight);
}
mirrorPopulation = new InitializedList <Candidate>(population.Count);
mirrorPopulation.ForEach((ii, vv) => mirrorPopulation[ii] = (Candidate)bestMemberEver_.Clone());
}
break;
case Strategy.CurrentToBest2Diffs:
{
population.Shuffle();
List <Candidate> shuffledPop1 = (List <Candidate>)population.Clone();
population.Shuffle();
for (int popIter = 0; popIter < population.Count; popIter++)
{
population[popIter].values = oldPopulation[popIter].values
+ configuration().stepsizeWeight
*(bestMemberEver_.values - oldPopulation[popIter].values)
+ configuration().stepsizeWeight
*(population[popIter].values - shuffledPop1[popIter].values);
}
mirrorPopulation = (List <Candidate>)shuffledPop1.Clone();
}
break;
case Strategy.Rand1DiffWithPerVectorDither:
{
population.Shuffle();
List <Candidate> shuffledPop1 = (List <Candidate>)population.Clone();
population.Shuffle();
List <Candidate> shuffledPop2 = (List <Candidate>)population.Clone();
population.Shuffle();
mirrorPopulation = (List <Candidate>)shuffledPop1.Clone();
Vector FWeight = new Vector(population.First().values.size(), 0.0);
for (int fwIter = 0; fwIter < FWeight.Count; fwIter++)
{
FWeight[fwIter] = (1.0 - configuration().stepsizeWeight)
* rng_.nextReal() + configuration().stepsizeWeight;
}
for (int popIter = 0; popIter < population.Count; popIter++)
{
population[popIter].values = population[popIter].values
+ Vector.DirectMultiply(FWeight,
shuffledPop1[popIter].values - shuffledPop2[popIter].values);
}
}
break;
case Strategy.Rand1DiffWithDither:
{
population.Shuffle();
List <Candidate> shuffledPop1 = (List <Candidate>)population.Clone();
population.Shuffle();
List <Candidate> shuffledPop2 = (List <Candidate>)population.Clone();
population.Shuffle();
mirrorPopulation = (List <Candidate>)shuffledPop1.Clone();
double FWeight = (1.0 - configuration().stepsizeWeight) * rng_.nextReal()
+ configuration().stepsizeWeight;
for (int popIter = 0; popIter < population.Count; popIter++)
{
population[popIter].values = population[popIter].values
+ FWeight * (shuffledPop1[popIter].values -
shuffledPop2[popIter].values);
}
}
break;
case Strategy.EitherOrWithOptimalRecombination:
{
population.Shuffle();
List <Candidate> shuffledPop1 = (List <Candidate>)population.Clone();
population.Shuffle();
List <Candidate> shuffledPop2 = (List <Candidate>)population.Clone();
population.Shuffle();
mirrorPopulation = (List <Candidate>)shuffledPop1.Clone();
double probFWeight = 0.5;
if (rng_.nextReal() < probFWeight)
{
for (int popIter = 0; popIter < population.Count; popIter++)
{
population[popIter].values = oldPopulation[popIter].values
+ configuration().stepsizeWeight
*(shuffledPop1[popIter].values - shuffledPop2[popIter].values);
}
}
else
{
double K = 0.5 * (configuration().stepsizeWeight + 1); // invariant with respect to probFWeight used
for (int popIter = 0; popIter < population.Count; popIter++)
{
population[popIter].values = oldPopulation[popIter].values
+ K
* (shuffledPop1[popIter].values - shuffledPop2[popIter].values
- 2.0 * population[popIter].values);
}
}
}
break;
case Strategy.Rand1SelfadaptiveWithRotation:
{
population.Shuffle();
List <Candidate> shuffledPop1 = (List <Candidate>)population.Clone();
population.Shuffle();
List <Candidate> shuffledPop2 = (List <Candidate>)population.Clone();
population.Shuffle();
mirrorPopulation = (List <Candidate>)shuffledPop1.Clone();
adaptSizeWeights();
for (int popIter = 0; popIter < population.Count; popIter++)
{
if (rng_.nextReal() < 0.1)
{
population[popIter].values = rotateArray(bestMemberEver_.values);
}
else
{
population[popIter].values = bestMemberEver_.values
+ currGenSizeWeights_[popIter]
* (shuffledPop1[popIter].values - shuffledPop2[popIter].values);
}
}
}
break;
default:
Utils.QL_FAIL("Unknown strategy ("
+ Convert.ToInt32(configuration().strategy) + ")");
break;
}
// in order to avoid unnecessary copying we use the same population object for mutants
crossover(oldPopulation, population, population, mirrorPopulation,
costFunction);
}
