MMDPExperiment.java

package mmdp;
import mmdp.instance.MMDPInstance;
import mmdp.instance.MMDPInstancesManager;
import mmdp.algorithm.MMDPRandomConstructive;

import java.io.IOException;

import jmh.algorithm.MultiStartSearch;
import jmh.algorithm.MultiStartSearch.MSConfig;
import jmh.algorithm.Objective;
import jmh.algorithm.Solution;
import jmh.experiment.ExperimentManager;
import jmh.instance.FormatException;

/**
 * Runs the MaxMin Diversity experiment: for every {@code GKD-Ia} instance it
 * compares the four multi-start configurations and reports the results. MMDP is
 * a maximization problem — the only difference from the CPH/CWP runners is
 * the {@link Objective}.
 */
public class MMDPExperiment {

    private static final long TIMEOUT_MILLIS = 10 * 1000;

    public static void main(String[] args) throws IOException, FormatException {

        var instancesManager = new MMDPInstancesManager();
        var instances = instancesManager.getInstancesByType("GKD-Ia");

        var manager = new ExperimentManager(Objective.MAXIMIZE);

        for (var instanceFile : instances) {

            MMDPInstance instance = instanceFile.loadInstance();
            Solution solution;

            solution = search(instance, MSConfig.RANDOM);
            manager.addSolution(instance, "Random", solution);

            solution = search(instance, MSConfig.FIRST_IMPROVEMENT_RANDOM);
            manager.addSolution(instance, "First Improvement Random", solution);

            solution = search(instance, MSConfig.FIRST_IMPROVEMENT_LEXICOGRAPHICAL);
            manager.addSolution(instance, "First Improvement Lexicographical", solution);

            solution = search(instance, MSConfig.BEST_IMPROVEMENT);
            manager.addSolution(instance, "Best Improvement", solution);
        }

        manager.experimentFinished();
    }

    private static Solution search(MMDPInstance instance, MSConfig config) {
        return new MultiStartSearch(new MMDPRandomConstructive(instance), config, Objective.MAXIMIZE)
                .calculateSolution(TIMEOUT_MILLIS);
    }

}