visibility determination
[vrlugburz] / libs / imago / libpng / pngwutil.c
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 */