halfdan/eva2
1
0
Fork 0
Code Issues Pull Requests Projects Releases Wiki Activity

Support seed populations in the matlab interface problem

Browse Source
This commit is contained in:
Marcel Kronfeld 2011-03-14 15:03:52 +00:00
parent 3cd5bcf73d
commit 5f24abf73a
1 changed files with 63 additions and 0 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

63
src/eva2/server/go/problems/MatlabProblem.java
View File

@ -44,6 +44,7 @@ public class MatlabProblem extends AbstractOptimizationProblem implements Interf
public static boolean TRACE = false; public static boolean TRACE = false;
transient OptimizerRunnable runnable = null; transient OptimizerRunnable runnable = null;
protected boolean allowSingleRunnable = true; protected boolean allowSingleRunnable = true;
private transient Population seedPopulation=null;
protected int problemDimension = 10; protected int problemDimension = 10;
transient PrintStream dos = null; transient PrintStream dos = null;
private double range[][] = null; private double range[][] = null;
@ -320,6 +321,60 @@ public class MatlabProblem extends AbstractOptimizationProblem implements Interf
optimize(optType, outputFilePrefix, null, null); optimize(optType, outputFilePrefix, null, null);
} }
public void clearSeedPopulation() {
seedPopulation=null;
}
public void setSeedPopulation(double[][] seedData, double[][] seedDataFit) {
if (seedData==null) seedPopulation=null;
else {
if ((seedData.length!=seedDataFit.length)
|| (seedData[0].length!=getProblemDimension())) {
System.err.println("Error, unable to set seed population due to mismatching dimensions");
seedPopulation=null;
} else {
log("Setting seed population of size " + seedData.length + "\n");
seedPopulation=new Population(seedData.length);
for (int i=0; i<seedData.length; i++) {
AbstractEAIndividual indy = (AbstractEAIndividual)m_Template.clone();
setIndyGenotype(indy, seedData[i]);
indy.SetFitness(seedDataFit[i]);
seedPopulation.add(indy);
}
}
}
}
private void setIndyGenotype(AbstractEAIndividual indy, double[] ds) {
switch (dataType) {
case typeDouble:
((InterfaceDataTypeDouble)indy).SetDoubleGenotype(ds);
break;
case typeBinary:
((InterfaceDataTypeBinary)indy).SetBinaryGenotype(toBinary(ds));
break;
case typeInteger:
((InterfaceDataTypeInteger)indy).SetIntGenotype(toInteger(ds));
break;
}
}
private int[] toInteger(double[] ds) {
int[] a=new int[ds.length];
for (int i=0; i<ds.length; i++) {
a[i]=(int)Math.round(ds[i]);
}
return a;
}
private BitSet toBinary(double[] ds) {
BitSet bs=new BitSet(ds.length);
for (int i=0; i<ds.length; i++) {
bs.set(i, (ds[i]>0.5));
}
return bs;
}
/** /**
* Start an optimization using the MatlabProblem. The optType references the standard * Start an optimization using the MatlabProblem. The optType references the standard
* optimizer as in OptimizerFactory. An output file prefix is optional. In two arrays, * optimizer as in OptimizerFactory. An output file prefix is optional. In two arrays,
@ -348,6 +403,14 @@ public class MatlabProblem extends AbstractOptimizationProblem implements Interf
else runnable.setOutputTo(1); // only window else runnable.setOutputTo(1); // only window
runnable.setOutputFullStatsToText(outputAllStatsField); runnable.setOutputFullStatsToText(outputAllStatsField);
if (seedPopulation!=null) {
runnable.getGOParams().getOptimizer().setPopulation(seedPopulation);
runnable.setDoRestart(true);
log("Setting seed population of size " + seedPopulation.size() + ", target size " + seedPopulation.getTargetSize() + "\n");
log(BeanInspector.toString(seedPopulation.getStringRepresentation())+"\n");
log("Restart of optimization targetted.\n");
}
// log("in MP optimize C\n"); // log("in MP optimize C\n");
if ((specParams != null) && (specParams.length > 0)) { if ((specParams != null) && (specParams.length > 0)) {
if ((specValues == null) || (specValues.length != specParams.length)) { if ((specValues == null) || (specValues.length != specParams.length)) {