contrib/tools/makesRGB.c

   1 /* makesRGB.c -- build sRGB-to-linear and linear-to-sRGB conversion tables
   2  *
   3  * Last changed in libpng 1.6.0 [February 14, 2013]
   4  *
   5  * COPYRIGHT: Written by John Cunningham Bowler, 2013.
   6  * To the extent possible under law, the author has waived all copyright and
   7  * related or neighboring rights to this work.  This work is published from:
   8  * United States.
   9  *
  10  * Make a table to convert 8-bit sRGB encoding values into the closest 16-bit
  11  * linear value.
  12  *
  13  * Make two tables to take a linear value scaled to 255*65535 and return an
  14  * approximation to the 8-bit sRGB encoded value.  Calculate the error in these
  15  * tables and display it.
  16  */
  17 #define _C99_SOURCE 1
  18 #include <stdio.h>
  19 #include <math.h>
  20 #include <stdlib.h>
  21
  22 /* pngpriv.h includes the definition of 'PNG_sRGB_FROM_LINEAR' which is required
  23  * to verify the actual code.
  24  */
  25 #include "../../pngpriv.h"
  26
  27 #include "sRGB.h"
  28
  29 /* The tables are declared 'const' in pngpriv.h, so this redefines the tables to
  30  * be used.
  31  */
  32 #define png_sRGB_table sRGB_table
  33 #define png_sRGB_base sRGB_base
  34 #define png_sRGB_delta sRGB_delta
  35
  36 static png_uint_16 png_sRGB_table[256];
  37 static png_uint_16 png_sRGB_base[512];
  38 static png_byte png_sRGB_delta[512];
  39
  40 static const unsigned int max_input = 255*65535;
  41
  42 double
  43 fsRGB(double l)
  44 {
  45    return sRGB_from_linear(l/max_input);
  46 }
  47
  48 double
  49 sRGB(unsigned int i)
  50 {
  51    return fsRGB(i);
  52 }
  53
  54 double
  55 finvsRGB(unsigned int i)
  56 {
  57    return 65535 * linear_from_sRGB(i/255.);
  58 }
  59
  60 png_uint_16
  61 invsRGB(unsigned int i)
  62 {
  63    unsigned int x = nearbyint(finvsRGB(i));
  64
  65    if (x > 65535)
  66    {
  67       fprintf(stderr, "invsRGB(%u) overflows to %u\n", i, x);
  68       exit(1);
  69    }
  70
  71    return (png_uint_16)x;
  72 }
  73
  74 int
  75 main(int argc, char **argv)
  76 {
  77    unsigned int i, i16, ibase;
  78    double min_error = 0;
  79    double max_error = 0;
  80    double min_error16 = 0;
  81    double max_error16 = 0;
  82    double adjust;
  83    double adjust_lo = 0.4, adjust_hi = 0.6, adjust_mid = 0.5;
  84    unsigned int ec_lo = 0, ec_hi = 0, ec_mid = 0;
  85    unsigned int error_count = 0;
  86    unsigned int error_count16 = 0;
  87    int test_only = 0;
  88
  89    if (argc > 1)
  90       test_only = strcmp("--test", argv[1]) == 0;
  91
  92    /* Initialize the encoding table first. */
  93    for (i=0; i<256; ++i)
  94    {
  95       png_sRGB_table[i] = invsRGB(i);
  96    }
  97
  98    /* Now work out the decoding tables (this is where the error comes in because
  99     * there are 512 set points and 512 straight lines between them.)
 100     */
 101    for (;;)
 102    {
 103       if (ec_lo == 0)
 104          adjust = adjust_lo;
 105
 106       else if (ec_hi == 0)
 107          adjust = adjust_hi;
 108
 109       else if (ec_mid == 0)
 110          adjust = adjust_mid;
 111
 112       else if (ec_mid < ec_hi)
 113          adjust = (adjust_mid + adjust_hi)/2;
 114
 115       else if (ec_mid < ec_lo)
 116          adjust = (adjust_mid + adjust_lo)/2;
 117
 118       else
 119       {
 120          fprintf(stderr, "not reached: %u .. %u .. %u\n", ec_lo, ec_mid, ec_hi);
 121          exit(1);
 122       }
 123
 124       /* Calculate the table using the current 'adjust' */
 125       for (i=0; i<=511; ++i)
 126       {
 127          double lo = 255 * sRGB(i << 15);
 128          double hi = 255 * sRGB((i+1) << 15);
 129          unsigned int calc;
 130
 131          calc = nearbyint((lo+adjust) * 256);
 132          if (calc > 65535)
 133          {
 134             fprintf(stderr, "table[%d][0]: overflow %08x (%d)\n", i, calc,
 135                calc);
 136             exit(1);
 137          }
 138          png_sRGB_base[i] = calc;
 139
 140          calc = nearbyint((hi-lo) * 32);
 141          if (calc > 255)
 142          {
 143             fprintf(stderr, "table[%d][1]: overflow %08x (%d)\n", i, calc,
 144                calc);
 145             exit(1);
 146          }
 147          png_sRGB_delta[i] = calc;
 148       }
 149
 150       /* Check the 16-bit linear values alone: */
 151       error_count16 = 0;
 152       for (i16=0; i16 <= 65535; ++i16)
 153       {
 154          unsigned int i = 255*i16;
 155          unsigned int iexact = nearbyint(255*sRGB(i));
 156          unsigned int icalc = PNG_sRGB_FROM_LINEAR(i);
 157
 158          if (icalc != iexact)
 159             ++error_count16;
 160       }
 161
 162       /* Now try changing the adjustment. */
 163       if (ec_lo == 0)
 164          ec_lo = error_count16;
 165
 166       else if (ec_hi == 0)
 167          ec_hi = error_count16;
 168
 169       else if (ec_mid == 0)
 170       {
 171          ec_mid = error_count16;
 172          printf("/* initial error counts: %u .. %u .. %u */\n", ec_lo, ec_mid,
 173             ec_hi);
 174       }
 175
 176       else if (error_count16 < ec_mid)
 177       {
 178          printf("/* adjust (mid ): %f: %u -> %u */\n", adjust, ec_mid,
 179             error_count16);
 180          ec_mid = error_count16;
 181          adjust_mid = adjust;
 182       }
 183
 184       else if (adjust < adjust_mid && error_count16 < ec_lo)
 185       {
 186          printf("/* adjust (low ): %f: %u -> %u */\n", adjust, ec_lo,
 187             error_count16);
 188          ec_lo = error_count16;
 189          adjust_lo = adjust;
 190       }
 191
 192       else if (adjust > adjust_mid && error_count16 < ec_hi)
 193       {
 194          printf("/* adjust (high): %f: %u -> %u */\n", adjust, ec_hi,
 195             error_count16);
 196          ec_hi = error_count16;
 197          adjust_hi = adjust;
 198       }
 199
 200       else
 201       {
 202          adjust = adjust_mid;
 203          printf("/* adjust: %f: %u */\n", adjust, ec_mid);
 204          break;
 205       }
 206    }
 207
 208    /* For each entry in the table try to adjust it to minimize the error count
 209     * in that entry.  Each entry corresponds to 128 input values.
 210     */
 211    for (ibase=0; ibase<65536; ibase+=128)
 212    {
 213       png_uint_16 base = png_sRGB_base[ibase >> 7], trybase = base, ob=base;
 214       png_byte delta = png_sRGB_delta[ibase >> 7], trydelta = delta, od=delta;
 215       unsigned int ecbase = 0, eco;
 216
 217       for (;;)
 218       {
 219          png_sRGB_base[ibase >> 7] = trybase;
 220          png_sRGB_delta[ibase >> 7] = trydelta;
 221
 222          /* Check the 16-bit linear values alone: */
 223          error_count16 = 0;
 224          for (i16=ibase; i16 < ibase+128; ++i16)
 225          {
 226             unsigned int i = 255*i16;
 227             unsigned int iexact = nearbyint(255*sRGB(i));
 228             unsigned int icalc = PNG_sRGB_FROM_LINEAR(i);
 229
 230             if (icalc != iexact)
 231                ++error_count16;
 232          }
 233
 234          if (error_count16 == 0)
 235             break;
 236
 237          if (ecbase == 0)
 238          {
 239             eco = ecbase = error_count16;
 240             ++trybase; /* First test */
 241          }
 242
 243          else if (error_count16 < ecbase)
 244          {
 245             if (trybase > base)
 246             {
 247                base = trybase;
 248                ++trybase;
 249             }
 250             else if (trybase < base)
 251             {
 252                base = trybase;
 253                --trybase;
 254             }
 255             else if (trydelta > delta)
 256             {
 257                delta = trydelta;
 258                ++trydelta;
 259             }
 260             else if (trydelta < delta)
 261             {
 262                delta = trydelta;
 263                --trydelta;
 264             }
 265             else
 266             {
 267                fprintf(stderr, "makesRGB: impossible\n");
 268                exit(1);
 269             }
 270             ecbase = error_count16;
 271          }
 272
 273          else
 274          {
 275             if (trybase > base)
 276                trybase = base-1;
 277             else if (trybase < base)
 278             {
 279                trybase = base;
 280                ++trydelta;
 281             }
 282             else if (trydelta > delta)
 283                trydelta = delta-1;
 284             else if (trydelta < delta)
 285                break; /* end of tests */
 286          }
 287       }
 288
 289       png_sRGB_base[ibase >> 7] = base;
 290       png_sRGB_delta[ibase >> 7] = delta;
 291       if (base != ob || delta != od)
 292       {
 293          printf("/* table[%u]={%u,%u} -> {%u,%u} %u -> %u errors */\n",
 294             ibase>>7, ob, od, base, delta, eco, ecbase);
 295       }
 296       else if (0)
 297          printf("/* table[%u]={%u,%u} %u errors */\n", ibase>>7, ob, od,
 298             ecbase);
 299    }
 300
 301    /* Only do the full (slow) test at the end: */
 302    min_error = -.4999;
 303    max_error = .4999;
 304    error_count = 0;
 305
 306    for (i=0; i <= max_input; ++i)
 307    {
 308       unsigned int iexact = nearbyint(255*sRGB(i));
 309       unsigned int icalc = PNG_sRGB_FROM_LINEAR(i);
 310
 311       if (icalc != iexact)
 312       {
 313          double err = 255*sRGB(i) - icalc;
 314
 315          if (err > (max_error+.001) || err < (min_error-.001))
 316          {
 317             printf(
 318                "/* 0x%08x: exact: %3d, got: %3d [tables: %08x, %08x] (%f) */\n",
 319                i, iexact, icalc, png_sRGB_base[i>>15],
 320                png_sRGB_delta[i>>15], err);
 321          }
 322
 323          ++error_count;
 324          if (err > max_error)
 325             max_error = err;
 326          else if (err < min_error)
 327             min_error = err;
 328       }
 329    }
 330
 331    /* Re-check the 16-bit cases too, including the warning if there is an error
 332     * bigger than 1.
 333     */
 334    error_count16 = 0;
 335    max_error16 = 0;
 336    min_error16 = 0;
 337    for (i16=0; i16 <= 65535; ++i16)
 338    {
 339       unsigned int i = 255*i16;
 340       unsigned int iexact = nearbyint(255*sRGB(i));
 341       unsigned int icalc = PNG_sRGB_FROM_LINEAR(i);
 342
 343       if (icalc != iexact)
 344       {
 345          double err = 255*sRGB(i) - icalc;
 346
 347          ++error_count16;
 348          if (err > max_error16)
 349             max_error16 = err;
 350          else if (err < min_error16)
 351             min_error16 = err;
 352
 353          if (abs(icalc - iexact) > 1)
 354             printf(
 355                "/* 0x%04x: exact: %3d, got: %3d [tables: %08x, %08x] (%f) */\n",
 356                i16, iexact, icalc, png_sRGB_base[i>>15],
 357                png_sRGB_delta[i>>15], err);
 358       }
 359    }
 360
 361    /* Check the round trip for each 8-bit sRGB value. */
 362    for (i16=0; i16 <= 255; ++i16)
 363    {
 364       unsigned int i = 255 * png_sRGB_table[i16];
 365       unsigned int iexact = nearbyint(255*sRGB(i));
 366       unsigned int icalc = PNG_sRGB_FROM_LINEAR(i);
 367
 368       if (i16 != iexact)
 369       {
 370          fprintf(stderr, "8-bit rounding error: %d -> %d\n", i16, iexact);
 371          exit(1);
 372       }
 373
 374       if (icalc != i16)
 375       {
 376          double finv = finvsRGB(i16);
 377
 378          printf("/* 8-bit roundtrip error: %d -> %f -> %d(%f) */\n",
 379             i16, finv, icalc, fsRGB(255*finv));
 380       }
 381    }
 382
 383
 384    printf("/* error: %g - %g, %u (%g%%) of readings inexact */\n",
 385       min_error, max_error, error_count, (100.*error_count)/max_input);
 386    printf("/* 16-bit error: %g - %g, %u (%g%%) of readings inexact */\n",
 387       min_error16, max_error16, error_count16, (100.*error_count16)/65535);
 388
 389    if (!test_only)
 390    {
 391       printf("const png_uint_16 png_sRGB_table[256] =\n{\n   ");
 392       for (i=0; i<255; )
 393       {
 394          do
 395          {
 396             printf("%d,", png_sRGB_table[i++]);
 397          }
 398          while ((i & 0x7) != 0 && i<255);
 399          if (i<255) printf("\n   ");
 400       }
 401       printf("%d\n};\n\n", png_sRGB_table[i]);
 402
 403
 404       printf("const png_uint_16 png_sRGB_base[512] =\n{\n   ");
 405       for (i=0; i<511; )
 406       {
 407          do
 408          {
 409             printf("%d,", png_sRGB_base[i++]);
 410          }
 411          while ((i & 0x7) != 0 && i<511);
 412          if (i<511) printf("\n   ");
 413       }
 414       printf("%d\n};\n\n", png_sRGB_base[i]);
 415
 416       printf("const png_byte png_sRGB_delta[512] =\n{\n   ");
 417       for (i=0; i<511; )
 418       {
 419          do
 420          {
 421             printf("%d,", png_sRGB_delta[i++]);
 422          }
 423          while ((i & 0xf) != 0 && i<511);
 424          if (i<511) printf("\n   ");
 425       }
 426       printf("%d\n};\n\n", png_sRGB_delta[i]);
 427    }
 428
 429    return 0;
 430 }