/**********************************************************************
# Copyright (c) EWA Systems Inc. 1998 - 2005 All Rights Reserved
#               No part of this program may be photocopied, reproduced,
#               or translated to another programming language without
#               the prior written consent of EWA Systems.
**********************************************************************/

package com.ewasystems.opt.ga.test;

import com.ewasystems.opt.ga.GeneticAlleleSet;
import com.ewasystems.opt.ga.GeneticEngine;
import com.ewasystems.opt.ga.GeneticGene;
import com.ewasystems.opt.ga.GeneticProblem;
import com.ewasystems.opt.ga.GeneticRunTask;

/**
 * <p>Genetic Algorithm Test Case: Ackley Function</p>
 * <p>Description: Minimize f(x) = 20+e-20*exp(-0.2*(sqrt((1/n)*sum(x(i)^2))-exp((1/n)*sum(cos(2*Pi*x(i))) where X = [-30, 30] by .001</p>
 * <p>Copyright: EWASystems Copyright (c) 1998 - 2005</p>
 * <p>Company: EWASystems</p>
 * @author Lincoln Evans-Beauchamp
 * @version 1.0
 */
public class Ackley extends GeneticProblem
{
  /**
   * Problem Constructor
   */
  public Ackley()
  {
    super("Ackley Function", "Minimize f(x) = 20+e-20*exp(-0.2*(sqrt((1/n)*sum(x(i)^2))-exp((1/n)*sum(cos(2*Pi*x(i))) where X = [-30, 30] by .001");
  }

  /**
   * The Objective Function
   * @param gene The Genetic Gene to Test
   * @return double: The Genetic Gene's Fitness
   */
  public double f(GeneticGene gene)
  {
    double f1 = 0d, f2 = 0d;
    for (int i=0; i<gene.getGeneLength(); i++)
    {
      double x = ((Double)gene.getGeneValue(i)).doubleValue();
      f1 += x*x;
      f2 += Math.cos(2d*Math.PI*x);
    }
    return 20d + Math.exp(-20d)*Math.exp(-0.2d*Math.sqrt(.1d)*f1 - Math.exp(.1d)*f2);
  }

  /**
   * Gets the Allele Set for the Problem
   * @return GeneticAlleleSet: The Allele Set
   */
  public GeneticAlleleSet getAlleleSet()
  {
    GeneticAlleleSet gas = new GeneticAlleleSet();
    Double[] values = new Double[] {new Double(-30d), new Double(.001d), new Double(30d)};
    gas.addContinuousAllele("Double [-30, 30] by .001", values, 10);
    return gas;
  }

  /**
   * Gets the Objective Type
   * @return int: The Objective Type
   */
  public int getObjectiveType()
    {return OBJECTIVE_TYPE_MIN;}

  /**
   * Main for Running the Problem
   * @param args Not Used
   */
  public static final void main(String[] args)
  {
    GeneticEngine ge = new GeneticEngine();
    GeneticRunTask grt = ge.createGeneticTask(new Ackley());
    grt.getConfig().setShowGenerationResults(true);
    grt.getConfig().setShowFinalReport(true);
    grt.run();
    System.out.println("Solution Gene: Fitness = "+grt.getResultGeneFitness());
    for (int i=0; i<grt.getResultGeneLength(); i++)
      {System.out.println("X"+i+": "+grt.getResultGeneValue(i));}
  }
}