/* nag_tsa_auto_corr_part(g13acc) Example Program
 *
 * Copyright 1991 Numerical Algorithms Group.
 *
 * Mark 2, 1991.
 */

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

#define NLMAX   5
#define NKMAX   10

int main(void)
{

  Integer i, nk, nl, nvl;
  double ar[NLMAX], p[NLMAX], r[NKMAX], v[NLMAX];
  static NagError fail;

  Vprintf("g13acc Example Program Results\n");
  /* Skip heading in data file */
  Vscanf("%*[^\n]");
  Vscanf("%ld %ld", &nk, &nl);

  if (nl>0 && nl<=NLMAX && nk>0 && nk<=NKMAX)
    {
      for (i=0; i<nk; ++i)
        Vscanf("%lf", &r[i]);

      g13acc(r, nk, nl, p, v, ar, &nvl, &fail);

      if (fail.code != NE_NOERROR)
        {
          Vprintf("  ****** ERROR ******  \n");
          Vprintf("g13acc fails with failure message \n%s\n", fail.message);
        }
      if (fail.code == NE_CORR_NOT_POS_DEF)
        Vprintf("   Only %2ld valid sets were generated\n\n", nvl);
      if (fail.code == NE_NOERROR || fail.code == NE_CORR_NOT_POS_DEF)
        {
          Vprintf("\n Lag  Partial    Predictor error  Autoregressive\n");
          Vprintf("      autocorrn  variance ratio   parameter\n\n");
          for (i=0; i<nvl; ++i)
            Vprintf(" %2ld%9.3f%16.3f%14.3f\n", i+1, p[i], v[i], ar[i]);
        }
      if (fail.code != NE_NOERROR)
        return EXIT_FAILURE;
    }
  else
    {
      Vfprintf(stderr, "One or both of nl and nk are out of range:"
               " nl = %-3ld while  nk = %-3ld\n", nl, nk);

      return EXIT_FAILURE;
    }
  return EXIT_SUCCESS;
}