/* This program converts the data files into a format acceptable by Matlab */

#include <stdio.h>
#include <iostream.h>

typedef struct {
  long type; /*type*/
  long mrows; /* row dimension */
  long ncols; /* column dimensions*/
  long imagf; /* flag to indicate imaginary part */
  long namelen; /* name lenth + null character */
} Fmatrix;

int main () {

  double W;
  long enddata;
  short w,w1,w2,dum;
  char file[85];

  cout<<"Enter File Name Of Raw Data File :";
  cin>>file;
  cout<<"\n";
  cout<<"Output File Is sound_dat.mat\n";

  FILE *fpraw;

  if ((fpraw = fopen(file,"rb")) == NULL) {
    cout<<"Error:cannot open raw data file \n";
    exit(0);
  }
  long counter=0;

  while(!feof(fpraw)) {
    fread(&w1,sizeof(short),1,fpraw);

/*    if(counter > 30000 && counter <30100)
      cout<<"w1="<<w1<<"    w2="<<w2<<"\n";    */
    counter +=1;
  }
 
  enddata = counter;
  rewind(fpraw);

  FILE *fpmat;

  if ((fpmat = fopen("sound_dat.mat","w")) == NULL) {
    cout<<"Error:cannot open file for matlab header\n";
    exit(0);
  }


  Fmatrix x;
  x.type=1100;
  char *pname;
  pname = "W";
  x.namelen = 2;
  x.mrows = enddata;
  x.ncols = 1;
  x.imagf = 0;

  fwrite(&x,sizeof(Fmatrix),1,fpmat);
  fwrite(pname,sizeof(char),x.namelen,fpmat);

  for(long i=0;i<enddata;i++){
   fread(&w,sizeof(short),1,fpraw);
   W = (double)((double)1*w);
   fwrite(&W,sizeof(double),1,fpmat);
  }

  fclose(fpraw);
  fclose(fpmat);
  
return 0;
}