#include <fstream.h>
#include <iostream.h>
#include <iomanip.h>
#include "Vector.hh"
#include "Matrix.hh"
#include "Solve.hh"
#include "Clock.hh"
#include "Random.hh"
#include "PeriodicOrbit.hh"
#include <string.h>
#include "lorenz.hh"
#include "lorenz1.hh"
#include "lorenz2.hh"
#include "lorenz3.hh"
#include "lorenz4.hh"
#include <stdlib.h>
#include <stdio.h>

void table2a(){
 char fname[100];
 sprintf(fname, "DATA2/WINNERS");
 ifstream ifile(fname);
 sprintf(fname, "FIGS/table2a");
 ofstream ofile(fname);
 char buf[256];
 int count = 0;
 loadOrbits();
 while(ifile){
   ifile.getline(buf, 256);
   int index;
   char s[50];
   int match = sscanf(buf, "orbit %d: %s", &index, s);
   if(match == 2)
     count++;
   else
     continue;
   if(count > 10)
     break;
   Vector pt(3);
   cout<<s<<endl;
   findFirstPoint(s, pt);
   double T;
   double err;
   Matrix * po = returnPO(s, T, err);
   double lyap = returnLyap(*po, T);
   delete po;
   ofile.setf(ios::scientific, ios::floatfield);
   ofile.precision(13);
   ofile<<s<<" & "<<pt(0)<<" & "<<pt(1)<<" & ";
   ofile.precision(3);
   ofile<<lyap;
   ofile.precision(13);
   ofile<<" & "<<T<<endl;
 }
 ifile.close();
 ofile.close();
}

void table2b(){
 char fname[100];
 sprintf(fname, "DATA2/Lmin.dat");
 ifstream ifile1(fname);
 sprintf(fname, "DATA2/Lmax.dat");
 ifstream ifile2(fname);
 sprintf(fname, "FIGS/table2b");
 ofstream ofile(fname);
 for(int i=2; i <= 20; i++){
   char s[100];
   double lyap;
   ifile1>>s;
   ifile1>>lyap;
   ofile<<s<<" & "<<lyap<<" & ";
   ifile2>>s;
   ifile2>>lyap;
   ofile<<s<<" & "<<lyap<<endl;
 }
 ifile1.close();
 ifile2.close();
 ofile.close();
}

void table2c(){
 char fname[100];
 sprintf(fname, "DATA2/Tmin.dat");
 ifstream ifile1(fname);
 sprintf(fname, "DATA2/Tmax.dat");
 ifstream ifile2(fname);
 sprintf(fname, "FIGS/table2c");
 ofstream ofile(fname);
 for(int i=2; i <= 20; i++){
   char s[100];
   double lyap;
   ifile1>>s;
   ifile1>>lyap;
   ofile<<s<<" & "<<lyap<<" & ";
   ifile2>>s;
   ifile2>>lyap;
   ofile<<s<<" & "<<lyap<<endl;
 }
 ifile1.close();
 ifile2.close();
 ofile.close();
}

void table5a(){
  char fname[50];
  sprintf(fname, "FIGS/table5a");
  ofstream ofile(fname);
  int index;
  char s[50];
  double lyap, T;
  Vector vin(1);
  int m;
  int n;
  n = 1;
  for(m=20; m<=23; m++){
    index = m;
    sprintf(s, "$A^{%d}B$", m);
    sprintf(fname, "DATA4/lyap%d.dat", index);
    vin.input(fname);
    lyap = vin(0);
    sprintf(fname, "DATA4/T%d.dat", index);
    vin.input(fname);
    T = vin(0);
    ofile.setf(0, ios::floatfield);
    ofile.precision(3);
    ofile<<s<<" & "<<lyap<<" & ";
    ofile.precision(13);
    ofile.setf(ios::scientific, ios::floatfield);
    ofile<<T<<endl;
  }
  
  
  for(n=22; n<=25; n++){
    index = n;
    m = n;
    sprintf(s, "$A^{%d}B^{%d}$", m, n);
    sprintf(fname, "DATA3/lyap%d.dat", index);
    vin.input(fname);
    lyap = vin(0);
    sprintf(fname, "DATA3/T%d.dat", index);
    vin.input(fname);
    T = vin(0);
    ofile.setf(0, ios::floatfield);
    ofile.precision(3);
    ofile<<s<<" & "<<lyap<<" & ";
    ofile.precision(13);
    ofile.setf(ios::scientific, ios::floatfield);
    ofile<<T<<endl;
  }
 
  for(n=10; n>=7; n--){
    m = 25;
    index = 100*m + n;
    sprintf(s, "$A^{%d}B^{%d}$", m, n);
    sprintf(fname, "DATA3/lyap%d.dat", index);
    vin.input(fname);
    lyap = vin(0);
    sprintf(fname, "DATA3/T%d.dat", index);
    vin.input(fname);
    T = vin(0);
    ofile.setf(0, ios::floatfield);
    ofile.precision(3);
    ofile<<s<<" & "<<lyap<<" & ";
    ofile.precision(13);
    ofile.setf(ios::scientific, ios::floatfield);
    ofile<<T<<endl;
  }
  ofile.close();
}

void table6a(){
  char fname[50];
  sprintf(fname, "FIGS/table6a");
  ofstream ofile(fname);
  ofile.setf(ios::scientific, ios::floatfield);
  ofile.precision(13);
  for(int i=1; i <= 11; i++){
    sprintf(fname, "DATA6/pt%d.dat", i);
    Vector pt(3);
    pt.input(fname);
    ofile<<i<<" & "<<pt(0)<<" & "<<pt(1)<<" & "<<pt(2)<<endl;
  }
  ofile.close();
}

void table6b(){
  char fname[50];
  sprintf(fname, "FIGS/table6b");
  ofstream ofile(fname);
  ofile.setf(ios::scientific, ios::floatfield);
  ofile.precision(2);
  for(int i=1; i <=11; i++){
    char buf[256];
    sprintf(fname,"DATA6/trace%d.dat",i);
    ifstream ifile(fname);
    for(int j=2; j <=32; j++)
      ifile.getline(buf, 256);
    double d;
    char s[100];
    buf[7+32] = ' ';
    buf[7+32+4+2+5-1] = ' ';
    sscanf(buf, "strng: %s 32) dist %e", s, &d);
    d = atof(buf+7+32+4+2+5);
    cout<<strlen(s)<<endl;
    ofile<<i<<" & "<<s<<" & "<<d<<endl;
  }
  ofile.close();
}

void table7a(){
  char fname[50];
  sprintf(fname, "FIGS/table7a");
  ofstream ofile(fname);
  ofile.setf(ios::scientific, ios::floatfield);
  ofile.precision(2);
  for(int i=1; i <= 11; i++){
    char buf[256];
    sprintf(fname, "DATA6/TRACE%d.dat",i);
    ifstream ifile(fname);
    buf[0] = 'X';
    while(buf[0] != 'B')
      ifile.getline(buf, 256);
    for(int j=0; j < strlen(buf); j++)
      if(buf[j]==':')
	buf[j]=' ';
    char s[100], ss[50];
    double d;
    //cout<<buf<<endl;
    sscanf(buf, "Best fit %s Dist %s",s, ss);
    d = atof(ss);
    //cout<<d<<endl;
    ofile<<i<<" & "<<s<<" & "<<d<<endl;
  }
  ofile.close();
}

void table7b(){
  char fname[50];
  sprintf(fname, "FIGS/table7b");
  ofstream ofile(fname);
  for(int i=1; i <=11; i++){
    sprintf(fname, "DATA6/pt%d.dat", i);
    Vector pt(3);
    pt.input(fname);
    advancePoint2Section(pt, 1.0/128.0);
    ofile.setf(ios::scientific, ios::floatfield);
    ofile.precision(13);
    ofile<<i<<" & "<<pt(0)<<" & "<<pt(1)<<" & ";
    sprintf(fname, "DATA6/lyap%d.dat", i);
    Vector lyap(1);
    lyap.input(fname);
    ofile.setf(0, ios::floatfield);
    ofile.precision(3);
    ofile<<lyap(0)<<endl;
  }
}

int main(){
  table5a();
}