struct AACISError ff_aac_is_encoding_err(AACEncContext *s, ChannelElement *cpe,

                                         int start, int w, int g, float ener0,

                                         float ener1, float ener01,

                                         int use_pcoeffs, int phase)

{

    int i, w2;

    SingleChannelElement *sce0 = &cpe->ch[0];

    SingleChannelElement *sce1 = &cpe->ch[1];

    float *L = use_pcoeffs ? sce0->pcoeffs : sce0->coeffs;

    float *R = use_pcoeffs ? sce1->pcoeffs : sce1->coeffs;

    float *L34 = &s->scoefs[256*0], *R34 = &s->scoefs[256*1];

    float *IS  = &s->scoefs[256*2], *I34 = &s->scoefs[256*3];

    float dist1 = 0.0f, dist2 = 0.0f;

    struct AACISError is_error = {0};



    for (w2 = 0; w2 < sce0->ics.group_len[w]; w2++) {

        FFPsyBand *band0 = &s->psy.ch[s->cur_channel+0].psy_bands[(w+w2)*16+g];

        FFPsyBand *band1 = &s->psy.ch[s->cur_channel+1].psy_bands[(w+w2)*16+g];

        int is_band_type, is_sf_idx = FFMAX(1, sce0->sf_idx[(w+w2)*16+g]-4);

        float e01_34 = phase*pow(ener1/ener0, 3.0/4.0);

        float maxval, dist_spec_err = 0.0f;

        float minthr = FFMIN(band0->threshold, band1->threshold);

        for (i = 0; i < sce0->ics.swb_sizes[g]; i++)

            IS[i] = (L[start+(w+w2)*128+i] + phase*R[start+(w+w2)*128+i])*sqrt(ener0/ener01);

        abs_pow34_v(L34, &L[start+(w+w2)*128], sce0->ics.swb_sizes[g]);

        abs_pow34_v(R34, &R[start+(w+w2)*128], sce0->ics.swb_sizes[g]);

        abs_pow34_v(I34, IS,                   sce0->ics.swb_sizes[g]);

        maxval = find_max_val(1, sce0->ics.swb_sizes[g], I34);

        is_band_type = find_min_book(maxval, is_sf_idx);

        dist1 += quantize_band_cost(s, &L[start + (w+w2)*128], L34,

                                    sce0->ics.swb_sizes[g],

                                    sce0->sf_idx[(w+w2)*16+g],

                                    sce0->band_type[(w+w2)*16+g],

                                    s->lambda / band0->threshold, INFINITY, NULL, 0);

        dist1 += quantize_band_cost(s, &R[start + (w+w2)*128], R34,

                                    sce1->ics.swb_sizes[g],

                                    sce1->sf_idx[(w+w2)*16+g],

                                    sce1->band_type[(w+w2)*16+g],

                                    s->lambda / band1->threshold, INFINITY, NULL, 0);

        dist2 += quantize_band_cost(s, IS, I34, sce0->ics.swb_sizes[g],

                                    is_sf_idx, is_band_type,

                                    s->lambda / minthr, INFINITY, NULL, 0);

        for (i = 0; i < sce0->ics.swb_sizes[g]; i++) {

            dist_spec_err += (L34[i] - I34[i])*(L34[i] - I34[i]);

            dist_spec_err += (R34[i] - I34[i]*e01_34)*(R34[i] - I34[i]*e01_34);

        }

        dist_spec_err *= s->lambda / minthr;

        dist2 += dist_spec_err;

    }



    is_error.pass = dist2 <= dist1;

    is_error.phase = phase;

    is_error.error = fabsf(dist1 - dist2);

    is_error.dist1 = dist1;

    is_error.dist2 = dist2;



    return is_error;

}

struct AACISError ff_aac_is_encoding_err(AACEncContext *s, ChannelElement *cpe,

                                         int start, int w, int g, float ener0,

                                         float ener1, float ener01,

                                         int use_pcoeffs, int phase)

{

    int i, w2;

    SingleChannelElement *sce0 = &cpe->ch[0];

    SingleChannelElement *sce1 = &cpe->ch[1];

    float *L = use_pcoeffs ? sce0->pcoeffs : sce0->coeffs;

    float *R = use_pcoeffs ? sce1->pcoeffs : sce1->coeffs;

    float *L34 = &s->scoefs[256*0], *R34 = &s->scoefs[256*1];

    float *IS  = &s->scoefs[256*2], *I34 = &s->scoefs[256*3];

    float dist1 = 0.0f, dist2 = 0.0f;

    struct AACISError is_error = {0};



    for (w2 = 0; w2 < sce0->ics.group_len[w]; w2++) {

        FFPsyBand *band0 = &s->psy.ch[s->cur_channel+0].psy_bands[(w+w2)*16+g];

        FFPsyBand *band1 = &s->psy.ch[s->cur_channel+1].psy_bands[(w+w2)*16+g];

        int is_band_type, is_sf_idx = FFMAX(1, sce0->sf_idx[(w+w2)*16+g]-4);

        float e01_34 = phase*pow(ener1/ener0, 3.0/4.0);

        float maxval, dist_spec_err = 0.0f;

        float minthr = FFMIN(band0->threshold, band1->threshold);

        for (i = 0; i < sce0->ics.swb_sizes[g]; i++)

            IS[i] = (L[start+(w+w2)*128+i] + phase*R[start+(w+w2)*128+i])*sqrt(ener0/ener01);

        abs_pow34_v(L34, &L[start+(w+w2)*128], sce0->ics.swb_sizes[g]);

        abs_pow34_v(R34, &R[start+(w+w2)*128], sce0->ics.swb_sizes[g]);

        abs_pow34_v(I34, IS,                   sce0->ics.swb_sizes[g]);

        maxval = find_max_val(1, sce0->ics.swb_sizes[g], I34);

        is_band_type = find_min_book(maxval, is_sf_idx);

        dist1 += quantize_band_cost(s, &L[start + (w+w2)*128], L34,

                                    sce0->ics.swb_sizes[g],

                                    sce0->sf_idx[(w+w2)*16+g],

                                    sce0->band_type[(w+w2)*16+g],

                                    s->lambda / band0->threshold, INFINITY, NULL, 0);

        dist1 += quantize_band_cost(s, &R[start + (w+w2)*128], R34,

                                    sce1->ics.swb_sizes[g],

                                    sce1->sf_idx[(w+w2)*16+g],

                                    sce1->band_type[(w+w2)*16+g],

                                    s->lambda / band1->threshold, INFINITY, NULL, 0);

        dist2 += quantize_band_cost(s, IS, I34, sce0->ics.swb_sizes[g],

                                    is_sf_idx, is_band_type,

                                    s->lambda / minthr, INFINITY, NULL, 0);

        for (i = 0; i < sce0->ics.swb_sizes[g]; i++) {

            dist_spec_err += (L34[i] - I34[i])*(L34[i] - I34[i]);

            dist_spec_err += (R34[i] - I34[i]*e01_34)*(R34[i] - I34[i]*e01_34);

        }

        dist_spec_err *= s->lambda / minthr;

        dist2 += dist_spec_err;

    }



    is_error.pass = dist2 <= dist1;

    is_error.phase = phase;

    is_error.error = fabsf(dist1 - dist2);

    is_error.dist1 = dist1;

    is_error.dist2 = dist2;



    return is_error;

}

                                    s->lambda / band0->threshold, INFINITY, NULL, 0);                                     s->lambda / band1->threshold, INFINITY, NULL, 0);                                     s->lambda / minthr, INFINITY, NULL, 0);
struct AACISError FUNC_0(AACEncContext *VAR_0, ChannelElement *VAR_1,

                                         int VAR_2, int VAR_3, int VAR_4, float VAR_5,

                                         float VAR_6, float VAR_7,

                                         int VAR_8, int VAR_9)

{

    int VAR_10, VAR_11;

    SingleChannelElement *sce0 = &VAR_1->ch[0];

    SingleChannelElement *sce1 = &VAR_1->ch[1];

    float *VAR_12 = VAR_8 ? sce0->pcoeffs : sce0->coeffs;

    float *VAR_13 = VAR_8 ? sce1->pcoeffs : sce1->coeffs;

    float *VAR_14 = &VAR_0->scoefs[256*0], *VAR_15 = &VAR_0->scoefs[256*1];

    float *VAR_16  = &VAR_0->scoefs[256*2], *VAR_17 = &VAR_0->scoefs[256*3];

    float VAR_18 = 0.0f, VAR_19 = 0.0f;

    struct AACISError VAR_20 = {0};



    for (VAR_11 = 0; VAR_11 < sce0->ics.group_len[VAR_3]; VAR_11++) {

        FFPsyBand *band0 = &VAR_0->psy.ch[VAR_0->cur_channel+0].psy_bands[(VAR_3+VAR_11)*16+VAR_4];

        FFPsyBand *band1 = &VAR_0->psy.ch[VAR_0->cur_channel+1].psy_bands[(VAR_3+VAR_11)*16+VAR_4];

        int is_band_type, is_sf_idx = FFMAX(1, sce0->sf_idx[(VAR_3+VAR_11)*16+VAR_4]-4);

        float e01_34 = VAR_9*pow(VAR_6/VAR_5, 3.0/4.0);

        float maxval, dist_spec_err = 0.0f;

        float minthr = FFMIN(band0->threshold, band1->threshold);

        for (VAR_10 = 0; VAR_10 < sce0->ics.swb_sizes[VAR_4]; VAR_10++)

            VAR_16[VAR_10] = (VAR_12[VAR_2+(VAR_3+VAR_11)*128+VAR_10] + VAR_9*VAR_13[VAR_2+(VAR_3+VAR_11)*128+VAR_10])*sqrt(VAR_5/VAR_7);

        abs_pow34_v(VAR_14, &VAR_12[VAR_2+(VAR_3+VAR_11)*128], sce0->ics.swb_sizes[VAR_4]);

        abs_pow34_v(VAR_15, &VAR_13[VAR_2+(VAR_3+VAR_11)*128], sce0->ics.swb_sizes[VAR_4]);

        abs_pow34_v(VAR_17, VAR_16,                   sce0->ics.swb_sizes[VAR_4]);

        maxval = find_max_val(1, sce0->ics.swb_sizes[VAR_4], VAR_17);

        is_band_type = find_min_book(maxval, is_sf_idx);

        VAR_18 += quantize_band_cost(VAR_0, &VAR_12[VAR_2 + (VAR_3+VAR_11)*128], VAR_14,

                                    sce0->ics.swb_sizes[VAR_4],

                                    sce0->sf_idx[(VAR_3+VAR_11)*16+VAR_4],

                                    sce0->band_type[(VAR_3+VAR_11)*16+VAR_4],

                                    VAR_0->lambda / band0->threshold, INFINITY, NULL, 0);

        VAR_18 += quantize_band_cost(VAR_0, &VAR_13[VAR_2 + (VAR_3+VAR_11)*128], VAR_15,

                                    sce1->ics.swb_sizes[VAR_4],

                                    sce1->sf_idx[(VAR_3+VAR_11)*16+VAR_4],

                                    sce1->band_type[(VAR_3+VAR_11)*16+VAR_4],

                                    VAR_0->lambda / band1->threshold, INFINITY, NULL, 0);

        VAR_19 += quantize_band_cost(VAR_0, VAR_16, VAR_17, sce0->ics.swb_sizes[VAR_4],

                                    is_sf_idx, is_band_type,

                                    VAR_0->lambda / minthr, INFINITY, NULL, 0);

        for (VAR_10 = 0; VAR_10 < sce0->ics.swb_sizes[VAR_4]; VAR_10++) {

            dist_spec_err += (VAR_14[VAR_10] - VAR_17[VAR_10])*(VAR_14[VAR_10] - VAR_17[VAR_10]);

            dist_spec_err += (VAR_15[VAR_10] - VAR_17[VAR_10]*e01_34)*(VAR_15[VAR_10] - VAR_17[VAR_10]*e01_34);

        }

        dist_spec_err *= VAR_0->lambda / minthr;

        VAR_19 += dist_spec_err;

    }



    VAR_20.pass = VAR_19 <= VAR_18;

    VAR_20.VAR_9 = VAR_9;

    VAR_20.error = fabsf(VAR_18 - VAR_19);

    VAR_20.VAR_18 = VAR_18;

    VAR_20.VAR_19 = VAR_19;



    return VAR_20;

}

struct AACISError FUNC_0(AACEncContext *VAR_0, ChannelElement *VAR_1,
int VAR_2, int VAR_3, int VAR_4, float VAR_5,
float VAR_6, float VAR_7,
int VAR_8, int VAR_9)
{ int VAR_10, VAR_11; SingleChannelElement *sce0 = &VAR_1->ch[0]; SingleChannelElement *sce1 = &VAR_1->ch[1]; float *VAR_12 = VAR_8 ? sce0->pcoeffs : sce0->coeffs; float *VAR_13 = VAR_8 ? sce1->pcoeffs : sce1->coeffs; float *VAR_14 = &VAR_0->scoefs[256*0], *VAR_15 = &VAR_0->scoefs[256*1]; float *VAR_16  = &VAR_0->scoefs[256*2], *VAR_17 = &VAR_0->scoefs[256*3]; float VAR_18 = 0.0f, VAR_19 = 0.0f; struct AACISError VAR_20 = {0}; for (VAR_11 = 0; VAR_11 < sce0->ics.group_len[VAR_3]; VAR_11++) { FFPsyBand *band0 = &VAR_0->psy.ch[VAR_0->cur_channel+0].psy_bands[(VAR_3+VAR_11)*16+VAR_4]; FFPsyBand *band1 = &VAR_0->psy.ch[VAR_0->cur_channel+1].psy_bands[(VAR_3+VAR_11)*16+VAR_4]; int is_band_type, is_sf_idx = FFMAX(1, sce0->sf_idx[(VAR_3+VAR_11)*16+VAR_4]-4); float e01_34 = VAR_9*pow(VAR_6/VAR_5, 3.0/4.0); float maxval, dist_spec_err = 0.0f; float minthr = FFMIN(band0->threshold, band1->threshold); for (VAR_10 = 0; VAR_10 < sce0->ics.swb_sizes[VAR_4]; VAR_10++) VAR_16[VAR_10] = (VAR_12[VAR_2+(VAR_3+VAR_11)*128+VAR_10] + VAR_9*VAR_13[VAR_2+(VAR_3+VAR_11)*128+VAR_10])*sqrt(VAR_5/VAR_7); abs_pow34_v(VAR_14, &VAR_12[VAR_2+(VAR_3+VAR_11)*128], sce0->ics.swb_sizes[VAR_4]); abs_pow34_v(VAR_15, &VAR_13[VAR_2+(VAR_3+VAR_11)*128], sce0->ics.swb_sizes[VAR_4]); abs_pow34_v(VAR_17, VAR_16,                   sce0->ics.swb_sizes[VAR_4]); maxval = find_max_val(1, sce0->ics.swb_sizes[VAR_4], VAR_17); is_band_type = find_min_book(maxval, is_sf_idx); VAR_18 += quantize_band_cost(VAR_0, &VAR_12[VAR_2 + (VAR_3+VAR_11)*128], VAR_14,
sce0->ics.swb_sizes[VAR_4],
sce0->sf_idx[(VAR_3+VAR_11)*16+VAR_4],
sce0->band_type[(VAR_3+VAR_11)*16+VAR_4],
VAR_0->lambda / band0->threshold, INFINITY, NULL, 0); VAR_18 += quantize_band_cost(VAR_0, &VAR_13[VAR_2 + (VAR_3+VAR_11)*128], VAR_15,
sce1->ics.swb_sizes[VAR_4],
sce1->sf_idx[(VAR_3+VAR_11)*16+VAR_4],
sce1->band_type[(VAR_3+VAR_11)*16+VAR_4],
VAR_0->lambda / band1->threshold, INFINITY, NULL, 0); VAR_19 += quantize_band_cost(VAR_0, VAR_16, VAR_17, sce0->ics.swb_sizes[VAR_4],
is_sf_idx, is_band_type,
VAR_0->lambda / minthr, INFINITY, NULL, 0); for (VAR_10 = 0; VAR_10 < sce0->ics.swb_sizes[VAR_4]; VAR_10++) { dist_spec_err += (VAR_14[VAR_10] - VAR_17[VAR_10])*(VAR_14[VAR_10] - VAR_17[VAR_10]); dist_spec_err += (VAR_15[VAR_10] - VAR_17[VAR_10]*e01_34)*(VAR_15[VAR_10] - VAR_17[VAR_10]*e01_34); } dist_spec_err *= VAR_0->lambda / minthr; VAR_19 += dist_spec_err; } VAR_20.pass = VAR_19 <= VAR_18; VAR_20.VAR_9 = VAR_9; VAR_20.error = fabsf(VAR_18 - VAR_19); VAR_20.VAR_18 = VAR_18; VAR_20.VAR_19 = VAR_19; return VAR_20; }