//      g13ac Example Program Text
//      C# version, NAG Copyright 2008
using System;
using NagLibrary;
using System.Globalization;
namespace NagDotNetExamples
{
  public class G13ACE
  {
    static string datafile = "ExampleData/g13ace.d";
    static void Main(String[] args)
    {
      if (args.Length == 1)
      {
        datafile = args[0];
      }
      StartExample();
    }
    public static void StartExample()
    {
      try
      {
        DataReader sr = new DataReader(datafile);
        int i,   nk,   nl,   nvl; int ifail;
        Console.WriteLine("g13ac Example Program Results");
        //      Skip heading in data file
        sr.Reset();
        sr.Reset();
        nk = int.Parse(sr.Next());
        nl = int.Parse(sr.Next());
        double[] ar = new double[nl];
        double[] p = new double[nl];
        double[] r = new double[nk];
        double[] v = new double[nl];
        Console.WriteLine("");
        if (nl > 0 && nk > 0)
        {
          sr.Reset();
          for (i = 1; i <= nk; i++)
          {
            r[i - 1] = double.Parse(sr.Next(), CultureInfo.InvariantCulture);
          }
          //
          G13.g13ac(r, nk, nl, p, v, ar, out nvl, out ifail);
          //
          if (ifail < 0)
          {
            Console.WriteLine("");
            Console.WriteLine("** g13ac failed with ifail = {0,5}", ifail);
          }
          else
          {
            if (ifail != 0)
            {
              Console.WriteLine("  {0}{1,5}", "g13ac fails. ifail = ", ifail);
              Console.WriteLine("");
            }
            if (ifail == 3)
            {
              Console.WriteLine("  {0}{1,2}{2}", "   Only", nvl, "valid sets were generated");
              Console.WriteLine("");
            }
            if ((ifail == 0) || (ifail == 3))
            {
              Console.WriteLine(" {0}", "Lag  Partial    Predictor error  Autoregressive");
              Console.WriteLine(" {0}", ("     autocorrn  variance ratio") + ("   parameter"));
              Console.WriteLine("");
              for (i = 1; i <= nvl; i++)
              {
                Console.Write(" {0}{1,10:f3}{2,15:f3}{3,15:f3}\r\n", i, p[i - 1], v[i - 1], ar[i - 1]);
              }
              Console.WriteLine("");
            }
          }
        }
        //
      }
      catch (Exception e)
      {
        Console.WriteLine(e.Message);
        Console.WriteLine("Exception Raised");
      }
    }
  }
}