libpng-1.2.8/pngwutil.c

   1
   2 /* pngwutil.c - utilities to write a PNG file
   3  *
   4  * libpng version 1.2.8 - December 3, 2004
   5  * For conditions of distribution and use, see copyright notice in png.h
   6  * Copyright (c) 1998-2004 Glenn Randers-Pehrson
   7  * (Version 0.96 Copyright (c) 1996, 1997 Andreas Dilger)
   8  * (Version 0.88 Copyright (c) 1995, 1996 Guy Eric Schalnat, Group 42, Inc.)
   9  */
  10
  11 #define PNG_INTERNAL
  12 #include "png.h"
  13 #ifdef PNG_WRITE_SUPPORTED
  14
  15 /* Place a 32-bit number into a buffer in PNG byte order.  We work
  16  * with unsigned numbers for convenience, although one supported
  17  * ancillary chunk uses signed (two's complement) numbers.
  18  */
  19 void /* PRIVATE */
  20 png_save_uint_32(png_bytep buf, png_uint_32 i)
  21 {
  22    buf[0] = (png_byte)((i >> 24) & 0xff);
  23    buf[1] = (png_byte)((i >> 16) & 0xff);
  24    buf[2] = (png_byte)((i >> 8) & 0xff);
  25    buf[3] = (png_byte)(i & 0xff);
  26 }
  27
  28 #if defined(PNG_WRITE_pCAL_SUPPORTED) || defined(PNG_WRITE_oFFs_SUPPORTED)
  29 /* The png_save_int_32 function assumes integers are stored in two's
  30  * complement format.  If this isn't the case, then this routine needs to
  31  * be modified to write data in two's complement format.
  32  */
  33 void /* PRIVATE */
  34 png_save_int_32(png_bytep buf, png_int_32 i)
  35 {
  36    buf[0] = (png_byte)((i >> 24) & 0xff);
  37    buf[1] = (png_byte)((i >> 16) & 0xff);
  38    buf[2] = (png_byte)((i >> 8) & 0xff);
  39    buf[3] = (png_byte)(i & 0xff);
  40 }
  41 #endif
  42
  43 /* Place a 16-bit number into a buffer in PNG byte order.
  44  * The parameter is declared unsigned int, not png_uint_16,
  45  * just to avoid potential problems on pre-ANSI C compilers.
  46  */
  47 void /* PRIVATE */
  48 png_save_uint_16(png_bytep buf, unsigned int i)
  49 {
  50    buf[0] = (png_byte)((i >> 8) & 0xff);
  51    buf[1] = (png_byte)(i & 0xff);
  52 }
  53
  54 /* Write a PNG chunk all at once.  The type is an array of ASCII characters
  55  * representing the chunk name.  The array must be at least 4 bytes in
  56  * length, and does not need to be null terminated.  To be safe, pass the
  57  * pre-defined chunk names here, and if you need a new one, define it
  58  * where the others are defined.  The length is the length of the data.
  59  * All the data must be present.  If that is not possible, use the
  60  * png_write_chunk_start(), png_write_chunk_data(), and png_write_chunk_end()
  61  * functions instead.
  62  */
  63 void PNGAPI
  64 png_write_chunk(png_structp png_ptr, png_bytep chunk_name,
  65    png_bytep data, png_size_t length)
  66 {
  67    png_write_chunk_start(png_ptr, chunk_name, (png_uint_32)length);
  68    png_write_chunk_data(png_ptr, data, length);
  69    png_write_chunk_end(png_ptr);
  70 }
  71
  72 /* Write the start of a PNG chunk.  The type is the chunk type.
  73  * The total_length is the sum of the lengths of all the data you will be
  74  * passing in png_write_chunk_data().
  75  */
  76 void PNGAPI
  77 png_write_chunk_start(png_structp png_ptr, png_bytep chunk_name,
  78    png_uint_32 length)
  79 {
  80    png_byte buf[4];
  81    png_debug2(0, "Writing %s chunk (%lu bytes)\n", chunk_name, length);
  82
  83    /* write the length */
  84    png_save_uint_32(buf, length);
  85    png_write_data(png_ptr, buf, (png_size_t)4);
  86
  87    /* write the chunk name */
  88    png_write_data(png_ptr, chunk_name, (png_size_t)4);
  89    /* reset the crc and run it over the chunk name */
  90    png_reset_crc(png_ptr);
  91    png_calculate_crc(png_ptr, chunk_name, (png_size_t)4);
  92 }
  93
  94 /* Write the data of a PNG chunk started with png_write_chunk_start().
  95  * Note that multiple calls to this function are allowed, and that the
  96  * sum of the lengths from these calls *must* add up to the total_length
  97  * given to png_write_chunk_start().
  98  */
  99 void PNGAPI
 100 png_write_chunk_data(png_structp png_ptr, png_bytep data, png_size_t length)
 101 {
 102    /* write the data, and run the CRC over it */
 103    if (data != NULL && length > 0)
 104    {
 105       png_calculate_crc(png_ptr, data, length);
 106       png_write_data(png_ptr, data, length);
 107    }
 108 }
 109
 110 /* Finish a chunk started with png_write_chunk_start(). */
 111 void PNGAPI
 112 png_write_chunk_end(png_structp png_ptr)
 113 {
 114    png_byte buf[4];
 115
 116    /* write the crc */
 117    png_save_uint_32(buf, png_ptr->crc);
 118
 119    png_write_data(png_ptr, buf, (png_size_t)4);
 120 }
 121
 122 /* Simple function to write the signature.  If we have already written
 123  * the magic bytes of the signature, or more likely, the PNG stream is
 124  * being embedded into another stream and doesn't need its own signature,
 125  * we should call png_set_sig_bytes() to tell libpng how many of the
 126  * bytes have already been written.
 127  */
 128 void /* PRIVATE */
 129 png_write_sig(png_structp png_ptr)
 130 {
 131    png_byte png_signature[8] = {137, 80, 78, 71, 13, 10, 26, 10};
 132    /* write the rest of the 8 byte signature */
 133    png_write_data(png_ptr, &png_signature[png_ptr->sig_bytes],
 134       (png_size_t)8 - png_ptr->sig_bytes);
 135    if(png_ptr->sig_bytes < 3)
 136       png_ptr->mode |= PNG_HAVE_PNG_SIGNATURE;
 137 }
 138
 139 #if defined(PNG_WRITE_TEXT_SUPPORTED) || defined(PNG_WRITE_iCCP_SUPPORTED)
 140 /*
 141  * This pair of functions encapsulates the operation of (a) compressing a
 142  * text string, and (b) issuing it later as a series of chunk data writes.
 143  * The compression_state structure is shared context for these functions
 144  * set up by the caller in order to make the whole mess thread-safe.
 145  */
 146
 147 typedef struct
 148 {
 149     char *input;   /* the uncompressed input data */
 150     int input_len;   /* its length */
 151     int num_output_ptr; /* number of output pointers used */
 152     int max_output_ptr; /* size of output_ptr */
 153     png_charpp output_ptr; /* array of pointers to output */
 154 } compression_state;
 155
 156 /* compress given text into storage in the png_ptr structure */
 157 static int /* PRIVATE */
 158 png_text_compress(png_structp png_ptr,
 159         png_charp text, png_size_t text_len, int compression,
 160         compression_state *comp)
 161 {
 162    int ret;
 163
 164    comp->num_output_ptr = comp->max_output_ptr = 0;
 165    comp->output_ptr = NULL;
 166    comp->input = NULL;
 167
 168    /* we may just want to pass the text right through */
 169    if (compression == PNG_TEXT_COMPRESSION_NONE)
 170    {
 171        comp->input = text;
 172        comp->input_len = text_len;
 173        return((int)text_len);
 174    }
 175
 176    if (compression >= PNG_TEXT_COMPRESSION_LAST)
 177    {
 178 #if !defined(PNG_NO_STDIO) && !defined(_WIN32_WCE)
 179       char msg[50];
 180       sprintf(msg, "Unknown compression type %d", compression);
 181       png_warning(png_ptr, msg);
 182 #else
 183       png_warning(png_ptr, "Unknown compression type");
 184 #endif
 185    }
 186
 187    /* We can't write the chunk until we find out how much data we have,
 188     * which means we need to run the compressor first and save the
 189     * output.  This shouldn't be a problem, as the vast majority of
 190     * comments should be reasonable, but we will set up an array of
 191     * malloc'd pointers to be sure.
 192     *
 193     * If we knew the application was well behaved, we could simplify this
 194     * greatly by assuming we can always malloc an output buffer large
 195     * enough to hold the compressed text ((1001 * text_len / 1000) + 12)
 196     * and malloc this directly.  The only time this would be a bad idea is
 197     * if we can't malloc more than 64K and we have 64K of random input
 198     * data, or if the input string is incredibly large (although this
 199     * wouldn't cause a failure, just a slowdown due to swapping).
 200     */
 201
 202    /* set up the compression buffers */
 203    png_ptr->zstream.avail_in = (uInt)text_len;
 204    png_ptr->zstream.next_in = (Bytef *)text;
 205    png_ptr->zstream.avail_out = (uInt)png_ptr->zbuf_size;
 206    png_ptr->zstream.next_out = (Bytef *)png_ptr->zbuf;
 207
 208    /* this is the same compression loop as in png_write_row() */
 209    do
 210    {
 211       /* compress the data */
 212       ret = deflate(&png_ptr->zstream, Z_NO_FLUSH);
 213       if (ret != Z_OK)
 214       {
 215          /* error */
 216          if (png_ptr->zstream.msg != NULL)
 217             png_error(png_ptr, png_ptr->zstream.msg);
 218          else
 219             png_error(png_ptr, "zlib error");
 220       }
 221       /* check to see if we need more room */
 222       if (!(png_ptr->zstream.avail_out))
 223       {
 224          /* make sure the output array has room */
 225          if (comp->num_output_ptr >= comp->max_output_ptr)
 226          {
 227             int old_max;
 228
 229             old_max = comp->max_output_ptr;
 230             comp->max_output_ptr = comp->num_output_ptr + 4;
 231             if (comp->output_ptr != NULL)
 232             {
 233                png_charpp old_ptr;
 234
 235                old_ptr = comp->output_ptr;
 236                comp->output_ptr = (png_charpp)png_malloc(png_ptr,
 237                   (png_uint_32)(comp->max_output_ptr *
 238                   png_sizeof (png_charpp)));
 239                png_memcpy(comp->output_ptr, old_ptr, old_max
 240                   * png_sizeof (png_charp));
 241                png_free(png_ptr, old_ptr);
 242             }
 243             else
 244                comp->output_ptr = (png_charpp)png_malloc(png_ptr,
 245                   (png_uint_32)(comp->max_output_ptr *
 246                   png_sizeof (png_charp)));
 247          }
 248
 249          /* save the data */
 250          comp->output_ptr[comp->num_output_ptr] = (png_charp)png_malloc(png_ptr,
 251             (png_uint_32)png_ptr->zbuf_size);
 252          png_memcpy(comp->output_ptr[comp->num_output_ptr], png_ptr->zbuf,
 253             png_ptr->zbuf_size);
 254          comp->num_output_ptr++;
 255
 256          /* and reset the buffer */
 257          png_ptr->zstream.avail_out = (uInt)png_ptr->zbuf_size;
 258          png_ptr->zstream.next_out = png_ptr->zbuf;
 259       }
 260    /* continue until we don't have any more to compress */
 261    } while (png_ptr->zstream.avail_in);
 262
 263    /* finish the compression */
 264    do
 265    {
 266       /* tell zlib we are finished */
 267       ret = deflate(&png_ptr->zstream, Z_FINISH);
 268
 269       if (ret == Z_OK)
 270       {
 271          /* check to see if we need more room */
 272          if (!(png_ptr->zstream.avail_out))
 273          {
 274             /* check to make sure our output array has room */
 275             if (comp->num_output_ptr >= comp->max_output_ptr)
 276             {
 277                int old_max;
 278
 279                old_max = comp->max_output_ptr;
 280                comp->max_output_ptr = comp->num_output_ptr + 4;
 281                if (comp->output_ptr != NULL)
 282                {
 283                   png_charpp old_ptr;
 284
 285                   old_ptr = comp->output_ptr;
 286                   /* This could be optimized to realloc() */
 287                   comp->output_ptr = (png_charpp)png_malloc(png_ptr,
 288                      (png_uint_32)(comp->max_output_ptr *
 289                      png_sizeof (png_charpp)));
 290                   png_memcpy(comp->output_ptr, old_ptr,
 291                      old_max * png_sizeof (png_charp));
 292                   png_free(png_ptr, old_ptr);
 293                }
 294                else
 295                   comp->output_ptr = (png_charpp)png_malloc(png_ptr,
 296                      (png_uint_32)(comp->max_output_ptr *
 297                      png_sizeof (png_charp)));
 298             }
 299
 300             /* save off the data */
 301             comp->output_ptr[comp->num_output_ptr] =
 302                (png_charp)png_malloc(png_ptr, (png_uint_32)png_ptr->zbuf_size);
 303             png_memcpy(comp->output_ptr[comp->num_output_ptr], png_ptr->zbuf,
 304                png_ptr->zbuf_size);
 305             comp->num_output_ptr++;
 306
 307             /* and reset the buffer pointers */
 308             png_ptr->zstream.avail_out = (uInt)png_ptr->zbuf_size;
 309             png_ptr->zstream.next_out = png_ptr->zbuf;
 310          }
 311       }
 312       else if (ret != Z_STREAM_END)
 313       {
 314          /* we got an error */
 315          if (png_ptr->zstream.msg != NULL)
 316             png_error(png_ptr, png_ptr->zstream.msg);
 317          else
 318             png_error(png_ptr, "zlib error");
 319       }
 320    } while (ret != Z_STREAM_END);
 321
 322    /* text length is number of buffers plus last buffer */
 323    text_len = png_ptr->zbuf_size * comp->num_output_ptr;
 324    if (png_ptr->zstream.avail_out < png_ptr->zbuf_size)
 325       text_len += png_ptr->zbuf_size - (png_size_t)png_ptr->zstream.avail_out;
 326
 327    return((int)text_len);
 328 }
 329
 330 /* ship the compressed text out via chunk writes */
 331 static void /* PRIVATE */
 332 png_write_compressed_data_out(png_structp png_ptr, compression_state *comp)
 333 {
 334    int i;
 335
 336    /* handle the no-compression case */
 337    if (comp->input)
 338    {
 339        png_write_chunk_data(png_ptr, (png_bytep)comp->input,
 340                             (png_size_t)comp->input_len);
 341        return;
 342    }
 343
 344    /* write saved output buffers, if any */
 345    for (i = 0; i < comp->num_output_ptr; i++)
 346    {
 347       png_write_chunk_data(png_ptr,(png_bytep)comp->output_ptr[i],
 348          png_ptr->zbuf_size);
 349       png_free(png_ptr, comp->output_ptr[i]);
 350       comp->output_ptr[i]=NULL;
 351    }
 352    if (comp->max_output_ptr != 0)
 353       png_free(png_ptr, comp->output_ptr);
 354       comp->output_ptr=NULL;
 355    /* write anything left in zbuf */
 356    if (png_ptr->zstream.avail_out < (png_uint_32)png_ptr->zbuf_size)
 357       png_write_chunk_data(png_ptr, png_ptr->zbuf,
 358          png_ptr->zbuf_size - png_ptr->zstream.avail_out);
 359
 360    /* reset zlib for another zTXt/iTXt or image data */
 361    deflateReset(&png_ptr->zstream);
 362    png_ptr->zstream.data_type = Z_BINARY;
 363 }
 364 #endif
 365
 366 /* Write the IHDR chunk, and update the png_struct with the necessary
 367  * information.  Note that the rest of this code depends upon this
 368  * information being correct.
 369  */
 370 void /* PRIVATE */
 371 png_write_IHDR(png_structp png_ptr, png_uint_32 width, png_uint_32 height,
 372    int bit_depth, int color_type, int compression_type, int filter_type,
 373    int interlace_type)
 374 {
 375 #ifdef PNG_USE_LOCAL_ARRAYS
 376    PNG_IHDR;
 377 #endif
 378    png_byte buf[13]; /* buffer to store the IHDR info */
 379
 380    png_debug(1, "in png_write_IHDR\n");
 381    /* Check that we have valid input data from the application info */
 382    switch (color_type)
 383    {
 384       case PNG_COLOR_TYPE_GRAY:
 385          switch (bit_depth)
 386          {
 387             case 1:
 388             case 2:
 389             case 4:
 390             case 8:
 391             case 16: png_ptr->channels = 1; break;
 392             default: png_error(png_ptr,"Invalid bit depth for grayscale image");
 393          }
 394          break;
 395       case PNG_COLOR_TYPE_RGB:
 396          if (bit_depth != 8 && bit_depth != 16)
 397             png_error(png_ptr, "Invalid bit depth for RGB image");
 398          png_ptr->channels = 3;
 399          break;
 400       case PNG_COLOR_TYPE_PALETTE:
 401          switch (bit_depth)
 402          {
 403             case 1:
 404             case 2:
 405             case 4:
 406             case 8: png_ptr->channels = 1; break;
 407             default: png_error(png_ptr, "Invalid bit depth for paletted image");
 408          }
 409          break;
 410       case PNG_COLOR_TYPE_GRAY_ALPHA:
 411          if (bit_depth != 8 && bit_depth != 16)
 412             png_error(png_ptr, "Invalid bit depth for grayscale+alpha image");
 413          png_ptr->channels = 2;
 414          break;
 415       case PNG_COLOR_TYPE_RGB_ALPHA:
 416          if (bit_depth != 8 && bit_depth != 16)
 417             png_error(png_ptr, "Invalid bit depth for RGBA image");
 418          png_ptr->channels = 4;
 419          break;
 420       default:
 421          png_error(png_ptr, "Invalid image color type specified");
 422    }
 423
 424    if (compression_type != PNG_COMPRESSION_TYPE_BASE)
 425    {
 426       png_warning(png_ptr, "Invalid compression type specified");
 427       compression_type = PNG_COMPRESSION_TYPE_BASE;
 428    }
 429
 430    /* Write filter_method 64 (intrapixel differencing) only if
 431     * 1. Libpng was compiled with PNG_MNG_FEATURES_SUPPORTED and
 432     * 2. Libpng did not write a PNG signature (this filter_method is only
 433     *    used in PNG datastreams that are embedded in MNG datastreams) and
 434     * 3. The application called png_permit_mng_features with a mask that
 435     *    included PNG_FLAG_MNG_FILTER_64 and
 436     * 4. The filter_method is 64 and
 437     * 5. The color_type is RGB or RGBA
 438     */
 439    if (
 440 #if defined(PNG_MNG_FEATURES_SUPPORTED)
 441       !((png_ptr->mng_features_permitted & PNG_FLAG_MNG_FILTER_64) &&
 442       ((png_ptr->mode&PNG_HAVE_PNG_SIGNATURE) == 0) &&
 443       (color_type == PNG_COLOR_TYPE_RGB ||
 444        color_type == PNG_COLOR_TYPE_RGB_ALPHA) &&
 445       (filter_type == PNG_INTRAPIXEL_DIFFERENCING)) &&
 446 #endif
 447       filter_type != PNG_FILTER_TYPE_BASE)
 448    {
 449       png_warning(png_ptr, "Invalid filter type specified");
 450       filter_type = PNG_FILTER_TYPE_BASE;
 451    }
 452
 453 #ifdef PNG_WRITE_INTERLACING_SUPPORTED
 454    if (interlace_type != PNG_INTERLACE_NONE &&
 455       interlace_type != PNG_INTERLACE_ADAM7)
 456    {
 457       png_warning(png_ptr, "Invalid interlace type specified");
 458       interlace_type = PNG_INTERLACE_ADAM7;
 459    }
 460 #else
 461    interlace_type=PNG_INTERLACE_NONE;
 462 #endif
 463
 464    /* save off the relevent information */
 465    png_ptr->bit_depth = (png_byte)bit_depth;
 466    png_ptr->color_type = (png_byte)color_type;
 467    png_ptr->interlaced = (png_byte)interlace_type;
 468 #if defined(PNG_MNG_FEATURES_SUPPORTED)
 469    png_ptr->filter_type = (png_byte)filter_type;
 470 #endif
 471    png_ptr->compression_type = (png_byte)compression_type;
 472    png_ptr->width = width;
 473    png_ptr->height = height;
 474
 475    png_ptr->pixel_depth = (png_byte)(bit_depth * png_ptr->channels);
 476    png_ptr->rowbytes = PNG_ROWBYTES(png_ptr->pixel_depth, width);
 477    /* set the usr info, so any transformations can modify it */
 478    png_ptr->usr_width = png_ptr->width;
 479    png_ptr->usr_bit_depth = png_ptr->bit_depth;
 480    png_ptr->usr_channels = png_ptr->channels;
 481
 482    /* pack the header information into the buffer */
 483    png_save_uint_32(buf, width);
 484    png_save_uint_32(buf + 4, height);
 485    buf[8] = (png_byte)bit_depth;
 486    buf[9] = (png_byte)color_type;
 487    buf[10] = (png_byte)compression_type;
 488    buf[11] = (png_byte)filter_type;
 489    buf[12] = (png_byte)interlace_type;
 490
 491    /* write the chunk */
 492    png_write_chunk(png_ptr, (png_bytep)png_IHDR, buf, (png_size_t)13);
 493
 494    /* initialize zlib with PNG info */
 495    png_ptr->zstream.zalloc = png_zalloc;
 496    png_ptr->zstream.zfree = png_zfree;
 497    png_ptr->zstream.opaque = (voidpf)png_ptr;
 498    if (!(png_ptr->do_filter))
 499    {
 500       if (png_ptr->color_type == PNG_COLOR_TYPE_PALETTE ||
 501          png_ptr->bit_depth < 8)
 502          png_ptr->do_filter = PNG_FILTER_NONE;
 503       else
 504          png_ptr->do_filter = PNG_ALL_FILTERS;
 505    }
 506    if (!(png_ptr->flags & PNG_FLAG_ZLIB_CUSTOM_STRATEGY))
 507    {
 508       if (png_ptr->do_filter != PNG_FILTER_NONE)
 509          png_ptr->zlib_strategy = Z_FILTERED;
 510       else
 511          png_ptr->zlib_strategy = Z_DEFAULT_STRATEGY;
 512    }
 513    if (!(png_ptr->flags & PNG_FLAG_ZLIB_CUSTOM_LEVEL))
 514       png_ptr->zlib_level = Z_DEFAULT_COMPRESSION;
 515    if (!(png_ptr->flags & PNG_FLAG_ZLIB_CUSTOM_MEM_LEVEL))
 516       png_ptr->zlib_mem_level = 8;
 517    if (!(png_ptr->flags & PNG_FLAG_ZLIB_CUSTOM_WINDOW_BITS))
 518       png_ptr->zlib_window_bits = 15;
 519    if (!(png_ptr->flags & PNG_FLAG_ZLIB_CUSTOM_METHOD))
 520       png_ptr->zlib_method = 8;
 521    deflateInit2(&png_ptr->zstream, png_ptr->zlib_level,
 522       png_ptr->zlib_method, png_ptr->zlib_window_bits,
 523       png_ptr->zlib_mem_level, png_ptr->zlib_strategy);
 524    png_ptr->zstream.next_out = png_ptr->zbuf;
 525    png_ptr->zstream.avail_out = (uInt)png_ptr->zbuf_size;
 526    /* libpng is not interested in zstream.data_type */
 527    /* set it to a predefined value, to avoid its evaluation inside zlib */
 528    png_ptr->zstream.data_type = Z_BINARY;
 529
 530    png_ptr->mode = PNG_HAVE_IHDR;
 531 }
 532
 533 /* write the palette.  We are careful not to trust png_color to be in the
 534  * correct order for PNG, so people can redefine it to any convenient
 535  * structure.
 536  */
 537 void /* PRIVATE */
 538 png_write_PLTE(png_structp png_ptr, png_colorp palette, png_uint_32 num_pal)
 539 {
 540 #ifdef PNG_USE_LOCAL_ARRAYS
 541    PNG_PLTE;
 542 #endif
 543    png_uint_32 i;
 544    png_colorp pal_ptr;
 545    png_byte buf[3];
 546
 547    png_debug(1, "in png_write_PLTE\n");
 548    if ((
 549 #if defined(PNG_MNG_FEATURES_SUPPORTED)
 550         !(png_ptr->mng_features_permitted & PNG_FLAG_MNG_EMPTY_PLTE) &&
 551 #endif
 552         num_pal == 0) || num_pal > 256)
 553    {
 554      if (png_ptr->color_type == PNG_COLOR_TYPE_PALETTE)
 555      {
 556         png_error(png_ptr, "Invalid number of colors in palette");
 557      }
 558      else
 559      {
 560         png_warning(png_ptr, "Invalid number of colors in palette");
 561         return;
 562      }
 563    }
 564
 565    if (!(png_ptr->color_type&PNG_COLOR_MASK_COLOR))
 566    {
 567       png_warning(png_ptr,
 568         "Ignoring request to write a PLTE chunk in grayscale PNG");
 569       return;
 570    }
 571
 572    png_ptr->num_palette = (png_uint_16)num_pal;
 573    png_debug1(3, "num_palette = %d\n", png_ptr->num_palette);
 574
 575    png_write_chunk_start(png_ptr, (png_bytep)png_PLTE, num_pal * 3);
 576 #ifndef PNG_NO_POINTER_INDEXING
 577    for (i = 0, pal_ptr = palette; i < num_pal; i++, pal_ptr++)
 578    {
 579       buf[0] = pal_ptr->red;
 580       buf[1] = pal_ptr->green;
 581       buf[2] = pal_ptr->blue;
 582       png_write_chunk_data(png_ptr, buf, (png_size_t)3);
 583    }
 584 #else
 585    /* This is a little slower but some buggy compilers need to do this instead */
 586    pal_ptr=palette;
 587    for (i = 0; i < num_pal; i++)
 588    {
 589       buf[0] = pal_ptr[i].red;
 590       buf[1] = pal_ptr[i].green;
 591       buf[2] = pal_ptr[i].blue;
 592       png_write_chunk_data(png_ptr, buf, (png_size_t)3);
 593    }
 594 #endif
 595    png_write_chunk_end(png_ptr);
 596    png_ptr->mode |= PNG_HAVE_PLTE;
 597 }
 598
 599 /* write an IDAT chunk */
 600 void /* PRIVATE */
 601 png_write_IDAT(png_structp png_ptr, png_bytep data, png_size_t length)
 602 {
 603 #ifdef PNG_USE_LOCAL_ARRAYS
 604    PNG_IDAT;
 605 #endif
 606    png_debug(1, "in png_write_IDAT\n");
 607
 608    /* Optimize the CMF field in the zlib stream. */
 609    /* This hack of the zlib stream is compliant to the stream specification. */
 610    if (!(png_ptr->mode & PNG_HAVE_IDAT) &&
 611        png_ptr->compression_type == PNG_COMPRESSION_TYPE_BASE)
 612    {
 613       unsigned int z_cmf = data[0];  /* zlib compression method and flags */
 614       if ((z_cmf & 0x0f) == 8 && (z_cmf & 0xf0) <= 0x70)
 615       {
 616          /* Avoid memory underflows and multiplication overflows. */
 617          /* The conditions below are practically always satisfied;
 618             however, they still must be checked. */
 619          if (length >= 2 &&
 620              png_ptr->height < 16384 && png_ptr->width < 16384)
 621          {
 622             png_uint_32 uncompressed_idat_size = png_ptr->height *
 623                ((png_ptr->width *
 624                png_ptr->channels * png_ptr->bit_depth + 15) >> 3);
 625             unsigned int z_cinfo = z_cmf >> 4;
 626             unsigned int half_z_window_size = 1 << (z_cinfo + 7);
 627             while (uncompressed_idat_size <= half_z_window_size &&
 628                    half_z_window_size >= 256)
 629             {
 630                z_cinfo--;
 631                half_z_window_size >>= 1;
 632             }
 633             z_cmf = (z_cmf & 0x0f) | (z_cinfo << 4);
 634             if (data[0] != (png_byte)z_cmf)
 635             {
 636                data[0] = (png_byte)z_cmf;
 637                data[1] &= 0xe0;
 638                data[1] += (png_byte)(0x1f - ((z_cmf << 8) + data[1]) % 0x1f);
 639             }
 640          }
 641       }
 642       else
 643          png_error(png_ptr,
 644             "Invalid zlib compression method or flags in IDAT");
 645    }
 646
 647    png_write_chunk(png_ptr, (png_bytep)png_IDAT, data, length);
 648    png_ptr->mode |= PNG_HAVE_IDAT;
 649 }
 650
 651 /* write an IEND chunk */
 652 void /* PRIVATE */
 653 png_write_IEND(png_structp png_ptr)
 654 {
 655 #ifdef PNG_USE_LOCAL_ARRAYS
 656    PNG_IEND;
 657 #endif
 658    png_debug(1, "in png_write_IEND\n");
 659    png_write_chunk(png_ptr, (png_bytep)png_IEND, png_bytep_NULL,
 660      (png_size_t)0);
 661    png_ptr->mode |= PNG_HAVE_IEND;
 662 }
 663
 664 #if defined(PNG_WRITE_gAMA_SUPPORTED)
 665 /* write a gAMA chunk */
 666 #ifdef PNG_FLOATING_POINT_SUPPORTED
 667 void /* PRIVATE */
 668 png_write_gAMA(png_structp png_ptr, double file_gamma)
 669 {
 670 #ifdef PNG_USE_LOCAL_ARRAYS
 671    PNG_gAMA;
 672 #endif
 673    png_uint_32 igamma;
 674    png_byte buf[4];
 675
 676    png_debug(1, "in png_write_gAMA\n");
 677    /* file_gamma is saved in 1/100,000ths */
 678    igamma = (png_uint_32)(file_gamma * 100000.0 + 0.5);
 679    png_save_uint_32(buf, igamma);
 680    png_write_chunk(png_ptr, (png_bytep)png_gAMA, buf, (png_size_t)4);
 681 }
 682 #endif
 683 #ifdef PNG_FIXED_POINT_SUPPORTED
 684 void /* PRIVATE */
 685 png_write_gAMA_fixed(png_structp png_ptr, png_fixed_point file_gamma)
 686 {
 687 #ifdef PNG_USE_LOCAL_ARRAYS
 688    PNG_gAMA;
 689 #endif
 690    png_byte buf[4];
 691
 692    png_debug(1, "in png_write_gAMA\n");
 693    /* file_gamma is saved in 1/100,000ths */
 694    png_save_uint_32(buf, (png_uint_32)file_gamma);
 695    png_write_chunk(png_ptr, (png_bytep)png_gAMA, buf, (png_size_t)4);
 696 }
 697 #endif
 698 #endif
 699
 700 #if defined(PNG_WRITE_sRGB_SUPPORTED)
 701 /* write a sRGB chunk */
 702 void /* PRIVATE */
 703 png_write_sRGB(png_structp png_ptr, int srgb_intent)
 704 {
 705 #ifdef PNG_USE_LOCAL_ARRAYS
 706    PNG_sRGB;
 707 #endif
 708    png_byte buf[1];
 709
 710    png_debug(1, "in png_write_sRGB\n");
 711    if(srgb_intent >= PNG_sRGB_INTENT_LAST)
 712          png_warning(png_ptr,
 713             "Invalid sRGB rendering intent specified");
 714    buf[0]=(png_byte)srgb_intent;
 715    png_write_chunk(png_ptr, (png_bytep)png_sRGB, buf, (png_size_t)1);
 716 }
 717 #endif
 718
 719 #if defined(PNG_WRITE_iCCP_SUPPORTED)
 720 /* write an iCCP chunk */
 721 void /* PRIVATE */
 722 png_write_iCCP(png_structp png_ptr, png_charp name, int compression_type,
 723    png_charp profile, int profile_len)
 724 {
 725 #ifdef PNG_USE_LOCAL_ARRAYS
 726    PNG_iCCP;
 727 #endif
 728    png_size_t name_len;
 729    png_charp new_name;
 730    compression_state comp;
 731
 732    png_debug(1, "in png_write_iCCP\n");
 733    if (name == NULL || (name_len = png_check_keyword(png_ptr, name,
 734       &new_name)) == 0)
 735    {
 736       png_warning(png_ptr, "Empty keyword in iCCP chunk");
 737       return;
 738    }
 739
 740    if (compression_type != PNG_COMPRESSION_TYPE_BASE)
 741       png_warning(png_ptr, "Unknown compression type in iCCP chunk");
 742
 743    if (profile == NULL)
 744       profile_len = 0;
 745
 746    if (profile_len)
 747        profile_len = png_text_compress(png_ptr, profile, (png_size_t)profile_len,
 748           PNG_COMPRESSION_TYPE_BASE, &comp);
 749
 750    /* make sure we include the NULL after the name and the compression type */
 751    png_write_chunk_start(png_ptr, (png_bytep)png_iCCP,
 752           (png_uint_32)name_len+profile_len+2);
 753    new_name[name_len+1]=0x00;
 754    png_write_chunk_data(png_ptr, (png_bytep)new_name, name_len + 2);
 755
 756    if (profile_len)
 757       png_write_compressed_data_out(png_ptr, &comp);
 758
 759    png_write_chunk_end(png_ptr);
 760    png_free(png_ptr, new_name);
 761 }
 762 #endif
 763
 764 #if defined(PNG_WRITE_sPLT_SUPPORTED)
 765 /* write a sPLT chunk */
 766 void /* PRIVATE */
 767 png_write_sPLT(png_structp png_ptr, png_sPLT_tp spalette)
 768 {
 769 #ifdef PNG_USE_LOCAL_ARRAYS
 770    PNG_sPLT;
 771 #endif
 772    png_size_t name_len;
 773    png_charp new_name;
 774    png_byte entrybuf[10];
 775    int entry_size = (spalette->depth == 8 ? 6 : 10);
 776    int palette_size = entry_size * spalette->nentries;
 777    png_sPLT_entryp ep;
 778 #ifdef PNG_NO_POINTER_INDEXING
 779    int i;
 780 #endif
 781
 782    png_debug(1, "in png_write_sPLT\n");
 783    if (spalette->name == NULL || (name_len = png_check_keyword(png_ptr,
 784       spalette->name, &new_name))==0)
 785    {
 786       png_warning(png_ptr, "Empty keyword in sPLT chunk");
 787       return;
 788    }
 789
 790    /* make sure we include the NULL after the name */
 791    png_write_chunk_start(png_ptr, (png_bytep)png_sPLT,
 792           (png_uint_32)(name_len + 2 + palette_size));
 793    png_write_chunk_data(png_ptr, (png_bytep)new_name, name_len + 1);
 794    png_write_chunk_data(png_ptr, (png_bytep)&spalette->depth, 1);
 795
 796    /* loop through each palette entry, writing appropriately */
 797 #ifndef PNG_NO_POINTER_INDEXING
 798    for (ep = spalette->entries; ep<spalette->entries+spalette->nentries; ep++)
 799    {
 800        if (spalette->depth == 8)
 801        {
 802            entrybuf[0] = (png_byte)ep->red;
 803            entrybuf[1] = (png_byte)ep->green;
 804            entrybuf[2] = (png_byte)ep->blue;
 805            entrybuf[3] = (png_byte)ep->alpha;
 806            png_save_uint_16(entrybuf + 4, ep->frequency);
 807        }
 808        else
 809        {
 810            png_save_uint_16(entrybuf + 0, ep->red);
 811            png_save_uint_16(entrybuf + 2, ep->green);
 812            png_save_uint_16(entrybuf + 4, ep->blue);
 813            png_save_uint_16(entrybuf + 6, ep->alpha);
 814            png_save_uint_16(entrybuf + 8, ep->frequency);
 815        }
 816        png_write_chunk_data(png_ptr, entrybuf, (png_size_t)entry_size);
 817    }
 818 #else
 819    ep=spalette->entries;
 820    for (i=0; i>spalette->nentries; i++)
 821    {
 822        if (spalette->depth == 8)
 823        {
 824            entrybuf[0] = (png_byte)ep[i].red;
 825            entrybuf[1] = (png_byte)ep[i].green;
 826            entrybuf[2] = (png_byte)ep[i].blue;
 827            entrybuf[3] = (png_byte)ep[i].alpha;
 828            png_save_uint_16(entrybuf + 4, ep[i].frequency);
 829        }
 830        else
 831        {
 832            png_save_uint_16(entrybuf + 0, ep[i].red);
 833            png_save_uint_16(entrybuf + 2, ep[i].green);
 834            png_save_uint_16(entrybuf + 4, ep[i].blue);
 835            png_save_uint_16(entrybuf + 6, ep[i].alpha);
 836            png_save_uint_16(entrybuf + 8, ep[i].frequency);
 837        }
 838        png_write_chunk_data(png_ptr, entrybuf, entry_size);
 839    }
 840 #endif
 841
 842    png_write_chunk_end(png_ptr);
 843    png_free(png_ptr, new_name);
 844 }
 845 #endif
 846
 847 #if defined(PNG_WRITE_sBIT_SUPPORTED)
 848 /* write the sBIT chunk */
 849 void /* PRIVATE */
 850 png_write_sBIT(png_structp png_ptr, png_color_8p sbit, int color_type)
 851 {
 852 #ifdef PNG_USE_LOCAL_ARRAYS
 853    PNG_sBIT;
 854 #endif
 855    png_byte buf[4];
 856    png_size_t size;
 857
 858    png_debug(1, "in png_write_sBIT\n");
 859    /* make sure we don't depend upon the order of PNG_COLOR_8 */
 860    if (color_type & PNG_COLOR_MASK_COLOR)
 861    {
 862       png_byte maxbits;
 863
 864       maxbits = (png_byte)(color_type==PNG_COLOR_TYPE_PALETTE ? 8 :
 865                 png_ptr->usr_bit_depth);
 866       if (sbit->red == 0 || sbit->red > maxbits ||
 867           sbit->green == 0 || sbit->green > maxbits ||
 868           sbit->blue == 0 || sbit->blue > maxbits)
 869       {
 870          png_warning(png_ptr, "Invalid sBIT depth specified");
 871          return;
 872       }
 873       buf[0] = sbit->red;
 874       buf[1] = sbit->green;
 875       buf[2] = sbit->blue;
 876       size = 3;
 877    }
 878    else
 879    {
 880       if (sbit->gray == 0 || sbit->gray > png_ptr->usr_bit_depth)
 881       {
 882          png_warning(png_ptr, "Invalid sBIT depth specified");
 883          return;
 884       }
 885       buf[0] = sbit->gray;
 886       size = 1;
 887    }
 888
 889    if (color_type & PNG_COLOR_MASK_ALPHA)
 890    {
 891       if (sbit->alpha == 0 || sbit->alpha > png_ptr->usr_bit_depth)
 892       {
 893          png_warning(png_ptr, "Invalid sBIT depth specified");
 894          return;
 895       }
 896       buf[size++] = sbit->alpha;
 897    }
 898
 899    png_write_chunk(png_ptr, (png_bytep)png_sBIT, buf, size);
 900 }
 901 #endif
 902
 903 #if defined(PNG_WRITE_cHRM_SUPPORTED)
 904 /* write the cHRM chunk */
 905 #ifdef PNG_FLOATING_POINT_SUPPORTED
 906 void /* PRIVATE */
 907 png_write_cHRM(png_structp png_ptr, double white_x, double white_y,
 908    double red_x, double red_y, double green_x, double green_y,
 909    double blue_x, double blue_y)
 910 {
 911 #ifdef PNG_USE_LOCAL_ARRAYS
 912    PNG_cHRM;
 913 #endif
 914    png_byte buf[32];
 915    png_uint_32 itemp;
 916
 917    png_debug(1, "in png_write_cHRM\n");
 918    /* each value is saved in 1/100,000ths */
 919    if (white_x < 0 || white_x > 0.8 || white_y < 0 || white_y > 0.8 ||
 920        white_x + white_y > 1.0)
 921    {
 922       png_warning(png_ptr, "Invalid cHRM white point specified");
 923 #if !defined(PNG_NO_CONSOLE_IO)
 924       fprintf(stderr,"white_x=%f, white_y=%f\n",white_x, white_y);
 925 #endif
 926       return;
 927    }
 928    itemp = (png_uint_32)(white_x * 100000.0 + 0.5);
 929    png_save_uint_32(buf, itemp);
 930    itemp = (png_uint_32)(white_y * 100000.0 + 0.5);
 931    png_save_uint_32(buf + 4, itemp);
 932
 933    if (red_x < 0 || red_x > 0.8 || red_y < 0 || red_y > 0.8 ||
 934        red_x + red_y > 1.0)
 935    {
 936       png_warning(png_ptr, "Invalid cHRM red point specified");
 937       return;
 938    }
 939    itemp = (png_uint_32)(red_x * 100000.0 + 0.5);
 940    png_save_uint_32(buf + 8, itemp);
 941    itemp = (png_uint_32)(red_y * 100000.0 + 0.5);
 942    png_save_uint_32(buf + 12, itemp);
 943
 944    if (green_x < 0 || green_x > 0.8 || green_y < 0 || green_y > 0.8 ||
 945        green_x + green_y > 1.0)
 946    {
 947       png_warning(png_ptr, "Invalid cHRM green point specified");
 948       return;
 949    }
 950    itemp = (png_uint_32)(green_x * 100000.0 + 0.5);
 951    png_save_uint_32(buf + 16, itemp);
 952    itemp = (png_uint_32)(green_y * 100000.0 + 0.5);
 953    png_save_uint_32(buf + 20, itemp);
 954
 955    if (blue_x < 0 || blue_x > 0.8 || blue_y < 0 || blue_y > 0.8 ||
 956        blue_x + blue_y > 1.0)
 957    {
 958       png_warning(png_ptr, "Invalid cHRM blue point specified");
 959       return;
 960    }
 961    itemp = (png_uint_32)(blue_x * 100000.0 + 0.5);
 962    png_save_uint_32(buf + 24, itemp);
 963    itemp = (png_uint_32)(blue_y * 100000.0 + 0.5);
 964    png_save_uint_32(buf + 28, itemp);
 965
 966    png_write_chunk(png_ptr, (png_bytep)png_cHRM, buf, (png_size_t)32);
 967 }
 968 #endif
 969 #ifdef PNG_FIXED_POINT_SUPPORTED
 970 void /* PRIVATE */
 971 png_write_cHRM_fixed(png_structp png_ptr, png_fixed_point white_x,
 972    png_fixed_point white_y, png_fixed_point red_x, png_fixed_point red_y,
 973    png_fixed_point green_x, png_fixed_point green_y, png_fixed_point blue_x,
 974    png_fixed_point blue_y)
 975 {
 976 #ifdef PNG_USE_LOCAL_ARRAYS
 977    PNG_cHRM;
 978 #endif
 979    png_byte buf[32];
 980
 981    png_debug(1, "in png_write_cHRM\n");
 982    /* each value is saved in 1/100,000ths */
 983    if (white_x > 80000L || white_y > 80000L || white_x + white_y > 100000L)
 984    {
 985       png_warning(png_ptr, "Invalid fixed cHRM white point specified");
 986 #if !defined(PNG_NO_CONSOLE_IO)
 987       fprintf(stderr,"white_x=%ld, white_y=%ld\n",white_x, white_y);
 988 #endif
 989       return;
 990    }
 991    png_save_uint_32(buf, (png_uint_32)white_x);
 992    png_save_uint_32(buf + 4, (png_uint_32)white_y);
 993
 994    if (red_x > 80000L || red_y > 80000L || red_x + red_y > 100000L)
 995    {
 996       png_warning(png_ptr, "Invalid cHRM fixed red point specified");
 997       return;
 998    }
 999    png_save_uint_32(buf + 8, (png_uint_32)red_x);
1000    png_save_uint_32(buf + 12, (png_uint_32)red_y);
1001
1002    if (green_x > 80000L || green_y > 80000L || green_x + green_y > 100000L)
1003    {
1004       png_warning(png_ptr, "Invalid fixed cHRM green point specified");
1005       return;
1006    }
1007    png_save_uint_32(buf + 16, (png_uint_32)green_x);
1008    png_save_uint_32(buf + 20, (png_uint_32)green_y);
1009
1010    if (blue_x > 80000L || blue_y > 80000L || blue_x + blue_y > 100000L)
1011    {
1012       png_warning(png_ptr, "Invalid fixed cHRM blue point specified");
1013       return;
1014    }
1015    png_save_uint_32(buf + 24, (png_uint_32)blue_x);
1016    png_save_uint_32(buf + 28, (png_uint_32)blue_y);
1017
1018    png_write_chunk(png_ptr, (png_bytep)png_cHRM, buf, (png_size_t)32);
1019 }
1020 #endif
1021 #endif
1022
1023 #if defined(PNG_WRITE_tRNS_SUPPORTED)
1024 /* write the tRNS chunk */
1025 void /* PRIVATE */
1026 png_write_tRNS(png_structp png_ptr, png_bytep trans, png_color_16p tran,
1027    int num_trans, int color_type)
1028 {
1029 #ifdef PNG_USE_LOCAL_ARRAYS
1030    PNG_tRNS;
1031 #endif
1032    png_byte buf[6];
1033
1034    png_debug(1, "in png_write_tRNS\n");
1035    if (color_type == PNG_COLOR_TYPE_PALETTE)
1036    {
1037       if (num_trans <= 0 || num_trans > (int)png_ptr->num_palette)
1038       {
1039          png_warning(png_ptr,"Invalid number of transparent colors specified");
1040          return;
1041       }
1042       /* write the chunk out as it is */
1043       png_write_chunk(png_ptr, (png_bytep)png_tRNS, trans, (png_size_t)num_trans);
1044    }
1045    else if (color_type == PNG_COLOR_TYPE_GRAY)
1046    {
1047       /* one 16 bit value */
1048       if(tran->gray >= (1 << png_ptr->bit_depth))
1049       {
1050          png_warning(png_ptr,
1051            "Ignoring attempt to write tRNS chunk out-of-range for bit_depth");
1052          return;
1053       }
1054       png_save_uint_16(buf, tran->gray);
1055       png_write_chunk(png_ptr, (png_bytep)png_tRNS, buf, (png_size_t)2);
1056    }
1057    else if (color_type == PNG_COLOR_TYPE_RGB)
1058    {
1059       /* three 16 bit values */
1060       png_save_uint_16(buf, tran->red);
1061       png_save_uint_16(buf + 2, tran->green);
1062       png_save_uint_16(buf + 4, tran->blue);
1063       if(png_ptr->bit_depth == 8 && (buf[0] | buf[2] | buf[4]))
1064          {
1065             png_warning(png_ptr,
1066               "Ignoring attempt to write 16-bit tRNS chunk when bit_depth is 8");
1067             return;
1068          }
1069       png_write_chunk(png_ptr, (png_bytep)png_tRNS, buf, (png_size_t)6);
1070    }
1071    else
1072    {
1073       png_warning(png_ptr, "Can't write tRNS with an alpha channel");
1074    }
1075 }
1076 #endif
1077
1078 #if defined(PNG_WRITE_bKGD_SUPPORTED)
1079 /* write the background chunk */
1080 void /* PRIVATE */
1081 png_write_bKGD(png_structp png_ptr, png_color_16p back, int color_type)
1082 {
1083 #ifdef PNG_USE_LOCAL_ARRAYS
1084    PNG_bKGD;
1085 #endif
1086    png_byte buf[6];
1087
1088    png_debug(1, "in png_write_bKGD\n");
1089    if (color_type == PNG_COLOR_TYPE_PALETTE)
1090    {
1091       if (
1092 #if defined(PNG_MNG_FEATURES_SUPPORTED)
1093           (png_ptr->num_palette ||
1094           (!(png_ptr->mng_features_permitted & PNG_FLAG_MNG_EMPTY_PLTE))) &&
1095 #endif
1096          back->index > png_ptr->num_palette)
1097       {
1098          png_warning(png_ptr, "Invalid background palette index");
1099          return;
1100       }
1101       buf[0] = back->index;
1102       png_write_chunk(png_ptr, (png_bytep)png_bKGD, buf, (png_size_t)1);
1103    }
1104    else if (color_type & PNG_COLOR_MASK_COLOR)
1105    {
1106       png_save_uint_16(buf, back->red);
1107       png_save_uint_16(buf + 2, back->green);
1108       png_save_uint_16(buf + 4, back->blue);
1109       if(png_ptr->bit_depth == 8 && (buf[0] | buf[2] | buf[4]))
1110          {
1111             png_warning(png_ptr,
1112               "Ignoring attempt to write 16-bit bKGD chunk when bit_depth is 8");
1113             return;
1114          }
1115       png_write_chunk(png_ptr, (png_bytep)png_bKGD, buf, (png_size_t)6);
1116    }
1117    else
1118    {
1119       if(back->gray >= (1 << png_ptr->bit_depth))
1120       {
1121          png_warning(png_ptr,
1122            "Ignoring attempt to write bKGD chunk out-of-range for bit_depth");
1123          return;
1124       }
1125       png_save_uint_16(buf, back->gray);
1126       png_write_chunk(png_ptr, (png_bytep)png_bKGD, buf, (png_size_t)2);
1127    }
1128 }
1129 #endif
1130
1131 #if defined(PNG_WRITE_hIST_SUPPORTED)
1132 /* write the histogram */
1133 void /* PRIVATE */
1134 png_write_hIST(png_structp png_ptr, png_uint_16p hist, int num_hist)
1135 {
1136 #ifdef PNG_USE_LOCAL_ARRAYS
1137    PNG_hIST;
1138 #endif
1139    int i;
1140    png_byte buf[3];
1141
1142    png_debug(1, "in png_write_hIST\n");
1143    if (num_hist > (int)png_ptr->num_palette)
1144    {
1145       png_debug2(3, "num_hist = %d, num_palette = %d\n", num_hist,
1146          png_ptr->num_palette);
1147       png_warning(png_ptr, "Invalid number of histogram entries specified");
1148       return;
1149    }
1150
1151    png_write_chunk_start(png_ptr, (png_bytep)png_hIST, (png_uint_32)(num_hist * 2));
1152    for (i = 0; i < num_hist; i++)
1153    {
1154       png_save_uint_16(buf, hist[i]);
1155       png_write_chunk_data(png_ptr, buf, (png_size_t)2);
1156    }
1157    png_write_chunk_end(png_ptr);
1158 }
1159 #endif
1160
1161 #if defined(PNG_WRITE_TEXT_SUPPORTED) || defined(PNG_WRITE_pCAL_SUPPORTED) || \
1162     defined(PNG_WRITE_iCCP_SUPPORTED) || defined(PNG_WRITE_sPLT_SUPPORTED)
1163 /* Check that the tEXt or zTXt keyword is valid per PNG 1.0 specification,
1164  * and if invalid, correct the keyword rather than discarding the entire
1165  * chunk.  The PNG 1.0 specification requires keywords 1-79 characters in
1166  * length, forbids leading or trailing whitespace, multiple internal spaces,
1167  * and the non-break space (0x80) from ISO 8859-1.  Returns keyword length.
1168  *
1169  * The new_key is allocated to hold the corrected keyword and must be freed
1170  * by the calling routine.  This avoids problems with trying to write to
1171  * static keywords without having to have duplicate copies of the strings.
1172  */
1173 png_size_t /* PRIVATE */
1174 png_check_keyword(png_structp png_ptr, png_charp key, png_charpp new_key)
1175 {
1176    png_size_t key_len;
1177    png_charp kp, dp;
1178    int kflag;
1179    int kwarn=0;
1180
1181    png_debug(1, "in png_check_keyword\n");
1182    *new_key = NULL;
1183
1184    if (key == NULL || (key_len = png_strlen(key)) == 0)
1185    {
1186       png_warning(png_ptr, "zero length keyword");
1187       return ((png_size_t)0);
1188    }
1189
1190    png_debug1(2, "Keyword to be checked is '%s'\n", key);
1191
1192    *new_key = (png_charp)png_malloc_warn(png_ptr, (png_uint_32)(key_len + 2));
1193    if (*new_key == NULL)
1194    {
1195       png_warning(png_ptr, "Out of memory while procesing keyword");
1196       return ((png_size_t)0);
1197    }
1198
1199    /* Replace non-printing characters with a blank and print a warning */
1200    for (kp = key, dp = *new_key; *kp != '\0'; kp++, dp++)
1201    {
1202       if (*kp < 0x20 || (*kp > 0x7E && (png_byte)*kp < 0xA1))
1203       {
1204 #if !defined(PNG_NO_STDIO) && !defined(_WIN32_WCE)
1205          char msg[40];
1206
1207          sprintf(msg, "invalid keyword character 0x%02X", *kp);
1208          png_warning(png_ptr, msg);
1209 #else
1210          png_warning(png_ptr, "invalid character in keyword");
1211 #endif
1212          *dp = ' ';
1213       }
1214       else
1215       {
1216          *dp = *kp;
1217       }
1218    }
1219    *dp = '\0';
1220
1221    /* Remove any trailing white space. */
1222    kp = *new_key + key_len - 1;
1223    if (*kp == ' ')
1224    {
1225       png_warning(png_ptr, "trailing spaces removed from keyword");
1226
1227       while (*kp == ' ')
1228       {
1229         *(kp--) = '\0';
1230         key_len--;
1231       }
1232    }
1233
1234    /* Remove any leading white space. */
1235    kp = *new_key;
1236    if (*kp == ' ')
1237    {
1238       png_warning(png_ptr, "leading spaces removed from keyword");
1239
1240       while (*kp == ' ')
1241       {
1242         kp++;
1243         key_len--;
1244       }
1245    }
1246
1247    png_debug1(2, "Checking for multiple internal spaces in '%s'\n", kp);
1248
1249    /* Remove multiple internal spaces. */
1250    for (kflag = 0, dp = *new_key; *kp != '\0'; kp++)
1251    {
1252       if (*kp == ' ' && kflag == 0)
1253       {
1254          *(dp++) = *kp;
1255          kflag = 1;
1256       }
1257       else if (*kp == ' ')
1258       {
1259          key_len--;
1260          kwarn=1;
1261       }
1262       else
1263       {
1264          *(dp++) = *kp;
1265          kflag = 0;
1266       }
1267    }
1268    *dp = '\0';
1269    if(kwarn)
1270       png_warning(png_ptr, "extra interior spaces removed from keyword");
1271
1272    if (key_len == 0)
1273    {
1274       png_free(png_ptr, *new_key);
1275       *new_key=NULL;
1276       png_warning(png_ptr, "Zero length keyword");
1277    }
1278
1279    if (key_len > 79)
1280    {
1281       png_warning(png_ptr, "keyword length must be 1 - 79 characters");
1282       new_key[79] = '\0';
1283       key_len = 79;
1284    }
1285
1286    return (key_len);
1287 }
1288 #endif
1289
1290 #if defined(PNG_WRITE_tEXt_SUPPORTED)
1291 /* write a tEXt chunk */
1292 void /* PRIVATE */
1293 png_write_tEXt(png_structp png_ptr, png_charp key, png_charp text,
1294    png_size_t text_len)
1295 {
1296 #ifdef PNG_USE_LOCAL_ARRAYS
1297    PNG_tEXt;
1298 #endif
1299    png_size_t key_len;
1300    png_charp new_key;
1301
1302    png_debug(1, "in png_write_tEXt\n");
1303    if (key == NULL || (key_len = png_check_keyword(png_ptr, key, &new_key))==0)
1304    {
1305       png_warning(png_ptr, "Empty keyword in tEXt chunk");
1306       return;
1307    }
1308
1309    if (text == NULL || *text == '\0')
1310       text_len = 0;
1311    else
1312       text_len = png_strlen(text);
1313
1314    /* make sure we include the 0 after the key */
1315    png_write_chunk_start(png_ptr, (png_bytep)png_tEXt, (png_uint_32)key_len+text_len+1);
1316    /*
1317     * We leave it to the application to meet PNG-1.0 requirements on the
1318     * contents of the text.  PNG-1.0 through PNG-1.2 discourage the use of
1319     * any non-Latin-1 characters except for NEWLINE.  ISO PNG will forbid them.
1320     * The NUL character is forbidden by PNG-1.0 through PNG-1.2 and ISO PNG.
1321     */
1322    png_write_chunk_data(png_ptr, (png_bytep)new_key, key_len + 1);
1323    if (text_len)
1324       png_write_chunk_data(png_ptr, (png_bytep)text, text_len);
1325
1326    png_write_chunk_end(png_ptr);
1327    png_free(png_ptr, new_key);
1328 }
1329 #endif
1330
1331 #if defined(PNG_WRITE_zTXt_SUPPORTED)
1332 /* write a compressed text chunk */
1333 void /* PRIVATE */
1334 png_write_zTXt(png_structp png_ptr, png_charp key, png_charp text,
1335    png_size_t text_len, int compression)
1336 {
1337 #ifdef PNG_USE_LOCAL_ARRAYS
1338    PNG_zTXt;
1339 #endif
1340    png_size_t key_len;
1341    char buf[1];
1342    png_charp new_key;
1343    compression_state comp;
1344
1345    png_debug(1, "in png_write_zTXt\n");
1346
1347    if (key == NULL || (key_len = png_check_keyword(png_ptr, key, &new_key))==0)
1348    {
1349       png_warning(png_ptr, "Empty keyword in zTXt chunk");
1350       return;
1351    }
1352
1353    if (text == NULL || *text == '\0' || compression==PNG_TEXT_COMPRESSION_NONE)
1354    {
1355       png_write_tEXt(png_ptr, new_key, text, (png_size_t)0);
1356       png_free(png_ptr, new_key);
1357       return;
1358    }
1359
1360    text_len = png_strlen(text);
1361
1362    png_free(png_ptr, new_key);
1363
1364    /* compute the compressed data; do it now for the length */
1365    text_len = png_text_compress(png_ptr, text, text_len, compression,
1366        &comp);
1367
1368    /* write start of chunk */
1369    png_write_chunk_start(png_ptr, (png_bytep)png_zTXt, (png_uint_32)
1370       (key_len+text_len+2));
1371    /* write key */
1372    png_write_chunk_data(png_ptr, (png_bytep)key, key_len + 1);
1373    buf[0] = (png_byte)compression;
1374    /* write compression */
1375    png_write_chunk_data(png_ptr, (png_bytep)buf, (png_size_t)1);
1376    /* write the compressed data */
1377    png_write_compressed_data_out(png_ptr, &comp);
1378
1379    /* close the chunk */
1380    png_write_chunk_end(png_ptr);
1381 }
1382 #endif
1383
1384 #if defined(PNG_WRITE_iTXt_SUPPORTED)
1385 /* write an iTXt chunk */
1386 void /* PRIVATE */
1387 png_write_iTXt(png_structp png_ptr, int compression, png_charp key,
1388     png_charp lang, png_charp lang_key, png_charp text)
1389 {
1390 #ifdef PNG_USE_LOCAL_ARRAYS
1391    PNG_iTXt;
1392 #endif
1393    png_size_t lang_len, key_len, lang_key_len, text_len;
1394    png_charp new_lang, new_key;
1395    png_byte cbuf[2];
1396    compression_state comp;
1397
1398    png_debug(1, "in png_write_iTXt\n");
1399
1400    if (key == NULL || (key_len = png_check_keyword(png_ptr, key, &new_key))==0)
1401    {
1402       png_warning(png_ptr, "Empty keyword in iTXt chunk");
1403       return;
1404    }
1405    if (lang == NULL || (lang_len = png_check_keyword(png_ptr, lang, &new_lang))==0)
1406    {
1407       png_warning(png_ptr, "Empty language field in iTXt chunk");
1408       new_lang = NULL;
1409       lang_len = 0;
1410    }
1411
1412    if (lang_key == NULL)
1413      lang_key_len = 0;
1414    else
1415      lang_key_len = png_strlen(lang_key);
1416
1417    if (text == NULL)
1418       text_len = 0;
1419    else
1420      text_len = png_strlen(text);
1421
1422    /* compute the compressed data; do it now for the length */
1423    text_len = png_text_compress(png_ptr, text, text_len, compression-2,
1424       &comp);
1425
1426
1427    /* make sure we include the compression flag, the compression byte,
1428     * and the NULs after the key, lang, and lang_key parts */
1429
1430    png_write_chunk_start(png_ptr, (png_bytep)png_iTXt,
1431           (png_uint_32)(
1432         5 /* comp byte, comp flag, terminators for key, lang and lang_key */
1433         + key_len
1434         + lang_len
1435         + lang_key_len
1436         + text_len));
1437
1438    /*
1439     * We leave it to the application to meet PNG-1.0 requirements on the
1440     * contents of the text.  PNG-1.0 through PNG-1.2 discourage the use of
1441     * any non-Latin-1 characters except for NEWLINE.  ISO PNG will forbid them.
1442     * The NUL character is forbidden by PNG-1.0 through PNG-1.2 and ISO PNG.
1443     */
1444    png_write_chunk_data(png_ptr, (png_bytep)new_key, key_len + 1);
1445
1446    /* set the compression flag */
1447    if (compression == PNG_ITXT_COMPRESSION_NONE || \
1448        compression == PNG_TEXT_COMPRESSION_NONE)
1449        cbuf[0] = 0;
1450    else /* compression == PNG_ITXT_COMPRESSION_zTXt */
1451        cbuf[0] = 1;
1452    /* set the compression method */
1453    cbuf[1] = 0;
1454    png_write_chunk_data(png_ptr, cbuf, 2);
1455
1456    cbuf[0] = 0;
1457    png_write_chunk_data(png_ptr, (new_lang ? (png_bytep)new_lang : cbuf), lang_len + 1);
1458    png_write_chunk_data(png_ptr, (lang_key ? (png_bytep)lang_key : cbuf), lang_key_len + 1);
1459    png_write_compressed_data_out(png_ptr, &comp);
1460
1461    png_write_chunk_end(png_ptr);
1462    png_free(png_ptr, new_key);
1463    if (new_lang)
1464      png_free(png_ptr, new_lang);
1465 }
1466 #endif
1467
1468 #if defined(PNG_WRITE_oFFs_SUPPORTED)
1469 /* write the oFFs chunk */
1470 void /* PRIVATE */
1471 png_write_oFFs(png_structp png_ptr, png_int_32 x_offset, png_int_32 y_offset,
1472    int unit_type)
1473 {
1474 #ifdef PNG_USE_LOCAL_ARRAYS
1475    PNG_oFFs;
1476 #endif
1477    png_byte buf[9];
1478
1479    png_debug(1, "in png_write_oFFs\n");
1480    if (unit_type >= PNG_OFFSET_LAST)
1481       png_warning(png_ptr, "Unrecognized unit type for oFFs chunk");
1482
1483    png_save_int_32(buf, x_offset);
1484    png_save_int_32(buf + 4, y_offset);
1485    buf[8] = (png_byte)unit_type;
1486
1487    png_write_chunk(png_ptr, (png_bytep)png_oFFs, buf, (png_size_t)9);
1488 }
1489 #endif
1490
1491 #if defined(PNG_WRITE_pCAL_SUPPORTED)
1492 /* write the pCAL chunk (described in the PNG extensions document) */
1493 void /* PRIVATE */
1494 png_write_pCAL(png_structp png_ptr, png_charp purpose, png_int_32 X0,
1495    png_int_32 X1, int type, int nparams, png_charp units, png_charpp params)
1496 {
1497 #ifdef PNG_USE_LOCAL_ARRAYS
1498    PNG_pCAL;
1499 #endif
1500    png_size_t purpose_len, units_len, total_len;
1501    png_uint_32p params_len;
1502    png_byte buf[10];
1503    png_charp new_purpose;
1504    int i;
1505
1506    png_debug1(1, "in png_write_pCAL (%d parameters)\n", nparams);
1507    if (type >= PNG_EQUATION_LAST)
1508       png_warning(png_ptr, "Unrecognized equation type for pCAL chunk");
1509
1510    purpose_len = png_check_keyword(png_ptr, purpose, &new_purpose) + 1;
1511    png_debug1(3, "pCAL purpose length = %d\n", (int)purpose_len);
1512    units_len = png_strlen(units) + (nparams == 0 ? 0 : 1);
1513    png_debug1(3, "pCAL units length = %d\n", (int)units_len);
1514    total_len = purpose_len + units_len + 10;
1515
1516    params_len = (png_uint_32p)png_malloc(png_ptr, (png_uint_32)(nparams
1517       *png_sizeof(png_uint_32)));
1518
1519    /* Find the length of each parameter, making sure we don't count the
1520       null terminator for the last parameter. */
1521    for (i = 0; i < nparams; i++)
1522    {
1523       params_len[i] = png_strlen(params[i]) + (i == nparams - 1 ? 0 : 1);
1524       png_debug2(3, "pCAL parameter %d length = %lu\n", i, params_len[i]);
1525       total_len += (png_size_t)params_len[i];
1526    }
1527
1528    png_debug1(3, "pCAL total length = %d\n", (int)total_len);
1529    png_write_chunk_start(png_ptr, (png_bytep)png_pCAL, (png_uint_32)total_len);
1530    png_write_chunk_data(png_ptr, (png_bytep)new_purpose, purpose_len);
1531    png_save_int_32(buf, X0);
1532    png_save_int_32(buf + 4, X1);
1533    buf[8] = (png_byte)type;
1534    buf[9] = (png_byte)nparams;
1535    png_write_chunk_data(png_ptr, buf, (png_size_t)10);
1536    png_write_chunk_data(png_ptr, (png_bytep)units, (png_size_t)units_len);
1537
1538    png_free(png_ptr, new_purpose);
1539
1540    for (i = 0; i < nparams; i++)
1541    {
1542       png_write_chunk_data(png_ptr, (png_bytep)params[i],
1543          (png_size_t)params_len[i]);
1544    }
1545
1546    png_free(png_ptr, params_len);
1547    png_write_chunk_end(png_ptr);
1548 }
1549 #endif
1550
1551 #if defined(PNG_WRITE_sCAL_SUPPORTED)
1552 /* write the sCAL chunk */
1553 #if defined(PNG_FLOATING_POINT_SUPPORTED) && !defined(PNG_NO_STDIO)
1554 void /* PRIVATE */
1555 png_write_sCAL(png_structp png_ptr, int unit, double width,double height)
1556 {
1557 #ifdef PNG_USE_LOCAL_ARRAYS
1558    PNG_sCAL;
1559 #endif
1560    png_size_t total_len;
1561    char wbuf[32], hbuf[32];
1562    png_byte bunit = unit;
1563
1564    png_debug(1, "in png_write_sCAL\n");
1565
1566 #if defined(_WIN32_WCE)
1567 /* sprintf() function is not supported on WindowsCE */
1568    {
1569       wchar_t wc_buf[32];
1570       swprintf(wc_buf, TEXT("%12.12e"), width);
1571       WideCharToMultiByte(CP_ACP, 0, wc_buf, -1, wbuf, 32, NULL, NULL);
1572       swprintf(wc_buf, TEXT("%12.12e"), height);
1573       WideCharToMultiByte(CP_ACP, 0, wc_buf, -1, hbuf, 32, NULL, NULL);
1574    }
1575 #else
1576    sprintf(wbuf, "%12.12e", width);
1577    sprintf(hbuf, "%12.12e", height);
1578 #endif
1579    total_len = 1 + png_strlen(wbuf)+1 + png_strlen(hbuf);
1580
1581    png_debug1(3, "sCAL total length = %d\n", (int)total_len);
1582    png_write_chunk_start(png_ptr, (png_bytep)png_sCAL, (png_uint_32)total_len);
1583    png_write_chunk_data(png_ptr, (png_bytep)&bunit, 1);
1584    png_write_chunk_data(png_ptr, (png_bytep)wbuf, png_strlen(wbuf)+1);
1585    png_write_chunk_data(png_ptr, (png_bytep)hbuf, png_strlen(hbuf));
1586
1587    png_write_chunk_end(png_ptr);
1588 }
1589 #else
1590 #ifdef PNG_FIXED_POINT_SUPPORTED
1591 void /* PRIVATE */
1592 png_write_sCAL_s(png_structp png_ptr, int unit, png_charp width,
1593    png_charp height)
1594 {
1595 #ifdef PNG_USE_LOCAL_ARRAYS
1596    PNG_sCAL;
1597 #endif
1598    png_size_t total_len;
1599    char wbuf[32], hbuf[32];
1600    png_byte bunit = unit;
1601
1602    png_debug(1, "in png_write_sCAL_s\n");
1603
1604    png_strcpy(wbuf,(const char *)width);
1605    png_strcpy(hbuf,(const char *)height);
1606    total_len = 1 + png_strlen(wbuf)+1 + png_strlen(hbuf);
1607
1608    png_debug1(3, "sCAL total length = %d\n", total_len);
1609    png_write_chunk_start(png_ptr, (png_bytep)png_sCAL, (png_uint_32)total_len);
1610    png_write_chunk_data(png_ptr, (png_bytep)&bunit, 1);
1611    png_write_chunk_data(png_ptr, (png_bytep)wbuf, png_strlen(wbuf)+1);
1612    png_write_chunk_data(png_ptr, (png_bytep)hbuf, png_strlen(hbuf));
1613
1614    png_write_chunk_end(png_ptr);
1615 }
1616 #endif
1617 #endif
1618 #endif
1619
1620 #if defined(PNG_WRITE_pHYs_SUPPORTED)
1621 /* write the pHYs chunk */
1622 void /* PRIVATE */
1623 png_write_pHYs(png_structp png_ptr, png_uint_32 x_pixels_per_unit,
1624    png_uint_32 y_pixels_per_unit,
1625    int unit_type)
1626 {
1627 #ifdef PNG_USE_LOCAL_ARRAYS
1628    PNG_pHYs;
1629 #endif
1630    png_byte buf[9];
1631
1632    png_debug(1, "in png_write_pHYs\n");
1633    if (unit_type >= PNG_RESOLUTION_LAST)
1634       png_warning(png_ptr, "Unrecognized unit type for pHYs chunk");
1635
1636    png_save_uint_32(buf, x_pixels_per_unit);
1637    png_save_uint_32(buf + 4, y_pixels_per_unit);
1638    buf[8] = (png_byte)unit_type;
1639
1640    png_write_chunk(png_ptr, (png_bytep)png_pHYs, buf, (png_size_t)9);
1641 }
1642 #endif
1643
1644 #if defined(PNG_WRITE_tIME_SUPPORTED)
1645 /* Write the tIME chunk.  Use either png_convert_from_struct_tm()
1646  * or png_convert_from_time_t(), or fill in the structure yourself.
1647  */
1648 void /* PRIVATE */
1649 png_write_tIME(png_structp png_ptr, png_timep mod_time)
1650 {
1651 #ifdef PNG_USE_LOCAL_ARRAYS
1652    PNG_tIME;
1653 #endif
1654    png_byte buf[7];
1655
1656    png_debug(1, "in png_write_tIME\n");
1657    if (mod_time->month  > 12 || mod_time->month  < 1 ||
1658        mod_time->day    > 31 || mod_time->day    < 1 ||
1659        mod_time->hour   > 23 || mod_time->second > 60)
1660    {
1661       png_warning(png_ptr, "Invalid time specified for tIME chunk");
1662       return;
1663    }
1664
1665    png_save_uint_16(buf, mod_time->year);
1666    buf[2] = mod_time->month;
1667    buf[3] = mod_time->day;
1668    buf[4] = mod_time->hour;
1669    buf[5] = mod_time->minute;
1670    buf[6] = mod_time->second;
1671
1672    png_write_chunk(png_ptr, (png_bytep)png_tIME, buf, (png_size_t)7);
1673 }
1674 #endif
1675
1676 /* initializes the row writing capability of libpng */
1677 void /* PRIVATE */
1678 png_write_start_row(png_structp png_ptr)
1679 {
1680 #ifdef PNG_USE_LOCAL_ARRAYS
1681    /* arrays to facilitate easy interlacing - use pass (0 - 6) as index */
1682
1683    /* start of interlace block */
1684    int png_pass_start[7] = {0, 4, 0, 2, 0, 1, 0};
1685
1686    /* offset to next interlace block */
1687    int png_pass_inc[7] = {8, 8, 4, 4, 2, 2, 1};
1688
1689    /* start of interlace block in the y direction */
1690    int png_pass_ystart[7] = {0, 0, 4, 0, 2, 0, 1};
1691
1692    /* offset to next interlace block in the y direction */
1693    int png_pass_yinc[7] = {8, 8, 8, 4, 4, 2, 2};
1694 #endif
1695
1696    png_size_t buf_size;
1697
1698    png_debug(1, "in png_write_start_row\n");
1699    buf_size = (png_size_t)(PNG_ROWBYTES(
1700       png_ptr->usr_channels*png_ptr->usr_bit_depth,png_ptr->width)+1);
1701
1702    /* set up row buffer */
1703    png_ptr->row_buf = (png_bytep)png_malloc(png_ptr, (png_uint_32)buf_size);
1704    png_ptr->row_buf[0] = PNG_FILTER_VALUE_NONE;
1705
1706    /* set up filtering buffer, if using this filter */
1707    if (png_ptr->do_filter & PNG_FILTER_SUB)
1708    {
1709       png_ptr->sub_row = (png_bytep)png_malloc(png_ptr,
1710          (png_ptr->rowbytes + 1));
1711       png_ptr->sub_row[0] = PNG_FILTER_VALUE_SUB;
1712    }
1713
1714    /* We only need to keep the previous row if we are using one of these. */
1715    if (png_ptr->do_filter & (PNG_FILTER_AVG | PNG_FILTER_UP | PNG_FILTER_PAETH))
1716    {
1717      /* set up previous row buffer */
1718       png_ptr->prev_row = (png_bytep)png_malloc(png_ptr, (png_uint_32)buf_size);
1719       png_memset(png_ptr->prev_row, 0, buf_size);
1720
1721       if (png_ptr->do_filter & PNG_FILTER_UP)
1722       {
1723          png_ptr->up_row = (png_bytep )png_malloc(png_ptr,
1724             (png_ptr->rowbytes + 1));
1725          png_ptr->up_row[0] = PNG_FILTER_VALUE_UP;
1726       }
1727
1728       if (png_ptr->do_filter & PNG_FILTER_AVG)
1729       {
1730          png_ptr->avg_row = (png_bytep)png_malloc(png_ptr,
1731             (png_ptr->rowbytes + 1));
1732          png_ptr->avg_row[0] = PNG_FILTER_VALUE_AVG;
1733       }
1734
1735       if (png_ptr->do_filter & PNG_FILTER_PAETH)
1736       {
1737          png_ptr->paeth_row = (png_bytep )png_malloc(png_ptr,
1738             (png_ptr->rowbytes + 1));
1739          png_ptr->paeth_row[0] = PNG_FILTER_VALUE_PAETH;
1740       }
1741    }
1742
1743 #ifdef PNG_WRITE_INTERLACING_SUPPORTED
1744    /* if interlaced, we need to set up width and height of pass */
1745    if (png_ptr->interlaced)
1746    {
1747       if (!(png_ptr->transformations & PNG_INTERLACE))
1748       {
1749          png_ptr->num_rows = (png_ptr->height + png_pass_yinc[0] - 1 -
1750             png_pass_ystart[0]) / png_pass_yinc[0];
1751          png_ptr->usr_width = (png_ptr->width + png_pass_inc[0] - 1 -
1752             png_pass_start[0]) / png_pass_inc[0];
1753       }
1754       else
1755       {
1756          png_ptr->num_rows = png_ptr->height;
1757          png_ptr->usr_width = png_ptr->width;
1758       }
1759    }
1760    else
1761 #endif
1762    {
1763       png_ptr->num_rows = png_ptr->height;
1764       png_ptr->usr_width = png_ptr->width;
1765    }
1766    png_ptr->zstream.avail_out = (uInt)png_ptr->zbuf_size;
1767    png_ptr->zstream.next_out = png_ptr->zbuf;
1768 }
1769
1770 /* Internal use only.  Called when finished processing a row of data. */
1771 void /* PRIVATE */
1772 png_write_finish_row(png_structp png_ptr)
1773 {
1774 #ifdef PNG_USE_LOCAL_ARRAYS
1775    /* arrays to facilitate easy interlacing - use pass (0 - 6) as index */
1776
1777    /* start of interlace block */
1778    int png_pass_start[7] = {0, 4, 0, 2, 0, 1, 0};
1779
1780    /* offset to next interlace block */
1781    int png_pass_inc[7] = {8, 8, 4, 4, 2, 2, 1};
1782
1783    /* start of interlace block in the y direction */
1784    int png_pass_ystart[7] = {0, 0, 4, 0, 2, 0, 1};
1785
1786    /* offset to next interlace block in the y direction */
1787    int png_pass_yinc[7] = {8, 8, 8, 4, 4, 2, 2};
1788 #endif
1789
1790    int ret;
1791
1792    png_debug(1, "in png_write_finish_row\n");
1793    /* next row */
1794    png_ptr->row_number++;
1795
1796    /* see if we are done */
1797    if (png_ptr->row_number < png_ptr->num_rows)
1798       return;
1799
1800 #ifdef PNG_WRITE_INTERLACING_SUPPORTED
1801    /* if interlaced, go to next pass */
1802    if (png_ptr->interlaced)
1803    {
1804       png_ptr->row_number = 0;
1805       if (png_ptr->transformations & PNG_INTERLACE)
1806       {
1807          png_ptr->pass++;
1808       }
1809       else
1810       {
1811          /* loop until we find a non-zero width or height pass */
1812          do
1813          {
1814             png_ptr->pass++;
1815             if (png_ptr->pass >= 7)
1816                break;
1817             png_ptr->usr_width = (png_ptr->width +
1818                png_pass_inc[png_ptr->pass] - 1 -
1819                png_pass_start[png_ptr->pass]) /
1820                png_pass_inc[png_ptr->pass];
1821             png_ptr->num_rows = (png_ptr->height +
1822                png_pass_yinc[png_ptr->pass] - 1 -
1823                png_pass_ystart[png_ptr->pass]) /
1824                png_pass_yinc[png_ptr->pass];
1825             if (png_ptr->transformations & PNG_INTERLACE)
1826                break;
1827          } while (png_ptr->usr_width == 0 || png_ptr->num_rows == 0);
1828
1829       }
1830
1831       /* reset the row above the image for the next pass */
1832       if (png_ptr->pass < 7)
1833       {
1834          if (png_ptr->prev_row != NULL)
1835             png_memset(png_ptr->prev_row, 0,
1836                (png_size_t)(PNG_ROWBYTES(png_ptr->usr_channels*
1837                png_ptr->usr_bit_depth,png_ptr->width))+1);
1838          return;
1839       }
1840    }
1841 #endif
1842
1843    /* if we get here, we've just written the last row, so we need
1844       to flush the compressor */
1845    do
1846    {
1847       /* tell the compressor we are done */
1848       ret = deflate(&png_ptr->zstream, Z_FINISH);
1849       /* check for an error */
1850       if (ret == Z_OK)
1851       {
1852          /* check to see if we need more room */
1853          if (!(png_ptr->zstream.avail_out))
1854          {
1855             png_write_IDAT(png_ptr, png_ptr->zbuf, png_ptr->zbuf_size);
1856             png_ptr->zstream.next_out = png_ptr->zbuf;
1857             png_ptr->zstream.avail_out = (uInt)png_ptr->zbuf_size;
1858          }
1859       }
1860       else if (ret != Z_STREAM_END)
1861       {
1862          if (png_ptr->zstream.msg != NULL)
1863             png_error(png_ptr, png_ptr->zstream.msg);
1864          else
1865             png_error(png_ptr, "zlib error");
1866       }
1867    } while (ret != Z_STREAM_END);
1868
1869    /* write any extra space */
1870    if (png_ptr->zstream.avail_out < png_ptr->zbuf_size)
1871    {
1872       png_write_IDAT(png_ptr, png_ptr->zbuf, png_ptr->zbuf_size -
1873          png_ptr->zstream.avail_out);
1874    }
1875
1876    deflateReset(&png_ptr->zstream);
1877    png_ptr->zstream.data_type = Z_BINARY;
1878 }
1879
1880 #if defined(PNG_WRITE_INTERLACING_SUPPORTED)
1881 /* Pick out the correct pixels for the interlace pass.
1882  * The basic idea here is to go through the row with a source
1883  * pointer and a destination pointer (sp and dp), and copy the
1884  * correct pixels for the pass.  As the row gets compacted,
1885  * sp will always be >= dp, so we should never overwrite anything.
1886  * See the default: case for the easiest code to understand.
1887  */
1888 void /* PRIVATE */
1889 png_do_write_interlace(png_row_infop row_info, png_bytep row, int pass)
1890 {
1891 #ifdef PNG_USE_LOCAL_ARRAYS
1892    /* arrays to facilitate easy interlacing - use pass (0 - 6) as index */
1893
1894    /* start of interlace block */
1895    int png_pass_start[7] = {0, 4, 0, 2, 0, 1, 0};
1896
1897    /* offset to next interlace block */
1898    int png_pass_inc[7] = {8, 8, 4, 4, 2, 2, 1};
1899 #endif
1900
1901    png_debug(1, "in png_do_write_interlace\n");
1902    /* we don't have to do anything on the last pass (6) */
1903 #if defined(PNG_USELESS_TESTS_SUPPORTED)
1904    if (row != NULL && row_info != NULL && pass < 6)
1905 #else
1906    if (pass < 6)
1907 #endif
1908    {
1909       /* each pixel depth is handled separately */
1910       switch (row_info->pixel_depth)
1911       {
1912          case 1:
1913          {
1914             png_bytep sp;
1915             png_bytep dp;
1916             int shift;
1917             int d;
1918             int value;
1919             png_uint_32 i;
1920             png_uint_32 row_width = row_info->width;
1921
1922             dp = row;
1923             d = 0;
1924             shift = 7;
1925             for (i = png_pass_start[pass]; i < row_width;
1926                i += png_pass_inc[pass])
1927             {
1928                sp = row + (png_size_t)(i >> 3);
1929                value = (int)(*sp >> (7 - (int)(i & 0x07))) & 0x01;
1930                d |= (value << shift);
1931
1932                if (shift == 0)
1933                {
1934                   shift = 7;
1935                   *dp++ = (png_byte)d;
1936                   d = 0;
1937                }
1938                else
1939                   shift--;
1940
1941             }
1942             if (shift != 7)
1943                *dp = (png_byte)d;
1944             break;
1945          }
1946          case 2:
1947          {
1948             png_bytep sp;
1949             png_bytep dp;
1950             int shift;
1951             int d;
1952             int value;
1953             png_uint_32 i;
1954             png_uint_32 row_width = row_info->width;
1955
1956             dp = row;
1957             shift = 6;
1958             d = 0;
1959             for (i = png_pass_start[pass]; i < row_width;
1960                i += png_pass_inc[pass])
1961             {
1962                sp = row + (png_size_t)(i >> 2);
1963                value = (*sp >> ((3 - (int)(i & 0x03)) << 1)) & 0x03;
1964                d |= (value << shift);
1965
1966                if (shift == 0)
1967                {
1968                   shift = 6;
1969                   *dp++ = (png_byte)d;
1970                   d = 0;
1971                }
1972                else
1973                   shift -= 2;
1974             }
1975             if (shift != 6)
1976                    *dp = (png_byte)d;
1977             break;
1978          }
1979          case 4:
1980          {
1981             png_bytep sp;
1982             png_bytep dp;
1983             int shift;
1984             int d;
1985             int value;
1986             png_uint_32 i;
1987             png_uint_32 row_width = row_info->width;
1988
1989             dp = row;
1990             shift = 4;
1991             d = 0;
1992             for (i = png_pass_start[pass]; i < row_width;
1993                i += png_pass_inc[pass])
1994             {
1995                sp = row + (png_size_t)(i >> 1);
1996                value = (*sp >> ((1 - (int)(i & 0x01)) << 2)) & 0x0f;
1997                d |= (value << shift);
1998
1999                if (shift == 0)
2000                {
2001                   shift = 4;
2002                   *dp++ = (png_byte)d;
2003                   d = 0;
2004                }
2005                else
2006                   shift -= 4;
2007             }
2008             if (shift != 4)
2009                *dp = (png_byte)d;
2010             break;
2011          }
2012          default:
2013          {
2014             png_bytep sp;
2015             png_bytep dp;
2016             png_uint_32 i;
2017             png_uint_32 row_width = row_info->width;
2018             png_size_t pixel_bytes;
2019
2020             /* start at the beginning */
2021             dp = row;
2022             /* find out how many bytes each pixel takes up */
2023             pixel_bytes = (row_info->pixel_depth >> 3);
2024             /* loop through the row, only looking at the pixels that
2025                matter */
2026             for (i = png_pass_start[pass]; i < row_width;
2027                i += png_pass_inc[pass])
2028             {
2029                /* find out where the original pixel is */
2030                sp = row + (png_size_t)i * pixel_bytes;
2031                /* move the pixel */
2032                if (dp != sp)
2033                   png_memcpy(dp, sp, pixel_bytes);
2034                /* next pixel */
2035                dp += pixel_bytes;
2036             }
2037             break;
2038          }
2039       }
2040       /* set new row width */
2041       row_info->width = (row_info->width +
2042          png_pass_inc[pass] - 1 -
2043          png_pass_start[pass]) /
2044          png_pass_inc[pass];
2045          row_info->rowbytes = PNG_ROWBYTES(row_info->pixel_depth,
2046             row_info->width);
2047    }
2048 }
2049 #endif
2050
2051 /* This filters the row, chooses which filter to use, if it has not already
2052  * been specified by the application, and then writes the row out with the
2053  * chosen filter.
2054  */
2055 #define PNG_MAXSUM (((png_uint_32)(-1)) >> 1)
2056 #define PNG_HISHIFT 10
2057 #define PNG_LOMASK ((png_uint_32)0xffffL)
2058 #define PNG_HIMASK ((png_uint_32)(~PNG_LOMASK >> PNG_HISHIFT))
2059 void /* PRIVATE */
2060 png_write_find_filter(png_structp png_ptr, png_row_infop row_info)
2061 {
2062    png_bytep prev_row, best_row, row_buf;
2063    png_uint_32 mins, bpp;
2064    png_byte filter_to_do = png_ptr->do_filter;
2065    png_uint_32 row_bytes = row_info->rowbytes;
2066 #if defined(PNG_WRITE_WEIGHTED_FILTER_SUPPORTED)
2067    int num_p_filters = (int)png_ptr->num_prev_filters;
2068 #endif
2069
2070    png_debug(1, "in png_write_find_filter\n");
2071    /* find out how many bytes offset each pixel is */
2072    bpp = (row_info->pixel_depth + 7) >> 3;
2073
2074    prev_row = png_ptr->prev_row;
2075    best_row = row_buf = png_ptr->row_buf;
2076    mins = PNG_MAXSUM;
2077
2078    /* The prediction method we use is to find which method provides the
2079     * smallest value when summing the absolute values of the distances
2080     * from zero, using anything >= 128 as negative numbers.  This is known
2081     * as the "minimum sum of absolute differences" heuristic.  Other
2082     * heuristics are the "weighted minimum sum of absolute differences"
2083     * (experimental and can in theory improve compression), and the "zlib
2084     * predictive" method (not implemented yet), which does test compressions
2085     * of lines using different filter methods, and then chooses the
2086     * (series of) filter(s) that give minimum compressed data size (VERY
2087     * computationally expensive).
2088     *
2089     * GRR 980525:  consider also
2090     *   (1) minimum sum of absolute differences from running average (i.e.,
2091     *       keep running sum of non-absolute differences & count of bytes)
2092     *       [track dispersion, too?  restart average if dispersion too large?]
2093     *  (1b) minimum sum of absolute differences from sliding average, probably
2094     *       with window size <= deflate window (usually 32K)
2095     *   (2) minimum sum of squared differences from zero or running average
2096     *       (i.e., ~ root-mean-square approach)
2097     */
2098
2099
2100    /* We don't need to test the 'no filter' case if this is the only filter
2101     * that has been chosen, as it doesn't actually do anything to the data.
2102     */
2103    if ((filter_to_do & PNG_FILTER_NONE) &&
2104        filter_to_do != PNG_FILTER_NONE)
2105    {
2106       png_bytep rp;
2107       png_uint_32 sum = 0;
2108       png_uint_32 i;
2109       int v;
2110
2111       for (i = 0, rp = row_buf + 1; i < row_bytes; i++, rp++)
2112       {
2113          v = *rp;
2114          sum += (v < 128) ? v : 256 - v;
2115       }
2116
2117 #if defined(PNG_WRITE_WEIGHTED_FILTER_SUPPORTED)
2118       if (png_ptr->heuristic_method == PNG_FILTER_HEURISTIC_WEIGHTED)
2119       {
2120          png_uint_32 sumhi, sumlo;
2121          int j;
2122          sumlo = sum & PNG_LOMASK;
2123          sumhi = (sum >> PNG_HISHIFT) & PNG_HIMASK; /* Gives us some footroom */
2124
2125          /* Reduce the sum if we match any of the previous rows */
2126          for (j = 0; j < num_p_filters; j++)
2127          {
2128             if (png_ptr->prev_filters[j] == PNG_FILTER_VALUE_NONE)
2129             {
2130                sumlo = (sumlo * png_ptr->filter_weights[j]) >>
2131                   PNG_WEIGHT_SHIFT;
2132                sumhi = (sumhi * png_ptr->filter_weights[j]) >>
2133                   PNG_WEIGHT_SHIFT;
2134             }
2135          }
2136
2137          /* Factor in the cost of this filter (this is here for completeness,
2138           * but it makes no sense to have a "cost" for the NONE filter, as
2139           * it has the minimum possible computational cost - none).
2140           */
2141          sumlo = (sumlo * png_ptr->filter_costs[PNG_FILTER_VALUE_NONE]) >>
2142             PNG_COST_SHIFT;
2143          sumhi = (sumhi * png_ptr->filter_costs[PNG_FILTER_VALUE_NONE]) >>
2144             PNG_COST_SHIFT;
2145
2146          if (sumhi > PNG_HIMASK)
2147             sum = PNG_MAXSUM;
2148          else
2149             sum = (sumhi << PNG_HISHIFT) + sumlo;
2150       }
2151 #endif
2152       mins = sum;
2153    }
2154
2155    /* sub filter */
2156    if (filter_to_do == PNG_FILTER_SUB)
2157    /* it's the only filter so no testing is needed */
2158    {
2159       png_bytep rp, lp, dp;
2160       png_uint_32 i;
2161       for (i = 0, rp = row_buf + 1, dp = png_ptr->sub_row + 1; i < bpp;
2162            i++, rp++, dp++)
2163       {
2164          *dp = *rp;
2165       }
2166       for (lp = row_buf + 1; i < row_bytes;
2167          i++, rp++, lp++, dp++)
2168       {
2169          *dp = (png_byte)(((int)*rp - (int)*lp) & 0xff);
2170       }
2171       best_row = png_ptr->sub_row;
2172    }
2173
2174    else if (filter_to_do & PNG_FILTER_SUB)
2175    {
2176       png_bytep rp, dp, lp;
2177       png_uint_32 sum = 0, lmins = mins;
2178       png_uint_32 i;
2179       int v;
2180
2181 #if defined(PNG_WRITE_WEIGHTED_FILTER_SUPPORTED)
2182       /* We temporarily increase the "minimum sum" by the factor we
2183        * would reduce the sum of this filter, so that we can do the
2184        * early exit comparison without scaling the sum each time.
2185        */
2186       if (png_ptr->heuristic_method == PNG_FILTER_HEURISTIC_WEIGHTED)
2187       {
2188          int j;
2189          png_uint_32 lmhi, lmlo;
2190          lmlo = lmins & PNG_LOMASK;
2191          lmhi = (lmins >> PNG_HISHIFT) & PNG_HIMASK;
2192
2193          for (j = 0; j < num_p_filters; j++)
2194          {
2195             if (png_ptr->prev_filters[j] == PNG_FILTER_VALUE_SUB)
2196             {
2197                lmlo = (lmlo * png_ptr->inv_filter_weights[j]) >>
2198                   PNG_WEIGHT_SHIFT;
2199                lmhi = (lmhi * png_ptr->inv_filter_weights[j]) >>
2200                   PNG_WEIGHT_SHIFT;
2201             }
2202          }
2203
2204          lmlo = (lmlo * png_ptr->inv_filter_costs[PNG_FILTER_VALUE_SUB]) >>
2205             PNG_COST_SHIFT;
2206          lmhi = (lmhi * png_ptr->inv_filter_costs[PNG_FILTER_VALUE_SUB]) >>
2207             PNG_COST_SHIFT;
2208
2209          if (lmhi > PNG_HIMASK)
2210             lmins = PNG_MAXSUM;
2211          else
2212             lmins = (lmhi << PNG_HISHIFT) + lmlo;
2213       }
2214 #endif
2215
2216       for (i = 0, rp = row_buf + 1, dp = png_ptr->sub_row + 1; i < bpp;
2217            i++, rp++, dp++)
2218       {
2219          v = *dp = *rp;
2220
2221          sum += (v < 128) ? v : 256 - v;
2222       }
2223       for (lp = row_buf + 1; i < row_bytes;
2224          i++, rp++, lp++, dp++)
2225       {
2226          v = *dp = (png_byte)(((int)*rp - (int)*lp) & 0xff);
2227
2228          sum += (v < 128) ? v : 256 - v;
2229
2230          if (sum > lmins)  /* We are already worse, don't continue. */
2231             break;
2232       }
2233
2234 #if defined(PNG_WRITE_WEIGHTED_FILTER_SUPPORTED)
2235       if (png_ptr->heuristic_method == PNG_FILTER_HEURISTIC_WEIGHTED)
2236       {
2237          int j;
2238          png_uint_32 sumhi, sumlo;
2239          sumlo = sum & PNG_LOMASK;
2240          sumhi = (sum >> PNG_HISHIFT) & PNG_HIMASK;
2241
2242          for (j = 0; j < num_p_filters; j++)
2243          {
2244             if (png_ptr->prev_filters[j] == PNG_FILTER_VALUE_SUB)
2245             {
2246                sumlo = (sumlo * png_ptr->inv_filter_weights[j]) >>
2247                   PNG_WEIGHT_SHIFT;
2248                sumhi = (sumhi * png_ptr->inv_filter_weights[j]) >>
2249                   PNG_WEIGHT_SHIFT;
2250             }
2251          }
2252
2253          sumlo = (sumlo * png_ptr->inv_filter_costs[PNG_FILTER_VALUE_SUB]) >>
2254             PNG_COST_SHIFT;
2255          sumhi = (sumhi * png_ptr->inv_filter_costs[PNG_FILTER_VALUE_SUB]) >>
2256             PNG_COST_SHIFT;
2257
2258          if (sumhi > PNG_HIMASK)
2259             sum = PNG_MAXSUM;
2260          else
2261             sum = (sumhi << PNG_HISHIFT) + sumlo;
2262       }
2263 #endif
2264
2265       if (sum < mins)
2266       {
2267          mins = sum;
2268          best_row = png_ptr->sub_row;
2269       }
2270    }
2271
2272    /* up filter */
2273    if (filter_to_do == PNG_FILTER_UP)
2274    {
2275       png_bytep rp, dp, pp;
2276       png_uint_32 i;
2277
2278       for (i = 0, rp = row_buf + 1, dp = png_ptr->up_row + 1,
2279            pp = prev_row + 1; i < row_bytes;
2280            i++, rp++, pp++, dp++)
2281       {
2282          *dp = (png_byte)(((int)*rp - (int)*pp) & 0xff);
2283       }
2284       best_row = png_ptr->up_row;
2285    }
2286
2287    else if (filter_to_do & PNG_FILTER_UP)
2288    {
2289       png_bytep rp, dp, pp;
2290       png_uint_32 sum = 0, lmins = mins;
2291       png_uint_32 i;
2292       int v;
2293
2294
2295 #if defined(PNG_WRITE_WEIGHTED_FILTER_SUPPORTED)
2296       if (png_ptr->heuristic_method == PNG_FILTER_HEURISTIC_WEIGHTED)
2297       {
2298          int j;
2299          png_uint_32 lmhi, lmlo;
2300          lmlo = lmins & PNG_LOMASK;
2301          lmhi = (lmins >> PNG_HISHIFT) & PNG_HIMASK;
2302
2303          for (j = 0; j < num_p_filters; j++)
2304          {
2305             if (png_ptr->prev_filters[j] == PNG_FILTER_VALUE_UP)
2306             {
2307                lmlo = (lmlo * png_ptr->inv_filter_weights[j]) >>
2308                   PNG_WEIGHT_SHIFT;
2309                lmhi = (lmhi * png_ptr->inv_filter_weights[j]) >>
2310                   PNG_WEIGHT_SHIFT;
2311             }
2312          }
2313
2314          lmlo = (lmlo * png_ptr->inv_filter_costs[PNG_FILTER_VALUE_UP]) >>
2315             PNG_COST_SHIFT;
2316          lmhi = (lmhi * png_ptr->inv_filter_costs[PNG_FILTER_VALUE_UP]) >>
2317             PNG_COST_SHIFT;
2318
2319          if (lmhi > PNG_HIMASK)
2320             lmins = PNG_MAXSUM;
2321          else
2322             lmins = (lmhi << PNG_HISHIFT) + lmlo;
2323       }
2324 #endif
2325
2326       for (i = 0, rp = row_buf + 1, dp = png_ptr->up_row + 1,
2327            pp = prev_row + 1; i < row_bytes; i++)
2328       {
2329          v = *dp++ = (png_byte)(((int)*rp++ - (int)*pp++) & 0xff);
2330
2331          sum += (v < 128) ? v : 256 - v;
2332
2333          if (sum > lmins)  /* We are already worse, don't continue. */
2334             break;
2335       }
2336
2337 #if defined(PNG_WRITE_WEIGHTED_FILTER_SUPPORTED)
2338       if (png_ptr->heuristic_method == PNG_FILTER_HEURISTIC_WEIGHTED)
2339       {
2340          int j;
2341          png_uint_32 sumhi, sumlo;
2342          sumlo = sum & PNG_LOMASK;
2343          sumhi = (sum >> PNG_HISHIFT) & PNG_HIMASK;
2344
2345          for (j = 0; j < num_p_filters; j++)
2346          {
2347             if (png_ptr->prev_filters[j] == PNG_FILTER_VALUE_UP)
2348             {
2349                sumlo = (sumlo * png_ptr->filter_weights[j]) >>
2350                   PNG_WEIGHT_SHIFT;
2351                sumhi = (sumhi * png_ptr->filter_weights[j]) >>
2352                   PNG_WEIGHT_SHIFT;
2353             }
2354          }
2355
2356          sumlo = (sumlo * png_ptr->filter_costs[PNG_FILTER_VALUE_UP]) >>
2357             PNG_COST_SHIFT;
2358          sumhi = (sumhi * png_ptr->filter_costs[PNG_FILTER_VALUE_UP]) >>
2359             PNG_COST_SHIFT;
2360
2361          if (sumhi > PNG_HIMASK)
2362             sum = PNG_MAXSUM;
2363          else
2364             sum = (sumhi << PNG_HISHIFT) + sumlo;
2365       }
2366 #endif
2367
2368       if (sum < mins)
2369       {
2370          mins = sum;
2371          best_row = png_ptr->up_row;
2372       }
2373    }
2374
2375    /* avg filter */
2376    if (filter_to_do == PNG_FILTER_AVG)
2377    {
2378       png_bytep rp, dp, pp, lp;
2379       png_uint_32 i;
2380       for (i = 0, rp = row_buf + 1, dp = png_ptr->avg_row + 1,
2381            pp = prev_row + 1; i < bpp; i++)
2382       {
2383          *dp++ = (png_byte)(((int)*rp++ - ((int)*pp++ / 2)) & 0xff);
2384       }
2385       for (lp = row_buf + 1; i < row_bytes; i++)
2386       {
2387          *dp++ = (png_byte)(((int)*rp++ - (((int)*pp++ + (int)*lp++) / 2))
2388                  & 0xff);
2389       }
2390       best_row = png_ptr->avg_row;
2391    }
2392
2393    else if (filter_to_do & PNG_FILTER_AVG)
2394    {
2395       png_bytep rp, dp, pp, lp;
2396       png_uint_32 sum = 0, lmins = mins;
2397       png_uint_32 i;
2398       int v;
2399
2400 #if defined(PNG_WRITE_WEIGHTED_FILTER_SUPPORTED)
2401       if (png_ptr->heuristic_method == PNG_FILTER_HEURISTIC_WEIGHTED)
2402       {
2403          int j;
2404          png_uint_32 lmhi, lmlo;
2405          lmlo = lmins & PNG_LOMASK;
2406          lmhi = (lmins >> PNG_HISHIFT) & PNG_HIMASK;
2407
2408          for (j = 0; j < num_p_filters; j++)
2409          {
2410             if (png_ptr->prev_filters[j] == PNG_FILTER_VALUE_AVG)
2411             {
2412                lmlo = (lmlo * png_ptr->inv_filter_weights[j]) >>
2413                   PNG_WEIGHT_SHIFT;
2414                lmhi = (lmhi * png_ptr->inv_filter_weights[j]) >>
2415                   PNG_WEIGHT_SHIFT;
2416             }
2417          }
2418
2419          lmlo = (lmlo * png_ptr->inv_filter_costs[PNG_FILTER_VALUE_AVG]) >>
2420             PNG_COST_SHIFT;
2421          lmhi = (lmhi * png_ptr->inv_filter_costs[PNG_FILTER_VALUE_AVG]) >>
2422             PNG_COST_SHIFT;
2423
2424          if (lmhi > PNG_HIMASK)
2425             lmins = PNG_MAXSUM;
2426          else
2427             lmins = (lmhi << PNG_HISHIFT) + lmlo;
2428       }
2429 #endif
2430
2431       for (i = 0, rp = row_buf + 1, dp = png_ptr->avg_row + 1,
2432            pp = prev_row + 1; i < bpp; i++)
2433       {
2434          v = *dp++ = (png_byte)(((int)*rp++ - ((int)*pp++ / 2)) & 0xff);
2435
2436          sum += (v < 128) ? v : 256 - v;
2437       }
2438       for (lp = row_buf + 1; i < row_bytes; i++)
2439       {
2440          v = *dp++ =
2441           (png_byte)(((int)*rp++ - (((int)*pp++ + (int)*lp++) / 2)) & 0xff);
2442
2443          sum += (v < 128) ? v : 256 - v;
2444
2445          if (sum > lmins)  /* We are already worse, don't continue. */
2446             break;
2447       }
2448
2449 #if defined(PNG_WRITE_WEIGHTED_FILTER_SUPPORTED)
2450       if (png_ptr->heuristic_method == PNG_FILTER_HEURISTIC_WEIGHTED)
2451       {
2452          int j;
2453          png_uint_32 sumhi, sumlo;
2454          sumlo = sum & PNG_LOMASK;
2455          sumhi = (sum >> PNG_HISHIFT) & PNG_HIMASK;
2456
2457          for (j = 0; j < num_p_filters; j++)
2458          {
2459             if (png_ptr->prev_filters[j] == PNG_FILTER_VALUE_NONE)
2460             {
2461                sumlo = (sumlo * png_ptr->filter_weights[j]) >>
2462                   PNG_WEIGHT_SHIFT;
2463                sumhi = (sumhi * png_ptr->filter_weights[j]) >>
2464                   PNG_WEIGHT_SHIFT;
2465             }
2466          }
2467
2468          sumlo = (sumlo * png_ptr->filter_costs[PNG_FILTER_VALUE_AVG]) >>
2469             PNG_COST_SHIFT;
2470          sumhi = (sumhi * png_ptr->filter_costs[PNG_FILTER_VALUE_AVG]) >>
2471             PNG_COST_SHIFT;
2472
2473          if (sumhi > PNG_HIMASK)
2474             sum = PNG_MAXSUM;
2475          else
2476             sum = (sumhi << PNG_HISHIFT) + sumlo;
2477       }
2478 #endif
2479
2480       if (sum < mins)
2481       {
2482          mins = sum;
2483          best_row = png_ptr->avg_row;
2484       }
2485    }
2486
2487    /* Paeth filter */
2488    if (filter_to_do == PNG_FILTER_PAETH)
2489    {
2490       png_bytep rp, dp, pp, cp, lp;
2491       png_uint_32 i;
2492       for (i = 0, rp = row_buf + 1, dp = png_ptr->paeth_row + 1,
2493            pp = prev_row + 1; i < bpp; i++)
2494       {
2495          *dp++ = (png_byte)(((int)*rp++ - (int)*pp++) & 0xff);
2496       }
2497
2498       for (lp = row_buf + 1, cp = prev_row + 1; i < row_bytes; i++)
2499       {
2500          int a, b, c, pa, pb, pc, p;
2501
2502          b = *pp++;
2503          c = *cp++;
2504          a = *lp++;
2505
2506          p = b - c;
2507          pc = a - c;
2508
2509 #ifdef PNG_USE_ABS
2510          pa = abs(p);
2511          pb = abs(pc);
2512          pc = abs(p + pc);
2513 #else
2514          pa = p < 0 ? -p : p;
2515          pb = pc < 0 ? -pc : pc;
2516          pc = (p + pc) < 0 ? -(p + pc) : p + pc;
2517 #endif
2518
2519          p = (pa <= pb && pa <=pc) ? a : (pb <= pc) ? b : c;
2520
2521          *dp++ = (png_byte)(((int)*rp++ - p) & 0xff);
2522       }
2523       best_row = png_ptr->paeth_row;
2524    }
2525
2526    else if (filter_to_do & PNG_FILTER_PAETH)
2527    {
2528       png_bytep rp, dp, pp, cp, lp;
2529       png_uint_32 sum = 0, lmins = mins;
2530       png_uint_32 i;
2531       int v;
2532
2533 #if defined(PNG_WRITE_WEIGHTED_FILTER_SUPPORTED)
2534       if (png_ptr->heuristic_method == PNG_FILTER_HEURISTIC_WEIGHTED)
2535       {
2536          int j;
2537          png_uint_32 lmhi, lmlo;
2538          lmlo = lmins & PNG_LOMASK;
2539          lmhi = (lmins >> PNG_HISHIFT) & PNG_HIMASK;
2540
2541          for (j = 0; j < num_p_filters; j++)
2542          {
2543             if (png_ptr->prev_filters[j] == PNG_FILTER_VALUE_PAETH)
2544             {
2545                lmlo = (lmlo * png_ptr->inv_filter_weights[j]) >>
2546                   PNG_WEIGHT_SHIFT;
2547                lmhi = (lmhi * png_ptr->inv_filter_weights[j]) >>
2548                   PNG_WEIGHT_SHIFT;
2549             }
2550          }
2551
2552          lmlo = (lmlo * png_ptr->inv_filter_costs[PNG_FILTER_VALUE_PAETH]) >>
2553             PNG_COST_SHIFT;
2554          lmhi = (lmhi * png_ptr->inv_filter_costs[PNG_FILTER_VALUE_PAETH]) >>
2555             PNG_COST_SHIFT;
2556
2557          if (lmhi > PNG_HIMASK)
2558             lmins = PNG_MAXSUM;
2559          else
2560             lmins = (lmhi << PNG_HISHIFT) + lmlo;
2561       }
2562 #endif
2563
2564       for (i = 0, rp = row_buf + 1, dp = png_ptr->paeth_row + 1,
2565            pp = prev_row + 1; i < bpp; i++)
2566       {
2567          v = *dp++ = (png_byte)(((int)*rp++ - (int)*pp++) & 0xff);
2568
2569          sum += (v < 128) ? v : 256 - v;
2570       }
2571
2572       for (lp = row_buf + 1, cp = prev_row + 1; i < row_bytes; i++)
2573       {
2574          int a, b, c, pa, pb, pc, p;
2575
2576          b = *pp++;
2577          c = *cp++;
2578          a = *lp++;
2579
2580 #ifndef PNG_SLOW_PAETH
2581          p = b - c;
2582          pc = a - c;
2583 #ifdef PNG_USE_ABS
2584          pa = abs(p);
2585          pb = abs(pc);
2586          pc = abs(p + pc);
2587 #else
2588          pa = p < 0 ? -p : p;
2589          pb = pc < 0 ? -pc : pc;
2590          pc = (p + pc) < 0 ? -(p + pc) : p + pc;
2591 #endif
2592          p = (pa <= pb && pa <=pc) ? a : (pb <= pc) ? b : c;
2593 #else /* PNG_SLOW_PAETH */
2594          p = a + b - c;
2595          pa = abs(p - a);
2596          pb = abs(p - b);
2597          pc = abs(p - c);
2598          if (pa <= pb && pa <= pc)
2599             p = a;
2600          else if (pb <= pc)
2601             p = b;
2602          else
2603             p = c;
2604 #endif /* PNG_SLOW_PAETH */
2605
2606          v = *dp++ = (png_byte)(((int)*rp++ - p) & 0xff);
2607
2608          sum += (v < 128) ? v : 256 - v;
2609
2610          if (sum > lmins)  /* We are already worse, don't continue. */
2611             break;
2612       }
2613
2614 #if defined(PNG_WRITE_WEIGHTED_FILTER_SUPPORTED)
2615       if (png_ptr->heuristic_method == PNG_FILTER_HEURISTIC_WEIGHTED)
2616       {
2617          int j;
2618          png_uint_32 sumhi, sumlo;
2619          sumlo = sum & PNG_LOMASK;
2620          sumhi = (sum >> PNG_HISHIFT) & PNG_HIMASK;
2621
2622          for (j = 0; j < num_p_filters; j++)
2623          {
2624             if (png_ptr->prev_filters[j] == PNG_FILTER_VALUE_PAETH)
2625             {
2626                sumlo = (sumlo * png_ptr->filter_weights[j]) >>
2627                   PNG_WEIGHT_SHIFT;
2628                sumhi = (sumhi * png_ptr->filter_weights[j]) >>
2629                   PNG_WEIGHT_SHIFT;
2630             }
2631          }
2632
2633          sumlo = (sumlo * png_ptr->filter_costs[PNG_FILTER_VALUE_PAETH]) >>
2634             PNG_COST_SHIFT;
2635          sumhi = (sumhi * png_ptr->filter_costs[PNG_FILTER_VALUE_PAETH]) >>
2636             PNG_COST_SHIFT;
2637
2638          if (sumhi > PNG_HIMASK)
2639             sum = PNG_MAXSUM;
2640          else
2641             sum = (sumhi << PNG_HISHIFT) + sumlo;
2642       }
2643 #endif
2644
2645       if (sum < mins)
2646       {
2647          best_row = png_ptr->paeth_row;
2648       }
2649    }
2650
2651    /* Do the actual writing of the filtered row data from the chosen filter. */
2652
2653    png_write_filtered_row(png_ptr, best_row);
2654
2655 #if defined(PNG_WRITE_WEIGHTED_FILTER_SUPPORTED)
2656    /* Save the type of filter we picked this time for future calculations */
2657    if (png_ptr->num_prev_filters > 0)
2658    {
2659       int j;
2660       for (j = 1; j < num_p_filters; j++)
2661       {
2662          png_ptr->prev_filters[j] = png_ptr->prev_filters[j - 1];
2663       }
2664       png_ptr->prev_filters[j] = best_row[0];
2665    }
2666 #endif
2667 }
2668
2669
2670 /* Do the actual writing of a previously filtered row. */
2671 void /* PRIVATE */
2672 png_write_filtered_row(png_structp png_ptr, png_bytep filtered_row)
2673 {
2674    png_debug(1, "in png_write_filtered_row\n");
2675    png_debug1(2, "filter = %d\n", filtered_row[0]);
2676    /* set up the zlib input buffer */
2677
2678    png_ptr->zstream.next_in = filtered_row;
2679    png_ptr->zstream.avail_in = (uInt)png_ptr->row_info.rowbytes + 1;
2680    /* repeat until we have compressed all the data */
2681    do
2682    {
2683       int ret; /* return of zlib */
2684
2685       /* compress the data */
2686       ret = deflate(&png_ptr->zstream, Z_NO_FLUSH);
2687       /* check for compression errors */
2688       if (ret != Z_OK)
2689       {
2690          if (png_ptr->zstream.msg != NULL)
2691             png_error(png_ptr, png_ptr->zstream.msg);
2692          else
2693             png_error(png_ptr, "zlib error");
2694       }
2695
2696       /* see if it is time to write another IDAT */
2697       if (!(png_ptr->zstream.avail_out))
2698       {
2699          /* write the IDAT and reset the zlib output buffer */
2700          png_write_IDAT(png_ptr, png_ptr->zbuf, png_ptr->zbuf_size);
2701          png_ptr->zstream.next_out = png_ptr->zbuf;
2702          png_ptr->zstream.avail_out = (uInt)png_ptr->zbuf_size;
2703       }
2704    /* repeat until all data has been compressed */
2705    } while (png_ptr->zstream.avail_in);
2706
2707    /* swap the current and previous rows */
2708    if (png_ptr->prev_row != NULL)
2709    {
2710       png_bytep tptr;
2711
2712       tptr = png_ptr->prev_row;
2713       png_ptr->prev_row = png_ptr->row_buf;
2714       png_ptr->row_buf = tptr;
2715    }
2716
2717    /* finish row - updates counters and flushes zlib if last row */
2718    png_write_finish_row(png_ptr);
2719
2720 #if defined(PNG_WRITE_FLUSH_SUPPORTED)
2721    png_ptr->flush_rows++;
2722
2723    if (png_ptr->flush_dist > 0 &&
2724        png_ptr->flush_rows >= png_ptr->flush_dist)
2725    {
2726       png_write_flush(png_ptr);
2727    }
2728 #endif
2729 }
2730 #endif /* PNG_WRITE_SUPPORTED */