/* nag_ken_spe_corr_coeff(g02brc) Example Program
 *
 * Copyright 1992 Numerical Algorithms Group.
 *
 * Mark 3, 1993.
 */

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

#define NMAX 20
#define MMAX 20

int main(void)
{
  Integer i, j, m, n, tdx, tdc;
  double x[NMAX][MMAX], corr[MMAX][MMAX];
  Integer svar[MMAX], sobs[NMAX];
  Integer *svarptr, *sobsptr;
  char s, w;

  tdx = MMAX;
  tdc = MMAX;

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

  /* Skip heading in data file */
  Vscanf(" %*[^\n]");

  /* Read data */
  Vscanf("%ld%ld\n",&m,&n);
  if (m>=1 && m<=MMAX && n>=1 && n<=NMAX)
    {
      for (i=0; i<n; i++)
        for (j=0; j<m; j++)
          Vscanf("%lf",&x[i][j]);

      /* Read flag specifying if svar is to be supplied */
      Vscanf(" %c",&s);
      if (s == 'S' || s == 's')
        {
          /* Assign pointer to svar and read in values for svar */
          svarptr = svar;
          for (i=0; i<m; i++)
            Vscanf("%ld",&svar[i]);
        }
      else
        {
          /* Assign pointer to NULL and discard rest of line */
          svarptr = (Integer *)0;
          Vscanf(" %*[^\n]"); /* skip rest of line */
        }

      /* Read flag specifying if sobs is to be supllied */
      Vscanf(" %c",&w);
      if (w == 'W' || w == 'w')
        {
          /* Assign pointer to sobs and read in values for sobs */
          sobsptr = sobs;
          for (i=0; i<n; i++)
            Vscanf("%ld",&sobs[i]);
        }
      else
        {
          /* Assign pointer to NULL and discard rest of line */
          sobsptr = (Integer *)0;
          Vscanf(" %*[^\n]"); /* skip rest of line */
        }

      /* Calculate the Kendall and Spearman coefficients */
      g02brc(n, m, &x[0][0], tdx, svarptr, sobsptr, &corr[0][0], tdc,
             NAGERR_DEFAULT);

      Vprintf("\nCorrelation coefficients:\n\n");
      for (i=0; i<m; i++)
        {
          for (j=0; j<m; j++)
            Vprintf("%8.5f ", corr[i][j]);
          Vprintf("\n");
        }
    }
  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;
}