//      e05jb Example Program Text
// 
//      C# version, NAG Copyright 2008
using System;
using NagLibrary;
using System.Globalization;
namespace NagDotNetExamples
{
  public class E05JBE
  {
    const int n=2;
    static string datafile = "ExampleData/e05jbe.d";
    static bool lplot;
    static void Main(String[] args)
    {
      if (args.Length == 1)
      {
        datafile = args[0];
      }
      StartExample();
    }
    public static void StartExample()
    {
      try
      {
        DataReader sr = new DataReader(datafile);
        //      values of optional parameters associated with e05jb
        // 
        double infbnd,   obj; int i,  ibound,  iinit,  j,  sdlist,  stclim;
        string lcsrch="";
        E05.E05JB_OBJFUN objfunE05JB = new E05.E05JB_OBJFUN(objfun);
        E05.E05JB_MONIT monitE05JB = new E05.E05JB_MONIT(monit);
        // 
        int ifail=0;
        Console.WriteLine("e05jb Example Program Results");
        // 
        //      Skip heading in data file
        // 
        sr.Reset();
        // 
        //      Read sdlist from data file
        // 
        sr.Reset();
        sdlist = int.Parse(sr.Next());
        double[] bl = new double[n];
        double[] bu = new double[n];
        double[,] list = new double[n, sdlist];
        double[] ruser = new double[1];
        double[] x = new double[n];
        int[] initpt = new int[n];
        int[] numpts = new int[n];
        // 
        // 
        //         Read ibound, bl, and bu from data file
        // 
        sr.Reset();
        ibound = int.Parse(sr.Next());
        // 
        if (ibound == 0)
        {
          // 
          //            Read in the whole of each bound
          // 
          sr.Reset();
          for (i = 1 ; i <= n ; i++)
          {
            bl[i - 1] = double.Parse(sr.Next(), CultureInfo.InvariantCulture);
          }
          sr.Reset();
          for (i = 1 ; i <= n ; i++)
          {
            bu[i - 1] = double.Parse(sr.Next(), CultureInfo.InvariantCulture);
          }
          // 
        }
        else if (ibound == 3)
        {
          // 
          //            Bounds are uniform: read in only the first entry of each
          // 
          sr.Reset();
          bl[0] = double.Parse(sr.Next(), CultureInfo.InvariantCulture);
          sr.Reset();
          bu[0] = double.Parse(sr.Next(), CultureInfo.InvariantCulture);
          // 
        }
        // 
        //         Read in iinit (and list, numpts and initpt if necessary)
        //         from data file
        // 
        sr.Reset();
        iinit = int.Parse(sr.Next());
        // 
        if (iinit == 3)
        {
          // 
          //            User is specifying the initialization list
          // 
          sr.Reset();
          for (i = 1 ; i <= n ; i++)
          {
            numpts[i - 1] = int.Parse(sr.Next());
          }
          sr.Reset();
          for (i = 1 ; i <= n ; i++)
          {
            for (j = 1 ; j <= numpts[i - 1] ; j++)
            {
              list[i - 1 , j - 1] = double.Parse(sr.Next(), CultureInfo.InvariantCulture);
            }
          }
          sr.Reset();
          for (i = 1 ; i <= n ; i++)
          {
            initpt[i - 1] = int.Parse(sr.Next());
          }
          // 
        }
        // 
        //         Read lplot. Its value determines whether monit displays
        //         information on the current search box
        // 
        sr.Reset();
        lplot = bool.Parse(sr.Next());
        // 
        //         Initialize e05jb
        // 
        E05.e05jbOptions options = new E05.e05jbOptions();
        // 
        if (ifail >= 0)
        {
          // 
          //            Use e05jb to read some options from the end of the data
          //            file
          // 
          options.Set(sr);
          Console.WriteLine("");
          // 
          //            Set 'Static Limit'
          // 
          stclim = 4 * n;
          options.Set("Static Limit", stclim);
          // 
          //            Set 'Infinite Bound Size' tenfold
          // 
          infbnd = 1.00e100;
          options.Set("Infinite Bound Size", infbnd);
          // 
          //            Set the option 'Local Searches Limit = 40'
          // 
          options.Set("Local Searches Limit = 40");
          // 
          //            Set the option 'Local Searches' to 'On'
          // 
          lcsrch = "On";
          options.Set("Local Searches", lcsrch);
          // 
          //            Solve the problem.
          // 
          E05.e05jb(n, objfunE05JB, ibound, iinit, bl, bu, list, numpts, initpt, monitE05JB, x,
                    out obj,
          options, out ifail);
          // 
          Console.WriteLine("");
          if (ifail >= 0)
          {
            Console.WriteLine("  On exit from e05jb, ifail ={0,5}",ifail);
            // 
            if (ifail == 0)
            {
              Console.WriteLine("  Final objective value ={0,11:f5}",obj);
              Console.Write("Global optimum x = ");
              for (i = 1 ; i <= n ; i++)
              {
                Console.Write(" {0, 9:f5}", x[i - 1]);
              }
              Console.WriteLine();
            }
          }
          else
          {
            Console.WriteLine(" ** e05jb returned with ifail = {0, 3}", ifail);
          }
          // 
        }
        // 
        // 
      }
      catch (Exception e)
      {
        Console.WriteLine(e.Message);
        Console.WriteLine("Exception Raised");
      }
    }
    // 
    public static void objfun (int n, double[] x, out double f, int nstate,
    out int inform)
    {
      // 
      //      Routine to evaluate objective function.
      // 
      double x1,   x2; 
      // 
      inform = 0;
      f = 0.0;
      if ((inform) >= 0)
      {
        // 
        //         If inform >= 0 then we're prepared to evaluate objfun
        //         at the current x
        // 
        if (nstate == 1)
        {
          // 
          //            This is the first call to objfun
          // 
          Console.WriteLine("");
          Console.WriteLine("  (OBJFUN was just called for the first time)");
        }
        // 
        x1 = x[0];
        x2 = x[1];
        // 
        f = 3.0e0*(1.0e0-x1)*(1.0e0-x1)*Math.Exp(-(x1*x1)-(x2+1)*(x2+1)) -
        1.0e1*(x1/5.0e0-Math.Pow(x1, 3)-Math.Pow(x2,5))*Math.Exp(-x1*x1-x2*x2) -
        1.0e0/3.0e0*Math.Exp(-(x1+1.0e0)*(x1+1.0e0)-x2*x2);
      }
      // 
      // 
    }
    // 
    public static void monit (int n, int ncall, double[] xbest, int[] icount,
    double[,] list, int[] numpts, int[] initpt, int nbaskt, double[,] xbaskt,
    double[] boxl, double[] boxu, int nstate, out int inform)
    {
      // 
      //      Monitoring method.
      // 
      int i,  j; //  inform;
      // 
      inform = 0;
      if ((inform) >= 0)
      {
        // 
        //         We are going to allow the iterations to continue.
        //         Extract plot from the integer communication-array
        // 
        // 
        if ((nstate == 0) || (nstate == 1))
        {
          // 
          //            When nstate == 1, monit is called for the first time. When
          //            nstate == 0, monit is called for the first AND last time.
          //            Display a welcome message
          // 
          Console.WriteLine("");
          Console.WriteLine("  *** Begin monitoring information ***");
          Console.WriteLine("");
          // 
          if (lplot && n == 2)
          {
            Console.WriteLine("  <Begin displaying search boxes>");
            Console.WriteLine("");
          }
          // 
        }
        // 
        if (lplot && n == 2)
        {
          // 
          //            Display the coordinates of the edges of the current search
          //            box
          // 
          outbox(ref n, boxl, boxu);
          // 
        }
        // 
        if (nstate <= 0)
        {
          // 
          //            monit is called for the last time
          // 
          if (lplot && n == 2)
          {
            Console.WriteLine("  <End displaying search boxes>");
            Console.WriteLine("");
          }
          // 
          Console.WriteLine("  Total sub-boxes ={0,5}",icount[0]);
          Console.WriteLine("  Total function evaluations ={0,5}",ncall);
          Console.WriteLine("  Total function evaluations used in local search ={0,5}",icount[1]);
          Console.WriteLine("  Total points used in local search ={0,5}",icount[2]);
          Console.WriteLine("  Total sweeps through levels ={0,5}",icount[3]);
          Console.WriteLine("  Total splits by init. list ={0,5}",icount[4]);
          Console.WriteLine("  Lowest level with nonsplit boxes ={0,5}",icount[5]);
          Console.WriteLine("  Number of candidate minima in the \"shopping basket\" ={0,5}",nbaskt);
          Console.WriteLine("  Shopping basket:");
          // 
          for (i = 1 ; i <= n ; i++)
          {
            Console.Write("xbaskt({0,3},:) = ",  i);
            for (j = 1 ; j <= nbaskt ; j++)
            {
              Console.Write(" {0, 9:f5}", xbaskt[i - 1 , j - 1]);
            }
            Console.WriteLine();
          }
          // 
          Console.WriteLine("");
          Console.WriteLine("  *** End monitoring information ***");
          Console.WriteLine("");
        }
        // 
      }
      // 
      // 
    }
    public static void outbox (ref int n, double[] boxl, double[] boxu)
    {
      // 
      //      Displays edges of box with bounds boxl and boxu in format suitable
      //      for plotting.
      // 
      // 
      // 
      Console.WriteLine("{0,20:f15}  {1,20:f15}",boxl[0],boxl[1]);
      Console.WriteLine("{0,20:f15}  {1,20:f15}",boxl[0],boxu[1]);
      Console.WriteLine("");
      Console.WriteLine("{0,20:f15}  {1,20:f15}",boxl[0],boxl[1]);
      Console.WriteLine("{0,20:f15}  {1,20:f15}",boxu[0],boxl[1]);
      Console.WriteLine("");
      Console.WriteLine("{0,20:f15}  {1,20:f15}",boxl[0],boxu[1]);
      Console.WriteLine("{0,20:f15}  {1,20:f15}",boxu[0],boxu[1]);
      Console.WriteLine("");
      Console.WriteLine("{0,20:f15}  {1,20:f15}",boxu[0],boxl[1]);
      Console.WriteLine("{0,20:f15}  {1,20:f15}",boxu[0],boxu[1]);
      Console.WriteLine("");
      // 
      // 
    }
  }
}