/* nag_summary_stats_1var(g01aac) Example Program
 *
 * Copyright 1990 Numerical Algorithms Group.
 *
 * Mark 1, 1990.
 *
 * Mark 5 revised, 1998.
 *
 */

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

#define NMAX 30

int main(void)
{
  double xsd, xskew, xkurt, wsum, xmean, xmax, xmin;
  double wt[NMAX], x[NMAX];
  Integer i, weight, j, n, nprob, nvalid;
  static NagError fail;

  /* Skip heading in data file */
  Vscanf("%*[^\n]");
  Vprintf("g01aac Example Program Results\n");
  Vscanf("%ld", &nprob);
  for (j=1; j<=nprob; j++)
    {
      Vscanf("%ld %ld", &n, &weight);
      Vprintf("Problem  %5ld\n", j);
      Vprintf("Number of cases %ld\n", n);
      if (n<=0 || n>NMAX)
        {
          Vfprintf(stderr, "Error: n is out of range: n = %5ld\n", n);
          return EXIT_FAILURE;
        }
      else
        {
          for (i=0; i<n; i++)
            Vscanf("%lf", &x[i]);
          Vprintf("Data as input -\n");
          for (i=0; i<n; i++)
            Vprintf("%12.1f%c", x[i], (i%5==4 || i==n-1) ? '\n' : ' ');
          if (weight)
            {
              Vprintf("Weights as input -\n");
              for (i=0; i<n; i++)
                Vscanf("%lf", &wt[i]);
              for (i=0; i<n; i++)
                Vprintf("%12.1f%c", wt[i], (i%5==4 || i==n-1) ? '\n' : ' ');
              g01aac(n, x, wt, &nvalid, &xmean, &xsd, &xskew, &xkurt, &xmin,
                     &xmax, &wsum, &fail);
            }
          else
            g01aac(n, x, (double *)0, &nvalid, &xmean, &xsd, &xskew, &xkurt,
                   &xmin, &xmax, &wsum, &fail);

          if (fail.code==NE_NOERROR)
            {
              Vprintf("\n");
              Vprintf("Successful call of g01aac\n");
              Vprintf("No. of valid cases %5ld\n", nvalid);
              Vprintf("Mean           %13.1f\n", xmean);
              Vprintf("Std devn       %13.1f\n", xsd);
              Vprintf("Skewness       %13.1f\n", xskew);
              Vprintf("Kurtosis       %13.1f\n", xkurt);
              Vprintf("Minimum        %13.1f\n", xmin);
              Vprintf("Maximum        %13.1f\n", xmax);
              Vprintf("Sum of weights %13.1f\n", wsum);
            }
          else
            {
              Vprintf("Unsuccessful call of g01aac\n");
              Vprintf("%s \n", fail.message);
              if (fail.code==NE_CASES_ONE)
                {
                  Vprintf("No. of valid cases %5ld\n", nvalid);
                  Vprintf("Mean           %13.1f\n", xmean);
                  Vprintf("Minimum        %13.1f\n", xmin);
                  Vprintf("Maximum        %13.1f\n", xmax);
                  Vprintf("Sum of weights %13.1f\n", wsum);
                  Vprintf("Std devn and coeffts of skewness\n");
                  Vprintf("and kurtosis not defined\n");
                }
              else
                Vprintf("%s \n", fail.message);
              return EXIT_FAILURE;
            }
        }
    }
  return EXIT_SUCCESS;
}