/* nag_running_median_smoother(g10cac) Example Program
 *
 * Copyright 1992 Numerical Algorithms Group.
 *
 * Mark 3, 1992.
 */

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

#define NMAX 50

int main(void)
{
  Integer i, isel, n;
  double  rough[NMAX], smooth[NMAX], y[NMAX];
  Nag_Smooth_Type smoother;

  Vprintf("g10cac Example Program Results\n");
  /* Skip heading in data file */
  Vscanf("%*[^\n]");

  Vscanf("%ld", &n);
  if (n>=1 && n<= NMAX)
    {
      for (i=0; i<n; ++i)
        Vscanf("%lf", &y[i]);

      for (isel = 1; isel <= 2; ++isel)
        {
          /* Select the smoothing method. */
          if (isel == 1)
            smoother = Nag_4253H;
          else if (isel == 2)
            smoother = Nag_3RSSH;

          if (smoother == Nag_4253H)
            Vprintf("\nExample 1: Data smoothing using  4253H method.\n\n");
          else if (smoother == Nag_3RSSH)
            Vprintf("\n\nExample 2: Data smoothing using  3RSSH method.\n\n");

          g10cac(smoother, n, y, smooth, rough, NAGERR_DEFAULT);

          Vprintf("     index        data          smooth          rough\n");
          for (i=0; i<n; ++i)
            Vprintf("%9ld%15.4f%15.4f%15.4f\n", i, y[i], smooth[i], rough[i]);
        }
    }
  else
    {
      Vfprintf(stderr, "n is out of range:"
               " n = %-3ld\n", n);
      return EXIT_FAILURE;
    }
  return EXIT_SUCCESS;
}