fix a calcuation bug in FLAC__lpc_compute_best_order()
[flac.git] / src / libFLAC / stream_encoder.c
1 /* libFLAC - Free Lossless Audio Codec library
2  * Copyright (C) 2000,2001,2002,2003,2004,2005,2006  Josh Coalson
3  *
4  * Redistribution and use in source and binary forms, with or without
5  * modification, are permitted provided that the following conditions
6  * are met:
7  *
8  * - Redistributions of source code must retain the above copyright
9  * notice, this list of conditions and the following disclaimer.
10  *
11  * - Redistributions in binary form must reproduce the above copyright
12  * notice, this list of conditions and the following disclaimer in the
13  * documentation and/or other materials provided with the distribution.
14  *
15  * - Neither the name of the Xiph.org Foundation nor the names of its
16  * contributors may be used to endorse or promote products derived from
17  * this software without specific prior written permission.
18  *
19  * THIS SOFTWARE IS PROVIDED BY THE COPYRIGHT HOLDERS AND CONTRIBUTORS
20  * ``AS IS'' AND ANY EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT
21  * LIMITED TO, THE IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS FOR
22  * A PARTICULAR PURPOSE ARE DISCLAIMED.  IN NO EVENT SHALL THE FOUNDATION OR
23  * CONTRIBUTORS BE LIABLE FOR ANY DIRECT, INDIRECT, INCIDENTAL, SPECIAL,
24  * EXEMPLARY, OR CONSEQUENTIAL DAMAGES (INCLUDING, BUT NOT LIMITED TO,
25  * PROCUREMENT OF SUBSTITUTE GOODS OR SERVICES; LOSS OF USE, DATA, OR
26  * PROFITS; OR BUSINESS INTERRUPTION) HOWEVER CAUSED AND ON ANY THEORY OF
27  * LIABILITY, WHETHER IN CONTRACT, STRICT LIABILITY, OR TORT (INCLUDING
28  * NEGLIGENCE OR OTHERWISE) ARISING IN ANY WAY OUT OF THE USE OF THIS
29  * SOFTWARE, EVEN IF ADVISED OF THE POSSIBILITY OF SUCH DAMAGE.
30  */
31
32 /*@@@@@@*/
33 #define WINDOW_DEBUG_OUTPUT
34
35 #include <limits.h>
36 #include <stdio.h>
37 #include <stdlib.h> /* for malloc() */
38 #include <string.h> /* for memcpy() */
39 #include "FLAC/assert.h"
40 #include "FLAC/stream_decoder.h"
41 #include "protected/stream_encoder.h"
42 #include "private/bitbuffer.h"
43 #include "private/bitmath.h"
44 #include "private/crc.h"
45 #include "private/cpu.h"
46 #include "private/fixed.h"
47 #include "private/format.h"
48 #include "private/lpc.h"
49 #include "private/md5.h"
50 #include "private/memory.h"
51 #include "private/stream_encoder_framing.h"
52 #include "private/window.h"
53
54 #ifdef HAVE_CONFIG_H
55 #include <config.h>
56 #endif
57
58 #ifdef min
59 #undef min
60 #endif
61 #define min(x,y) ((x)<(y)?(x):(y))
62
63 #ifdef max
64 #undef max
65 #endif
66 #define max(x,y) ((x)>(y)?(x):(y))
67
68 typedef struct {
69         FLAC__int32 *data[FLAC__MAX_CHANNELS];
70         unsigned size; /* of each data[] in samples */
71         unsigned tail;
72 } verify_input_fifo;
73
74 typedef struct {
75         const FLAC__byte *data;
76         unsigned capacity;
77         unsigned bytes;
78 } verify_output;
79
80 typedef enum {
81         ENCODER_IN_MAGIC = 0,
82         ENCODER_IN_METADATA = 1,
83         ENCODER_IN_AUDIO = 2
84 } EncoderStateHint;
85
86 /***********************************************************************
87  *
88  * Private class method prototypes
89  *
90  ***********************************************************************/
91
92 static void set_defaults_(FLAC__StreamEncoder *encoder);
93 static void free_(FLAC__StreamEncoder *encoder);
94 static FLAC__bool resize_buffers_(FLAC__StreamEncoder *encoder, unsigned new_size);
95 static FLAC__bool write_bitbuffer_(FLAC__StreamEncoder *encoder, unsigned samples);
96 static FLAC__bool process_frame_(FLAC__StreamEncoder *encoder, FLAC__bool is_last_frame);
97 static FLAC__bool process_subframes_(FLAC__StreamEncoder *encoder, FLAC__bool is_last_frame);
98
99 static FLAC__bool process_subframe_(
100         FLAC__StreamEncoder *encoder,
101         unsigned min_partition_order,
102         unsigned max_partition_order,
103         FLAC__bool precompute_partition_sums,
104         const FLAC__FrameHeader *frame_header,
105         unsigned subframe_bps,
106         const FLAC__int32 integer_signal[],
107 #ifndef FLAC__INTEGER_ONLY_LIBRARY
108         const FLAC__real real_signal[],
109 #endif
110         FLAC__Subframe *subframe[2],
111         FLAC__EntropyCodingMethod_PartitionedRiceContents *partitioned_rice_contents[2],
112         FLAC__int32 *residual[2],
113         unsigned *best_subframe,
114         unsigned *best_bits
115 #ifdef WINDOW_DEBUG_OUTPUT
116         ,unsigned subframe_number
117 #endif
118 );
119
120 static FLAC__bool add_subframe_(
121         FLAC__StreamEncoder *encoder,
122         const FLAC__FrameHeader *frame_header,
123         unsigned subframe_bps,
124         const FLAC__Subframe *subframe,
125         FLAC__BitBuffer *frame
126 #ifdef WINDOW_DEBUG_OUTPUT
127 ,unsigned subframe_bits
128 #endif
129 );
130
131 static unsigned evaluate_constant_subframe_(
132         const FLAC__int32 signal,
133         unsigned subframe_bps,
134         FLAC__Subframe *subframe
135 );
136
137 static unsigned evaluate_fixed_subframe_(
138         FLAC__StreamEncoder *encoder,
139         const FLAC__int32 signal[],
140         FLAC__int32 residual[],
141         FLAC__uint32 abs_residual[],
142         FLAC__uint64 abs_residual_partition_sums[],
143         unsigned raw_bits_per_partition[],
144         unsigned blocksize,
145         unsigned subframe_bps,
146         unsigned order,
147         unsigned rice_parameter,
148         unsigned min_partition_order,
149         unsigned max_partition_order,
150         FLAC__bool precompute_partition_sums,
151         FLAC__bool do_escape_coding,
152         unsigned rice_parameter_search_dist,
153         FLAC__Subframe *subframe,
154         FLAC__EntropyCodingMethod_PartitionedRiceContents *partitioned_rice_contents
155 );
156
157 #ifndef FLAC__INTEGER_ONLY_LIBRARY
158 static unsigned evaluate_lpc_subframe_(
159         FLAC__StreamEncoder *encoder,
160         const FLAC__int32 signal[],
161         FLAC__int32 residual[],
162         FLAC__uint32 abs_residual[],
163         FLAC__uint64 abs_residual_partition_sums[],
164         unsigned raw_bits_per_partition[],
165         const FLAC__real lp_coeff[],
166         unsigned blocksize,
167         unsigned subframe_bps,
168         unsigned order,
169         unsigned qlp_coeff_precision,
170         unsigned rice_parameter,
171         unsigned min_partition_order,
172         unsigned max_partition_order,
173         FLAC__bool precompute_partition_sums,
174         FLAC__bool do_escape_coding,
175         unsigned rice_parameter_search_dist,
176         FLAC__Subframe *subframe,
177         FLAC__EntropyCodingMethod_PartitionedRiceContents *partitioned_rice_contents
178 #ifdef WINDOW_DEBUG_OUTPUT
179         ,unsigned frame_number
180         ,unsigned subframe_number
181         ,FLAC__ApodizationSpecification aspec
182 #endif
183 );
184 #endif
185
186 static unsigned evaluate_verbatim_subframe_(
187         const FLAC__int32 signal[],
188         unsigned blocksize,
189         unsigned subframe_bps,
190         FLAC__Subframe *subframe
191 );
192
193 static unsigned find_best_partition_order_(
194         struct FLAC__StreamEncoderPrivate *private_,
195         const FLAC__int32 residual[],
196         FLAC__uint32 abs_residual[],
197         FLAC__uint64 abs_residual_partition_sums[],
198         unsigned raw_bits_per_partition[],
199         unsigned residual_samples,
200         unsigned predictor_order,
201         unsigned rice_parameter,
202         unsigned min_partition_order,
203         unsigned max_partition_order,
204         FLAC__bool precompute_partition_sums,
205         FLAC__bool do_escape_coding,
206         unsigned rice_parameter_search_dist,
207         FLAC__EntropyCodingMethod_PartitionedRice *best_partitioned_rice
208 );
209
210 static void precompute_partition_info_sums_(
211         const FLAC__uint32 abs_residual[],
212         FLAC__uint64 abs_residual_partition_sums[],
213         unsigned residual_samples,
214         unsigned predictor_order,
215         unsigned min_partition_order,
216         unsigned max_partition_order
217 );
218
219 static void precompute_partition_info_escapes_(
220         const FLAC__int32 residual[],
221         unsigned raw_bits_per_partition[],
222         unsigned residual_samples,
223         unsigned predictor_order,
224         unsigned min_partition_order,
225         unsigned max_partition_order
226 );
227
228 #ifdef DONT_ESTIMATE_RICE_BITS
229 static FLAC__bool set_partitioned_rice_(
230         const FLAC__uint32 abs_residual[],
231         const FLAC__int32 residual[],
232         const unsigned residual_samples,
233         const unsigned predictor_order,
234         const unsigned suggested_rice_parameter,
235         const unsigned rice_parameter_search_dist,
236         const unsigned partition_order,
237         FLAC__EntropyCodingMethod_PartitionedRiceContents *partitioned_rice_contents,
238         unsigned *bits
239 );
240
241 static FLAC__bool set_partitioned_rice_with_precompute_(
242         const FLAC__int32 residual[],
243         const FLAC__uint64 abs_residual_partition_sums[],
244         const unsigned raw_bits_per_partition[],
245         const unsigned residual_samples,
246         const unsigned predictor_order,
247         const unsigned suggested_rice_parameter,
248         const unsigned rice_parameter_search_dist,
249         const unsigned partition_order,
250         const FLAC__bool search_for_escapes,
251         FLAC__EntropyCodingMethod_PartitionedRiceContents *partitioned_rice_contents,
252         unsigned *bits
253 );
254 #else
255 static FLAC__bool set_partitioned_rice_(
256         const FLAC__uint32 abs_residual[],
257         const unsigned residual_samples,
258         const unsigned predictor_order,
259         const unsigned suggested_rice_parameter,
260         const unsigned rice_parameter_search_dist,
261         const unsigned partition_order,
262         FLAC__EntropyCodingMethod_PartitionedRiceContents *partitioned_rice_contents,
263         unsigned *bits
264 );
265
266 static FLAC__bool set_partitioned_rice_with_precompute_(
267         const FLAC__uint32 abs_residual[],
268         const FLAC__uint64 abs_residual_partition_sums[],
269         const unsigned raw_bits_per_partition[],
270         const unsigned residual_samples,
271         const unsigned predictor_order,
272         const unsigned suggested_rice_parameter,
273         const unsigned rice_parameter_search_dist,
274         const unsigned partition_order,
275         const FLAC__bool search_for_escapes,
276         FLAC__EntropyCodingMethod_PartitionedRiceContents *partitioned_rice_contents,
277         unsigned *bits
278 );
279 #endif
280
281 static unsigned get_wasted_bits_(FLAC__int32 signal[], unsigned samples);
282
283 /* verify-related routines: */
284 static void append_to_verify_fifo_(
285         verify_input_fifo *fifo,
286         const FLAC__int32 * const input[],
287         unsigned input_offset,
288         unsigned channels,
289         unsigned wide_samples
290 );
291
292 static void append_to_verify_fifo_interleaved_(
293         verify_input_fifo *fifo,
294         const FLAC__int32 input[],
295         unsigned input_offset,
296         unsigned channels,
297         unsigned wide_samples
298 );
299
300 static FLAC__StreamDecoderReadStatus verify_read_callback_(
301         const FLAC__StreamDecoder *decoder,
302         FLAC__byte buffer[],
303         unsigned *bytes,
304         void *client_data
305 );
306
307 static FLAC__StreamDecoderWriteStatus verify_write_callback_(
308         const FLAC__StreamDecoder *decoder,
309         const FLAC__Frame *frame,
310         const FLAC__int32 * const buffer[],
311         void *client_data
312 );
313
314 static void verify_metadata_callback_(
315         const FLAC__StreamDecoder *decoder,
316         const FLAC__StreamMetadata *metadata,
317         void *client_data
318 );
319
320 static void verify_error_callback_(
321         const FLAC__StreamDecoder *decoder,
322         FLAC__StreamDecoderErrorStatus status,
323         void *client_data
324 );
325
326 #ifdef WINDOW_DEBUG_OUTPUT
327 static const char * const winstr[] = {
328         "bartlett",
329         "bartlett_hann",
330         "blackman",
331         "blackman_harris_4term_92db_sidelobe",
332         "connes",
333         "flattop",
334         "gauss",
335         "hamming",
336         "hann",
337         "kaiser_bessel",
338         "nuttall",
339         "rectangular",
340         "triangle",
341         "tukey",
342         "welch"
343 };
344 #endif
345
346 /***********************************************************************
347  *
348  * Private class data
349  *
350  ***********************************************************************/
351
352 typedef struct FLAC__StreamEncoderPrivate {
353         unsigned input_capacity;                          /* current size (in samples) of the signal and residual buffers */
354         FLAC__int32 *integer_signal[FLAC__MAX_CHANNELS];  /* the integer version of the input signal */
355         FLAC__int32 *integer_signal_mid_side[2];          /* the integer version of the mid-side input signal (stereo only) */
356 #ifndef FLAC__INTEGER_ONLY_LIBRARY
357         FLAC__real *real_signal[FLAC__MAX_CHANNELS];      /* the floating-point version of the input signal */
358         FLAC__real *real_signal_mid_side[2];              /* the floating-point version of the mid-side input signal (stereo only) */
359         FLAC__real *window[FLAC__MAX_APODIZATION_FUNCTIONS]; /* the pre-computed floating-point window for each apodization function */
360         FLAC__real *windowed_signal;                      /* the real_signal[] * current window[] */
361 #endif
362         unsigned subframe_bps[FLAC__MAX_CHANNELS];        /* the effective bits per sample of the input signal (stream bps - wasted bits) */
363         unsigned subframe_bps_mid_side[2];                /* the effective bits per sample of the mid-side input signal (stream bps - wasted bits + 0/1) */
364         FLAC__int32 *residual_workspace[FLAC__MAX_CHANNELS][2]; /* each channel has a candidate and best workspace where the subframe residual signals will be stored */
365         FLAC__int32 *residual_workspace_mid_side[2][2];
366         FLAC__Subframe subframe_workspace[FLAC__MAX_CHANNELS][2];
367         FLAC__Subframe subframe_workspace_mid_side[2][2];
368         FLAC__Subframe *subframe_workspace_ptr[FLAC__MAX_CHANNELS][2];
369         FLAC__Subframe *subframe_workspace_ptr_mid_side[2][2];
370         FLAC__EntropyCodingMethod_PartitionedRiceContents partitioned_rice_contents_workspace[FLAC__MAX_CHANNELS][2];
371         FLAC__EntropyCodingMethod_PartitionedRiceContents partitioned_rice_contents_workspace_mid_side[FLAC__MAX_CHANNELS][2];
372         FLAC__EntropyCodingMethod_PartitionedRiceContents *partitioned_rice_contents_workspace_ptr[FLAC__MAX_CHANNELS][2];
373         FLAC__EntropyCodingMethod_PartitionedRiceContents *partitioned_rice_contents_workspace_ptr_mid_side[FLAC__MAX_CHANNELS][2];
374         unsigned best_subframe[FLAC__MAX_CHANNELS];       /* index into the above workspaces */
375         unsigned best_subframe_mid_side[2];
376         unsigned best_subframe_bits[FLAC__MAX_CHANNELS];  /* size in bits of the best subframe for each channel */
377         unsigned best_subframe_bits_mid_side[2];
378         FLAC__uint32 *abs_residual;                       /* workspace where abs(candidate residual) is stored */
379         FLAC__uint64 *abs_residual_partition_sums;        /* workspace where the sum of abs(candidate residual) for each partition is stored */
380         unsigned *raw_bits_per_partition;                 /* workspace where the sum of silog2(candidate residual) for each partition is stored */
381         FLAC__BitBuffer *frame;                           /* the current frame being worked on */
382         unsigned loose_mid_side_stereo_frames;            /* rounded number of frames the encoder will use before trying both independent and mid/side frames again */
383         unsigned loose_mid_side_stereo_frame_count;       /* number of frames using the current channel assignment */
384         FLAC__ChannelAssignment last_channel_assignment;
385         FLAC__StreamMetadata metadata;
386         unsigned current_sample_number;
387         unsigned current_frame_number;
388         struct FLAC__MD5Context md5context;
389         FLAC__CPUInfo cpuinfo;
390 #ifndef FLAC__INTEGER_ONLY_LIBRARY
391         unsigned (*local_fixed_compute_best_predictor)(const FLAC__int32 data[], unsigned data_len, FLAC__float residual_bits_per_sample[FLAC__MAX_FIXED_ORDER+1]);
392 #else
393         unsigned (*local_fixed_compute_best_predictor)(const FLAC__int32 data[], unsigned data_len, FLAC__fixedpoint residual_bits_per_sample[FLAC__MAX_FIXED_ORDER+1]);
394 #endif
395 #ifndef FLAC__INTEGER_ONLY_LIBRARY
396         void (*local_lpc_compute_autocorrelation)(const FLAC__real data[], unsigned data_len, unsigned lag, FLAC__real autoc[]);
397         void (*local_lpc_compute_residual_from_qlp_coefficients)(const FLAC__int32 *data, unsigned data_len, const FLAC__int32 qlp_coeff[], unsigned order, int lp_quantization, FLAC__int32 residual[]);
398         void (*local_lpc_compute_residual_from_qlp_coefficients_64bit)(const FLAC__int32 *data, unsigned data_len, const FLAC__int32 qlp_coeff[], unsigned order, int lp_quantization, FLAC__int32 residual[]);
399         void (*local_lpc_compute_residual_from_qlp_coefficients_16bit)(const FLAC__int32 *data, unsigned data_len, const FLAC__int32 qlp_coeff[], unsigned order, int lp_quantization, FLAC__int32 residual[]);
400 #endif
401         FLAC__bool use_wide_by_block;          /* use slow 64-bit versions of some functions because of the block size */
402         FLAC__bool use_wide_by_partition;      /* use slow 64-bit versions of some functions because of the min partition order and blocksize */
403         FLAC__bool use_wide_by_order;          /* use slow 64-bit versions of some functions because of the lpc order */
404         FLAC__bool precompute_partition_sums;  /* our initial guess as to whether precomputing the partitions sums will be a speed improvement */
405         FLAC__bool disable_constant_subframes;
406         FLAC__bool disable_fixed_subframes;
407         FLAC__bool disable_verbatim_subframes;
408         FLAC__StreamEncoderWriteCallback write_callback;
409         FLAC__StreamEncoderMetadataCallback metadata_callback;
410         void *client_data;
411         /* unaligned (original) pointers to allocated data */
412         FLAC__int32 *integer_signal_unaligned[FLAC__MAX_CHANNELS];
413         FLAC__int32 *integer_signal_mid_side_unaligned[2];
414 #ifndef FLAC__INTEGER_ONLY_LIBRARY
415         FLAC__real *real_signal_unaligned[FLAC__MAX_CHANNELS];
416         FLAC__real *real_signal_mid_side_unaligned[2];
417         FLAC__real *window_unaligned[FLAC__MAX_APODIZATION_FUNCTIONS];
418         FLAC__real *windowed_signal_unaligned;
419 #endif
420         FLAC__int32 *residual_workspace_unaligned[FLAC__MAX_CHANNELS][2];
421         FLAC__int32 *residual_workspace_mid_side_unaligned[2][2];
422         FLAC__uint32 *abs_residual_unaligned;
423         FLAC__uint64 *abs_residual_partition_sums_unaligned;
424         unsigned *raw_bits_per_partition_unaligned;
425         /*
426          * These fields have been moved here from private function local
427          * declarations merely to save stack space during encoding.
428          */
429 #ifndef FLAC__INTEGER_ONLY_LIBRARY
430         FLAC__real lp_coeff[FLAC__MAX_LPC_ORDER][FLAC__MAX_LPC_ORDER]; /* from process_subframe_() */
431 #endif
432         FLAC__EntropyCodingMethod_PartitionedRiceContents partitioned_rice_contents_extra[2]; /* from find_best_partition_order_() */
433         /*
434          * The data for the verify section
435          */
436         struct {
437                 FLAC__StreamDecoder *decoder;
438                 EncoderStateHint state_hint;
439                 FLAC__bool needs_magic_hack;
440                 verify_input_fifo input_fifo;
441                 verify_output output;
442                 struct {
443                         FLAC__uint64 absolute_sample;
444                         unsigned frame_number;
445                         unsigned channel;
446                         unsigned sample;
447                         FLAC__int32 expected;
448                         FLAC__int32 got;
449                 } error_stats;
450         } verify;
451         FLAC__bool is_being_deleted; /* if true, call to ..._finish() from ..._delete() will not call the callbacks */
452 } FLAC__StreamEncoderPrivate;
453
454 /***********************************************************************
455  *
456  * Public static class data
457  *
458  ***********************************************************************/
459
460 FLAC_API const char * const FLAC__StreamEncoderStateString[] = {
461         "FLAC__STREAM_ENCODER_OK",
462         "FLAC__STREAM_ENCODER_VERIFY_DECODER_ERROR",
463         "FLAC__STREAM_ENCODER_VERIFY_MISMATCH_IN_AUDIO_DATA",
464         "FLAC__STREAM_ENCODER_INVALID_CALLBACK",
465         "FLAC__STREAM_ENCODER_INVALID_NUMBER_OF_CHANNELS",
466         "FLAC__STREAM_ENCODER_INVALID_BITS_PER_SAMPLE",
467         "FLAC__STREAM_ENCODER_INVALID_SAMPLE_RATE",
468         "FLAC__STREAM_ENCODER_INVALID_BLOCK_SIZE",
469         "FLAC__STREAM_ENCODER_INVALID_MAX_LPC_ORDER",
470         "FLAC__STREAM_ENCODER_INVALID_QLP_COEFF_PRECISION",
471         "FLAC__STREAM_ENCODER_MID_SIDE_CHANNELS_MISMATCH",
472         "FLAC__STREAM_ENCODER_MID_SIDE_SAMPLE_SIZE_MISMATCH",
473         "FLAC__STREAM_ENCODER_ILLEGAL_MID_SIDE_FORCE",
474         "FLAC__STREAM_ENCODER_BLOCK_SIZE_TOO_SMALL_FOR_LPC_ORDER",
475         "FLAC__STREAM_ENCODER_NOT_STREAMABLE",
476         "FLAC__STREAM_ENCODER_FRAMING_ERROR",
477         "FLAC__STREAM_ENCODER_INVALID_METADATA",
478         "FLAC__STREAM_ENCODER_FATAL_ERROR_WHILE_ENCODING",
479         "FLAC__STREAM_ENCODER_FATAL_ERROR_WHILE_WRITING",
480         "FLAC__STREAM_ENCODER_MEMORY_ALLOCATION_ERROR",
481         "FLAC__STREAM_ENCODER_ALREADY_INITIALIZED",
482         "FLAC__STREAM_ENCODER_UNINITIALIZED"
483 };
484
485 FLAC_API const char * const FLAC__StreamEncoderWriteStatusString[] = {
486         "FLAC__STREAM_ENCODER_WRITE_STATUS_OK",
487         "FLAC__STREAM_ENCODER_WRITE_STATUS_FATAL_ERROR"
488 };
489
490 /***********************************************************************
491  *
492  * Class constructor/destructor
493  *
494  */
495 FLAC_API FLAC__StreamEncoder *FLAC__stream_encoder_new()
496 {
497         FLAC__StreamEncoder *encoder;
498         unsigned i;
499
500         FLAC__ASSERT(sizeof(int) >= 4); /* we want to die right away if this is not true */
501
502         encoder = (FLAC__StreamEncoder*)calloc(1, sizeof(FLAC__StreamEncoder));
503         if(encoder == 0) {
504                 return 0;
505         }
506
507         encoder->protected_ = (FLAC__StreamEncoderProtected*)calloc(1, sizeof(FLAC__StreamEncoderProtected));
508         if(encoder->protected_ == 0) {
509                 free(encoder);
510                 return 0;
511         }
512
513         encoder->private_ = (FLAC__StreamEncoderPrivate*)calloc(1, sizeof(FLAC__StreamEncoderPrivate));
514         if(encoder->private_ == 0) {
515                 free(encoder->protected_);
516                 free(encoder);
517                 return 0;
518         }
519
520         encoder->private_->frame = FLAC__bitbuffer_new();
521         if(encoder->private_->frame == 0) {
522                 free(encoder->private_);
523                 free(encoder->protected_);
524                 free(encoder);
525                 return 0;
526         }
527
528         set_defaults_(encoder);
529
530         encoder->private_->is_being_deleted = false;
531
532         for(i = 0; i < FLAC__MAX_CHANNELS; i++) {
533                 encoder->private_->subframe_workspace_ptr[i][0] = &encoder->private_->subframe_workspace[i][0];
534                 encoder->private_->subframe_workspace_ptr[i][1] = &encoder->private_->subframe_workspace[i][1];
535         }
536         for(i = 0; i < 2; i++) {
537                 encoder->private_->subframe_workspace_ptr_mid_side[i][0] = &encoder->private_->subframe_workspace_mid_side[i][0];
538                 encoder->private_->subframe_workspace_ptr_mid_side[i][1] = &encoder->private_->subframe_workspace_mid_side[i][1];
539         }
540         for(i = 0; i < FLAC__MAX_CHANNELS; i++) {
541                 encoder->private_->partitioned_rice_contents_workspace_ptr[i][0] = &encoder->private_->partitioned_rice_contents_workspace[i][0];
542                 encoder->private_->partitioned_rice_contents_workspace_ptr[i][1] = &encoder->private_->partitioned_rice_contents_workspace[i][1];
543         }
544         for(i = 0; i < 2; i++) {
545                 encoder->private_->partitioned_rice_contents_workspace_ptr_mid_side[i][0] = &encoder->private_->partitioned_rice_contents_workspace_mid_side[i][0];
546                 encoder->private_->partitioned_rice_contents_workspace_ptr_mid_side[i][1] = &encoder->private_->partitioned_rice_contents_workspace_mid_side[i][1];
547         }
548
549         for(i = 0; i < FLAC__MAX_CHANNELS; i++) {
550                 FLAC__format_entropy_coding_method_partitioned_rice_contents_init(&encoder->private_->partitioned_rice_contents_workspace[i][0]);
551                 FLAC__format_entropy_coding_method_partitioned_rice_contents_init(&encoder->private_->partitioned_rice_contents_workspace[i][1]);
552         }
553         for(i = 0; i < 2; i++) {
554                 FLAC__format_entropy_coding_method_partitioned_rice_contents_init(&encoder->private_->partitioned_rice_contents_workspace_mid_side[i][0]);
555                 FLAC__format_entropy_coding_method_partitioned_rice_contents_init(&encoder->private_->partitioned_rice_contents_workspace_mid_side[i][1]);
556         }
557         for(i = 0; i < 2; i++)
558                 FLAC__format_entropy_coding_method_partitioned_rice_contents_init(&encoder->private_->partitioned_rice_contents_extra[i]);
559
560         encoder->protected_->state = FLAC__STREAM_ENCODER_UNINITIALIZED;
561
562         return encoder;
563 }
564
565 FLAC_API void FLAC__stream_encoder_delete(FLAC__StreamEncoder *encoder)
566 {
567         unsigned i;
568
569         FLAC__ASSERT(0 != encoder);
570         FLAC__ASSERT(0 != encoder->protected_);
571         FLAC__ASSERT(0 != encoder->private_);
572         FLAC__ASSERT(0 != encoder->private_->frame);
573
574         encoder->private_->is_being_deleted = true;
575
576         FLAC__stream_encoder_finish(encoder);
577
578         if(0 != encoder->private_->verify.decoder)
579                 FLAC__stream_decoder_delete(encoder->private_->verify.decoder);
580
581         for(i = 0; i < FLAC__MAX_CHANNELS; i++) {
582                 FLAC__format_entropy_coding_method_partitioned_rice_contents_clear(&encoder->private_->partitioned_rice_contents_workspace[i][0]);
583                 FLAC__format_entropy_coding_method_partitioned_rice_contents_clear(&encoder->private_->partitioned_rice_contents_workspace[i][1]);
584         }
585         for(i = 0; i < 2; i++) {
586                 FLAC__format_entropy_coding_method_partitioned_rice_contents_clear(&encoder->private_->partitioned_rice_contents_workspace_mid_side[i][0]);
587                 FLAC__format_entropy_coding_method_partitioned_rice_contents_clear(&encoder->private_->partitioned_rice_contents_workspace_mid_side[i][1]);
588         }
589         for(i = 0; i < 2; i++)
590                 FLAC__format_entropy_coding_method_partitioned_rice_contents_clear(&encoder->private_->partitioned_rice_contents_extra[i]);
591
592         FLAC__bitbuffer_delete(encoder->private_->frame);
593         free(encoder->private_);
594         free(encoder->protected_);
595         free(encoder);
596 }
597
598 /***********************************************************************
599  *
600  * Public class methods
601  *
602  ***********************************************************************/
603
604 FLAC_API FLAC__StreamEncoderState FLAC__stream_encoder_init(FLAC__StreamEncoder *encoder)
605 {
606         unsigned i;
607         FLAC__bool metadata_has_seektable, metadata_has_vorbis_comment;
608
609         FLAC__ASSERT(0 != encoder);
610
611         if(encoder->protected_->state != FLAC__STREAM_ENCODER_UNINITIALIZED)
612                 return encoder->protected_->state = FLAC__STREAM_ENCODER_ALREADY_INITIALIZED;
613
614         encoder->protected_->state = FLAC__STREAM_ENCODER_OK;
615
616         if(0 == encoder->private_->write_callback || 0 == encoder->private_->metadata_callback)
617                 return encoder->protected_->state = FLAC__STREAM_ENCODER_INVALID_CALLBACK;
618
619         if(encoder->protected_->channels == 0 || encoder->protected_->channels > FLAC__MAX_CHANNELS)
620                 return encoder->protected_->state = FLAC__STREAM_ENCODER_INVALID_NUMBER_OF_CHANNELS;
621
622         if(encoder->protected_->do_mid_side_stereo && encoder->protected_->channels != 2)
623                 return encoder->protected_->state = FLAC__STREAM_ENCODER_MID_SIDE_CHANNELS_MISMATCH;
624
625         if(encoder->protected_->loose_mid_side_stereo && !encoder->protected_->do_mid_side_stereo)
626                 return encoder->protected_->state = FLAC__STREAM_ENCODER_ILLEGAL_MID_SIDE_FORCE;
627
628         if(encoder->protected_->bits_per_sample >= 32)
629                 encoder->protected_->do_mid_side_stereo = false; /* since we do 32-bit math, the side channel would have 33 bps and overflow */
630
631         if(encoder->protected_->bits_per_sample < FLAC__MIN_BITS_PER_SAMPLE || encoder->protected_->bits_per_sample > FLAC__REFERENCE_CODEC_MAX_BITS_PER_SAMPLE)
632                 return encoder->protected_->state = FLAC__STREAM_ENCODER_INVALID_BITS_PER_SAMPLE;
633
634         if(!FLAC__format_sample_rate_is_valid(encoder->protected_->sample_rate))
635                 return encoder->protected_->state = FLAC__STREAM_ENCODER_INVALID_SAMPLE_RATE;
636
637         if(encoder->protected_->blocksize < FLAC__MIN_BLOCK_SIZE || encoder->protected_->blocksize > FLAC__MAX_BLOCK_SIZE)
638                 return encoder->protected_->state = FLAC__STREAM_ENCODER_INVALID_BLOCK_SIZE;
639
640         if(encoder->protected_->max_lpc_order > FLAC__MAX_LPC_ORDER)
641                 return encoder->protected_->state = FLAC__STREAM_ENCODER_INVALID_MAX_LPC_ORDER;
642
643         if(encoder->protected_->blocksize < encoder->protected_->max_lpc_order)
644                 return encoder->protected_->state = FLAC__STREAM_ENCODER_BLOCK_SIZE_TOO_SMALL_FOR_LPC_ORDER;
645
646         if(encoder->protected_->qlp_coeff_precision == 0) {
647                 if(encoder->protected_->bits_per_sample < 16) {
648                         /* @@@ need some data about how to set this here w.r.t. blocksize and sample rate */
649                         /* @@@ until then we'll make a guess */
650                         encoder->protected_->qlp_coeff_precision = max(FLAC__MIN_QLP_COEFF_PRECISION, 2 + encoder->protected_->bits_per_sample / 2);
651                 }
652                 else if(encoder->protected_->bits_per_sample == 16) {
653                         if(encoder->protected_->blocksize <= 192)
654                                 encoder->protected_->qlp_coeff_precision = 7;
655                         else if(encoder->protected_->blocksize <= 384)
656                                 encoder->protected_->qlp_coeff_precision = 8;
657                         else if(encoder->protected_->blocksize <= 576)
658                                 encoder->protected_->qlp_coeff_precision = 9;
659                         else if(encoder->protected_->blocksize <= 1152)
660                                 encoder->protected_->qlp_coeff_precision = 10;
661                         else if(encoder->protected_->blocksize <= 2304)
662                                 encoder->protected_->qlp_coeff_precision = 11;
663                         else if(encoder->protected_->blocksize <= 4608)
664                                 encoder->protected_->qlp_coeff_precision = 12;
665                         else
666                                 encoder->protected_->qlp_coeff_precision = 13;
667                 }
668                 else {
669                         if(encoder->protected_->blocksize <= 384)
670                                 encoder->protected_->qlp_coeff_precision = FLAC__MAX_QLP_COEFF_PRECISION-2;
671                         else if(encoder->protected_->blocksize <= 1152)
672                                 encoder->protected_->qlp_coeff_precision = FLAC__MAX_QLP_COEFF_PRECISION-1;
673                         else
674                                 encoder->protected_->qlp_coeff_precision = FLAC__MAX_QLP_COEFF_PRECISION;
675                 }
676                 FLAC__ASSERT(encoder->protected_->qlp_coeff_precision <= FLAC__MAX_QLP_COEFF_PRECISION);
677         }
678         else if(encoder->protected_->qlp_coeff_precision < FLAC__MIN_QLP_COEFF_PRECISION || encoder->protected_->qlp_coeff_precision > FLAC__MAX_QLP_COEFF_PRECISION)
679                 return encoder->protected_->state = FLAC__STREAM_ENCODER_INVALID_QLP_COEFF_PRECISION;
680
681         if(encoder->protected_->streamable_subset) {
682                 if(
683                         encoder->protected_->blocksize != 192 &&
684                         encoder->protected_->blocksize != 576 &&
685                         encoder->protected_->blocksize != 1152 &&
686                         encoder->protected_->blocksize != 2304 &&
687                         encoder->protected_->blocksize != 4608 &&
688                         encoder->protected_->blocksize != 256 &&
689                         encoder->protected_->blocksize != 512 &&
690                         encoder->protected_->blocksize != 1024 &&
691                         encoder->protected_->blocksize != 2048 &&
692                         encoder->protected_->blocksize != 4096 &&
693                         encoder->protected_->blocksize != 8192 &&
694                         encoder->protected_->blocksize != 16384
695                 )
696                         return encoder->protected_->state = FLAC__STREAM_ENCODER_NOT_STREAMABLE;
697                 if(
698                         encoder->protected_->sample_rate != 8000 &&
699                         encoder->protected_->sample_rate != 16000 &&
700                         encoder->protected_->sample_rate != 22050 &&
701                         encoder->protected_->sample_rate != 24000 &&
702                         encoder->protected_->sample_rate != 32000 &&
703                         encoder->protected_->sample_rate != 44100 &&
704                         encoder->protected_->sample_rate != 48000 &&
705                         encoder->protected_->sample_rate != 96000
706                 )
707                         return encoder->protected_->state = FLAC__STREAM_ENCODER_NOT_STREAMABLE;
708                 if(
709                         encoder->protected_->bits_per_sample != 8 &&
710                         encoder->protected_->bits_per_sample != 12 &&
711                         encoder->protected_->bits_per_sample != 16 &&
712                         encoder->protected_->bits_per_sample != 20 &&
713                         encoder->protected_->bits_per_sample != 24
714                 )
715                         return encoder->protected_->state = FLAC__STREAM_ENCODER_NOT_STREAMABLE;
716                 if(encoder->protected_->max_residual_partition_order > FLAC__SUBSET_MAX_RICE_PARTITION_ORDER)
717                         return encoder->protected_->state = FLAC__STREAM_ENCODER_NOT_STREAMABLE;
718         }
719
720         if(encoder->protected_->max_residual_partition_order >= (1u << FLAC__ENTROPY_CODING_METHOD_PARTITIONED_RICE_ORDER_LEN))
721                 encoder->protected_->max_residual_partition_order = (1u << FLAC__ENTROPY_CODING_METHOD_PARTITIONED_RICE_ORDER_LEN) - 1;
722         if(encoder->protected_->min_residual_partition_order >= encoder->protected_->max_residual_partition_order)
723                 encoder->protected_->min_residual_partition_order = encoder->protected_->max_residual_partition_order;
724
725         /* validate metadata */
726         if(0 == encoder->protected_->metadata && encoder->protected_->num_metadata_blocks > 0)
727                 return encoder->protected_->state = FLAC__STREAM_ENCODER_INVALID_METADATA;
728         metadata_has_seektable = false;
729         metadata_has_vorbis_comment = false;
730         for(i = 0; i < encoder->protected_->num_metadata_blocks; i++) {
731                 if(encoder->protected_->metadata[i]->type == FLAC__METADATA_TYPE_STREAMINFO)
732                         return encoder->protected_->state = FLAC__STREAM_ENCODER_INVALID_METADATA;
733                 else if(encoder->protected_->metadata[i]->type == FLAC__METADATA_TYPE_SEEKTABLE) {
734                         if(metadata_has_seektable) /* only one is allowed */
735                                 return encoder->protected_->state = FLAC__STREAM_ENCODER_INVALID_METADATA;
736                         metadata_has_seektable = true;
737                         if(!FLAC__format_seektable_is_legal(&encoder->protected_->metadata[i]->data.seek_table))
738                                 return encoder->protected_->state = FLAC__STREAM_ENCODER_INVALID_METADATA;
739                 }
740                 else if(encoder->protected_->metadata[i]->type == FLAC__METADATA_TYPE_VORBIS_COMMENT) {
741                         if(metadata_has_vorbis_comment) /* only one is allowed */
742                                 return encoder->protected_->state = FLAC__STREAM_ENCODER_INVALID_METADATA;
743                         metadata_has_vorbis_comment = true;
744                 }
745                 else if(encoder->protected_->metadata[i]->type == FLAC__METADATA_TYPE_CUESHEET) {
746                         if(!FLAC__format_cuesheet_is_legal(&encoder->protected_->metadata[i]->data.cue_sheet, encoder->protected_->metadata[i]->data.cue_sheet.is_cd, /*violation=*/0))
747                                 return encoder->protected_->state = FLAC__STREAM_ENCODER_INVALID_METADATA;
748                 }
749         }
750
751         encoder->private_->input_capacity = 0;
752         for(i = 0; i < encoder->protected_->channels; i++) {
753                 encoder->private_->integer_signal_unaligned[i] = encoder->private_->integer_signal[i] = 0;
754 #ifndef FLAC__INTEGER_ONLY_LIBRARY
755                 encoder->private_->real_signal_unaligned[i] = encoder->private_->real_signal[i] = 0;
756 #endif
757         }
758         for(i = 0; i < 2; i++) {
759                 encoder->private_->integer_signal_mid_side_unaligned[i] = encoder->private_->integer_signal_mid_side[i] = 0;
760 #ifndef FLAC__INTEGER_ONLY_LIBRARY
761                 encoder->private_->real_signal_mid_side_unaligned[i] = encoder->private_->real_signal_mid_side[i] = 0;
762 #endif
763         }
764 #ifndef FLAC__INTEGER_ONLY_LIBRARY
765         for(i = 0; i < encoder->protected_->num_apodizations; i++)
766                 encoder->private_->window_unaligned[i] = encoder->private_->window[i] = 0;
767         encoder->private_->windowed_signal_unaligned = encoder->private_->windowed_signal = 0;
768 #endif
769         for(i = 0; i < encoder->protected_->channels; i++) {
770                 encoder->private_->residual_workspace_unaligned[i][0] = encoder->private_->residual_workspace[i][0] = 0;
771                 encoder->private_->residual_workspace_unaligned[i][1] = encoder->private_->residual_workspace[i][1] = 0;
772                 encoder->private_->best_subframe[i] = 0;
773         }
774         for(i = 0; i < 2; i++) {
775                 encoder->private_->residual_workspace_mid_side_unaligned[i][0] = encoder->private_->residual_workspace_mid_side[i][0] = 0;
776                 encoder->private_->residual_workspace_mid_side_unaligned[i][1] = encoder->private_->residual_workspace_mid_side[i][1] = 0;
777                 encoder->private_->best_subframe_mid_side[i] = 0;
778         }
779         encoder->private_->abs_residual_unaligned = encoder->private_->abs_residual = 0;
780         encoder->private_->abs_residual_partition_sums_unaligned = encoder->private_->abs_residual_partition_sums = 0;
781         encoder->private_->raw_bits_per_partition_unaligned = encoder->private_->raw_bits_per_partition = 0;
782 #ifndef FLAC__INTEGER_ONLY_LIBRARY
783         encoder->private_->loose_mid_side_stereo_frames = (unsigned)((FLAC__double)encoder->protected_->sample_rate * 0.4 / (FLAC__double)encoder->protected_->blocksize + 0.5);
784 #else
785         /* 26214 is the approximate fixed-point equivalent to 0.4 (0.4 * 2^16) */
786         /* sample rate can be up to 655350 Hz, and thus use 20 bits, so we do the multiply&divide by hand */
787         FLAC__ASSERT(FLAC__MAX_SAMPLE_RATE <= 655350);
788         FLAC__ASSERT(FLAC__MAX_BLOCK_SIZE <= 65535);
789         FLAC__ASSERT(encoder->protected_->sample_rate <= 655350);
790         FLAC__ASSERT(encoder->protected_->blocksize <= 65535);
791         encoder->private_->loose_mid_side_stereo_frames = (unsigned)FLAC__fixedpoint_trunc((((FLAC__uint64)(encoder->protected_->sample_rate) * (FLAC__uint64)(26214)) << 16) / (encoder->protected_->blocksize<<16) + FLAC__FP_ONE_HALF);
792 #endif
793         if(encoder->private_->loose_mid_side_stereo_frames == 0)
794                 encoder->private_->loose_mid_side_stereo_frames = 1;
795         encoder->private_->loose_mid_side_stereo_frame_count = 0;
796         encoder->private_->current_sample_number = 0;
797         encoder->private_->current_frame_number = 0;
798
799         encoder->private_->use_wide_by_block = (encoder->protected_->bits_per_sample + FLAC__bitmath_ilog2(encoder->protected_->blocksize)+1 > 30);
800         encoder->private_->use_wide_by_order = (encoder->protected_->bits_per_sample + FLAC__bitmath_ilog2(max(encoder->protected_->max_lpc_order, FLAC__MAX_FIXED_ORDER))+1 > 30); /*@@@ need to use this? */
801         encoder->private_->use_wide_by_partition = (false); /*@@@ need to set this */
802
803         /*
804          * get the CPU info and set the function pointers
805          */
806         FLAC__cpu_info(&encoder->private_->cpuinfo);
807         /* first default to the non-asm routines */
808 #ifndef FLAC__INTEGER_ONLY_LIBRARY
809         encoder->private_->local_lpc_compute_autocorrelation = FLAC__lpc_compute_autocorrelation;
810 #endif
811         encoder->private_->local_fixed_compute_best_predictor = FLAC__fixed_compute_best_predictor;
812 #ifndef FLAC__INTEGER_ONLY_LIBRARY
813         encoder->private_->local_lpc_compute_residual_from_qlp_coefficients = FLAC__lpc_compute_residual_from_qlp_coefficients;
814         encoder->private_->local_lpc_compute_residual_from_qlp_coefficients_64bit = FLAC__lpc_compute_residual_from_qlp_coefficients_wide;
815         encoder->private_->local_lpc_compute_residual_from_qlp_coefficients_16bit = FLAC__lpc_compute_residual_from_qlp_coefficients;
816 #endif
817         /* now override with asm where appropriate */
818 #ifndef FLAC__INTEGER_ONLY_LIBRARY
819 # ifndef FLAC__NO_ASM
820         if(encoder->private_->cpuinfo.use_asm) {
821 #  ifdef FLAC__CPU_IA32
822                 FLAC__ASSERT(encoder->private_->cpuinfo.type == FLAC__CPUINFO_TYPE_IA32);
823 #   ifdef FLAC__HAS_NASM
824 #    ifdef FLAC__SSE_OS
825                 if(encoder->private_->cpuinfo.data.ia32.sse) {
826                         if(encoder->protected_->max_lpc_order < 4)
827                                 encoder->private_->local_lpc_compute_autocorrelation = FLAC__lpc_compute_autocorrelation_asm_ia32_sse_lag_4;
828                         else if(encoder->protected_->max_lpc_order < 8)
829                                 encoder->private_->local_lpc_compute_autocorrelation = FLAC__lpc_compute_autocorrelation_asm_ia32_sse_lag_8;
830                         else if(encoder->protected_->max_lpc_order < 12)
831                                 encoder->private_->local_lpc_compute_autocorrelation = FLAC__lpc_compute_autocorrelation_asm_ia32_sse_lag_12;
832                         else
833                                 encoder->private_->local_lpc_compute_autocorrelation = FLAC__lpc_compute_autocorrelation_asm_ia32;
834                 }
835                 else
836 #    endif /* FLAC__SSE_OS */
837                 if(encoder->private_->cpuinfo.data.ia32._3dnow)
838                         encoder->private_->local_lpc_compute_autocorrelation = FLAC__lpc_compute_autocorrelation_asm_ia32_3dnow;
839                 else
840                         encoder->private_->local_lpc_compute_autocorrelation = FLAC__lpc_compute_autocorrelation_asm_ia32;
841                 if(encoder->private_->cpuinfo.data.ia32.mmx) {
842                         encoder->private_->local_lpc_compute_residual_from_qlp_coefficients = FLAC__lpc_compute_residual_from_qlp_coefficients_asm_ia32;
843                         encoder->private_->local_lpc_compute_residual_from_qlp_coefficients_16bit = FLAC__lpc_compute_residual_from_qlp_coefficients_asm_ia32_mmx;
844                 }
845                 else {
846                         encoder->private_->local_lpc_compute_residual_from_qlp_coefficients = FLAC__lpc_compute_residual_from_qlp_coefficients_asm_ia32;
847                         encoder->private_->local_lpc_compute_residual_from_qlp_coefficients_16bit = FLAC__lpc_compute_residual_from_qlp_coefficients_asm_ia32;
848                 }
849                 if(encoder->private_->cpuinfo.data.ia32.mmx && encoder->private_->cpuinfo.data.ia32.cmov)
850                         encoder->private_->local_fixed_compute_best_predictor = FLAC__fixed_compute_best_predictor_asm_ia32_mmx_cmov;
851 #   endif /* FLAC__HAS_NASM */
852 #  endif /* FLAC__CPU_IA32 */
853         }
854 # endif /* !FLAC__NO_ASM */
855 #endif /* !FLAC__INTEGER_ONLY_LIBRARY */
856         /* finally override based on wide-ness if necessary */
857         if(encoder->private_->use_wide_by_block) {
858                 encoder->private_->local_fixed_compute_best_predictor = FLAC__fixed_compute_best_predictor_wide;
859         }
860
861         /* we require precompute_partition_sums if do_escape_coding because of their intertwined nature */
862         encoder->private_->precompute_partition_sums = (encoder->protected_->max_residual_partition_order > encoder->protected_->min_residual_partition_order) || encoder->protected_->do_escape_coding;
863
864         if(!resize_buffers_(encoder, encoder->protected_->blocksize)) {
865                 /* the above function sets the state for us in case of an error */
866                 return encoder->protected_->state;
867         }
868
869         if(!FLAC__bitbuffer_init(encoder->private_->frame))
870                 return encoder->protected_->state = FLAC__STREAM_ENCODER_MEMORY_ALLOCATION_ERROR;
871
872         /*
873          * Set up the verify stuff if necessary
874          */
875         if(encoder->protected_->verify) {
876                 /*
877                  * First, set up the fifo which will hold the
878                  * original signal to compare against
879                  */
880                 encoder->private_->verify.input_fifo.size = encoder->protected_->blocksize;
881                 for(i = 0; i < encoder->protected_->channels; i++) {
882                         if(0 == (encoder->private_->verify.input_fifo.data[i] = (FLAC__int32*)malloc(sizeof(FLAC__int32) * encoder->private_->verify.input_fifo.size)))
883                                 return encoder->protected_->state = FLAC__STREAM_ENCODER_MEMORY_ALLOCATION_ERROR;
884                 }
885                 encoder->private_->verify.input_fifo.tail = 0;
886
887                 /*
888                  * Now set up a stream decoder for verification
889                  */
890                 encoder->private_->verify.decoder = FLAC__stream_decoder_new();
891                 if(0 == encoder->private_->verify.decoder)
892                         return encoder->protected_->state = FLAC__STREAM_ENCODER_VERIFY_DECODER_ERROR;
893
894                 FLAC__stream_decoder_set_read_callback(encoder->private_->verify.decoder, verify_read_callback_);
895                 FLAC__stream_decoder_set_write_callback(encoder->private_->verify.decoder, verify_write_callback_);
896                 FLAC__stream_decoder_set_metadata_callback(encoder->private_->verify.decoder, verify_metadata_callback_);
897                 FLAC__stream_decoder_set_error_callback(encoder->private_->verify.decoder, verify_error_callback_);
898                 FLAC__stream_decoder_set_client_data(encoder->private_->verify.decoder, encoder);
899                 if(FLAC__stream_decoder_init(encoder->private_->verify.decoder) != FLAC__STREAM_DECODER_SEARCH_FOR_METADATA)
900                         return encoder->protected_->state = FLAC__STREAM_ENCODER_VERIFY_DECODER_ERROR;
901         }
902         encoder->private_->verify.error_stats.absolute_sample = 0;
903         encoder->private_->verify.error_stats.frame_number = 0;
904         encoder->private_->verify.error_stats.channel = 0;
905         encoder->private_->verify.error_stats.sample = 0;
906         encoder->private_->verify.error_stats.expected = 0;
907         encoder->private_->verify.error_stats.got = 0;
908
909         /*
910          * write the stream header
911          */
912         if(encoder->protected_->verify)
913                 encoder->private_->verify.state_hint = ENCODER_IN_MAGIC;
914         if(!FLAC__bitbuffer_write_raw_uint32(encoder->private_->frame, FLAC__STREAM_SYNC, FLAC__STREAM_SYNC_LEN))
915                 return encoder->protected_->state = FLAC__STREAM_ENCODER_FRAMING_ERROR;
916         if(!write_bitbuffer_(encoder, 0)) {
917                 /* the above function sets the state for us in case of an error */
918                 return encoder->protected_->state;
919         }
920
921         /*
922          * write the STREAMINFO metadata block
923          */
924         if(encoder->protected_->verify)
925                 encoder->private_->verify.state_hint = ENCODER_IN_METADATA;
926         encoder->private_->metadata.type = FLAC__METADATA_TYPE_STREAMINFO;
927         encoder->private_->metadata.is_last = false; /* we will have at a minimum a VORBIS_COMMENT afterwards */
928         encoder->private_->metadata.length = FLAC__STREAM_METADATA_STREAMINFO_LENGTH;
929         encoder->private_->metadata.data.stream_info.min_blocksize = encoder->protected_->blocksize; /* this encoder uses the same blocksize for the whole stream */
930         encoder->private_->metadata.data.stream_info.max_blocksize = encoder->protected_->blocksize;
931         encoder->private_->metadata.data.stream_info.min_framesize = 0; /* we don't know this yet; have to fill it in later */
932         encoder->private_->metadata.data.stream_info.max_framesize = 0; /* we don't know this yet; have to fill it in later */
933         encoder->private_->metadata.data.stream_info.sample_rate = encoder->protected_->sample_rate;
934         encoder->private_->metadata.data.stream_info.channels = encoder->protected_->channels;
935         encoder->private_->metadata.data.stream_info.bits_per_sample = encoder->protected_->bits_per_sample;
936         encoder->private_->metadata.data.stream_info.total_samples = encoder->protected_->total_samples_estimate; /* we will replace this later with the real total */
937         memset(encoder->private_->metadata.data.stream_info.md5sum, 0, 16); /* we don't know this yet; have to fill it in later */
938         FLAC__MD5Init(&encoder->private_->md5context);
939         if(!FLAC__bitbuffer_clear(encoder->private_->frame))
940                 return encoder->protected_->state = FLAC__STREAM_ENCODER_MEMORY_ALLOCATION_ERROR;
941         if(!FLAC__add_metadata_block(&encoder->private_->metadata, encoder->private_->frame))
942                 return encoder->protected_->state = FLAC__STREAM_ENCODER_FRAMING_ERROR;
943         if(!write_bitbuffer_(encoder, 0)) {
944                 /* the above function sets the state for us in case of an error */
945                 return encoder->protected_->state;
946         }
947
948         /*
949          * Now that the STREAMINFO block is written, we can init this to an
950          * absurdly-high value...
951          */
952         encoder->private_->metadata.data.stream_info.min_framesize = (1u << FLAC__STREAM_METADATA_STREAMINFO_MIN_FRAME_SIZE_LEN) - 1;
953         /* ... and clear this to 0 */
954         encoder->private_->metadata.data.stream_info.total_samples = 0;
955
956         /*
957          * Check to see if the supplied metadata contains a VORBIS_COMMENT;
958          * if not, we will write an empty one (FLAC__add_metadata_block()
959          * automatically supplies the vendor string).
960          *
961          * WATCHOUT: libOggFLAC depends on us to write this block after the
962          * STREAMINFO since that's what the mapping requires.  (In the case
963          * that metadata_has_vorbis_comment is true it will have already
964          * insured that the metadata list is properly ordered.)
965          */
966         if(!metadata_has_vorbis_comment) {
967                 FLAC__StreamMetadata vorbis_comment;
968                 vorbis_comment.type = FLAC__METADATA_TYPE_VORBIS_COMMENT;
969                 vorbis_comment.is_last = (encoder->protected_->num_metadata_blocks == 0);
970                 vorbis_comment.length = 4 + 4; /* MAGIC NUMBER */
971                 vorbis_comment.data.vorbis_comment.vendor_string.length = 0;
972                 vorbis_comment.data.vorbis_comment.vendor_string.entry = 0;
973                 vorbis_comment.data.vorbis_comment.num_comments = 0;
974                 vorbis_comment.data.vorbis_comment.comments = 0;
975                 if(!FLAC__bitbuffer_clear(encoder->private_->frame))
976                         return encoder->protected_->state = FLAC__STREAM_ENCODER_MEMORY_ALLOCATION_ERROR;
977                 if(!FLAC__add_metadata_block(&vorbis_comment, encoder->private_->frame))
978                         return encoder->protected_->state = FLAC__STREAM_ENCODER_FRAMING_ERROR;
979                 if(!write_bitbuffer_(encoder, 0)) {
980                         /* the above function sets the state for us in case of an error */
981                         return encoder->protected_->state;
982                 }
983         }
984
985         /*
986          * write the user's metadata blocks
987          */
988         for(i = 0; i < encoder->protected_->num_metadata_blocks; i++) {
989                 encoder->protected_->metadata[i]->is_last = (i == encoder->protected_->num_metadata_blocks - 1);
990                 if(!FLAC__bitbuffer_clear(encoder->private_->frame))
991                         return encoder->protected_->state = FLAC__STREAM_ENCODER_MEMORY_ALLOCATION_ERROR;
992                 if(!FLAC__add_metadata_block(encoder->protected_->metadata[i], encoder->private_->frame))
993                         return encoder->protected_->state = FLAC__STREAM_ENCODER_FRAMING_ERROR;
994                 if(!write_bitbuffer_(encoder, 0)) {
995                         /* the above function sets the state for us in case of an error */
996                         return encoder->protected_->state;
997                 }
998         }
999
1000         if(encoder->protected_->verify)
1001                 encoder->private_->verify.state_hint = ENCODER_IN_AUDIO;
1002
1003         return encoder->protected_->state;
1004 }
1005
1006 FLAC_API void FLAC__stream_encoder_finish(FLAC__StreamEncoder *encoder)
1007 {
1008         FLAC__ASSERT(0 != encoder);
1009
1010         if(encoder->protected_->state == FLAC__STREAM_ENCODER_UNINITIALIZED)
1011                 return;
1012
1013         if(encoder->protected_->state == FLAC__STREAM_ENCODER_OK && !encoder->private_->is_being_deleted) {
1014                 if(encoder->private_->current_sample_number != 0) {
1015                         encoder->protected_->blocksize = encoder->private_->current_sample_number;
1016                         process_frame_(encoder, true); /* true => is last frame */
1017                 }
1018         }
1019
1020         FLAC__MD5Final(encoder->private_->metadata.data.stream_info.md5sum, &encoder->private_->md5context);
1021
1022         if(encoder->protected_->state == FLAC__STREAM_ENCODER_OK && !encoder->private_->is_being_deleted) {
1023                 encoder->private_->metadata_callback(encoder, &encoder->private_->metadata, encoder->private_->client_data);
1024         }
1025
1026         if(encoder->protected_->verify && 0 != encoder->private_->verify.decoder)
1027                 FLAC__stream_decoder_finish(encoder->private_->verify.decoder);
1028
1029         free_(encoder);
1030         set_defaults_(encoder);
1031
1032         encoder->protected_->state = FLAC__STREAM_ENCODER_UNINITIALIZED;
1033 }
1034
1035 FLAC_API FLAC__bool FLAC__stream_encoder_set_verify(FLAC__StreamEncoder *encoder, FLAC__bool value)
1036 {
1037         FLAC__ASSERT(0 != encoder);
1038         if(encoder->protected_->state != FLAC__STREAM_ENCODER_UNINITIALIZED)
1039                 return false;
1040 #ifndef FLAC__MANDATORY_VERIFY_WHILE_ENCODING
1041         encoder->protected_->verify = value;
1042 #endif
1043         return true;
1044 }
1045
1046 FLAC_API FLAC__bool FLAC__stream_encoder_set_streamable_subset(FLAC__StreamEncoder *encoder, FLAC__bool value)
1047 {
1048         FLAC__ASSERT(0 != encoder);
1049         if(encoder->protected_->state != FLAC__STREAM_ENCODER_UNINITIALIZED)
1050                 return false;
1051         encoder->protected_->streamable_subset = value;
1052         return true;
1053 }
1054
1055 FLAC_API FLAC__bool FLAC__stream_encoder_set_do_mid_side_stereo(FLAC__StreamEncoder *encoder, FLAC__bool value)
1056 {
1057         FLAC__ASSERT(0 != encoder);
1058         if(encoder->protected_->state != FLAC__STREAM_ENCODER_UNINITIALIZED)
1059                 return false;
1060         encoder->protected_->do_mid_side_stereo = value;
1061         return true;
1062 }
1063
1064 FLAC_API FLAC__bool FLAC__stream_encoder_set_loose_mid_side_stereo(FLAC__StreamEncoder *encoder, FLAC__bool value)
1065 {
1066         FLAC__ASSERT(0 != encoder);
1067         if(encoder->protected_->state != FLAC__STREAM_ENCODER_UNINITIALIZED)
1068                 return false;
1069         encoder->protected_->loose_mid_side_stereo = value;
1070         return true;
1071 }
1072
1073 FLAC_API FLAC__bool FLAC__stream_encoder_set_channels(FLAC__StreamEncoder *encoder, unsigned value)
1074 {
1075         FLAC__ASSERT(0 != encoder);
1076         if(encoder->protected_->state != FLAC__STREAM_ENCODER_UNINITIALIZED)
1077                 return false;
1078         encoder->protected_->channels = value;
1079         return true;
1080 }
1081
1082 FLAC_API FLAC__bool FLAC__stream_encoder_set_bits_per_sample(FLAC__StreamEncoder *encoder, unsigned value)
1083 {
1084         FLAC__ASSERT(0 != encoder);
1085         if(encoder->protected_->state != FLAC__STREAM_ENCODER_UNINITIALIZED)
1086                 return false;
1087         encoder->protected_->bits_per_sample = value;
1088         return true;
1089 }
1090
1091 FLAC_API FLAC__bool FLAC__stream_encoder_set_sample_rate(FLAC__StreamEncoder *encoder, unsigned value)
1092 {
1093         FLAC__ASSERT(0 != encoder);
1094         if(encoder->protected_->state != FLAC__STREAM_ENCODER_UNINITIALIZED)
1095                 return false;
1096         encoder->protected_->sample_rate = value;
1097         return true;
1098 }
1099
1100 FLAC_API FLAC__bool FLAC__stream_encoder_set_blocksize(FLAC__StreamEncoder *encoder, unsigned value)
1101 {
1102         FLAC__ASSERT(0 != encoder);
1103         if(encoder->protected_->state != FLAC__STREAM_ENCODER_UNINITIALIZED)
1104                 return false;
1105         encoder->protected_->blocksize = value;
1106         return true;
1107 }
1108
1109 FLAC_API FLAC__bool FLAC__stream_encoder_set_apodization(FLAC__StreamEncoder *encoder, const char *specification)
1110 {
1111         FLAC__ASSERT(0 != encoder);
1112         FLAC__ASSERT(0 != specification);
1113         if(encoder->protected_->state != FLAC__STREAM_ENCODER_UNINITIALIZED)
1114                 return false;
1115 #ifdef FLAC__INTEGER_ONLY_LIBRARY
1116         (void)specification; /* silently ignore since we haven't integerized; will always use a rectangular window */
1117 #else
1118         encoder->protected_->num_apodizations = 0;
1119         while(1) {
1120                 const char *s = strchr(specification, ';');
1121                 const size_t n = s? (size_t)(s - specification) : strlen(specification);
1122                 if     (n==8  && 0 == strncmp("bartlett"     , specification, n))
1123                         encoder->protected_->apodizations[encoder->protected_->num_apodizations++].type = FLAC__APODIZATION_BARTLETT;
1124                 else if(n==13 && 0 == strncmp("bartlett_hann", specification, n))
1125                         encoder->protected_->apodizations[encoder->protected_->num_apodizations++].type = FLAC__APODIZATION_BARTLETT_HANN;
1126                 else if(n==8  && 0 == strncmp("blackman"     , specification, n))
1127                         encoder->protected_->apodizations[encoder->protected_->num_apodizations++].type = FLAC__APODIZATION_BLACKMAN;
1128                 else if(n==26 && 0 == strncmp("blackman_harris_4term_92db", specification, n))
1129                         encoder->protected_->apodizations[encoder->protected_->num_apodizations++].type = FLAC__APODIZATION_BLACKMAN_HARRIS_4TERM_92DB_SIDELOBE;
1130                 else if(n==6  && 0 == strncmp("connes"       , specification, n))
1131                         encoder->protected_->apodizations[encoder->protected_->num_apodizations++].type = FLAC__APODIZATION_CONNES;
1132                 else if(n==7  && 0 == strncmp("flattop"      , specification, n))
1133                         encoder->protected_->apodizations[encoder->protected_->num_apodizations++].type = FLAC__APODIZATION_FLATTOP;
1134                 else if(n>7   && 0 == strncmp("gauss("       , specification, 6)) {
1135                         FLAC__real stddev = (FLAC__real)strtod(specification+6, 0);
1136                         if (stddev > 0.0 && stddev <= 0.5) {
1137                                 encoder->protected_->apodizations[encoder->protected_->num_apodizations].parameters.gauss.stddev = stddev;
1138                                 encoder->protected_->apodizations[encoder->protected_->num_apodizations++].type = FLAC__APODIZATION_GAUSS;
1139                         }
1140                 }
1141                 else if(n==7  && 0 == strncmp("hamming"      , specification, n))
1142                         encoder->protected_->apodizations[encoder->protected_->num_apodizations++].type = FLAC__APODIZATION_HAMMING;
1143                 else if(n==4  && 0 == strncmp("hann"         , specification, n))
1144                         encoder->protected_->apodizations[encoder->protected_->num_apodizations++].type = FLAC__APODIZATION_HANN;
1145                 else if(n==13 && 0 == strncmp("kaiser_bessel", specification, n))
1146                         encoder->protected_->apodizations[encoder->protected_->num_apodizations++].type = FLAC__APODIZATION_KAISER_BESSEL;
1147                 else if(n==7  && 0 == strncmp("nuttall"      , specification, n))
1148                         encoder->protected_->apodizations[encoder->protected_->num_apodizations++].type = FLAC__APODIZATION_NUTTALL;
1149                 else if(n==9  && 0 == strncmp("rectangle"    , specification, n))
1150                         encoder->protected_->apodizations[encoder->protected_->num_apodizations++].type = FLAC__APODIZATION_RECTANGLE;
1151                 else if(n==8  && 0 == strncmp("triangle"     , specification, n))
1152                         encoder->protected_->apodizations[encoder->protected_->num_apodizations++].type = FLAC__APODIZATION_TRIANGLE;
1153                 else if(n>7   && 0 == strncmp("tukey("       , specification, 6)) {
1154                         FLAC__real p = (FLAC__real)strtod(specification+6, 0);
1155                         if (p >= 0.0 && p <= 1.0) {
1156                                 encoder->protected_->apodizations[encoder->protected_->num_apodizations].parameters.tukey.p = p;
1157                                 encoder->protected_->apodizations[encoder->protected_->num_apodizations++].type = FLAC__APODIZATION_TUKEY;
1158                         }
1159                 }
1160                 else if(n==5  && 0 == strncmp("welch"        , specification, n))
1161                         encoder->protected_->apodizations[encoder->protected_->num_apodizations++].type = FLAC__APODIZATION_WELCH;
1162                 if (encoder->protected_->num_apodizations == 32)
1163                         break;
1164                 if (s)
1165                         specification = s+1;
1166                 else
1167                         break;
1168         }
1169         if(encoder->protected_->num_apodizations == 0) {
1170                 encoder->protected_->num_apodizations = 1;
1171                 encoder->protected_->apodizations[0].type = FLAC__APODIZATION_HANN;
1172         }
1173 #ifdef WINDOW_DEBUG_OUTPUT
1174 {unsigned n;for(n=0;n<encoder->protected_->num_apodizations;n++)fprintf(stderr,"@@@@@@ parsed apodization[%zu]: %s\n",n,winstr[encoder->protected_->apodizations[n].type]);}
1175 #endif
1176 #endif
1177         return true;
1178 }
1179
1180 FLAC_API FLAC__bool FLAC__stream_encoder_set_max_lpc_order(FLAC__StreamEncoder *encoder, unsigned value)
1181 {
1182         FLAC__ASSERT(0 != encoder);
1183         if(encoder->protected_->state != FLAC__STREAM_ENCODER_UNINITIALIZED)
1184                 return false;
1185         encoder->protected_->max_lpc_order = value;
1186         return true;
1187 }
1188
1189 FLAC_API FLAC__bool FLAC__stream_encoder_set_qlp_coeff_precision(FLAC__StreamEncoder *encoder, unsigned value)
1190 {
1191         FLAC__ASSERT(0 != encoder);
1192         if(encoder->protected_->state != FLAC__STREAM_ENCODER_UNINITIALIZED)
1193                 return false;
1194         encoder->protected_->qlp_coeff_precision = value;
1195         return true;
1196 }
1197
1198 FLAC_API FLAC__bool FLAC__stream_encoder_set_do_qlp_coeff_prec_search(FLAC__StreamEncoder *encoder, FLAC__bool value)
1199 {
1200         FLAC__ASSERT(0 != encoder);
1201         if(encoder->protected_->state != FLAC__STREAM_ENCODER_UNINITIALIZED)
1202                 return false;
1203         encoder->protected_->do_qlp_coeff_prec_search = value;
1204         return true;
1205 }
1206
1207 FLAC_API FLAC__bool FLAC__stream_encoder_set_do_escape_coding(FLAC__StreamEncoder *encoder, FLAC__bool value)
1208 {
1209         FLAC__ASSERT(0 != encoder);
1210         if(encoder->protected_->state != FLAC__STREAM_ENCODER_UNINITIALIZED)
1211                 return false;
1212 #if 0
1213         /*@@@ deprecated: */
1214         encoder->protected_->do_escape_coding = value;
1215 #else
1216         (void)value;
1217 #endif
1218         return true;
1219 }
1220
1221 FLAC_API FLAC__bool FLAC__stream_encoder_set_do_exhaustive_model_search(FLAC__StreamEncoder *encoder, FLAC__bool value)
1222 {
1223         FLAC__ASSERT(0 != encoder);
1224         if(encoder->protected_->state != FLAC__STREAM_ENCODER_UNINITIALIZED)
1225                 return false;
1226         encoder->protected_->do_exhaustive_model_search = value;
1227         return true;
1228 }
1229
1230 FLAC_API FLAC__bool FLAC__stream_encoder_set_min_residual_partition_order(FLAC__StreamEncoder *encoder, unsigned value)
1231 {
1232         FLAC__ASSERT(0 != encoder);
1233         if(encoder->protected_->state != FLAC__STREAM_ENCODER_UNINITIALIZED)
1234                 return false;
1235         encoder->protected_->min_residual_partition_order = value;
1236         return true;
1237 }
1238
1239 FLAC_API FLAC__bool FLAC__stream_encoder_set_max_residual_partition_order(FLAC__StreamEncoder *encoder, unsigned value)
1240 {
1241         FLAC__ASSERT(0 != encoder);
1242         if(encoder->protected_->state != FLAC__STREAM_ENCODER_UNINITIALIZED)
1243                 return false;
1244         encoder->protected_->max_residual_partition_order = value;
1245         return true;
1246 }
1247
1248 FLAC_API FLAC__bool FLAC__stream_encoder_set_rice_parameter_search_dist(FLAC__StreamEncoder *encoder, unsigned value)
1249 {
1250         FLAC__ASSERT(0 != encoder);
1251         if(encoder->protected_->state != FLAC__STREAM_ENCODER_UNINITIALIZED)
1252                 return false;
1253 #if 0
1254         /*@@@ deprecated: */
1255         encoder->protected_->rice_parameter_search_dist = value;
1256 #else
1257         (void)value;
1258 #endif
1259         return true;
1260 }
1261
1262 FLAC_API FLAC__bool FLAC__stream_encoder_set_total_samples_estimate(FLAC__StreamEncoder *encoder, FLAC__uint64 value)
1263 {
1264         FLAC__ASSERT(0 != encoder);
1265         if(encoder->protected_->state != FLAC__STREAM_ENCODER_UNINITIALIZED)
1266                 return false;
1267         encoder->protected_->total_samples_estimate = value;
1268         return true;
1269 }
1270
1271 FLAC_API FLAC__bool FLAC__stream_encoder_set_metadata(FLAC__StreamEncoder *encoder, FLAC__StreamMetadata **metadata, unsigned num_blocks)
1272 {
1273         FLAC__ASSERT(0 != encoder);
1274         if(encoder->protected_->state != FLAC__STREAM_ENCODER_UNINITIALIZED)
1275                 return false;
1276         encoder->protected_->metadata = metadata;
1277         encoder->protected_->num_metadata_blocks = num_blocks;
1278         return true;
1279 }
1280
1281 FLAC_API FLAC__bool FLAC__stream_encoder_set_write_callback(FLAC__StreamEncoder *encoder, FLAC__StreamEncoderWriteCallback value)
1282 {
1283         FLAC__ASSERT(0 != encoder);
1284         FLAC__ASSERT(0 != value);
1285         if(encoder->protected_->state != FLAC__STREAM_ENCODER_UNINITIALIZED)
1286                 return false;
1287         encoder->private_->write_callback = value;
1288         return true;
1289 }
1290
1291 FLAC_API FLAC__bool FLAC__stream_encoder_set_metadata_callback(FLAC__StreamEncoder *encoder, FLAC__StreamEncoderMetadataCallback value)
1292 {
1293         FLAC__ASSERT(0 != encoder);
1294         FLAC__ASSERT(0 != value);
1295         if(encoder->protected_->state != FLAC__STREAM_ENCODER_UNINITIALIZED)
1296                 return false;
1297         encoder->private_->metadata_callback = value;
1298         return true;
1299 }
1300
1301 FLAC_API FLAC__bool FLAC__stream_encoder_set_client_data(FLAC__StreamEncoder *encoder, void *value)
1302 {
1303         FLAC__ASSERT(0 != encoder);
1304         if(encoder->protected_->state != FLAC__STREAM_ENCODER_UNINITIALIZED)
1305                 return false;
1306         encoder->private_->client_data = value;
1307         return true;
1308 }
1309
1310 /*
1311  * These three functions are not static, but not publically exposed in
1312  * include/FLAC/ either.  They are used by the test suite.
1313  */
1314 FLAC_API FLAC__bool FLAC__stream_encoder_disable_constant_subframes(FLAC__StreamEncoder *encoder, FLAC__bool value)
1315 {
1316         FLAC__ASSERT(0 != encoder);
1317         if(encoder->protected_->state != FLAC__STREAM_ENCODER_UNINITIALIZED)
1318                 return false;
1319         encoder->private_->disable_constant_subframes = value;
1320         return true;
1321 }
1322
1323 FLAC_API FLAC__bool FLAC__stream_encoder_disable_fixed_subframes(FLAC__StreamEncoder *encoder, FLAC__bool value)
1324 {
1325         FLAC__ASSERT(0 != encoder);
1326         if(encoder->protected_->state != FLAC__STREAM_ENCODER_UNINITIALIZED)
1327                 return false;
1328         encoder->private_->disable_fixed_subframes = value;
1329         return true;
1330 }
1331
1332 FLAC_API FLAC__bool FLAC__stream_encoder_disable_verbatim_subframes(FLAC__StreamEncoder *encoder, FLAC__bool value)
1333 {
1334         FLAC__ASSERT(0 != encoder);
1335         if(encoder->protected_->state != FLAC__STREAM_ENCODER_UNINITIALIZED)
1336                 return false;
1337         encoder->private_->disable_verbatim_subframes = value;
1338         return true;
1339 }
1340
1341 FLAC_API FLAC__StreamEncoderState FLAC__stream_encoder_get_state(const FLAC__StreamEncoder *encoder)
1342 {
1343         FLAC__ASSERT(0 != encoder);
1344         return encoder->protected_->state;
1345 }
1346
1347 FLAC_API FLAC__StreamDecoderState FLAC__stream_encoder_get_verify_decoder_state(const FLAC__StreamEncoder *encoder)
1348 {
1349         FLAC__ASSERT(0 != encoder);
1350         if(encoder->protected_->verify)
1351                 return FLAC__stream_decoder_get_state(encoder->private_->verify.decoder);
1352         else
1353                 return FLAC__STREAM_DECODER_UNINITIALIZED;
1354 }
1355
1356 FLAC_API const char *FLAC__stream_encoder_get_resolved_state_string(const FLAC__StreamEncoder *encoder)
1357 {
1358         if(encoder->protected_->state != FLAC__STREAM_ENCODER_VERIFY_DECODER_ERROR)
1359                 return FLAC__StreamEncoderStateString[encoder->protected_->state];
1360         else
1361                 return FLAC__stream_decoder_get_resolved_state_string(encoder->private_->verify.decoder);
1362 }
1363
1364 FLAC_API void FLAC__stream_encoder_get_verify_decoder_error_stats(const FLAC__StreamEncoder *encoder, FLAC__uint64 *absolute_sample, unsigned *frame_number, unsigned *channel, unsigned *sample, FLAC__int32 *expected, FLAC__int32 *got)
1365 {
1366         FLAC__ASSERT(0 != encoder);
1367         if(0 != absolute_sample)
1368                 *absolute_sample = encoder->private_->verify.error_stats.absolute_sample;
1369         if(0 != frame_number)
1370                 *frame_number = encoder->private_->verify.error_stats.frame_number;
1371         if(0 != channel)
1372                 *channel = encoder->private_->verify.error_stats.channel;
1373         if(0 != sample)
1374                 *sample = encoder->private_->verify.error_stats.sample;
1375         if(0 != expected)
1376                 *expected = encoder->private_->verify.error_stats.expected;
1377         if(0 != got)
1378                 *got = encoder->private_->verify.error_stats.got;
1379 }
1380
1381 FLAC_API FLAC__bool FLAC__stream_encoder_get_verify(const FLAC__StreamEncoder *encoder)
1382 {
1383         FLAC__ASSERT(0 != encoder);
1384         return encoder->protected_->verify;
1385 }
1386
1387 FLAC_API FLAC__bool FLAC__stream_encoder_get_streamable_subset(const FLAC__StreamEncoder *encoder)
1388 {
1389         FLAC__ASSERT(0 != encoder);
1390         return encoder->protected_->streamable_subset;
1391 }
1392
1393 FLAC_API FLAC__bool FLAC__stream_encoder_get_do_mid_side_stereo(const FLAC__StreamEncoder *encoder)
1394 {
1395         FLAC__ASSERT(0 != encoder);
1396         return encoder->protected_->do_mid_side_stereo;
1397 }
1398
1399 FLAC_API FLAC__bool FLAC__stream_encoder_get_loose_mid_side_stereo(const FLAC__StreamEncoder *encoder)
1400 {
1401         FLAC__ASSERT(0 != encoder);
1402         return encoder->protected_->loose_mid_side_stereo;
1403 }
1404
1405 FLAC_API unsigned FLAC__stream_encoder_get_channels(const FLAC__StreamEncoder *encoder)
1406 {
1407         FLAC__ASSERT(0 != encoder);
1408         return encoder->protected_->channels;
1409 }
1410
1411 FLAC_API unsigned FLAC__stream_encoder_get_bits_per_sample(const FLAC__StreamEncoder *encoder)
1412 {
1413         FLAC__ASSERT(0 != encoder);
1414         return encoder->protected_->bits_per_sample;
1415 }
1416
1417 FLAC_API unsigned FLAC__stream_encoder_get_sample_rate(const FLAC__StreamEncoder *encoder)
1418 {
1419         FLAC__ASSERT(0 != encoder);
1420         return encoder->protected_->sample_rate;
1421 }
1422
1423 FLAC_API unsigned FLAC__stream_encoder_get_blocksize(const FLAC__StreamEncoder *encoder)
1424 {
1425         FLAC__ASSERT(0 != encoder);
1426         return encoder->protected_->blocksize;
1427 }
1428
1429 FLAC_API unsigned FLAC__stream_encoder_get_max_lpc_order(const FLAC__StreamEncoder *encoder)
1430 {
1431         FLAC__ASSERT(0 != encoder);
1432         return encoder->protected_->max_lpc_order;
1433 }
1434
1435 FLAC_API unsigned FLAC__stream_encoder_get_qlp_coeff_precision(const FLAC__StreamEncoder *encoder)
1436 {
1437         FLAC__ASSERT(0 != encoder);
1438         return encoder->protected_->qlp_coeff_precision;
1439 }
1440
1441 FLAC_API FLAC__bool FLAC__stream_encoder_get_do_qlp_coeff_prec_search(const FLAC__StreamEncoder *encoder)
1442 {
1443         FLAC__ASSERT(0 != encoder);
1444         return encoder->protected_->do_qlp_coeff_prec_search;
1445 }
1446
1447 FLAC_API FLAC__bool FLAC__stream_encoder_get_do_escape_coding(const FLAC__StreamEncoder *encoder)
1448 {
1449         FLAC__ASSERT(0 != encoder);
1450         return encoder->protected_->do_escape_coding;
1451 }
1452
1453 FLAC_API FLAC__bool FLAC__stream_encoder_get_do_exhaustive_model_search(const FLAC__StreamEncoder *encoder)
1454 {
1455         FLAC__ASSERT(0 != encoder);
1456         return encoder->protected_->do_exhaustive_model_search;
1457 }
1458
1459 FLAC_API unsigned FLAC__stream_encoder_get_min_residual_partition_order(const FLAC__StreamEncoder *encoder)
1460 {
1461         FLAC__ASSERT(0 != encoder);
1462         return encoder->protected_->min_residual_partition_order;
1463 }
1464
1465 FLAC_API unsigned FLAC__stream_encoder_get_max_residual_partition_order(const FLAC__StreamEncoder *encoder)
1466 {
1467         FLAC__ASSERT(0 != encoder);
1468         return encoder->protected_->max_residual_partition_order;
1469 }
1470
1471 FLAC_API unsigned FLAC__stream_encoder_get_rice_parameter_search_dist(const FLAC__StreamEncoder *encoder)
1472 {
1473         FLAC__ASSERT(0 != encoder);
1474         return encoder->protected_->rice_parameter_search_dist;
1475 }
1476
1477 FLAC_API FLAC__uint64 FLAC__stream_encoder_get_total_samples_estimate(const FLAC__StreamEncoder *encoder)
1478 {
1479         FLAC__ASSERT(0 != encoder);
1480         return encoder->protected_->total_samples_estimate;
1481 }
1482
1483 FLAC_API FLAC__bool FLAC__stream_encoder_process(FLAC__StreamEncoder *encoder, const FLAC__int32 * const buffer[], unsigned samples)
1484 {
1485         unsigned i, j, channel;
1486         FLAC__int32 x, mid, side;
1487         const unsigned channels = encoder->protected_->channels, blocksize = encoder->protected_->blocksize;
1488
1489         FLAC__ASSERT(0 != encoder);
1490         FLAC__ASSERT(encoder->protected_->state == FLAC__STREAM_ENCODER_OK);
1491
1492         j = 0;
1493         /*
1494          * we have several flavors of the same basic loop, optimized for
1495          * different conditions:
1496          */
1497         if(encoder->protected_->max_lpc_order > 0) {
1498                 if(encoder->protected_->do_mid_side_stereo && channels == 2) {
1499                         /*
1500                          * stereo coding: unroll channel loop
1501                          * with LPC: calculate floating point version of signal
1502                          */
1503                         do {
1504                                 if(encoder->protected_->verify)
1505                                         append_to_verify_fifo_(&encoder->private_->verify.input_fifo, buffer, j, channels, min(blocksize-encoder->private_->current_sample_number, samples-j));
1506
1507                                 for(i = encoder->private_->current_sample_number; i < blocksize && j < samples; i++, j++) {
1508                                         x = mid = side = buffer[0][j];
1509                                         encoder->private_->integer_signal[0][i] = x;
1510 #ifndef FLAC__INTEGER_ONLY_LIBRARY
1511                                         encoder->private_->real_signal[0][i] = (FLAC__real)x;
1512 #endif
1513                                         x = buffer[1][j];
1514                                         encoder->private_->integer_signal[1][i] = x;
1515 #ifndef FLAC__INTEGER_ONLY_LIBRARY
1516                                         encoder->private_->real_signal[1][i] = (FLAC__real)x;
1517 #endif
1518                                         mid += x;
1519                                         side -= x;
1520                                         mid >>= 1; /* NOTE: not the same as 'mid = (buffer[0][j] + buffer[1][j]) / 2' ! */
1521                                         encoder->private_->integer_signal_mid_side[1][i] = side;
1522                                         encoder->private_->integer_signal_mid_side[0][i] = mid;
1523 #ifndef FLAC__INTEGER_ONLY_LIBRARY
1524                                         encoder->private_->real_signal_mid_side[1][i] = (FLAC__real)side;
1525                                         encoder->private_->real_signal_mid_side[0][i] = (FLAC__real)mid;
1526 #endif
1527                                         encoder->private_->current_sample_number++;
1528                                 }
1529                                 if(i == blocksize) {
1530                                         if(!process_frame_(encoder, false)) /* false => not last frame */
1531                                                 return false;
1532                                 }
1533                         } while(j < samples);
1534                 }
1535                 else {
1536                         /*
1537                          * independent channel coding: buffer each channel in inner loop
1538                          * with LPC: calculate floating point version of signal
1539                          */
1540                         do {
1541                                 if(encoder->protected_->verify)
1542                                         append_to_verify_fifo_(&encoder->private_->verify.input_fifo, buffer, j, channels, min(blocksize-encoder->private_->current_sample_number, samples-j));
1543
1544                                 for(i = encoder->private_->current_sample_number; i < blocksize && j < samples; i++, j++) {
1545                                         for(channel = 0; channel < channels; channel++) {
1546                                                 x = buffer[channel][j];
1547                                                 encoder->private_->integer_signal[channel][i] = x;
1548 #ifndef FLAC__INTEGER_ONLY_LIBRARY
1549                                                 encoder->private_->real_signal[channel][i] = (FLAC__real)x;
1550 #endif
1551                                         }
1552                                         encoder->private_->current_sample_number++;
1553                                 }
1554                                 if(i == blocksize) {
1555                                         if(!process_frame_(encoder, false)) /* false => not last frame */
1556                                                 return false;
1557                                 }
1558                         } while(j < samples);
1559                 }
1560         }
1561         else {
1562                 if(encoder->protected_->do_mid_side_stereo && channels == 2) {
1563                         /*
1564                          * stereo coding: unroll channel loop
1565                          * without LPC: no need to calculate floating point version of signal
1566                          */
1567                         do {
1568                                 if(encoder->protected_->verify)
1569                                         append_to_verify_fifo_(&encoder->private_->verify.input_fifo, buffer, j, channels, min(blocksize-encoder->private_->current_sample_number, samples-j));
1570
1571                                 for(i = encoder->private_->current_sample_number; i < blocksize && j < samples; i++, j++) {
1572                                         encoder->private_->integer_signal[0][i] = mid = side = buffer[0][j];
1573                                         x = buffer[1][j];
1574                                         encoder->private_->integer_signal[1][i] = x;
1575                                         mid += x;
1576                                         side -= x;
1577                                         mid >>= 1; /* NOTE: not the same as 'mid = (buffer[0][j] + buffer[1][j]) / 2' ! */
1578                                         encoder->private_->integer_signal_mid_side[1][i] = side;
1579                                         encoder->private_->integer_signal_mid_side[0][i] = mid;
1580                                         encoder->private_->current_sample_number++;
1581                                 }
1582                                 if(i == blocksize) {
1583                                         if(!process_frame_(encoder, false)) /* false => not last frame */
1584                                                 return false;
1585                                 }
1586                         } while(j < samples);
1587                 }
1588                 else {
1589                         /*
1590                          * independent channel coding: buffer each channel in inner loop
1591                          * without LPC: no need to calculate floating point version of signal
1592                          */
1593                         do {
1594                                 if(encoder->protected_->verify)
1595                                         append_to_verify_fifo_(&encoder->private_->verify.input_fifo, buffer, j, channels, min(blocksize-encoder->private_->current_sample_number, samples-j));
1596
1597                                 for(i = encoder->private_->current_sample_number; i < blocksize && j < samples; i++, j++) {
1598                                         for(channel = 0; channel < channels; channel++)
1599                                                 encoder->private_->integer_signal[channel][i] = buffer[channel][j];
1600                                         encoder->private_->current_sample_number++;
1601                                 }
1602                                 if(i == blocksize) {
1603                                         if(!process_frame_(encoder, false)) /* false => not last frame */
1604                                                 return false;
1605                                 }
1606                         } while(j < samples);
1607                 }
1608         }
1609
1610         return true;
1611 }
1612
1613 FLAC_API FLAC__bool FLAC__stream_encoder_process_interleaved(FLAC__StreamEncoder *encoder, const FLAC__int32 buffer[], unsigned samples)
1614 {
1615         unsigned i, j, k, channel;
1616         FLAC__int32 x, mid, side;
1617         const unsigned channels = encoder->protected_->channels, blocksize = encoder->protected_->blocksize;
1618
1619         FLAC__ASSERT(0 != encoder);
1620         FLAC__ASSERT(encoder->protected_->state == FLAC__STREAM_ENCODER_OK);
1621
1622         j = k = 0;
1623         /*
1624          * we have several flavors of the same basic loop, optimized for
1625          * different conditions:
1626          */
1627         if(encoder->protected_->max_lpc_order > 0) {
1628                 if(encoder->protected_->do_mid_side_stereo && channels == 2) {
1629                         /*
1630                          * stereo coding: unroll channel loop
1631                          * with LPC: calculate floating point version of signal
1632                          */
1633                         do {
1634                                 if(encoder->protected_->verify)
1635                                         append_to_verify_fifo_interleaved_(&encoder->private_->verify.input_fifo, buffer, j, channels, min(blocksize-encoder->private_->current_sample_number, samples-j));
1636
1637                                 for(i = encoder->private_->current_sample_number; i < blocksize && j < samples; i++, j++) {
1638                                         x = mid = side = buffer[k++];
1639                                         encoder->private_->integer_signal[0][i] = x;
1640 #ifndef FLAC__INTEGER_ONLY_LIBRARY
1641                                         encoder->private_->real_signal[0][i] = (FLAC__real)x;
1642 #endif
1643                                         x = buffer[k++];
1644                                         encoder->private_->integer_signal[1][i] = x;
1645 #ifndef FLAC__INTEGER_ONLY_LIBRARY
1646                                         encoder->private_->real_signal[1][i] = (FLAC__real)x;
1647 #endif
1648                                         mid += x;
1649                                         side -= x;
1650                                         mid >>= 1; /* NOTE: not the same as 'mid = (left + right) / 2' ! */
1651                                         encoder->private_->integer_signal_mid_side[1][i] = side;
1652                                         encoder->private_->integer_signal_mid_side[0][i] = mid;
1653 #ifndef FLAC__INTEGER_ONLY_LIBRARY
1654                                         encoder->private_->real_signal_mid_side[1][i] = (FLAC__real)side;
1655                                         encoder->private_->real_signal_mid_side[0][i] = (FLAC__real)mid;
1656 #endif
1657                                         encoder->private_->current_sample_number++;
1658                                 }
1659                                 if(i == blocksize) {
1660                                         if(!process_frame_(encoder, false)) /* false => not last frame */
1661                                                 return false;
1662                                 }
1663                         } while(j < samples);
1664                 }
1665                 else {
1666                         /*
1667                          * independent channel coding: buffer each channel in inner loop
1668                          * with LPC: calculate floating point version of signal
1669                          */
1670                         do {
1671                                 if(encoder->protected_->verify)
1672                                         append_to_verify_fifo_interleaved_(&encoder->private_->verify.input_fifo, buffer, j, channels, min(blocksize-encoder->private_->current_sample_number, samples-j));
1673
1674                                 for(i = encoder->private_->current_sample_number; i < blocksize && j < samples; i++, j++) {
1675                                         for(channel = 0; channel < channels; channel++) {
1676                                                 x = buffer[k++];
1677                                                 encoder->private_->integer_signal[channel][i] = x;
1678 #ifndef FLAC__INTEGER_ONLY_LIBRARY
1679                                                 encoder->private_->real_signal[channel][i] = (FLAC__real)x;
1680 #endif
1681                                         }
1682                                         encoder->private_->current_sample_number++;
1683                                 }
1684                                 if(i == blocksize) {
1685                                         if(!process_frame_(encoder, false)) /* false => not last frame */
1686                                                 return false;
1687                                 }
1688                         } while(j < samples);
1689                 }
1690         }
1691         else {
1692                 if(encoder->protected_->do_mid_side_stereo && channels == 2) {
1693                         /*
1694                          * stereo coding: unroll channel loop
1695                          * without LPC: no need to calculate floating point version of signal
1696                          */
1697                         do {
1698                                 if(encoder->protected_->verify)
1699                                         append_to_verify_fifo_interleaved_(&encoder->private_->verify.input_fifo, buffer, j, channels, min(blocksize-encoder->private_->current_sample_number, samples-j));
1700
1701                                 for(i = encoder->private_->current_sample_number; i < blocksize && j < samples; i++, j++) {
1702                                         encoder->private_->integer_signal[0][i] = mid = side = buffer[k++];
1703                                         x = buffer[k++];
1704                                         encoder->private_->integer_signal[1][i] = x;
1705                                         mid += x;
1706                                         side -= x;
1707                                         mid >>= 1; /* NOTE: not the same as 'mid = (left + right) / 2' ! */
1708                                         encoder->private_->integer_signal_mid_side[1][i] = side;
1709                                         encoder->private_->integer_signal_mid_side[0][i] = mid;
1710                                         encoder->private_->current_sample_number++;
1711                                 }
1712                                 if(i == blocksize) {
1713                                         if(!process_frame_(encoder, false)) /* false => not last frame */
1714                                                 return false;
1715                                 }
1716                         } while(j < samples);
1717                 }
1718                 else {
1719                         /*
1720                          * independent channel coding: buffer each channel in inner loop
1721                          * without LPC: no need to calculate floating point version of signal
1722                          */
1723                         do {
1724                                 if(encoder->protected_->verify)
1725                                         append_to_verify_fifo_interleaved_(&encoder->private_->verify.input_fifo, buffer, j, channels, min(blocksize-encoder->private_->current_sample_number, samples-j));
1726
1727                                 for(i = encoder->private_->current_sample_number; i < blocksize && j < samples; i++, j++) {
1728                                         for(channel = 0; channel < channels; channel++)
1729                                                 encoder->private_->integer_signal[channel][i] = buffer[k++];
1730                                         encoder->private_->current_sample_number++;
1731                                 }
1732                                 if(i == blocksize) {
1733                                         if(!process_frame_(encoder, false)) /* false => not last frame */
1734                                                 return false;
1735                                 }
1736                         } while(j < samples);
1737                 }
1738         }
1739
1740         return true;
1741 }
1742
1743 /***********************************************************************
1744  *
1745  * Private class methods
1746  *
1747  ***********************************************************************/
1748
1749 void set_defaults_(FLAC__StreamEncoder *encoder)
1750 {
1751         FLAC__ASSERT(0 != encoder);
1752
1753 #ifdef FLAC__MANDATORY_VERIFY_WHILE_ENCODING
1754         encoder->protected_->verify = true;
1755 #else
1756         encoder->protected_->verify = false;
1757 #endif
1758         encoder->protected_->streamable_subset = true;
1759         encoder->protected_->do_mid_side_stereo = false;
1760         encoder->protected_->loose_mid_side_stereo = false;
1761         encoder->protected_->channels = 2;
1762         encoder->protected_->bits_per_sample = 16;
1763         encoder->protected_->sample_rate = 44100;
1764         encoder->protected_->blocksize = 1152;
1765 #ifndef FLAC__INTEGER_ONLY_LIBRARY
1766         encoder->protected_->num_apodizations = 1;
1767         encoder->protected_->apodizations[0].type = FLAC__APODIZATION_HANN;
1768 #endif
1769         encoder->protected_->max_lpc_order = 0;
1770         encoder->protected_->qlp_coeff_precision = 0;
1771         encoder->protected_->do_qlp_coeff_prec_search = false;
1772         encoder->protected_->do_exhaustive_model_search = false;
1773         encoder->protected_->do_escape_coding = false;
1774         encoder->protected_->min_residual_partition_order = 0;
1775         encoder->protected_->max_residual_partition_order = 0;
1776         encoder->protected_->rice_parameter_search_dist = 0;
1777         encoder->protected_->total_samples_estimate = 0;
1778         encoder->protected_->metadata = 0;
1779         encoder->protected_->num_metadata_blocks = 0;
1780
1781         encoder->private_->disable_constant_subframes = false;
1782         encoder->private_->disable_fixed_subframes = false;
1783         encoder->private_->disable_verbatim_subframes = false;
1784         encoder->private_->write_callback = 0;
1785         encoder->private_->metadata_callback = 0;
1786         encoder->private_->client_data = 0;
1787 }
1788
1789 void free_(FLAC__StreamEncoder *encoder)
1790 {
1791         unsigned i, channel;
1792
1793         FLAC__ASSERT(0 != encoder);
1794         for(i = 0; i < encoder->protected_->channels; i++) {
1795                 if(0 != encoder->private_->integer_signal_unaligned[i]) {
1796                         free(encoder->private_->integer_signal_unaligned[i]);
1797                         encoder->private_->integer_signal_unaligned[i] = 0;
1798                 }
1799 #ifndef FLAC__INTEGER_ONLY_LIBRARY
1800                 if(0 != encoder->private_->real_signal_unaligned[i]) {
1801                         free(encoder->private_->real_signal_unaligned[i]);
1802                         encoder->private_->real_signal_unaligned[i] = 0;
1803                 }
1804 #endif
1805         }
1806         for(i = 0; i < 2; i++) {
1807                 if(0 != encoder->private_->integer_signal_mid_side_unaligned[i]) {
1808                         free(encoder->private_->integer_signal_mid_side_unaligned[i]);
1809                         encoder->private_->integer_signal_mid_side_unaligned[i] = 0;
1810                 }
1811 #ifndef FLAC__INTEGER_ONLY_LIBRARY
1812                 if(0 != encoder->private_->real_signal_mid_side_unaligned[i]) {
1813                         free(encoder->private_->real_signal_mid_side_unaligned[i]);
1814                         encoder->private_->real_signal_mid_side_unaligned[i] = 0;
1815                 }
1816 #endif
1817         }
1818 #ifndef FLAC__INTEGER_ONLY_LIBRARY
1819         for(i = 0; i < encoder->protected_->num_apodizations; i++) {
1820                 if(0 != encoder->private_->window_unaligned[i]) {
1821                         free(encoder->private_->window_unaligned[i]);
1822                         encoder->private_->window_unaligned[i] = 0;
1823                 }
1824         }
1825         if(0 != encoder->private_->windowed_signal_unaligned) {
1826                 free(encoder->private_->windowed_signal_unaligned);
1827                 encoder->private_->windowed_signal_unaligned = 0;
1828         }
1829 #endif
1830         for(channel = 0; channel < encoder->protected_->channels; channel++) {
1831                 for(i = 0; i < 2; i++) {
1832                         if(0 != encoder->private_->residual_workspace_unaligned[channel][i]) {
1833                                 free(encoder->private_->residual_workspace_unaligned[channel][i]);
1834                                 encoder->private_->residual_workspace_unaligned[channel][i] = 0;
1835                         }
1836                 }
1837         }
1838         for(channel = 0; channel < 2; channel++) {
1839                 for(i = 0; i < 2; i++) {
1840                         if(0 != encoder->private_->residual_workspace_mid_side_unaligned[channel][i]) {
1841                                 free(encoder->private_->residual_workspace_mid_side_unaligned[channel][i]);
1842                                 encoder->private_->residual_workspace_mid_side_unaligned[channel][i] = 0;
1843                         }
1844                 }
1845         }
1846         if(0 != encoder->private_->abs_residual_unaligned) {
1847                 free(encoder->private_->abs_residual_unaligned);
1848                 encoder->private_->abs_residual_unaligned = 0;
1849         }
1850         if(0 != encoder->private_->abs_residual_partition_sums_unaligned) {
1851                 free(encoder->private_->abs_residual_partition_sums_unaligned);
1852                 encoder->private_->abs_residual_partition_sums_unaligned = 0;
1853         }
1854         if(0 != encoder->private_->raw_bits_per_partition_unaligned) {
1855                 free(encoder->private_->raw_bits_per_partition_unaligned);
1856                 encoder->private_->raw_bits_per_partition_unaligned = 0;
1857         }
1858         if(encoder->protected_->verify) {
1859                 for(i = 0; i < encoder->protected_->channels; i++) {
1860                         if(0 != encoder->private_->verify.input_fifo.data[i]) {
1861                                 free(encoder->private_->verify.input_fifo.data[i]);
1862                                 encoder->private_->verify.input_fifo.data[i] = 0;
1863                         }
1864                 }
1865         }
1866         FLAC__bitbuffer_free(encoder->private_->frame);
1867 }
1868
1869 FLAC__bool resize_buffers_(FLAC__StreamEncoder *encoder, unsigned new_size)
1870 {
1871         FLAC__bool ok;
1872         unsigned i, channel;
1873
1874         FLAC__ASSERT(new_size > 0);
1875         FLAC__ASSERT(encoder->protected_->state == FLAC__STREAM_ENCODER_OK);
1876         FLAC__ASSERT(encoder->private_->current_sample_number == 0);
1877
1878         /* To avoid excessive malloc'ing, we only grow the buffer; no shrinking. */
1879         if(new_size <= encoder->private_->input_capacity)
1880                 return true;
1881
1882         ok = true;
1883
1884         /* WATCHOUT: FLAC__lpc_compute_residual_from_qlp_coefficients_asm_ia32_mmx()
1885          * requires that the input arrays (in our case the integer signals)
1886          * have a buffer of up to 3 zeroes in front (at negative indices) for
1887          * alignment purposes; we use 4 to keep the data well-aligned.
1888          */
1889
1890         for(i = 0; ok && i < encoder->protected_->channels; i++) {
1891                 ok = ok && FLAC__memory_alloc_aligned_int32_array(new_size+4, &encoder->private_->integer_signal_unaligned[i], &encoder->private_->integer_signal[i]);
1892 #ifndef FLAC__INTEGER_ONLY_LIBRARY
1893                 if(encoder->protected_->max_lpc_order > 0)
1894                         ok = ok && FLAC__memory_alloc_aligned_real_array(new_size, &encoder->private_->real_signal_unaligned[i], &encoder->private_->real_signal[i]);
1895 #endif
1896                 memset(encoder->private_->integer_signal[i], 0, sizeof(FLAC__int32)*4);
1897                 encoder->private_->integer_signal[i] += 4;
1898         }
1899         for(i = 0; ok && i < 2; i++) {
1900                 ok = ok && FLAC__memory_alloc_aligned_int32_array(new_size+4, &encoder->private_->integer_signal_mid_side_unaligned[i], &encoder->private_->integer_signal_mid_side[i]);
1901 #ifndef FLAC__INTEGER_ONLY_LIBRARY
1902                 if(encoder->protected_->max_lpc_order > 0)
1903                         ok = ok && FLAC__memory_alloc_aligned_real_array(new_size, &encoder->private_->real_signal_mid_side_unaligned[i], &encoder->private_->real_signal_mid_side[i]);
1904 #endif
1905                 memset(encoder->private_->integer_signal_mid_side[i], 0, sizeof(FLAC__int32)*4);
1906                 encoder->private_->integer_signal_mid_side[i] += 4;
1907         }
1908 #ifndef FLAC__INTEGER_ONLY_LIBRARY
1909         if(ok && encoder->protected_->max_lpc_order > 0) {
1910                 for(i = 0; ok && i < encoder->protected_->num_apodizations; i++)
1911                         ok = ok && FLAC__memory_alloc_aligned_real_array(new_size, &encoder->private_->window_unaligned[i], &encoder->private_->window[i]);
1912                 ok = ok && FLAC__memory_alloc_aligned_real_array(new_size, &encoder->private_->windowed_signal_unaligned, &encoder->private_->windowed_signal);
1913         }
1914 #endif
1915         for(channel = 0; ok && channel < encoder->protected_->channels; channel++) {
1916                 for(i = 0; ok && i < 2; i++) {
1917                         ok = ok && FLAC__memory_alloc_aligned_int32_array(new_size, &encoder->private_->residual_workspace_unaligned[channel][i], &encoder->private_->residual_workspace[channel][i]);
1918                 }
1919         }
1920         for(channel = 0; ok && channel < 2; channel++) {
1921                 for(i = 0; ok && i < 2; i++) {
1922                         ok = ok && FLAC__memory_alloc_aligned_int32_array(new_size, &encoder->private_->residual_workspace_mid_side_unaligned[channel][i], &encoder->private_->residual_workspace_mid_side[channel][i]);
1923                 }
1924         }
1925         ok = ok && FLAC__memory_alloc_aligned_uint32_array(new_size, &encoder->private_->abs_residual_unaligned, &encoder->private_->abs_residual);
1926         if(encoder->private_->precompute_partition_sums || encoder->protected_->do_escape_coding) /* we require precompute_partition_sums if do_escape_coding because of their intertwined nature */
1927                 ok = ok && FLAC__memory_alloc_aligned_uint64_array(new_size * 2, &encoder->private_->abs_residual_partition_sums_unaligned, &encoder->private_->abs_residual_partition_sums);
1928         if(encoder->protected_->do_escape_coding)
1929                 ok = ok && FLAC__memory_alloc_aligned_unsigned_array(new_size * 2, &encoder->private_->raw_bits_per_partition_unaligned, &encoder->private_->raw_bits_per_partition);
1930
1931         if(ok)
1932                 encoder->private_->input_capacity = new_size;
1933         else
1934                 encoder->protected_->state = FLAC__STREAM_ENCODER_MEMORY_ALLOCATION_ERROR;
1935
1936 #ifndef FLAC__INTEGER_ONLY_LIBRARY
1937         if(ok && encoder->protected_->max_lpc_order > 0) {
1938                 for(i = 0; ok && i < encoder->protected_->num_apodizations; i++) {
1939                         switch(encoder->protected_->apodizations[i].type) {
1940                                 case FLAC__APODIZATION_BARTLETT:
1941                                         FLAC__window_bartlett(encoder->private_->window[i], new_size);
1942                                         break;
1943                                 case FLAC__APODIZATION_BARTLETT_HANN:
1944                                         FLAC__window_bartlett_hann(encoder->private_->window[i], new_size);
1945                                         break;
1946                                 case FLAC__APODIZATION_BLACKMAN:
1947                                         FLAC__window_blackman(encoder->private_->window[i], new_size);
1948                                         break;
1949                                 case FLAC__APODIZATION_BLACKMAN_HARRIS_4TERM_92DB_SIDELOBE:
1950                                         FLAC__window_blackman_harris_4term_92db_sidelobe(encoder->private_->window[i], new_size);
1951                                         break;
1952                                 case FLAC__APODIZATION_CONNES:
1953                                         FLAC__window_connes(encoder->private_->window[i], new_size);
1954                                         break;
1955                                 case FLAC__APODIZATION_FLATTOP:
1956                                         FLAC__window_flattop(encoder->private_->window[i], new_size);
1957                                         break;
1958                                 case FLAC__APODIZATION_GAUSS:
1959                                         FLAC__window_gauss(encoder->private_->window[i], new_size, encoder->protected_->apodizations[i].parameters.gauss.stddev);
1960                                         break;
1961                                 case FLAC__APODIZATION_HAMMING:
1962                                         FLAC__window_hamming(encoder->private_->window[i], new_size);
1963                                         break;
1964                                 case FLAC__APODIZATION_HANN:
1965                                         FLAC__window_hann(encoder->private_->window[i], new_size);
1966                                         break;
1967                                 case FLAC__APODIZATION_KAISER_BESSEL:
1968                                         FLAC__window_kaiser_bessel(encoder->private_->window[i], new_size);
1969                                         break;
1970                                 case FLAC__APODIZATION_NUTTALL:
1971                                         FLAC__window_nuttall(encoder->private_->window[i], new_size);
1972                                         break;
1973                                 case FLAC__APODIZATION_RECTANGLE:
1974                                         FLAC__window_rectangle(encoder->private_->window[i], new_size);
1975                                         break;
1976                                 case FLAC__APODIZATION_TRIANGLE:
1977                                         FLAC__window_triangle(encoder->private_->window[i], new_size);
1978                                         break;
1979                                 case FLAC__APODIZATION_TUKEY:
1980                                         FLAC__window_tukey(encoder->private_->window[i], new_size, encoder->protected_->apodizations[i].parameters.tukey.p);
1981                                         break;
1982                                 case FLAC__APODIZATION_WELCH:
1983                                         FLAC__window_welch(encoder->private_->window[i], new_size);
1984                                         break;
1985                                 default:
1986                                         FLAC__ASSERT(0);
1987                                         /* double protection */
1988                                         FLAC__window_rectangle(encoder->private_->window[i], new_size);
1989                                         break;
1990                         }
1991                 }
1992         }
1993 #endif
1994
1995         return ok;
1996 }
1997
1998 FLAC__bool write_bitbuffer_(FLAC__StreamEncoder *encoder, unsigned samples)
1999 {
2000         const FLAC__byte *buffer;
2001         unsigned bytes;
2002
2003         FLAC__ASSERT(FLAC__bitbuffer_is_byte_aligned(encoder->private_->frame));
2004
2005         FLAC__bitbuffer_get_buffer(encoder->private_->frame, &buffer, &bytes);
2006
2007         if(encoder->protected_->verify) {
2008                 encoder->private_->verify.output.data = buffer;
2009                 encoder->private_->verify.output.bytes = bytes;
2010                 if(encoder->private_->verify.state_hint == ENCODER_IN_MAGIC) {
2011                         encoder->private_->verify.needs_magic_hack = true;
2012                 }
2013                 else {
2014                         if(!FLAC__stream_decoder_process_single(encoder->private_->verify.decoder)) {
2015                                 FLAC__bitbuffer_release_buffer(encoder->private_->frame);
2016                                 if(encoder->protected_->state != FLAC__STREAM_ENCODER_VERIFY_MISMATCH_IN_AUDIO_DATA)
2017                                         encoder->protected_->state = FLAC__STREAM_ENCODER_VERIFY_DECODER_ERROR;
2018                                 return false;
2019                         }
2020                 }
2021         }
2022
2023         if(encoder->private_->write_callback(encoder, buffer, bytes, samples, encoder->private_->current_frame_number, encoder->private_->client_data) != FLAC__STREAM_ENCODER_WRITE_STATUS_OK) {
2024                 FLAC__bitbuffer_release_buffer(encoder->private_->frame);
2025                 encoder->protected_->state = FLAC__STREAM_ENCODER_FATAL_ERROR_WHILE_WRITING;
2026                 return false;
2027         }
2028
2029         FLAC__bitbuffer_release_buffer(encoder->private_->frame);
2030
2031         if(samples > 0) {
2032                 encoder->private_->metadata.data.stream_info.min_framesize = min(bytes, encoder->private_->metadata.data.stream_info.min_framesize);
2033                 encoder->private_->metadata.data.stream_info.max_framesize = max(bytes, encoder->private_->metadata.data.stream_info.max_framesize);
2034         }
2035
2036         return true;
2037 }
2038
2039 FLAC__bool process_frame_(FLAC__StreamEncoder *encoder, FLAC__bool is_last_frame)
2040 {
2041         FLAC__ASSERT(encoder->protected_->state == FLAC__STREAM_ENCODER_OK);
2042
2043         /*
2044          * Accumulate raw signal to the MD5 signature
2045          */
2046         if(!FLAC__MD5Accumulate(&encoder->private_->md5context, (const FLAC__int32 * const *)encoder->private_->integer_signal, encoder->protected_->channels, encoder->protected_->blocksize, (encoder->protected_->bits_per_sample+7) / 8)) {
2047                 encoder->protected_->state = FLAC__STREAM_ENCODER_MEMORY_ALLOCATION_ERROR;
2048                 return false;
2049         }
2050
2051         /*
2052          * Process the frame header and subframes into the frame bitbuffer
2053          */
2054         if(!process_subframes_(encoder, is_last_frame)) {
2055                 /* the above function sets the state for us in case of an error */
2056                 return false;
2057         }
2058
2059         /*
2060          * Zero-pad the frame to a byte_boundary
2061          */
2062         if(!FLAC__bitbuffer_zero_pad_to_byte_boundary(encoder->private_->frame)) {
2063                 encoder->protected_->state = FLAC__STREAM_ENCODER_MEMORY_ALLOCATION_ERROR;
2064                 return false;
2065         }
2066
2067         /*
2068          * CRC-16 the whole thing
2069          */
2070         FLAC__ASSERT(FLAC__bitbuffer_is_byte_aligned(encoder->private_->frame));
2071         FLAC__bitbuffer_write_raw_uint32(encoder->private_->frame, FLAC__bitbuffer_get_write_crc16(encoder->private_->frame), FLAC__FRAME_FOOTER_CRC_LEN);
2072
2073         /*
2074          * Write it
2075          */
2076         if(!write_bitbuffer_(encoder, encoder->protected_->blocksize)) {
2077                 /* the above function sets the state for us in case of an error */
2078                 return false;
2079         }
2080
2081         /*
2082          * Get ready for the next frame
2083          */
2084         encoder->private_->current_sample_number = 0;
2085         encoder->private_->current_frame_number++;
2086         encoder->private_->metadata.data.stream_info.total_samples += (FLAC__uint64)encoder->protected_->blocksize;
2087
2088         return true;
2089 }
2090
2091 FLAC__bool process_subframes_(FLAC__StreamEncoder *encoder, FLAC__bool is_last_frame)
2092 {
2093         FLAC__FrameHeader frame_header;
2094         unsigned channel, min_partition_order = encoder->protected_->min_residual_partition_order, max_partition_order;
2095         FLAC__bool do_independent, do_mid_side, precompute_partition_sums;
2096
2097         /*
2098          * Calculate the min,max Rice partition orders
2099          */
2100         if(is_last_frame) {
2101                 max_partition_order = 0;
2102         }
2103         else {
2104                 max_partition_order = FLAC__format_get_max_rice_partition_order_from_blocksize(encoder->protected_->blocksize);
2105                 max_partition_order = min(max_partition_order, encoder->protected_->max_residual_partition_order);
2106         }
2107         min_partition_order = min(min_partition_order, max_partition_order);
2108
2109         precompute_partition_sums = encoder->private_->precompute_partition_sums && ((max_partition_order > min_partition_order) || encoder->protected_->do_escape_coding);
2110
2111         /*
2112          * Setup the frame
2113          */
2114         if(!FLAC__bitbuffer_clear(encoder->private_->frame)) {
2115                 encoder->protected_->state = FLAC__STREAM_ENCODER_MEMORY_ALLOCATION_ERROR;
2116                 return false;
2117         }
2118         frame_header.blocksize = encoder->protected_->blocksize;
2119         frame_header.sample_rate = encoder->protected_->sample_rate;
2120         frame_header.channels = encoder->protected_->channels;
2121         frame_header.channel_assignment = FLAC__CHANNEL_ASSIGNMENT_INDEPENDENT; /* the default unless the encoder determines otherwise */
2122         frame_header.bits_per_sample = encoder->protected_->bits_per_sample;
2123         frame_header.number_type = FLAC__FRAME_NUMBER_TYPE_FRAME_NUMBER;
2124         frame_header.number.frame_number = encoder->private_->current_frame_number;
2125
2126         /*
2127          * Figure out what channel assignments to try
2128          */
2129         if(encoder->protected_->do_mid_side_stereo) {
2130                 if(encoder->protected_->loose_mid_side_stereo) {
2131                         if(encoder->private_->loose_mid_side_stereo_frame_count == 0) {
2132                                 do_independent = true;
2133                                 do_mid_side = true;
2134                         }
2135                         else {
2136                                 do_independent = (encoder->private_->last_channel_assignment == FLAC__CHANNEL_ASSIGNMENT_INDEPENDENT);
2137                                 do_mid_side = !do_independent;
2138                         }
2139                 }
2140                 else {
2141                         do_independent = true;
2142                         do_mid_side = true;
2143                 }
2144         }
2145         else {
2146                 do_independent = true;
2147                 do_mid_side = false;
2148         }
2149
2150         FLAC__ASSERT(do_independent || do_mid_side);
2151
2152         /*
2153          * Check for wasted bits; set effective bps for each subframe
2154          */
2155         if(do_independent) {
2156                 for(channel = 0; channel < encoder->protected_->channels; channel++) {
2157                         const unsigned w = get_wasted_bits_(encoder->private_->integer_signal[channel], encoder->protected_->blocksize);
2158                         encoder->private_->subframe_workspace[channel][0].wasted_bits = encoder->private_->subframe_workspace[channel][1].wasted_bits = w;
2159                         encoder->private_->subframe_bps[channel] = encoder->protected_->bits_per_sample - w;
2160                 }
2161         }
2162         if(do_mid_side) {
2163                 FLAC__ASSERT(encoder->protected_->channels == 2);
2164                 for(channel = 0; channel < 2; channel++) {
2165                         const unsigned w = get_wasted_bits_(encoder->private_->integer_signal_mid_side[channel], encoder->protected_->blocksize);
2166                         encoder->private_->subframe_workspace_mid_side[channel][0].wasted_bits = encoder->private_->subframe_workspace_mid_side[channel][1].wasted_bits = w;
2167                         encoder->private_->subframe_bps_mid_side[channel] = encoder->protected_->bits_per_sample - w + (channel==0? 0:1);
2168                 }
2169         }
2170
2171         /*
2172          * First do a normal encoding pass of each independent channel
2173          */
2174         if(do_independent) {
2175                 for(channel = 0; channel < encoder->protected_->channels; channel++) {
2176                         if(!
2177                                 process_subframe_(
2178                                         encoder,
2179                                         min_partition_order,
2180                                         max_partition_order,
2181                                         precompute_partition_sums,
2182                                         &frame_header,
2183                                         encoder->private_->subframe_bps[channel],
2184                                         encoder->private_->integer_signal[channel],
2185 #ifndef FLAC__INTEGER_ONLY_LIBRARY
2186                                         encoder->private_->real_signal[channel],
2187 #endif
2188                                         encoder->private_->subframe_workspace_ptr[channel],
2189                                         encoder->private_->partitioned_rice_contents_workspace_ptr[channel],
2190                                         encoder->private_->residual_workspace[channel],
2191                                         encoder->private_->best_subframe+channel,
2192                                         encoder->private_->best_subframe_bits+channel
2193 #ifdef WINDOW_DEBUG_OUTPUT
2194                                         ,channel
2195 #endif
2196                                 )
2197                         )
2198                                 return false;
2199                 }
2200         }
2201
2202         /*
2203          * Now do mid and side channels if requested
2204          */
2205         if(do_mid_side) {
2206                 FLAC__ASSERT(encoder->protected_->channels == 2);
2207
2208                 for(channel = 0; channel < 2; channel++) {
2209                         if(!
2210                                 process_subframe_(
2211                                         encoder,
2212                                         min_partition_order,
2213                                         max_partition_order,
2214                                         precompute_partition_sums,
2215                                         &frame_header,
2216                                         encoder->private_->subframe_bps_mid_side[channel],
2217                                         encoder->private_->integer_signal_mid_side[channel],
2218 #ifndef FLAC__INTEGER_ONLY_LIBRARY
2219                                         encoder->private_->real_signal_mid_side[channel],
2220 #endif
2221                                         encoder->private_->subframe_workspace_ptr_mid_side[channel],
2222                                         encoder->private_->partitioned_rice_contents_workspace_ptr_mid_side[channel],
2223                                         encoder->private_->residual_workspace_mid_side[channel],
2224                                         encoder->private_->best_subframe_mid_side+channel,
2225                                         encoder->private_->best_subframe_bits_mid_side+channel
2226 #ifdef WINDOW_DEBUG_OUTPUT
2227                                         ,channel
2228 #endif
2229                                 )
2230                         )
2231                                 return false;
2232                 }
2233         }
2234
2235         /*
2236          * Compose the frame bitbuffer
2237          */
2238         if(do_mid_side) {
2239                 unsigned left_bps = 0, right_bps = 0; /* initialized only to prevent superfluous compiler warning */
2240                 FLAC__Subframe *left_subframe = 0, *right_subframe = 0; /* initialized only to prevent superfluous compiler warning */
2241                 FLAC__ChannelAssignment channel_assignment;
2242 #ifdef WINDOW_DEBUG_OUTPUT
2243                 unsigned left_bits = 0, right_bits = 0;
2244 #endif
2245
2246                 FLAC__ASSERT(encoder->protected_->channels == 2);
2247
2248                 if(encoder->protected_->loose_mid_side_stereo && encoder->private_->loose_mid_side_stereo_frame_count > 0) {
2249                         channel_assignment = (encoder->private_->last_channel_assignment == FLAC__CHANNEL_ASSIGNMENT_INDEPENDENT? FLAC__CHANNEL_ASSIGNMENT_INDEPENDENT : FLAC__CHANNEL_ASSIGNMENT_MID_SIDE);
2250                 }
2251                 else {
2252                         unsigned bits[4]; /* WATCHOUT - indexed by FLAC__ChannelAssignment */
2253                         unsigned min_bits;
2254                         FLAC__ChannelAssignment ca;
2255
2256                         FLAC__ASSERT(do_independent && do_mid_side);
2257
2258                         /* We have to figure out which channel assignent results in the smallest frame */
2259                         bits[FLAC__CHANNEL_ASSIGNMENT_INDEPENDENT] = encoder->private_->best_subframe_bits         [0] + encoder->private_->best_subframe_bits         [1];
2260                         bits[FLAC__CHANNEL_ASSIGNMENT_LEFT_SIDE  ] = encoder->private_->best_subframe_bits         [0] + encoder->private_->best_subframe_bits_mid_side[1];
2261                         bits[FLAC__CHANNEL_ASSIGNMENT_RIGHT_SIDE ] = encoder->private_->best_subframe_bits         [1] + encoder->private_->best_subframe_bits_mid_side[1];
2262                         bits[FLAC__CHANNEL_ASSIGNMENT_MID_SIDE   ] = encoder->private_->best_subframe_bits_mid_side[0] + encoder->private_->best_subframe_bits_mid_side[1];
2263
2264                         for(channel_assignment = (FLAC__ChannelAssignment)0, min_bits = bits[0], ca = (FLAC__ChannelAssignment)1; (int)ca <= 3; ca = (FLAC__ChannelAssignment)((int)ca + 1)) {
2265                                 if(bits[ca] < min_bits) {
2266                                         min_bits = bits[ca];
2267                                         channel_assignment = ca;
2268                                 }
2269                         }
2270                 }
2271
2272                 frame_header.channel_assignment = channel_assignment;
2273
2274                 if(!FLAC__frame_add_header(&frame_header, encoder->protected_->streamable_subset, encoder->private_->frame)) {
2275                         encoder->protected_->state = FLAC__STREAM_ENCODER_FRAMING_ERROR;
2276                         return false;
2277                 }
2278
2279                 switch(channel_assignment) {
2280                         case FLAC__CHANNEL_ASSIGNMENT_INDEPENDENT:
2281                                 left_subframe  = &encoder->private_->subframe_workspace         [0][encoder->private_->best_subframe         [0]];
2282                                 right_subframe = &encoder->private_->subframe_workspace         [1][encoder->private_->best_subframe         [1]];
2283 #ifdef WINDOW_DEBUG_OUTPUT
2284                                 left_bits      = encoder->private_->best_subframe_bits          [0];
2285                                 right_bits     = encoder->private_->best_subframe_bits          [1];
2286 #endif
2287                                 break;
2288                         case FLAC__CHANNEL_ASSIGNMENT_LEFT_SIDE:
2289                                 left_subframe  = &encoder->private_->subframe_workspace         [0][encoder->private_->best_subframe         [0]];
2290                                 right_subframe = &encoder->private_->subframe_workspace_mid_side[1][encoder->private_->best_subframe_mid_side[1]];
2291 #ifdef WINDOW_DEBUG_OUTPUT
2292                                 left_bits      = encoder->private_->best_subframe_bits          [0];
2293                                 right_bits     = encoder->private_->best_subframe_bits_mid_side [1];
2294 #endif
2295                                 break;
2296                         case FLAC__CHANNEL_ASSIGNMENT_RIGHT_SIDE:
2297                                 left_subframe  = &encoder->private_->subframe_workspace_mid_side[1][encoder->private_->best_subframe_mid_side[1]];
2298                                 right_subframe = &encoder->private_->subframe_workspace         [1][encoder->private_->best_subframe         [1]];
2299 #ifdef WINDOW_DEBUG_OUTPUT
2300                                 left_bits      = encoder->private_->best_subframe_bits_mid_side [1];
2301                                 right_bits     = encoder->private_->best_subframe_bits          [1];
2302 #endif
2303                                 break;
2304                         case FLAC__CHANNEL_ASSIGNMENT_MID_SIDE:
2305                                 left_subframe  = &encoder->private_->subframe_workspace_mid_side[0][encoder->private_->best_subframe_mid_side[0]];
2306                                 right_subframe = &encoder->private_->subframe_workspace_mid_side[1][encoder->private_->best_subframe_mid_side[1]];
2307 #ifdef WINDOW_DEBUG_OUTPUT
2308                                 left_bits      = encoder->private_->best_subframe_bits_mid_side [0];
2309                                 right_bits     = encoder->private_->best_subframe_bits_mid_side [1];
2310 #endif
2311                                 break;
2312                         default:
2313                                 FLAC__ASSERT(0);
2314                 }
2315
2316                 switch(channel_assignment) {
2317                         case FLAC__CHANNEL_ASSIGNMENT_INDEPENDENT:
2318                                 left_bps  = encoder->private_->subframe_bps         [0];
2319                                 right_bps = encoder->private_->subframe_bps         [1];
2320                                 break;
2321                         case FLAC__CHANNEL_ASSIGNMENT_LEFT_SIDE:
2322                                 left_bps  = encoder->private_->subframe_bps         [0];
2323                                 right_bps = encoder->private_->subframe_bps_mid_side[1];
2324                                 break;
2325                         case FLAC__CHANNEL_ASSIGNMENT_RIGHT_SIDE:
2326                                 left_bps  = encoder->private_->subframe_bps_mid_side[1];
2327                                 right_bps = encoder->private_->subframe_bps         [1];
2328                                 break;
2329                         case FLAC__CHANNEL_ASSIGNMENT_MID_SIDE:
2330                                 left_bps  = encoder->private_->subframe_bps_mid_side[0];
2331                                 right_bps = encoder->private_->subframe_bps_mid_side[1];
2332                                 break;
2333                         default:
2334                                 FLAC__ASSERT(0);
2335                 }
2336
2337                 /* note that encoder_add_subframe_ sets the state for us in case of an error */
2338 #ifdef WINDOW_DEBUG_OUTPUT
2339                 if(!add_subframe_(encoder, &frame_header, left_bps , left_subframe , encoder->private_->frame, left_bits))
2340                         return false;
2341                 if(!add_subframe_(encoder, &frame_header, right_bps, right_subframe, encoder->private_->frame, right_bits))
2342                         return false;
2343 #else
2344                 if(!add_subframe_(encoder, &frame_header, left_bps , left_subframe , encoder->private_->frame))
2345                         return false;
2346                 if(!add_subframe_(encoder, &frame_header, right_bps, right_subframe, encoder->private_->frame))
2347                         return false;
2348 #endif
2349         }
2350         else {
2351                 if(!FLAC__frame_add_header(&frame_header, encoder->protected_->streamable_subset, encoder->private_->frame)) {
2352                         encoder->protected_->state = FLAC__STREAM_ENCODER_FRAMING_ERROR;
2353                         return false;
2354                 }
2355
2356                 for(channel = 0; channel < encoder->protected_->channels; channel++) {
2357                         if(!add_subframe_(encoder, &frame_header, encoder->private_->subframe_bps[channel], &encoder->private_->subframe_workspace[channel][encoder->private_->best_subframe[channel]], encoder->private_->frame, encoder->private_->best_subframe_bits[channel])) {
2358                                 /* the above function sets the state for us in case of an error */
2359                                 return false;
2360                         }
2361                 }
2362         }
2363
2364         if(encoder->protected_->loose_mid_side_stereo) {
2365                 encoder->private_->loose_mid_side_stereo_frame_count++;
2366                 if(encoder->private_->loose_mid_side_stereo_frame_count >= encoder->private_->loose_mid_side_stereo_frames)
2367                         encoder->private_->loose_mid_side_stereo_frame_count = 0;
2368         }
2369
2370         encoder->private_->last_channel_assignment = frame_header.channel_assignment;
2371
2372         return true;
2373 }
2374
2375 FLAC__bool process_subframe_(
2376         FLAC__StreamEncoder *encoder,
2377         unsigned min_partition_order,
2378         unsigned max_partition_order,
2379         FLAC__bool precompute_partition_sums,
2380         const FLAC__FrameHeader *frame_header,
2381         unsigned subframe_bps,
2382         const FLAC__int32 integer_signal[],
2383 #ifndef FLAC__INTEGER_ONLY_LIBRARY
2384         const FLAC__real real_signal[],
2385 #endif
2386         FLAC__Subframe *subframe[2],
2387         FLAC__EntropyCodingMethod_PartitionedRiceContents *partitioned_rice_contents[2],
2388         FLAC__int32 *residual[2],
2389         unsigned *best_subframe,
2390         unsigned *best_bits
2391 #ifdef WINDOW_DEBUG_OUTPUT
2392         ,unsigned subframe_number
2393 #endif
2394 )
2395 {
2396 #ifndef FLAC__INTEGER_ONLY_LIBRARY
2397         FLAC__float fixed_residual_bits_per_sample[FLAC__MAX_FIXED_ORDER+1];
2398 #else
2399         FLAC__fixedpoint fixed_residual_bits_per_sample[FLAC__MAX_FIXED_ORDER+1];
2400 #endif
2401 #ifndef FLAC__INTEGER_ONLY_LIBRARY
2402         FLAC__double lpc_residual_bits_per_sample;
2403         FLAC__real autoc[FLAC__MAX_LPC_ORDER+1]; /* WATCHOUT: the size is important even though encoder->protected_->max_lpc_order might be less; some asm routines need all the space */
2404         FLAC__double lpc_error[FLAC__MAX_LPC_ORDER];
2405         unsigned min_lpc_order, max_lpc_order, lpc_order;
2406         unsigned min_qlp_coeff_precision, max_qlp_coeff_precision, qlp_coeff_precision;
2407 #endif
2408         unsigned min_fixed_order, max_fixed_order, guess_fixed_order, fixed_order;
2409         unsigned rice_parameter;
2410         unsigned _candidate_bits, _best_bits;
2411         unsigned _best_subframe;
2412
2413         /* verbatim subframe is the baseline against which we measure other compressed subframes */
2414         _best_subframe = 0;
2415         if(encoder->private_->disable_verbatim_subframes && frame_header->blocksize >= FLAC__MAX_FIXED_ORDER)
2416                 _best_bits = UINT_MAX;
2417         else
2418                 _best_bits = evaluate_verbatim_subframe_(integer_signal, frame_header->blocksize, subframe_bps, subframe[_best_subframe]);
2419
2420         if(frame_header->blocksize >= FLAC__MAX_FIXED_ORDER) {
2421                 unsigned signal_is_constant = false;
2422                 guess_fixed_order = encoder->private_->local_fixed_compute_best_predictor(integer_signal+FLAC__MAX_FIXED_ORDER, frame_header->blocksize-FLAC__MAX_FIXED_ORDER, fixed_residual_bits_per_sample);
2423                 /* check for constant subframe */
2424                 if(
2425                         !encoder->private_->disable_constant_subframes &&
2426 #ifndef FLAC__INTEGER_ONLY_LIBRARY
2427                         fixed_residual_bits_per_sample[1] == 0.0
2428 #else
2429                         fixed_residual_bits_per_sample[1] == FLAC__FP_ZERO
2430 #endif
2431                 ) {
2432                         /* the above means it's possible all samples are the same value; now double-check it: */
2433                         unsigned i;
2434                         signal_is_constant = true;
2435                         for(i = 1; i < frame_header->blocksize; i++) {
2436                                 if(integer_signal[0] != integer_signal[i]) {
2437                                         signal_is_constant = false;
2438                                         break;
2439                                 }
2440                         }
2441                 }
2442                 if(signal_is_constant) {
2443                         _candidate_bits = evaluate_constant_subframe_(integer_signal[0], subframe_bps, subframe[!_best_subframe]);
2444                         if(_candidate_bits < _best_bits) {
2445                                 _best_subframe = !_best_subframe;
2446                                 _best_bits = _candidate_bits;
2447                         }
2448                 }
2449                 else {
2450                         if(!encoder->private_->disable_fixed_subframes || (encoder->protected_->max_lpc_order == 0 && _best_bits == UINT_MAX)) {
2451                                 /* encode fixed */
2452                                 if(encoder->protected_->do_exhaustive_model_search) {
2453                                         min_fixed_order = 0;
2454                                         max_fixed_order = FLAC__MAX_FIXED_ORDER;
2455                                 }
2456                                 else {
2457                                         min_fixed_order = max_fixed_order = guess_fixed_order;
2458                                 }
2459                                 for(fixed_order = min_fixed_order; fixed_order <= max_fixed_order; fixed_order++) {
2460 #ifndef FLAC__INTEGER_ONLY_LIBRARY
2461                                         if(fixed_residual_bits_per_sample[fixed_order] >= (FLAC__float)subframe_bps)
2462                                                 continue; /* don't even try */
2463                                         rice_parameter = (fixed_residual_bits_per_sample[fixed_order] > 0.0)? (unsigned)(fixed_residual_bits_per_sample[fixed_order]+0.5) : 0; /* 0.5 is for rounding */
2464 #else
2465                                         if(FLAC__fixedpoint_trunc(fixed_residual_bits_per_sample[fixed_order]) >= (int)subframe_bps)
2466                                                 continue; /* don't even try */
2467                                         rice_parameter = (fixed_residual_bits_per_sample[fixed_order] > FLAC__FP_ZERO)? (unsigned)FLAC__fixedpoint_trunc(fixed_residual_bits_per_sample[fixed_order]+FLAC__FP_ONE_HALF) : 0; /* 0.5 is for rounding */
2468 #endif
2469                                         rice_parameter++; /* to account for the signed->unsigned conversion during rice coding */
2470                                         if(rice_parameter >= FLAC__ENTROPY_CODING_METHOD_PARTITIONED_RICE_ESCAPE_PARAMETER) {
2471 #ifdef DEBUG_VERBOSE
2472                                                 fprintf(stderr, "clipping rice_parameter (%u -> %u) @0\n", rice_parameter, FLAC__ENTROPY_CODING_METHOD_PARTITIONED_RICE_ESCAPE_PARAMETER - 1);
2473 #endif
2474                                                 rice_parameter = FLAC__ENTROPY_CODING_METHOD_PARTITIONED_RICE_ESCAPE_PARAMETER - 1;
2475                                         }
2476                                         _candidate_bits =
2477                                                 evaluate_fixed_subframe_(
2478                                                         encoder,
2479                                                         integer_signal,
2480                                                         residual[!_best_subframe],
2481                                                         encoder->private_->abs_residual,
2482                                                         encoder->private_->abs_residual_partition_sums,
2483                                                         encoder->private_->raw_bits_per_partition,
2484                                                         frame_header->blocksize,
2485                                                         subframe_bps,
2486                                                         fixed_order,
2487                                                         rice_parameter,
2488                                                         min_partition_order,
2489                                                         max_partition_order,
2490                                                         precompute_partition_sums,
2491                                                         encoder->protected_->do_escape_coding,
2492                                                         encoder->protected_->rice_parameter_search_dist,
2493                                                         subframe[!_best_subframe],
2494                                                         partitioned_rice_contents[!_best_subframe]
2495                                                 );
2496                                         if(_candidate_bits < _best_bits) {
2497                                                 _best_subframe = !_best_subframe;
2498                                                 _best_bits = _candidate_bits;
2499                                         }
2500                                 }
2501                         }
2502
2503 #ifndef FLAC__INTEGER_ONLY_LIBRARY
2504                         /* encode lpc */
2505                         if(encoder->protected_->max_lpc_order > 0) {
2506                                 if(encoder->protected_->max_lpc_order >= frame_header->blocksize)
2507                                         max_lpc_order = frame_header->blocksize-1;
2508                                 else
2509                                         max_lpc_order = encoder->protected_->max_lpc_order;
2510                                 if(max_lpc_order > 0) {
2511                                         unsigned a;
2512                                         for (a = 0; a < encoder->protected_->num_apodizations; a++) {
2513                                                 FLAC__lpc_apply_apodization(real_signal, encoder->private_->window[a], encoder->private_->windowed_signal, frame_header->blocksize);
2514                                                 encoder->private_->local_lpc_compute_autocorrelation(encoder->private_->windowed_signal, frame_header->blocksize, max_lpc_order+1, autoc);
2515                                                 /* if autoc[0] == 0.0, the signal is constant and we usually won't get here, but it can happen */
2516                                                 if(autoc[0] != 0.0) {
2517                                                         FLAC__lpc_compute_lp_coefficients(autoc, max_lpc_order, encoder->private_->lp_coeff, lpc_error);
2518                                                         if(encoder->protected_->do_exhaustive_model_search) {
2519                                                                 min_lpc_order = 1;
2520                                                         }
2521                                                         else {
2522                                                                 unsigned guess_lpc_order = FLAC__lpc_compute_best_order(lpc_error, max_lpc_order, frame_header->blocksize, subframe_bps);
2523                                                                 min_lpc_order = max_lpc_order = guess_lpc_order;
2524                                                         }
2525                                                         for(lpc_order = min_lpc_order; lpc_order <= max_lpc_order; lpc_order++) {
2526                                                                 lpc_residual_bits_per_sample = FLAC__lpc_compute_expected_bits_per_residual_sample(lpc_error[lpc_order-1], frame_header->blocksize-lpc_order);
2527                                                                 if(lpc_residual_bits_per_sample >= (FLAC__double)subframe_bps)
2528                                                                         continue; /* don't even try */
2529                                                                 rice_parameter = (lpc_residual_bits_per_sample > 0.0)? (unsigned)(lpc_residual_bits_per_sample+0.5) : 0; /* 0.5 is for rounding */
2530                                                                 rice_parameter++; /* to account for the signed->unsigned conversion during rice coding */
2531                                                                 if(rice_parameter >= FLAC__ENTROPY_CODING_METHOD_PARTITIONED_RICE_ESCAPE_PARAMETER) {
2532         #ifdef DEBUG_VERBOSE
2533                                                                         fprintf(stderr, "clipping rice_parameter (%u -> %u) @1\n", rice_parameter, FLAC__ENTROPY_CODING_METHOD_PARTITIONED_RICE_ESCAPE_PARAMETER - 1);
2534         #endif
2535                                                                         rice_parameter = FLAC__ENTROPY_CODING_METHOD_PARTITIONED_RICE_ESCAPE_PARAMETER - 1;
2536                                                                 }
2537                                                                 if(encoder->protected_->do_qlp_coeff_prec_search) {
2538                                                                         min_qlp_coeff_precision = FLAC__MIN_QLP_COEFF_PRECISION;
2539                                                                         /* ensure a 32-bit datapath throughout for 16bps or less */
2540                                                                         if(subframe_bps <= 16)
2541                                                                                 max_qlp_coeff_precision = min(32 - subframe_bps - lpc_order, FLAC__MAX_QLP_COEFF_PRECISION);
2542                                                                         else
2543                                                                                 max_qlp_coeff_precision = FLAC__MAX_QLP_COEFF_PRECISION;
2544                                                                 }
2545                                                                 else {
2546                                                                         min_qlp_coeff_precision = max_qlp_coeff_precision = encoder->protected_->qlp_coeff_precision;
2547                                                                 }
2548                                                                 for(qlp_coeff_precision = min_qlp_coeff_precision; qlp_coeff_precision <= max_qlp_coeff_precision; qlp_coeff_precision++) {
2549                                                                         _candidate_bits =
2550                                                                                 evaluate_lpc_subframe_(
2551                                                                                         encoder,
2552                                                                                         integer_signal,
2553                                                                                         residual[!_best_subframe],
2554                                                                                         encoder->private_->abs_residual,
2555                                                                                         encoder->private_->abs_residual_partition_sums,
2556                                                                                         encoder->private_->raw_bits_per_partition,
2557                                                                                         encoder->private_->lp_coeff[lpc_order-1],
2558                                                                                         frame_header->blocksize,
2559                                                                                         subframe_bps,
2560                                                                                         lpc_order,
2561                                                                                         qlp_coeff_precision,
2562                                                                                         rice_parameter,
2563                                                                                         min_partition_order,
2564                                                                                         max_partition_order,
2565                                                                                         precompute_partition_sums,
2566                                                                                         encoder->protected_->do_escape_coding,
2567                                                                                         encoder->protected_->rice_parameter_search_dist,
2568                                                                                         subframe[!_best_subframe],
2569                                                                                         partitioned_rice_contents[!_best_subframe]
2570 #ifdef WINDOW_DEBUG_OUTPUT
2571                                                                                         ,frame_header->number.frame_number
2572                                                                                         ,subframe_number
2573                                                                                         ,encoder->protected_->apodizations[a]
2574 #endif
2575                                                                                 );
2576                                                                         if(_candidate_bits > 0) { /* if == 0, there was a problem quantizing the lpcoeffs */
2577                                                                                 if(_candidate_bits < _best_bits) {
2578                                                                                         _best_subframe = !_best_subframe;
2579                                                                                         _best_bits = _candidate_bits;
2580                                                                                 }
2581                                                                         }
2582                                                                 }
2583                                                         }
2584                                                 }
2585                                         }
2586                                 }
2587                         }
2588 #endif /* !defined FLAC__INTEGER_ONLY_LIBRARY */
2589                 }
2590         }
2591
2592         /* under rare circumstances this can happen when all but lpc subframe types are disabled: */
2593         if(_best_bits == UINT_MAX) {
2594                 FLAC__ASSERT(_best_subframe == 0);
2595                 _best_bits = evaluate_verbatim_subframe_(integer_signal, frame_header->blocksize, subframe_bps, subframe[_best_subframe]);
2596         }
2597
2598         *best_subframe = _best_subframe;
2599         *best_bits = _best_bits;
2600
2601         return true;
2602 }
2603
2604 FLAC__bool add_subframe_(
2605         FLAC__StreamEncoder *encoder,
2606         const FLAC__FrameHeader *frame_header,
2607         unsigned subframe_bps,
2608         const FLAC__Subframe *subframe,
2609         FLAC__BitBuffer *frame
2610 #ifdef WINDOW_DEBUG_OUTPUT
2611 ,unsigned subframe_bits
2612 #endif
2613 )
2614 {
2615         switch(subframe->type) {
2616                 case FLAC__SUBFRAME_TYPE_CONSTANT:
2617                         if(!FLAC__subframe_add_constant(&(subframe->data.constant), subframe_bps, subframe->wasted_bits, frame)) {
2618                                 encoder->protected_->state = FLAC__STREAM_ENCODER_FATAL_ERROR_WHILE_ENCODING;
2619                                 return false;
2620                         }
2621                         break;
2622                 case FLAC__SUBFRAME_TYPE_FIXED:
2623                         if(!FLAC__subframe_add_fixed(&(subframe->data.fixed), frame_header->blocksize - subframe->data.fixed.order, subframe_bps, subframe->wasted_bits, frame)) {
2624                                 encoder->protected_->state = FLAC__STREAM_ENCODER_FATAL_ERROR_WHILE_ENCODING;
2625                                 return false;
2626                         }
2627                         break;
2628                 case FLAC__SUBFRAME_TYPE_LPC:
2629 #ifdef WINDOW_DEBUG_OUTPUT
2630                         fprintf(stderr, "WIN:\tframe=%u\tsubframe=?\torder=%u\twindow=%s\tbits=%u\n", frame_header->number.frame_number, subframe->data.lpc.order, subframe->data.lpc.window_type, subframe_bits);
2631 #endif
2632                         if(!FLAC__subframe_add_lpc(&(subframe->data.lpc), frame_header->blocksize - subframe->data.lpc.order, subframe_bps, subframe->wasted_bits, frame)) {
2633                                 encoder->protected_->state = FLAC__STREAM_ENCODER_FATAL_ERROR_WHILE_ENCODING;
2634                                 return false;
2635                         }
2636                         break;
2637                 case FLAC__SUBFRAME_TYPE_VERBATIM:
2638                         if(!FLAC__subframe_add_verbatim(&(subframe->data.verbatim), frame_header->blocksize, subframe_bps, subframe->wasted_bits, frame)) {
2639                                 encoder->protected_->state = FLAC__STREAM_ENCODER_FATAL_ERROR_WHILE_ENCODING;
2640                                 return false;
2641                         }
2642                         break;
2643                 default:
2644                         FLAC__ASSERT(0);
2645         }
2646
2647         return true;
2648 }
2649
2650 unsigned evaluate_constant_subframe_(
2651         const FLAC__int32 signal,
2652         unsigned subframe_bps,
2653         FLAC__Subframe *subframe
2654 )
2655 {
2656         subframe->type = FLAC__SUBFRAME_TYPE_CONSTANT;
2657         subframe->data.constant.value = signal;
2658
2659         return FLAC__SUBFRAME_ZERO_PAD_LEN + FLAC__SUBFRAME_TYPE_LEN + FLAC__SUBFRAME_WASTED_BITS_FLAG_LEN + subframe_bps;
2660 }
2661
2662 unsigned evaluate_fixed_subframe_(
2663         FLAC__StreamEncoder *encoder,
2664         const FLAC__int32 signal[],
2665         FLAC__int32 residual[],
2666         FLAC__uint32 abs_residual[],
2667         FLAC__uint64 abs_residual_partition_sums[],
2668         unsigned raw_bits_per_partition[],
2669         unsigned blocksize,
2670         unsigned subframe_bps,
2671         unsigned order,
2672         unsigned rice_parameter,
2673         unsigned min_partition_order,
2674         unsigned max_partition_order,
2675         FLAC__bool precompute_partition_sums,
2676         FLAC__bool do_escape_coding,
2677         unsigned rice_parameter_search_dist,
2678         FLAC__Subframe *subframe,
2679         FLAC__EntropyCodingMethod_PartitionedRiceContents *partitioned_rice_contents
2680 )
2681 {
2682         unsigned i, residual_bits;
2683         const unsigned residual_samples = blocksize - order;
2684
2685         FLAC__fixed_compute_residual(signal+order, residual_samples, order, residual);
2686
2687         subframe->type = FLAC__SUBFRAME_TYPE_FIXED;
2688
2689         subframe->data.fixed.entropy_coding_method.type = FLAC__ENTROPY_CODING_METHOD_PARTITIONED_RICE;
2690         subframe->data.fixed.entropy_coding_method.data.partitioned_rice.contents = partitioned_rice_contents;
2691         subframe->data.fixed.residual = residual;
2692
2693         residual_bits =
2694                 find_best_partition_order_(
2695                         encoder->private_,
2696                         residual,
2697                         abs_residual,
2698                         abs_residual_partition_sums,
2699                         raw_bits_per_partition,
2700                         residual_samples,
2701                         order,
2702                         rice_parameter,
2703                         min_partition_order,
2704                         max_partition_order,
2705                         precompute_partition_sums,
2706                         do_escape_coding,
2707                         rice_parameter_search_dist,
2708                         &subframe->data.fixed.entropy_coding_method.data.partitioned_rice
2709                 );
2710
2711         subframe->data.fixed.order = order;
2712         for(i = 0; i < order; i++)
2713                 subframe->data.fixed.warmup[i] = signal[i];
2714
2715         return FLAC__SUBFRAME_ZERO_PAD_LEN + FLAC__SUBFRAME_TYPE_LEN + FLAC__SUBFRAME_WASTED_BITS_FLAG_LEN + (order * subframe_bps) + residual_bits;
2716 }
2717
2718 #ifndef FLAC__INTEGER_ONLY_LIBRARY
2719 unsigned evaluate_lpc_subframe_(
2720         FLAC__StreamEncoder *encoder,
2721         const FLAC__int32 signal[],
2722         FLAC__int32 residual[],
2723         FLAC__uint32 abs_residual[],
2724         FLAC__uint64 abs_residual_partition_sums[],
2725         unsigned raw_bits_per_partition[],
2726         const FLAC__real lp_coeff[],
2727         unsigned blocksize,
2728         unsigned subframe_bps,
2729         unsigned order,
2730         unsigned qlp_coeff_precision,
2731         unsigned rice_parameter,
2732         unsigned min_partition_order,
2733         unsigned max_partition_order,
2734         FLAC__bool precompute_partition_sums,
2735         FLAC__bool do_escape_coding,
2736         unsigned rice_parameter_search_dist,
2737         FLAC__Subframe *subframe,
2738         FLAC__EntropyCodingMethod_PartitionedRiceContents *partitioned_rice_contents
2739 #ifdef WINDOW_DEBUG_OUTPUT
2740         ,unsigned frame_number
2741         ,unsigned subframe_number
2742         ,FLAC__ApodizationSpecification aspec
2743 #endif
2744 )
2745 {
2746         FLAC__int32 qlp_coeff[FLAC__MAX_LPC_ORDER];
2747         unsigned i, residual_bits;
2748         int quantization, ret;
2749         const unsigned residual_samples = blocksize - order;
2750
2751         /* try to keep qlp coeff precision such that only 32-bit math is required for decode of <=16bps streams */
2752         if(subframe_bps <= 16) {
2753                 FLAC__ASSERT(order > 0);
2754                 FLAC__ASSERT(order <= FLAC__MAX_LPC_ORDER);
2755                 qlp_coeff_precision = min(qlp_coeff_precision, 32 - subframe_bps - FLAC__bitmath_ilog2(order));
2756         }
2757
2758 #ifdef WINDOW_DEBUG_OUTPUT
2759         if (aspec.type == FLAC__APODIZATION_GAUSS)
2760                 snprintf(subframe->data.lpc.window_type, sizeof subframe->data.lpc.window_type, "%s(%0.5f)", winstr[aspec.type], aspec.parameters.gauss.stddev);
2761         else if (aspec.type == FLAC__APODIZATION_TUKEY)
2762                 snprintf(subframe->data.lpc.window_type, sizeof subframe->data.lpc.window_type, "%s(%0.5f)", winstr[aspec.type], aspec.parameters.tukey.p);
2763         else
2764                 strncpy(subframe->data.lpc.window_type, winstr[aspec.type], sizeof subframe->data.lpc.window_type);
2765 #endif
2766
2767         ret = FLAC__lpc_quantize_coefficients(lp_coeff, order, qlp_coeff_precision, qlp_coeff, &quantization);
2768         if(ret != 0)
2769                 return 0; /* this is a hack to indicate to the caller that we can't do lp at this order on this subframe */
2770
2771         if(subframe_bps + qlp_coeff_precision + FLAC__bitmath_ilog2(order) <= 32)
2772                 if(subframe_bps <= 16 && qlp_coeff_precision <= 16)
2773                         encoder->private_->local_lpc_compute_residual_from_qlp_coefficients_16bit(signal+order, residual_samples, qlp_coeff, order, quantization, residual);
2774                 else
2775                         encoder->private_->local_lpc_compute_residual_from_qlp_coefficients(signal+order, residual_samples, qlp_coeff, order, quantization, residual);
2776         else
2777                 encoder->private_->local_lpc_compute_residual_from_qlp_coefficients_64bit(signal+order, residual_samples, qlp_coeff, order, quantization, residual);
2778
2779         subframe->type = FLAC__SUBFRAME_TYPE_LPC;
2780
2781         subframe->data.lpc.entropy_coding_method.type = FLAC__ENTROPY_CODING_METHOD_PARTITIONED_RICE;
2782         subframe->data.lpc.entropy_coding_method.data.partitioned_rice.contents = partitioned_rice_contents;
2783         subframe->data.lpc.residual = residual;
2784
2785         residual_bits =
2786                 find_best_partition_order_(
2787                         encoder->private_,
2788                         residual,
2789                         abs_residual,
2790                         abs_residual_partition_sums,
2791                         raw_bits_per_partition,
2792                         residual_samples,
2793                         order,
2794                         rice_parameter,
2795                         min_partition_order,
2796                         max_partition_order,
2797                         precompute_partition_sums,
2798                         do_escape_coding,
2799                         rice_parameter_search_dist,
2800                         &subframe->data.fixed.entropy_coding_method.data.partitioned_rice
2801                 );
2802
2803         subframe->data.lpc.order = order;
2804         subframe->data.lpc.qlp_coeff_precision = qlp_coeff_precision;
2805         subframe->data.lpc.quantization_level = quantization;
2806         memcpy(subframe->data.lpc.qlp_coeff, qlp_coeff, sizeof(FLAC__int32)*FLAC__MAX_LPC_ORDER);
2807         for(i = 0; i < order; i++)
2808                 subframe->data.lpc.warmup[i] = signal[i];
2809
2810 #ifdef WINDOW_DEBUG_OUTPUT
2811         fprintf(stderr, "SWIN:\tframe=%u\tsubframe=%u\torder=%u\twindow=%s\tbits=%u\n", frame_number, subframe_number, order, subframe->data.lpc.window_type, FLAC__SUBFRAME_ZERO_PAD_LEN + FLAC__SUBFRAME_TYPE_LEN + FLAC__SUBFRAME_WASTED_BITS_FLAG_LEN + FLAC__SUBFRAME_LPC_QLP_COEFF_PRECISION_LEN + FLAC__SUBFRAME_LPC_QLP_SHIFT_LEN + (order * (qlp_coeff_precision + subframe_bps)) + residual_bits);
2812 #endif
2813         return FLAC__SUBFRAME_ZERO_PAD_LEN + FLAC__SUBFRAME_TYPE_LEN + FLAC__SUBFRAME_WASTED_BITS_FLAG_LEN + FLAC__SUBFRAME_LPC_QLP_COEFF_PRECISION_LEN + FLAC__SUBFRAME_LPC_QLP_SHIFT_LEN + (order * (qlp_coeff_precision + subframe_bps)) + residual_bits;
2814 }
2815 #endif
2816
2817 unsigned evaluate_verbatim_subframe_(
2818         const FLAC__int32 signal[],
2819         unsigned blocksize,
2820         unsigned subframe_bps,
2821         FLAC__Subframe *subframe
2822 )
2823 {
2824         subframe->type = FLAC__SUBFRAME_TYPE_VERBATIM;
2825
2826         subframe->data.verbatim.data = signal;
2827
2828         return FLAC__SUBFRAME_ZERO_PAD_LEN + FLAC__SUBFRAME_TYPE_LEN + FLAC__SUBFRAME_WASTED_BITS_FLAG_LEN + (blocksize * subframe_bps);
2829 }
2830
2831 unsigned find_best_partition_order_(
2832         FLAC__StreamEncoderPrivate *private_,
2833         const FLAC__int32 residual[],
2834         FLAC__uint32 abs_residual[],
2835         FLAC__uint64 abs_residual_partition_sums[],
2836         unsigned raw_bits_per_partition[],
2837         unsigned residual_samples,
2838         unsigned predictor_order,
2839         unsigned rice_parameter,
2840         unsigned min_partition_order,
2841         unsigned max_partition_order,
2842         FLAC__bool precompute_partition_sums,
2843         FLAC__bool do_escape_coding,
2844         unsigned rice_parameter_search_dist,
2845         FLAC__EntropyCodingMethod_PartitionedRice *best_partitioned_rice
2846 )
2847 {
2848         FLAC__int32 r;
2849         unsigned residual_bits, best_residual_bits = 0;
2850         unsigned residual_sample;
2851         unsigned best_parameters_index = 0;
2852         const unsigned blocksize = residual_samples + predictor_order;
2853
2854         /* compute abs(residual) for use later */
2855         for(residual_sample = 0; residual_sample < residual_samples; residual_sample++) {
2856                 r = residual[residual_sample];
2857                 abs_residual[residual_sample] = (FLAC__uint32)(r<0? -r : r);
2858         }
2859
2860         max_partition_order = FLAC__format_get_max_rice_partition_order_from_blocksize_limited_max_and_predictor_order(max_partition_order, blocksize, predictor_order);
2861         min_partition_order = min(min_partition_order, max_partition_order);
2862
2863         if(precompute_partition_sums) {
2864                 int partition_order;
2865                 unsigned sum;
2866
2867                 precompute_partition_info_sums_(abs_residual, abs_residual_partition_sums, residual_samples, predictor_order, min_partition_order, max_partition_order);
2868
2869                 if(do_escape_coding)
2870                         precompute_partition_info_escapes_(residual, raw_bits_per_partition, residual_samples, predictor_order, min_partition_order, max_partition_order);
2871
2872                 for(partition_order = (int)max_partition_order, sum = 0; partition_order >= (int)min_partition_order; partition_order--) {
2873 #ifdef DONT_ESTIMATE_RICE_BITS
2874                         if(!
2875                                 set_partitioned_rice_with_precompute_(
2876                                         residual,
2877                                         abs_residual_partition_sums+sum,
2878                                         raw_bits_per_partition+sum,
2879                                         residual_samples,
2880                                         predictor_order,
2881                                         rice_parameter,
2882                                         rice_parameter_search_dist,
2883                                         (unsigned)partition_order,
2884                                         do_escape_coding,
2885                                         &private_->partitioned_rice_contents_extra[!best_parameters_index],
2886                                         &residual_bits
2887                                 )
2888                         )
2889 #else
2890                         if(!
2891                                 set_partitioned_rice_with_precompute_(
2892                                         abs_residual,
2893                                         abs_residual_partition_sums+sum,
2894                                         raw_bits_per_partition+sum,
2895                                         residual_samples,
2896                                         predictor_order,
2897                                         rice_parameter,
2898                                         rice_parameter_search_dist,
2899                                         (unsigned)partition_order,
2900                                         do_escape_coding,
2901                                         &private_->partitioned_rice_contents_extra[!best_parameters_index],
2902                                         &residual_bits
2903                                 )
2904                         )
2905 #endif
2906                         {
2907                                 FLAC__ASSERT(best_residual_bits != 0);
2908                                 break;
2909                         }
2910                         sum += 1u << partition_order;
2911                         if(best_residual_bits == 0 || residual_bits < best_residual_bits) {
2912                                 best_residual_bits = residual_bits;
2913                                 best_parameters_index = !best_parameters_index;
2914                                 best_partitioned_rice->order = partition_order;
2915                         }
2916                 }
2917         }
2918         else {
2919                 unsigned partition_order;
2920                 for(partition_order = min_partition_order; partition_order <= max_partition_order; partition_order++) {
2921 #ifdef DONT_ESTIMATE_RICE_BITS
2922                         if(!
2923                                 set_partitioned_rice_(
2924                                         abs_residual,
2925                                         residual,
2926                                         residual_samples,
2927                                         predictor_order,
2928                                         rice_parameter,
2929                                         rice_parameter_search_dist,
2930                                         partition_order,
2931                                         &private_->partitioned_rice_contents_extra[!best_parameters_index],
2932                                         &residual_bits
2933                                 )
2934                         )
2935 #else
2936                         if(!
2937                                 set_partitioned_rice_(
2938                                         abs_residual,
2939                                         residual_samples,
2940                                         predictor_order,
2941                                         rice_parameter,
2942                                         rice_parameter_search_dist,
2943                                         partition_order,
2944                                         &private_->partitioned_rice_contents_extra[!best_parameters_index],
2945                                         &residual_bits
2946                                 )
2947                         )
2948 #endif
2949                         {
2950                                 FLAC__ASSERT(best_residual_bits != 0);
2951                                 break;
2952                         }
2953                         if(best_residual_bits == 0 || residual_bits < best_residual_bits) {
2954                                 best_residual_bits = residual_bits;
2955                                 best_parameters_index = !best_parameters_index;
2956                                 best_partitioned_rice->order = partition_order;
2957                         }
2958                 }
2959         }
2960
2961         /*
2962          * We are allowed to de-const the pointer based on our special knowledge;
2963          * it is const to the outside world.
2964          */
2965         {
2966                 FLAC__EntropyCodingMethod_PartitionedRiceContents* best_partitioned_rice_contents = (FLAC__EntropyCodingMethod_PartitionedRiceContents*)best_partitioned_rice->contents;
2967                 FLAC__format_entropy_coding_method_partitioned_rice_contents_ensure_size(best_partitioned_rice_contents, max(6, best_partitioned_rice->order));
2968                 memcpy(best_partitioned_rice_contents->parameters, private_->partitioned_rice_contents_extra[best_parameters_index].parameters, sizeof(unsigned)*(1<<(best_partitioned_rice->order)));
2969                 memcpy(best_partitioned_rice_contents->raw_bits, private_->partitioned_rice_contents_extra[best_parameters_index].raw_bits, sizeof(unsigned)*(1<<(best_partitioned_rice->order)));
2970         }
2971
2972         return best_residual_bits;
2973 }
2974
2975 void precompute_partition_info_sums_(
2976         const FLAC__uint32 abs_residual[],
2977         FLAC__uint64 abs_residual_partition_sums[],
2978         unsigned residual_samples,
2979         unsigned predictor_order,
2980         unsigned min_partition_order,
2981         unsigned max_partition_order
2982 )
2983 {
2984         int partition_order;
2985         unsigned from_partition, to_partition = 0;
2986         const unsigned blocksize = residual_samples + predictor_order;
2987
2988         /* first do max_partition_order */
2989         for(partition_order = (int)max_partition_order; partition_order >= 0; partition_order--) {
2990                 FLAC__uint64 abs_residual_partition_sum;
2991                 FLAC__uint32 abs_r;
2992                 unsigned partition, partition_sample, partition_samples, residual_sample;
2993                 const unsigned partitions = 1u << partition_or