compfile.c

Reference Implementation of G.711 standard and other voice codecs

		CR;
		PRINT_RULE;
		printf("\n Short compare of %s and %s\n", File1, File2);
		PRINT_RULE;
	      }
	      printf("\nBLOCK: %6ld, SAMPLE: %4ld ===>", i + 1, j + 1);
	      if (KindOfDump == 'H')
		printf("\t%6.4X\t\t%6.4X", (short unsigned) a[j],
		       (short unsigned) b[j]);
	      else
		printf("\t%6d\t\t%6d", a[j], b[j]);
	      if (NrDiffs % 22 == 0)
		if ((c = toupper(getchar())) == 'Q' || c == 'X' ||
		    c == 27 || c == 'S')
		  break;
	    }
	    else
	      NrDiffs++;
	  }
	  j++;
	}
      else
      {
	if (l < 0)
	{
	  KILL(File1, 5);
	}
	else if (k<0)
	{
	  KILL(File2, 6);
	}
	else
	  break;
      }
    }
  }
  if (NrDiffs > 0 && !quiet)
    CR;

  /* Release allocated memory */
  free(a); free(b);

  /* Return the number of different samples */
  return(NrDiffs);
}
/* ...................... end of compare_shorts() ..........................*/


/*
  -------------------------------------------------------------------------
  long compare_longs(char *File1, char *File2, int fh1, int fh2,
		      long N, long N1, long N2, char KindOfDump, char quiet);

  Compare long data from 2 files and print in decimal or hex format the
  different samples, indicating block/sample no.

  Return the number of different samples OR -1 on error

  30.Dec.93  v1.0  Simao
  --------------------------------------------------------------------------
*/
long compare_longs(File1,File2,fh1,fh2,N,N1,N2,KindOfDump, quiet)
char *File1, *File2,KindOfDump;
int fh1, fh2;
long N,N1,N2;
char quiet;
{
  long i,j,l,k,NrDiffs;
  char c;
  long *a, *b;

  /* Allocate memory for data vectors */
  if ((a=(long *)calloc(N,sizeof(long)))==NULL)
    return -1;
  if ((b=(long *)calloc(N,sizeof(long)))==NULL) {
    free(a);
    return -1;
  }

  /* Start loop */
  for (c=0, NrDiffs= i = j = 0; i < N2; i++, j = 0)
  {
    if (!(c == 'Q' || c == 'X' || c == 27 || c == 'S'))
    {
      if ((l = read(fh1, a, sizeof(long) * N) / sizeof(long)) > 0 &&
	  (k = read(fh2, b, sizeof(long) * N) / sizeof(long)) > 0)
	while (j < l && j < k)
	{
	  if (a[j] != b[j])
	  {
	    if (!quiet)
	    {
	      if (NrDiffs++ % 22 == 0)
	      {
		CR;
		PRINT_RULE;
		printf("\n Long compare of %s and %s\n", File1, File2);
		PRINT_RULE;
	      }
	      printf("\nBLOCK: %6ld, SAMPLE: %4ld ===>", i + 1, j + 1);
	      if (KindOfDump == 'H')
		printf("\t%12.4lX\t\t%12.4lX", (long unsigned) a[j],
		       (long unsigned) b[j]);
	      else
		printf("\t%12ld\t\t%12ld", a[j], b[j]);
	      if (NrDiffs % 22 == 0)
		if ((c = toupper(getchar())) == 'Q' || c == 'X' ||
		    c == 27 || c == 'S')
		  break;
	    }
	    else
	      NrDiffs++;
	  }
	  j++;
	}
      else
      {
	if (l < 0)
	{
	  KILL(File1, 5);
	}
	else if (k<0)
	{
	  KILL(File2, 6);
	}
	else
	  break;
      }
    }
  }
  if (NrDiffs > 0 && !quiet)
    CR;

  /* Release allocated memory */
  free(a); free(b);

  /* Return the number of different samples */
  return(NrDiffs);
}
/* ...................... end of compare_longs() ..........................*/


/*============================== */
int main(argc, argv)
  int             argc;
  char           *argv[];
/*============================== */
{
  char            C[1];
  int             fh1, fh2;

  long            N, N1, N2, NrDiffs = 0, tot_smp;
  long            delay=0, start_byte1, start_byte2, samplesize;
  char            File1[50], File2[50];
  char            KindOfDump = 'D', TypeOfData = 'I', quiet=0;
  FILE           *f1, *f2;
#ifdef VMS
  char            mrs[15] = "mrs=";
#endif


  /* ......... GET PARAMETERS ......... */

  /* Check options */
  if (argc < 2)
    display_usage();
  else
  {
    while (argc > 1 && argv[1][0] == '-')
      if (strcmp(argv[1],"-delay")==0)
      {
	/* Get skip length */
	delay = atol(argv[2]);

	/* Move arg{c,v} over the option to the next argument */
	argc -= 2;
	argv += 2;
      }
      else if (strcmp(argv[1], "-h") == 0)
      {
	/* Set dump as hex type */
	KindOfDump = 'H';

	/* Move arg{c,v} over the option to the next argument */
	argc--;
	argv++;
      }
      else if (strcmp(argv[1], "-d") == 0)
      {
	/* Set dump in decimal format */
	KindOfDump = 'D';

	/* Move arg{c,v} over the option to the next argument */
	argc--;
	argv++;
      }
      else if (strcmp(argv[1], "-float") == 0)
      {
	/* Set data type as real */
	TypeOfData = 'R';

	/* Move arg{c,v} over the option to the next argument */
	argc--;
	argv++;
      }
      else if (strcmp(argv[1], "-double") == 0)
      {
	/* Set data type as double real */
	TypeOfData = 'D';

	/* Move arg{c,v} over the option to the next argument */
	argc--;
	argv++;
      }
      else if (strcmp(argv[1], "-short") == 0)
      {
	/* Set data type as real */
	TypeOfData = 'I';

	/* Move arg{c,v} over the option to the next argument */
	argc--;
	argv++;
      }
      else if (strcmp(argv[1], "-long") == 0)
      {
	/* Set data type as real */
	TypeOfData = 'L';

	/* Move arg{c,v} over the option to the next argument */
	argc--;
	argv++;
      }
      else if (strcmp(argv[1], "-q") == 0)
      {
	/* Set quiet compare - only log the total differences */
	quiet = 1;

	/* Move arg{c,v} over the option to the next argument */
	argc--;
	argv++;
      }
      else
      {
	fprintf(stderr, "ERROR! Invalid option \"%s\" in command line\n\n",
		argv[1]);
	display_usage();
      }
  }

  /* Read parameters for processing */
  GET_PAR_S(1,  "First file? ........................... ", File1);
  GET_PAR_S(2,  "Second file? .......................... ", File2);
  FIND_PAR_L(3, "Record Length? ........................ ", N, 256);
  FIND_PAR_L(4, "Starting Record? ...................... ", N1, 1);
  FIND_PAR_L(5, "Number of Records? .................... ", N2, 0);
  FIND_PAR_C(6, "Decimal or Hex dump? (D,H) [D] ........ ", C[0], KindOfDump);
  KindOfDump = (C[0] != 'H') ? 'D' : 'H';
  FIND_PAR_C(7, "Short(I), long(L), float (R) or double (D) data? .. ", 
                C[0], TypeOfData);
  TypeOfData = C[0];


  /* ......... SOME INITIALIZATIONS ......... */
  --N1;
  /* Define sample size */
  switch(TypeOfData)
  {
    case 'R':
      samplesize = sizeof(float);
      break;
    case 'D':
      samplesize = sizeof(double);
      break;
    case 'I':
      samplesize = sizeof(short);
      break;
    case 'L':
      samplesize = sizeof(long);
      break;
    default:
      HARAKIRI("++==++==++ UNSUPPORTED DATA TYPE ++==++==++\007\n",7);
  }
  
  /* Define 1st sample to compare */
  start_byte1 = start_byte2 = samplesize;
  if (delay>=0)
  {
    start_byte1 *= (N1*N + delay);
    start_byte2 *= (N1*N);
  }
  else
  {
    start_byte1 *= (N1*N);
    start_byte2 *= (N1*N - delay);
  }

  /* Check if is to process the whole file */
  if (N2 == 0)
  {
    struct stat     st;
    long            k, l, s1, s2;

    /* ... find the shortest of the 2 files and the number of blks from it */
    /* ... hey, need to skip the delayed samples! ... */
    stat(File1, &st);
    s1=st.st_size - start_byte1;
    k = ceil((st.st_size - start_byte1) / (double)(N * samplesize));
    stat(File2, &st);
    s2=st.st_size - start_byte2;
    l = ceil((st.st_size - start_byte2) / (double)(N * samplesize));
    N2 = k < l ? k : l;
    tot_smp = (s1 > s2? s1 : s2)/samplesize;
    if (k != l)
      fprintf(stderr, "%%CMP-W-DIFSIZ: Files have different sizes!\n");
  }

  /* Opening test file; abort if there's any problem */
#ifdef VMS
  sprintf(mrs, "mrs=%d", sizeof(short) * N);
#endif

  /* Open input files */
  if ((f1=fopen(File1,RB))==NULL) KILL(File1,3);
  if ((f2=fopen(File2,RB))==NULL) KILL(File2,4);
  fh1= fileno(f1); fh2= fileno(f2);


  /* Move pointer to 1st block of interest */
  if (lseek(fh1, start_byte1, 0l) < 0l)
    KILL(File1, 3);
  if (lseek(fh2, start_byte2, 0l) < 0l)
    KILL(File2, 4);

  /* Some preliminaries */
  N1++;

  /* Dumps the file to the screen */
  switch (TypeOfData)
  {
  case 'I': /* short data */
    NrDiffs = compare_shorts(File1,File2,fh1,fh2,N,N1,N2,KindOfDump,quiet);
    break;

  case 'L': /* short data */
    NrDiffs = compare_longs(File1,File2,fh1,fh2,N,N1,N2,KindOfDump,quiet);
    break;

  case 'R': /* float data */
    NrDiffs = compare_floats(File1,File2,fh1,fh2,N,N1,N2,quiet);
    break;

  case 'D': /* double data */
    NrDiffs = compare_doubles(File1,File2,fh1,fh2,N,N1,N2,quiet);
    break;
  }


  /* Final make-ups */
  if (NrDiffs > 0 && !quiet)
  {
    PRINT_RULE2;
    CR;
  }
  if (quiet)
    printf("%s and %s: %ld different samples found out of %ld.\n",
	   File1, File2, NrDiffs, tot_smp);
  else
    printf("%%CMP-I-NROFDIFFS, %ld different samples found out of %ld.\n",
	   NrDiffs, tot_smp);

  if (NrDiffs > 0 && !quiet)
  {
    PRINT_RULE2;
    CR;
  }

  /* Finalizations */
  close(fh1);
  close(fh2);
#ifndef VMS
  return(0);
#endif
}