Merge commit '7a00653be6b13131ce1b2cdeca696429f57caaf8'
authorDerek Buitenhuis <derek.buitenhuis@gmail.com>
Tue, 16 Feb 2016 18:11:05 +0000 (18:11 +0000)
committerDerek Buitenhuis <derek.buitenhuis@gmail.com>
Tue, 16 Feb 2016 18:11:05 +0000 (18:11 +0000)
* commit '7a00653be6b13131ce1b2cdeca696429f57caaf8':
  tiny_psnr: Support large files

Merged-by: Derek Buitenhuis <derek.buitenhuis@gmail.com>
1  2 
tests/tiny_psnr.c

@@@ -130,47 -115,79 +130,47 @@@ static double get_f64l(uint8_t *p
      return av_int2double(AV_RL64(p));
  }
  
 -int main(int argc, char *argv[])
 +static int run_psnr(FILE *f[2], int len, int shift, int skip_bytes)
  {
-     int i, j;
+     uint64_t i, j;
      uint64_t sse = 0;
      double sse_d = 0.0;
 -    FILE *f[2];
      uint8_t buf[2][SIZE];
 -    int len = 1;
 -    int64_t max;
 -    int shift      = argc < 5 ? 0 : atoi(argv[4]);
 -    int skip_bytes = argc < 6 ? 0 : atoi(argv[5]);
 +    int64_t max    = (1LL << (8 * len)) - 1;
-     int size0      = 0;
-     int size1      = 0;
+     uint64_t size0   = 0;
+     uint64_t size1   = 0;
      uint64_t maxdist = 0;
      double maxdist_d = 0.0;
 +    int noseek;
  
 -    if (argc < 3) {
 -        printf("tiny_psnr <file1> <file2> [<elem size> [<shift> [<skip bytes>]]]\n");
 -        printf("WAV headers are skipped automatically.\n");
 -        return 1;
 -    }
 -
 -    if (argc > 3) {
 -        if (!strcmp(argv[3], "u8")) {
 -            len = 1;
 -        } else if (!strcmp(argv[3], "s16")) {
 -            len = 2;
 -        } else if (!strcmp(argv[3], "f32")) {
 -            len = 4;
 -        } else if (!strcmp(argv[3], "f64")) {
 -            len = 8;
 -        } else {
 -            char *end;
 -            len = strtol(argv[3], &end, 0);
 -            if (*end || len < 1 || len > 2) {
 -                fprintf(stderr, "Unsupported sample format: %s\n", argv[3]);
 -                return 1;
 -            }
 -        }
 -    }
 -
 -    max = (1LL << (8 * len)) - 1;
 +    noseek = fseek(f[0], 0, SEEK_SET) ||
 +             fseek(f[1], 0, SEEK_SET);
  
 -    f[0] = fopen(argv[1], "rb");
 -    f[1] = fopen(argv[2], "rb");
 -    if (!f[0] || !f[1]) {
 -        fprintf(stderr, "Could not open input files.\n");
 -        return 1;
 -    }
 -
 -    for (i = 0; i < 2; i++) {
 -        uint8_t *p = buf[i];
 -        if (fread(p, 1, 12, f[i]) != 12)
 -            return 1;
 -        if (!memcmp(p, "RIFF", 4) &&
 -            !memcmp(p + 8, "WAVE", 4)) {
 -            if (fread(p, 1, 8, f[i]) != 8)
 -                return 1;
 -            while (memcmp(p, "data", 4)) {
 -                int s = p[4] | p[5] << 8 | p[6] << 16 | p[7] << 24;
 -                fseek(f[i], s, SEEK_CUR);
 +    if (!noseek) {
 +        for (i = 0; i < 2; i++) {
 +            uint8_t *p = buf[i];
 +            if (fread(p, 1, 12, f[i]) != 12)
 +                return -1;
 +            if (!memcmp(p, "RIFF", 4) &&
 +                !memcmp(p + 8, "WAVE", 4)) {
                  if (fread(p, 1, 8, f[i]) != 8)
 -                    return 1;
 +                    return -1;
 +                while (memcmp(p, "data", 4)) {
 +                    int s = p[4] | p[5] << 8 | p[6] << 16 | p[7] << 24;
 +                    fseek(f[i], s, SEEK_CUR);
 +                    if (fread(p, 1, 8, f[i]) != 8)
 +                        return -1;
 +                }
 +            } else {
 +                fseek(f[i], -12, SEEK_CUR);
              }
 -        } else {
 -            fseek(f[i], -12, SEEK_CUR);
          }
 -    }
  
 -    fseek(f[shift < 0], abs(shift), SEEK_CUR);
 +        fseek(f[shift < 0], abs(shift), SEEK_CUR);
  
 -    fseek(f[0], skip_bytes, SEEK_CUR);
 -    fseek(f[1], skip_bytes, SEEK_CUR);
 +        fseek(f[0], skip_bytes, SEEK_CUR);
 +        fseek(f[1], skip_bytes, SEEK_CUR);
 +    }
  
      for (;;) {
          int s0 = fread(buf[0], 1, SIZE, f[0]);
  
          maxdist = maxdist_d * scale;
  
-         printf("stddev:%10.2f PSNR:%s MAXDIFF:%10"PRIu64" bytes:%9d/%9d\n",
+         printf("stddev:%10.2f PSNR:%s MAXDIFF:%10"PRIu64" bytes:%9"PRIu64"/%9"PRIu64"\n",
                 dev * scale, psnr_str, maxdist, size0, size1);
 -        break;
 +        return psnr;
 +    }
 +    }
 +    return -1;
 +}
 +
 +int main(int argc, char *argv[])
 +{
 +    FILE *f[2];
 +    int len = 1;
 +    int shift_first= argc < 5 ? 0 : atoi(argv[4]);
 +    int skip_bytes = argc < 6 ? 0 : atoi(argv[5]);
 +    int shift_last = shift_first + (argc < 7 ? 0 : atoi(argv[6]));
 +    int shift;
 +    int max_psnr   = -1;
 +    int max_psnr_shift = 0;
 +
 +    if (shift_last > shift_first)
 +        shift_first -= shift_last - shift_first;
 +
 +    if (argc > 3) {
 +        if (!strcmp(argv[3], "u8")) {
 +            len = 1;
 +        } else if (!strcmp(argv[3], "s16")) {
 +            len = 2;
 +        } else if (!strcmp(argv[3], "f32")) {
 +            len = 4;
 +        } else if (!strcmp(argv[3], "f64")) {
 +            len = 8;
 +        } else {
 +            char *end;
 +            len = strtol(argv[3], &end, 0);
 +            if (*end || len < 1 || len > 2) {
 +                fprintf(stderr, "Unsupported sample format: %s\nSupported: u8, s16, f32, f64\n", argv[3]);
 +                return 1;
 +            }
 +        }
 +    }
 +
 +    if (argc < 3) {
 +        printf("tiny_psnr <file1> <file2> [<elem size>|u8|s16|f32|f64 [<shift> [<skip bytes> [<shift search range>]]]]\n");
 +        printf("WAV headers are skipped automatically.\n");
 +        return 1;
      }
 +
 +    f[0] = fopen(argv[1], "rb");
 +    f[1] = fopen(argv[2], "rb");
 +    if (!f[0] || !f[1]) {
 +        fprintf(stderr, "Could not open input files.\n");
 +        return 1;
      }
 +
 +    for (shift = shift_first; shift <= shift_last; shift++) {
 +        int psnr = run_psnr(f, len, shift, skip_bytes);
 +        if (psnr > max_psnr || (shift < 0 && psnr == max_psnr)) {
 +            max_psnr = psnr;
 +            max_psnr_shift = shift;
 +        }
 +    }
 +    if (max_psnr < 0)
 +        return 2;
 +
 +    if (shift_last > shift_first)
 +        printf("Best PSNR is %3d.%02d for shift %i\n", (int)(max_psnr / F), (int)(max_psnr % F), max_psnr_shift);
      return 0;
  }