//
// Purpose:
//    to implement Yao's algorithm
//
// Authors:
//    Frederico R. B. Cruz
//    Rosangela H. Loschi
//    Departamento de Estatistica
//    Universidade Federal de Minas Gerais
//    E-mail: {fcruz,loschi}@est.ufmg.br
//
// Version:
//    1.00
//
// Date:
//    May/2003
//
//

// inclusion of standard libraries
#include <stdio.h>
#include <math.h>

// prior parameters of mu (normal case)
#define mPrior(i,j) (0.0)
#define vPrior(i,j) (1.0)

// prior parameters of sigma^2 (normal case)
#define aPrior(i,j) (0.01)
#define dPrior(i,j) (4.0)

// inclusion of change point library
#include "chpyao.hpp"

// main program
int main(int argc, char *argv[]) {
   // checking input
   if (argc < 2) {
      fprintf(stderr,"Usage: %s <operation> <script file name>\n", argv[0]);
      fprintf(stderr,"Operation\n");
      fprintf(stderr,"\t1 -> Yao's\t\t[p=k]\n");
      exit(0);
   }
   if (argc < 3) {
      fprintf(stderr, "Usage: %s %s <script_file>\n", argv[0],argv[1]);
      exit(0);
   }
   int aux = 0;
   while ((argv[2][aux]!='\n')&&(argv[2][aux]!='\0')) aux++;
   argv[2][aux] = '\0';
   FILE *inputFile = fopen(argv[2],"r");
   if (inputFile == NULL) {
      fprintf(stderr, "%s: No such file\n", argv[2]);
      exit(0);
   }
   fclose(inputFile);
//
// create object
//
   ChangePoints myChPoint;
//
// perform operation
//
   if (argv[1][0]=='1') {
      // read data
      myChPoint.ReadInputScript(argv[2]);
      myChPoint.Yao();
   }
   return 0;
}