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