//      g02ce Example Program Text
//      C# version, NAG Copyright 2008
using System;
using NagLibrary;
using System.Globalization;
using System.IO;
namespace NagDotNetExamples
{
  public class G02CEE
  {
    static string datafile = "ExampleData/g02cee.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,   j;
        const int n = 4;
        const int m = 3;
        double[,] r = new double[n, n];
        double[,] r2 = new double[m, m];
        double[,] ssp = new double[n, n];
        double[,] ssp2 = new double[m, m];
        double[] std = new double[n];
        double[] std2 = new double[m];
        double[] xbar = new double[n];
        double[] xbar2 = new double[m];
        int[] korder = new int[m];
        int ifail;
        Console.WriteLine("g02ce Example Program Results");
        Console.WriteLine(" ");
        Console.WriteLine(" ");
        //      Skip heading in data file
        sr.Reset();
        sr.Reset();
        for (i = 1; i <= n; i++)
        {
          xbar[i - 1] = double.Parse(sr.Next(), CultureInfo.InvariantCulture);
        }
        for (i = 1; i <= n; i++)
        {
          std[i - 1] = double.Parse(sr.Next(), CultureInfo.InvariantCulture);
        }
        for (i = 1; i <= n; i++)
        {
          for (j = 1; j <= n; j++)
          {
            ssp[i - 1, j - 1] = double.Parse(sr.Next(), CultureInfo.InvariantCulture);
          }
        }
        for (i = 1; i <= n; i++)
        {
          for (j = 1; j <= n; j++)
          {
            r[i - 1, j - 1] = double.Parse(sr.Next(), CultureInfo.InvariantCulture);
          }
        }
        Console.WriteLine(" ");
        Console.Write(" {0}", "Original vector XBAR :   ");
        for (i = 1; i <= n; i++)
        {
          Console.Write(" {0,11:f4}", xbar[i - 1]);
        }
        Console.WriteLine(" ");
        Console.WriteLine(" ");
        Console.Write(" {0}", "Original vector STD  :   ");
        for (i = 1; i <= n; i++)
        {
          Console.Write(" {0,11:f4}", std[i - 1]);
        }
        Console.WriteLine(" ");
        Console.WriteLine(" ");
        Console.Write(" {0}", "Original matrix SSP  :");
        Console.WriteLine(" ");
        for (i = 1; i <= n; i++)
        {
          for (j = 1; j <= n; j++)
          {
            Console.Write(" {0,11:f4}", ssp[i - 1, j - 1]);
          }
          Console.WriteLine(" ");
        }
        Console.WriteLine(" ");
        Console.WriteLine(" ");
        Console.Write(" {0}", "Original matrix R    :");
        Console.WriteLine(" ");
        for (i = 1; i <= n; i++)
        {
          for (j = 1; j <= n; j++)
          {
            Console.Write(" {0,11:f4}", r[i - 1, j - 1]);
          }
          Console.WriteLine(" ");
        }
        Console.WriteLine(" ");
        Console.WriteLine(" ");
        korder[0] = 4;
        korder[1] = 1;
        korder[2] = 2;
        // 
        G02.g02ce(n, xbar, std, ssp, r, m, korder, xbar2, std2, ssp2, r2, out ifail);
        // 
        if (ifail != 0)
        {
          Console.WriteLine(" ");
          Console.WriteLine("** g02ce failed with ifail = {0,5}", ifail);
        }
        else
        {
          Console.Write(" {0}", "New vector XBAR2 :   ");
          for (i = 1; i <= m; i++)
          {
            Console.Write(" {0,11:f4}", xbar2[i - 1]);
          }
          Console.WriteLine(" ");
          Console.WriteLine(" ");
          Console.Write(" {0}", "New vector STD2  :   ");
          for (i = 1; i <= m; i++)
          {
            Console.Write(" {0,11:f4}", std2[i - 1]);
          }
          Console.WriteLine(" ");
          Console.WriteLine(" ");
          Console.Write(" {0}", "New matrix SSP2  :");
          Console.WriteLine(" ");
          for (i = 1; i <= m; i++)
          {
            for (j = 1; j <= m; j++)
            {
              Console.Write(" {0,11:f4}", ssp2[i - 1, j - 1]);
            }
            Console.WriteLine(" ");
          }
          Console.WriteLine(" ");
          Console.WriteLine(" ");
          Console.Write(" {0,11:f4}", "New matrix R2    :");
          Console.WriteLine(" ");
          for (i = 1; i <= m; i++)
          {
            for (j = 1; j <= m; j++)
            {
              Console.Write(" {0,11:f4}", r2[i - 1, j - 1]);
            }
            Console.WriteLine(" ");
          }
          Console.WriteLine(" ");
        }
        // 
      }
      catch (Exception e)
      {
        Console.WriteLine(e.Message);
        Console.Write( "Exception Raised");
      }
    }
  }
}