source: trunk/EventBenchConsole/src/de/ugoe/cs/eventbench/commands/CMDgenerateFixedLengthSequences.java @ 211

package de.ugoe.cs.eventbench.commands;

import java.security.InvalidParameterException;
import java.util.ArrayList;
import java.util.Collection;
import java.util.HashSet;
import java.util.LinkedHashSet;
import java.util.LinkedList;
import java.util.List;
import java.util.Random;
import java.util.Set;

import de.ugoe.cs.eventbench.data.Event;
import de.ugoe.cs.eventbench.data.GlobalDataContainer;
import de.ugoe.cs.eventbench.models.IStochasticProcess;
import de.ugoe.cs.util.console.Command;
import de.ugoe.cs.util.console.Console;

/**
 * <p>
 * Command to generate all sequences of a given length.
 * </p>
 * 
 * @author Steffen Herbold
 * @version 1.0
 */
public class CMDgenerateFixedLengthSequences implements Command {

        /*
         * (non-Javadoc)
         * 
         * @see de.ugoe.cs.util.console.Command#run(java.util.List)
         */
        @Override
        public void run(List<Object> parameters) {
                String modelname;
                String sequencesName;
                int minLength;
                int maxLength;
                boolean all = true;
                int numSequences = -1;
                try {
                        modelname = (String) parameters.get(0);
                        sequencesName = (String) parameters.get(1);
                        minLength = Integer.parseInt((String) parameters.get(2));
                        maxLength = Integer.parseInt((String) parameters.get(3));
                        if (parameters.size() >= 5) {
                                all = Boolean.parseBoolean((String) parameters.get(4));
                        }
                        if (parameters.size() == 6) {
                                numSequences = Integer.parseInt((String) parameters.get(5));
                        }
                } catch (Exception e) {
                        throw new InvalidParameterException();
                }

                IStochasticProcess model = null;
                Object dataObject = GlobalDataContainer.getInstance()
                                .getData(modelname);
                if (dataObject == null) {
                        Console.println("Model " + modelname + " not found in storage.");
                        return;
                } else if (!(dataObject instanceof IStochasticProcess)) {
                        Console.println("Object " + modelname + " not of type MarkovModel!");
                        return;
                }
                model = (IStochasticProcess) dataObject;
                Collection<List<? extends Event<?>>> sequences = new LinkedHashSet<List<? extends Event<?>>>();
                for (int length = minLength; length <= maxLength; length++) {
                        sequences.addAll(model.generateValidSequences(length + 2));
                }
                Console.traceln("" + sequences.size() + " possible");
                if (!all && numSequences < sequences.size()) {
                        List<Double> probabilities = new ArrayList<Double>(sequences.size());
                        double probSum = 0.0;
                        for (List<? extends Event<?>> sequence : sequences) {
                                double prob = model.getProbability(sequence);
                                probabilities.add(prob);
                                probSum += prob;
                        }
                        Set<Integer> drawnSequences = new HashSet<Integer>(numSequences);
                        Random r = new Random();
                        while (drawnSequences.size() < numSequences) {
                                double randVal = r.nextDouble() * probSum;
                                double sum = 0.0d;
                                int index = -1;
                                while (sum < randVal) {
                                        index++;
                                        double currentProb = probabilities.get(index);
                                        sum += currentProb;
                                }
                                if (!drawnSequences.contains(index)) {
                                        drawnSequences.add(index);
                                        probSum -= probabilities.get(index);
                                        probabilities.set(index, 0.0d);
                                }
                        }
                        Collection<List<? extends Event<?>>> retainedSequences = new LinkedList<List<? extends Event<?>>>();
                        int index = 0;
                        for (List<? extends Event<?>> sequence : sequences) {
                                if (drawnSequences.contains(index)) {
                                        retainedSequences.add(sequence);
                                }
                                index++;
                        }
                        sequences = retainedSequences;
                }
                if (GlobalDataContainer.getInstance().addData(sequencesName, sequences)) {
                        Console.traceln("Old data \"" + sequencesName + "\" overwritten");
                }
                Console.println("" + sequences.size() + " sequences generated");
        }

        /*
         * (non-Javadoc)
         * 
         * @see de.ugoe.cs.util.console.Command#help()
         */
        @Override
        public void help() {
                Console.println("Usage: generateFixedLengthSequences <modelname> <sequencesName> <length> {<all>} {<numSequences>}");
        }

}