/* nag_corr_cov(g02bxc) Example Program
 *
 * Copyright 1992 Numerical Algorithms Group.
 *
 * Mark 3, 1992.
 */

#include <nag.h>
#include <stdio.h>
#include <nag_stdlib.h>
#include <nagg02.h>

#define NMAX 5
#define MMAX 5
#define TDX MMAX
#define TDV MMAX
#define TDR MMAX
int main(void)
{
  double x[NMAX][TDX], r[MMAX][TDR], v[MMAX][TDV];
  double wt[NMAX], *wtptr;
  double sw, wmean[MMAX], std[MMAX];
  Integer i, j, n, m;
  char  w;
  Integer tdx, tdr, tdv;
  Integer test;

  Vprintf("g02bxc Example Program Results\n");

  /* Skip heading in data file */
  Vscanf("%*[^\n]");
  tdx = TDX;
  tdr = TDR;
  tdv = TDV;
  test = 0;

  while ((scanf("%ld%ld %c",&m,&n, &w) != EOF))
    {
      if (m>=1 && m<=MMAX && n>=1 && n<=NMAX)
        {
          for(i=0; i<n; i++)
            Vscanf("%lf",&wt[i]);
          for(i=0; i<n; i++)
            for(j=0; j<m; j++)
              Vscanf("%lf",&x[i][j]);

          if (w == 'w')
            wtptr = wt;
          else
            wtptr = (double *)0;

          g02bxc(n, m, &x[0][0], tdx, (Integer *)0, wtptr, &sw, wmean, std,
                 &r[0][0], tdr, &v[0][0], tdv, NAGERR_DEFAULT);

          if (wtptr)
            Vprintf("\nCase %ld --- Using weights\n", ++test);
          else
            Vprintf("\nCase %ld --- Not using weights\n", ++test);

          Vprintf("\nInput data\n");
          for(i=0; i<n; i++)
            Vprintf("%6.1f%6.1f%6.1f%6.1f\n",x[i][0],x[i][1],x[i][2],wt[i]);

          Vprintf("\n");
          Vprintf("Sample means.\n");
          for(i=0; i<m; i++)
            Vprintf("%6.1f\n",wmean[i]);
          Vprintf("\nStandard deviation.\n");
          for(i=0; i<m; i++)
            Vprintf("%6.1f\n",std[i]);

          Vprintf("\nCorrelation matrix.\n");
          for(i=0; i<m; i++)
            {
              for(j=0; j<m; j++)
                Vprintf("  %7.4f  ",r[i][j]);
              Vprintf("\n");
            }

          Vprintf("\nVariance matrix.\n");
          for(i=0; i<m; i++)
            {
              for(j=0; j<m; j++)
                Vprintf("  %7.3f  ",v[i][j]);
              Vprintf("\n");
            }
          Vprintf("\nSum of weights %6.1f\n", sw);
        }
      else
        {
          Vfprintf(stderr, "One or both of m and n are out of range:"
                   " m = %-3ld while  n = %-3ld\n", m, n);

          return EXIT_FAILURE;
        }
    }
  return EXIT_SUCCESS;
}