//      g02cf Example Program Text
//      C# version, NAG Copyright 2008
using System;
using NagLibrary;
using System.Globalization;
using System.IO;
namespace NagDotNetExamples
{
  public class G02CFE
  {
    static string datafile = "ExampleData/g02cfe.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 = 3;
        double[,] r = new double[n, n];
        double[,] ssp = new double[n, n];
        double[] std = new double[n];
        double[] xbar = new double[n];
        int[] korder = new int[n];
        int ifail;
        Console.WriteLine("g02cf Example Program Results");
        //      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,11:f4}", "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] = 1;
        korder[1] = 3;
        korder[2] = 2;
        // 
        G02.g02cf(n, korder, xbar, std, ssp, r, out ifail);
        // 
        if (ifail != 0)
        {
          Console.WriteLine(" ");
          Console.WriteLine("** g02cf failed with ifail = {0,5}", ifail);
        }
        else
        {
          Console.Write(" {0}", "New vector XBAR :   ");
          for (i = 1; i <= n; i++)
          {
            Console.Write(" {0,11:f4}", xbar[i - 1]);
          }
          Console.WriteLine(" ");
          Console.WriteLine(" ");
          Console.Write(" {0}", "New vector STD  :   ");
          for (i = 1; i <= n; i++)
          {
            Console.Write(" {0,11:f4}", std[i - 1]);
          }
          Console.WriteLine(" ");
          Console.WriteLine(" ");
          Console.Write(" {0}", "New 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}", "New 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(" ");
        }
        // 
      }
      catch (Exception e)
      {
        Console.WriteLine(e.Message);
        Console.Write("Exception Raised");
      }
    }
  }
}