Merge commit 'dd33637c18629c3e554ebb146bbeb45c9745a5cf'
authorMichael Niedermayer <michaelni@gmx.at>
Wed, 16 Oct 2013 10:02:38 +0000 (12:02 +0200)
committerMichael Niedermayer <michaelni@gmx.at>
Wed, 16 Oct 2013 10:07:02 +0000 (12:07 +0200)
* commit 'dd33637c18629c3e554ebb146bbeb45c9745a5cf':
  tiny_psnr: switch f32 handling to floating point

Conflicts:
tests/tiny_psnr.c

Merged-by: Michael Niedermayer <michaelni@gmx.at>
1  2 
tests/tiny_psnr.c

  #include <stdlib.h>
  #include <string.h>
  #include <inttypes.h>
 -#include <assert.h>
+ #include <math.h>
+ #include <float.h>
++#include <limits.h>
  
  #include "libavutil/intfloat.h"
+ #include "libavutil/intreadwrite.h"
  
  #define FFMIN(a, b) ((a) > (b) ? (b) : (a))
  #define F 100
@@@ -121,47 -110,84 +125,52 @@@ static float get_f32l(uint8_t *p
      return v.f;
  }
  
 -int main(int argc, char *argv[])
+ static double get_f64l(uint8_t *p)
+ {
+     return av_int2double(AV_RL64(p));
+ }
 +static int run_psnr(FILE *f[2], int len, int shift, int skip_bytes)
  {
      int i, j;
      uint64_t sse = 0;
-     uint64_t dev;
+     double sse_d = 0.0;
 -    FILE *f[2];
      uint8_t buf[2][SIZE];
-     uint64_t psnr;
 -    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;
-     int maxdist    = 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)
 +    if (!noseek) {
 +        for (i = 0; i < 2; i++) {
 +            uint8_t *p = buf[i];
 +            if (fread(p, 1, 12, f[i]) != 12)
                  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 (!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 (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]);
      i = FFMIN(size0, size1) / len;
      if (!i)
          i = 1;
-     dev = int_sqrt(((sse / i) * F * F) + (((sse % i) * F * F) + i / 2) / i);
-     if (sse)
-         psnr = ((2 * log16(max << 16) + log16(i) - log16(sse)) *
-                 284619LL * F + (1LL << 31)) / (1LL << 32);
-     else
-         psnr = 1000 * F - 1; // floating point free infinity :)
-     printf("stddev:%5d.%02d PSNR:%3d.%02d MAXDIFF:%5d bytes:%9d/%9d\n",
-            (int)(dev / F), (int)(dev % F),
-            (int)(psnr / F), (int)(psnr % F),
-            maxdist, size0, size1);
-     return psnr;
+     switch (len) {
+     case 1:
+     case 2: {
+         uint64_t psnr;
+         uint64_t dev = int_sqrt(((sse / i) * F * F) + (((sse % i) * F * F) + i / 2) / i);
+         if (sse)
+             psnr = ((2 * log16(max << 16) + log16(i) - log16(sse)) *
+                     284619LL * F + (1LL << 31)) / (1LL << 32);
+         else
+             psnr = 1000 * F - 1; // floating point free infinity :)
+         printf("stddev:%5d.%02d PSNR:%3d.%02d MAXDIFF:%5"PRIu64" bytes:%9d/%9d\n",
+                (int)(dev / F), (int)(dev % F),
+                (int)(psnr / F), (int)(psnr % F),
+                maxdist, size0, size1);
 -        break;
++        return psnr;
+         }
+     case 4:
+     case 8: {
+         char psnr_str[64];
++        double psnr = INT_MAX;
+         double dev = sqrt(sse_d / i);
+         uint64_t scale = (len == 4) ? (1ULL << 24) : (1ULL << 32);
+         if (sse_d) {
 -            double psnr = 2 * log(DBL_MAX) - log(i / sse_d);
++            psnr = 2 * log(DBL_MAX) - log(i / sse_d);
+             snprintf(psnr_str, sizeof(psnr_str), "%5.02f", psnr);
+         } else
+             snprintf(psnr_str, sizeof(psnr_str), "inf");
+         maxdist = maxdist_d * scale;
+         printf("stddev:%10.2f PSNR:%s MAXDIFF:%10"PRIu64" bytes:%9d/%9d\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 (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;
 +            }
 +        }
 +    }
 +
 +    if (argc < 3) {
 +        printf("tiny_psnr <file1> <file2> [<elem size> [<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 (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;
  }