int ff_jpeg2000_init_component(Jpeg2000Component *comp,

                               Jpeg2000CodingStyle *codsty,

                               Jpeg2000QuantStyle *qntsty,

                               int cbps, int dx, int dy,

                               AVCodecContext *avctx)

{

    uint8_t log2_band_prec_width, log2_band_prec_height;

    int reslevelno, bandno, gbandno = 0, ret, i, j;

    uint32_t csize = 1;



    if (!codsty->nreslevels2decode) {

        av_log(avctx, AV_LOG_ERROR, "nreslevels2decode uninitialized\n");

        return AVERROR_INVALIDDATA;

    }



    if (ret = ff_jpeg2000_dwt_init(&comp->dwt, comp->coord,

                                   codsty->nreslevels2decode - 1,

                                   codsty->transform))

        return ret;

    // component size comp->coord is uint16_t so ir cannot overflow

    csize = (comp->coord[0][1] - comp->coord[0][0]) *

            (comp->coord[1][1] - comp->coord[1][0]);



    comp->data = av_malloc_array(csize, sizeof(*comp->data));

    if (!comp->data)

        return AVERROR(ENOMEM);

    comp->reslevel = av_malloc_array(codsty->nreslevels, sizeof(*comp->reslevel));

    if (!comp->reslevel)

        return AVERROR(ENOMEM);

    /* LOOP on resolution levels */

    for (reslevelno = 0; reslevelno < codsty->nreslevels; reslevelno++) {

        int declvl = codsty->nreslevels - reslevelno;    // N_L -r see  ISO/IEC 15444-1:2002 B.5

        Jpeg2000ResLevel *reslevel = comp->reslevel + reslevelno;



        /* Compute borders for each resolution level.

         * Computation of trx_0, trx_1, try_0 and try_1.

         * see ISO/IEC 15444-1:2002 eq. B.5 and B-14 */

        for (i = 0; i < 2; i++)

            for (j = 0; j < 2; j++)

                reslevel->coord[i][j] =

                    ff_jpeg2000_ceildivpow2(comp->coord_o[i][j], declvl - 1);

        // update precincts size: 2^n value

        reslevel->log2_prec_width  = codsty->log2_prec_widths[reslevelno];

        reslevel->log2_prec_height = codsty->log2_prec_heights[reslevelno];



        /* Number of bands for each resolution level */

        if (reslevelno == 0)

            reslevel->nbands = 1;

        else

            reslevel->nbands = 3;



        /* Number of precincts wich span the tile for resolution level reslevelno

         * see B.6 in ISO/IEC 15444-1:2002 eq. B-16

         * num_precincts_x = |- trx_1 / 2 ^ log2_prec_width) -| - (trx_0 / 2 ^ log2_prec_width)

         * num_precincts_y = |- try_1 / 2 ^ log2_prec_width) -| - (try_0 / 2 ^ log2_prec_width)

         * for Dcinema profiles in JPEG 2000

         * num_precincts_x = |- trx_1 / 2 ^ log2_prec_width) -|

         * num_precincts_y = |- try_1 / 2 ^ log2_prec_width) -| */

        if (reslevel->coord[0][1] == reslevel->coord[0][0])

            reslevel->num_precincts_x = 0;

        else

            reslevel->num_precincts_x =

                ff_jpeg2000_ceildivpow2(reslevel->coord[0][1],

                                        reslevel->log2_prec_width) -

                (reslevel->coord[0][0] >> reslevel->log2_prec_width);



        if (reslevel->coord[1][1] == reslevel->coord[1][0])

            reslevel->num_precincts_y = 0;

        else

            reslevel->num_precincts_y =

                ff_jpeg2000_ceildivpow2(reslevel->coord[1][1],

                                        reslevel->log2_prec_height) -

                (reslevel->coord[1][0] >> reslevel->log2_prec_height);



        reslevel->band = av_malloc_array(reslevel->nbands, sizeof(*reslevel->band));

        if (!reslevel->band)

            return AVERROR(ENOMEM);



        for (bandno = 0; bandno < reslevel->nbands; bandno++, gbandno++) {

            Jpeg2000Band *band = reslevel->band + bandno;

            int cblkno, precno;

            int nb_precincts;



            /* TODO: Implementation of quantization step not finished,

             * see ISO/IEC 15444-1:2002 E.1 and A.6.4. */

            switch (qntsty->quantsty) {

                uint8_t gain;

                int numbps;

            case JPEG2000_QSTY_NONE:

                /* TODO: to verify. No quantization in this case */

                numbps = cbps +

                         lut_gain[codsty->transform][bandno + reslevelno > 0];

                band->stepsize = (float)SHL(2048 + qntsty->mant[gbandno],

                                            2 + numbps - qntsty->expn[gbandno]);

                break;

            case JPEG2000_QSTY_SI:

                /*TODO: Compute formula to implement. */

                band->stepsize = (float) (1 << 13);

                break;

            case JPEG2000_QSTY_SE:

                /* Exponent quantization step.

                 * Formula:

                 * delta_b = 2 ^ (R_b - expn_b) * (1 + (mant_b / 2 ^ 11))

                 * R_b = R_I + log2 (gain_b )

                 * see ISO/IEC 15444-1:2002 E.1.1 eqn. E-3 and E-4 */

                /* TODO/WARN: value of log2 (gain_b ) not taken into account

                 * but it works (compared to OpenJPEG). Why?

                 * Further investigation needed. */

                gain            = cbps;

                band->stepsize  = pow(2.0, gain - qntsty->expn[gbandno]);

                band->stepsize *= (float)qntsty->mant[gbandno] / 2048.0 + 1.0;

                /* FIXME: In openjepg code stespize = stepsize * 0.5. Why?

                 * If not set output of entropic decoder is not correct. */

                band->stepsize *= 0.5;

                break;

            default:

                band->stepsize = 0;

                av_log(avctx, AV_LOG_ERROR, "Unknown quantization format\n");

                break;

            }

            /* BITEXACT computing case --> convert to int */

            if (avctx->flags & CODEC_FLAG_BITEXACT)

                band->stepsize = (int32_t)(band->stepsize * (1 << 16));



            /* computation of tbx_0, tbx_1, tby_0, tby_1

             * see ISO/IEC 15444-1:2002 B.5 eq. B-15 and tbl B.1

             * codeblock width and height is computed for

             * DCI JPEG 2000 codeblock_width = codeblock_width = 32 = 2 ^ 5 */

            if (reslevelno == 0) {

                /* for reslevelno = 0, only one band, x0_b = y0_b = 0 */

                for (i = 0; i < 2; i++)

                    for (j = 0; j < 2; j++)

                        band->coord[i][j] =

                            ff_jpeg2000_ceildivpow2(comp->coord_o[i][j],

                                                    declvl - 1);



                log2_band_prec_width  = reslevel->log2_prec_width;

                log2_band_prec_height = reslevel->log2_prec_height;

                /* see ISO/IEC 15444-1:2002 eq. B-17 and eq. B-15 */

                band->log2_cblk_width  = FFMIN(codsty->log2_cblk_width,

                                               reslevel->log2_prec_width);

                band->log2_cblk_height = FFMIN(codsty->log2_cblk_height,

                                               reslevel->log2_prec_height);

            } else {

                /* 3 bands x0_b = 1 y0_b = 0; x0_b = 0 y0_b = 1; x0_b = y0_b = 1 */

                /* x0_b and y0_b are computed with ((bandno + 1 >> i) & 1) */

                for (i = 0; i < 2; i++)

                    for (j = 0; j < 2; j++)

                        /* Formula example for tbx_0 = ceildiv((tcx_0 - 2 ^ (declvl - 1) * x0_b) / declvl) */

                        band->coord[i][j] =

                            ff_jpeg2000_ceildivpow2(comp->coord_o[i][j] -

                                                    (((bandno + 1 >> i) & 1) << declvl - 1),

                                                    declvl);

                /* TODO: Manage case of 3 band offsets here or

                 * in coding/decoding function? */



                /* see ISO/IEC 15444-1:2002 eq. B-17 and eq. B-15 */

                band->log2_cblk_width  = FFMIN(codsty->log2_cblk_width,

                                               reslevel->log2_prec_width - 1);

                band->log2_cblk_height = FFMIN(codsty->log2_cblk_height,

                                               reslevel->log2_prec_height - 1);



                log2_band_prec_width  = reslevel->log2_prec_width  - 1;

                log2_band_prec_height = reslevel->log2_prec_height - 1;

            }



            band->prec = av_malloc_array(reslevel->num_precincts_x *

                                         reslevel->num_precincts_y,

                                         sizeof(*band->prec));

            if (!band->prec)

                return AVERROR(ENOMEM);



            nb_precincts = reslevel->num_precincts_x * reslevel->num_precincts_y;



            for (precno = 0; precno < nb_precincts; precno++) {

                Jpeg2000Prec *prec = band->prec + precno;



                /* TODO: Explain formula for JPEG200 DCINEMA. */

                /* TODO: Verify with previous count of codeblocks per band */



                /* Compute P_x0 */

                prec->coord[0][0] = (precno % reslevel->num_precincts_x) *

                                    (1 << log2_band_prec_width);

                prec->coord[0][0] = FFMAX(prec->coord[0][0], band->coord[0][0]);



                /* Compute P_y0 */

                prec->coord[1][0] = (precno / reslevel->num_precincts_x) *

                                    (1 << log2_band_prec_height);

                prec->coord[1][0] = FFMAX(prec->coord[1][0], band->coord[1][0]);



                /* Compute P_x1 */

                prec->coord[0][1] = prec->coord[0][0] +

                                    (1 << log2_band_prec_width);

                prec->coord[0][1] = FFMIN(prec->coord[0][1], band->coord[0][1]);



                /* Compute P_y1 */

                prec->coord[1][1] = prec->coord[1][0] +

                                    (1 << log2_band_prec_height);

                prec->coord[1][1] = FFMIN(prec->coord[1][1], band->coord[1][1]);



                prec->nb_codeblocks_width =

                    ff_jpeg2000_ceildivpow2(prec->coord[0][1] -

                                            prec->coord[0][0],

                                            band->log2_cblk_width);

                prec->nb_codeblocks_height =

                    ff_jpeg2000_ceildivpow2(prec->coord[1][1] -

                                            prec->coord[1][0],

                                            band->log2_cblk_height);



                /* Tag trees initialization */

                prec->cblkincl =

                    ff_jpeg2000_tag_tree_init(prec->nb_codeblocks_width,

                                              prec->nb_codeblocks_height);

                if (!prec->cblkincl)

                    return AVERROR(ENOMEM);



                prec->zerobits =

                    ff_jpeg2000_tag_tree_init(prec->nb_codeblocks_width,

                                              prec->nb_codeblocks_height);

                if (!prec->zerobits)

                    return AVERROR(ENOMEM);



                prec->cblk = av_malloc_array(prec->nb_codeblocks_width *

                                             prec->nb_codeblocks_height,

                                             sizeof(*prec->cblk));

                if (!prec->cblk)

                    return AVERROR(ENOMEM);

                for (cblkno = 0; cblkno < prec->nb_codeblocks_width * prec->nb_codeblocks_height; cblkno++) {

                    Jpeg2000Cblk *cblk = prec->cblk + cblkno;

                    uint16_t Cx0, Cy0;



                    /* Compute coordinates of codeblocks */

                    /* Compute Cx0*/

                    Cx0 = (prec->coord[0][0] >> band->log2_cblk_width) << band->log2_cblk_width;

                    Cx0 = Cx0 + ((cblkno % prec->nb_codeblocks_width)  << band->log2_cblk_width);

                    cblk->coord[0][0] = FFMAX(Cx0, prec->coord[0][0]);



                    /* Compute Cy0*/

                    Cy0 = (prec->coord[1][0] >> band->log2_cblk_height) << band->log2_cblk_height;

                    Cy0 = Cy0 + ((cblkno / prec->nb_codeblocks_width)   << band->log2_cblk_height);

                    cblk->coord[1][0] = FFMAX(Cy0, prec->coord[1][0]);



                    /* Compute Cx1 */

                    cblk->coord[0][1] = FFMIN(Cx0 + (1 << band->log2_cblk_width),

                                              prec->coord[0][1]);



                    /* Compute Cy1 */

                    cblk->coord[1][1] = FFMIN(Cy0 + (1 << band->log2_cblk_height),

                                              prec->coord[1][1]);

                    cblk->zero      = 0;

                    cblk->lblock    = 3;

                    cblk->length    = 0;

                    cblk->lengthinc = 0;

                    cblk->npasses   = 0;

                }

            }

        }

    }

    return 0;

}

int ff_jpeg2000_init_component(Jpeg2000Component *comp,

                               Jpeg2000CodingStyle *codsty,

                               Jpeg2000QuantStyle *qntsty,

                               int cbps, int dx, int dy,

                               AVCodecContext *avctx)

{

    uint8_t log2_band_prec_width, log2_band_prec_height;

    int reslevelno, bandno, gbandno = 0, ret, i, j;

    uint32_t csize = 1;



    if (!codsty->nreslevels2decode) {

        av_log(avctx, AV_LOG_ERROR, "nreslevels2decode uninitialized\n");

        return AVERROR_INVALIDDATA;

    }



    if (ret = ff_jpeg2000_dwt_init(&comp->dwt, comp->coord,

                                   codsty->nreslevels2decode - 1,

                                   codsty->transform))

        return ret;

    

    csize = (comp->coord[0][1] - comp->coord[0][0]) *

            (comp->coord[1][1] - comp->coord[1][0]);



    comp->data = av_malloc_array(csize, sizeof(*comp->data));

    if (!comp->data)

        return AVERROR(ENOMEM);

    comp->reslevel = av_malloc_array(codsty->nreslevels, sizeof(*comp->reslevel));

    if (!comp->reslevel)

        return AVERROR(ENOMEM);

    

    for (reslevelno = 0; reslevelno < codsty->nreslevels; reslevelno++) {

        int declvl = codsty->nreslevels - reslevelno;    

        Jpeg2000ResLevel *reslevel = comp->reslevel + reslevelno;



        





        for (i = 0; i < 2; i++)

            for (j = 0; j < 2; j++)

                reslevel->coord[i][j] =

                    ff_jpeg2000_ceildivpow2(comp->coord_o[i][j], declvl - 1);

        

        reslevel->log2_prec_width  = codsty->log2_prec_widths[reslevelno];

        reslevel->log2_prec_height = codsty->log2_prec_heights[reslevelno];



        

        if (reslevelno == 0)

            reslevel->nbands = 1;

        else

            reslevel->nbands = 3;



        













        if (reslevel->coord[0][1] == reslevel->coord[0][0])

            reslevel->num_precincts_x = 0;

        else

            reslevel->num_precincts_x =

                ff_jpeg2000_ceildivpow2(reslevel->coord[0][1],

                                        reslevel->log2_prec_width) -

                (reslevel->coord[0][0] >> reslevel->log2_prec_width);



        if (reslevel->coord[1][1] == reslevel->coord[1][0])

            reslevel->num_precincts_y = 0;

        else

            reslevel->num_precincts_y =

                ff_jpeg2000_ceildivpow2(reslevel->coord[1][1],

                                        reslevel->log2_prec_height) -

                (reslevel->coord[1][0] >> reslevel->log2_prec_height);



        reslevel->band = av_malloc_array(reslevel->nbands, sizeof(*reslevel->band));

        if (!reslevel->band)

            return AVERROR(ENOMEM);



        for (bandno = 0; bandno < reslevel->nbands; bandno++, gbandno++) {

            Jpeg2000Band *band = reslevel->band + bandno;

            int cblkno, precno;

            int nb_precincts;



            



            switch (qntsty->quantsty) {

                uint8_t gain;

                int numbps;

            case JPEG2000_QSTY_NONE:

                

                numbps = cbps +

                         lut_gain[codsty->transform][bandno + reslevelno > 0];

                band->stepsize = (float)SHL(2048 + qntsty->mant[gbandno],

                                            2 + numbps - qntsty->expn[gbandno]);

                break;

            case JPEG2000_QSTY_SI:

                

                band->stepsize = (float) (1 << 13);

                break;

            case JPEG2000_QSTY_SE:

                









                





                gain            = cbps;

                band->stepsize  = pow(2.0, gain - qntsty->expn[gbandno]);

                band->stepsize *= (float)qntsty->mant[gbandno] / 2048.0 + 1.0;

                



                band->stepsize *= 0.5;

                break;

            default:

                band->stepsize = 0;

                av_log(avctx, AV_LOG_ERROR, "Unknown quantization format\n");

                break;

            }

            

            if (avctx->flags & CODEC_FLAG_BITEXACT)

                band->stepsize = (int32_t)(band->stepsize * (1 << 16));



            







            if (reslevelno == 0) {

                

                for (i = 0; i < 2; i++)

                    for (j = 0; j < 2; j++)

                        band->coord[i][j] =

                            ff_jpeg2000_ceildivpow2(comp->coord_o[i][j],

                                                    declvl - 1);



                log2_band_prec_width  = reslevel->log2_prec_width;

                log2_band_prec_height = reslevel->log2_prec_height;

                

                band->log2_cblk_width  = FFMIN(codsty->log2_cblk_width,

                                               reslevel->log2_prec_width);

                band->log2_cblk_height = FFMIN(codsty->log2_cblk_height,

                                               reslevel->log2_prec_height);

            } else {

                

                

                for (i = 0; i < 2; i++)

                    for (j = 0; j < 2; j++)

                        

                        band->coord[i][j] =

                            ff_jpeg2000_ceildivpow2(comp->coord_o[i][j] -

                                                    (((bandno + 1 >> i) & 1) << declvl - 1),

                                                    declvl);

                





                

                band->log2_cblk_width  = FFMIN(codsty->log2_cblk_width,

                                               reslevel->log2_prec_width - 1);

                band->log2_cblk_height = FFMIN(codsty->log2_cblk_height,

                                               reslevel->log2_prec_height - 1);



                log2_band_prec_width  = reslevel->log2_prec_width  - 1;

                log2_band_prec_height = reslevel->log2_prec_height - 1;

            }



            band->prec = av_malloc_array(reslevel->num_precincts_x *

                                         reslevel->num_precincts_y,

                                         sizeof(*band->prec));

            if (!band->prec)

                return AVERROR(ENOMEM);



            nb_precincts = reslevel->num_precincts_x * reslevel->num_precincts_y;



            for (precno = 0; precno < nb_precincts; precno++) {

                Jpeg2000Prec *prec = band->prec + precno;



                

                



                

                prec->coord[0][0] = (precno % reslevel->num_precincts_x) *

                                    (1 << log2_band_prec_width);

                prec->coord[0][0] = FFMAX(prec->coord[0][0], band->coord[0][0]);



                

                prec->coord[1][0] = (precno / reslevel->num_precincts_x) *

                                    (1 << log2_band_prec_height);

                prec->coord[1][0] = FFMAX(prec->coord[1][0], band->coord[1][0]);



                

                prec->coord[0][1] = prec->coord[0][0] +

                                    (1 << log2_band_prec_width);

                prec->coord[0][1] = FFMIN(prec->coord[0][1], band->coord[0][1]);



                

                prec->coord[1][1] = prec->coord[1][0] +

                                    (1 << log2_band_prec_height);

                prec->coord[1][1] = FFMIN(prec->coord[1][1], band->coord[1][1]);



                prec->nb_codeblocks_width =

                    ff_jpeg2000_ceildivpow2(prec->coord[0][1] -

                                            prec->coord[0][0],

                                            band->log2_cblk_width);

                prec->nb_codeblocks_height =

                    ff_jpeg2000_ceildivpow2(prec->coord[1][1] -

                                            prec->coord[1][0],

                                            band->log2_cblk_height);



                

                prec->cblkincl =

                    ff_jpeg2000_tag_tree_init(prec->nb_codeblocks_width,

                                              prec->nb_codeblocks_height);

                if (!prec->cblkincl)

                    return AVERROR(ENOMEM);



                prec->zerobits =

                    ff_jpeg2000_tag_tree_init(prec->nb_codeblocks_width,

                                              prec->nb_codeblocks_height);

                if (!prec->zerobits)

                    return AVERROR(ENOMEM);



                prec->cblk = av_malloc_array(prec->nb_codeblocks_width *

                                             prec->nb_codeblocks_height,

                                             sizeof(*prec->cblk));

                if (!prec->cblk)

                    return AVERROR(ENOMEM);

                for (cblkno = 0; cblkno < prec->nb_codeblocks_width * prec->nb_codeblocks_height; cblkno++) {

                    Jpeg2000Cblk *cblk = prec->cblk + cblkno;

                    uint16_t Cx0, Cy0;



                    

                    

                    Cx0 = (prec->coord[0][0] >> band->log2_cblk_width) << band->log2_cblk_width;

                    Cx0 = Cx0 + ((cblkno % prec->nb_codeblocks_width)  << band->log2_cblk_width);

                    cblk->coord[0][0] = FFMAX(Cx0, prec->coord[0][0]);



                    

                    Cy0 = (prec->coord[1][0] >> band->log2_cblk_height) << band->log2_cblk_height;

                    Cy0 = Cy0 + ((cblkno / prec->nb_codeblocks_width)   << band->log2_cblk_height);

                    cblk->coord[1][0] = FFMAX(Cy0, prec->coord[1][0]);



                    

                    cblk->coord[0][1] = FFMIN(Cx0 + (1 << band->log2_cblk_width),

                                              prec->coord[0][1]);



                    

                    cblk->coord[1][1] = FFMIN(Cy0 + (1 << band->log2_cblk_height),

                                              prec->coord[1][1]);

                    cblk->zero      = 0;

                    cblk->lblock    = 3;

                    cblk->length    = 0;

                    cblk->lengthinc = 0;

                    cblk->npasses   = 0;

                }

            }

        }

    }

    return 0;

}

                prec->cblk = av_malloc_array(prec->nb_codeblocks_width *                                              prec->nb_codeblocks_height,                                              sizeof(*prec->cblk));
int FUNC_0(Jpeg2000Component *VAR_0,

                               Jpeg2000CodingStyle *VAR_1,

                               Jpeg2000QuantStyle *VAR_2,

                               int VAR_3, int VAR_4, int VAR_5,

                               AVCodecContext *VAR_6)

{

    uint8_t log2_band_prec_width, log2_band_prec_height;

    int VAR_7, VAR_8, VAR_9 = 0, VAR_10, VAR_11, VAR_12;

    uint32_t csize = 1;



    if (!VAR_1->nreslevels2decode) {

        av_log(VAR_6, AV_LOG_ERROR, "nreslevels2decode uninitialized\n");

        return AVERROR_INVALIDDATA;

    }



    if (VAR_10 = ff_jpeg2000_dwt_init(&VAR_0->dwt, VAR_0->coord,

                                   VAR_1->nreslevels2decode - 1,

                                   VAR_1->transform))

        return VAR_10;

    

    csize = (VAR_0->coord[0][1] - VAR_0->coord[0][0]) *

            (VAR_0->coord[1][1] - VAR_0->coord[1][0]);



    VAR_0->data = av_malloc_array(csize, sizeof(*VAR_0->data));

    if (!VAR_0->data)

        return AVERROR(ENOMEM);

    VAR_0->reslevel = av_malloc_array(VAR_1->nreslevels, sizeof(*VAR_0->reslevel));

    if (!VAR_0->reslevel)

        return AVERROR(ENOMEM);

    

    for (VAR_7 = 0; VAR_7 < VAR_1->nreslevels; VAR_7++) {

        int declvl = VAR_1->nreslevels - VAR_7;    

        Jpeg2000ResLevel *reslevel = VAR_0->reslevel + VAR_7;



        





        for (VAR_11 = 0; VAR_11 < 2; VAR_11++)

            for (VAR_12 = 0; VAR_12 < 2; VAR_12++)

                reslevel->coord[VAR_11][VAR_12] =

                    ff_jpeg2000_ceildivpow2(VAR_0->coord_o[VAR_11][VAR_12], declvl - 1);

        

        reslevel->log2_prec_width  = VAR_1->log2_prec_widths[VAR_7];

        reslevel->log2_prec_height = VAR_1->log2_prec_heights[VAR_7];



        

        if (VAR_7 == 0)

            reslevel->nbands = 1;

        else

            reslevel->nbands = 3;



        













        if (reslevel->coord[0][1] == reslevel->coord[0][0])

            reslevel->num_precincts_x = 0;

        else

            reslevel->num_precincts_x =

                ff_jpeg2000_ceildivpow2(reslevel->coord[0][1],

                                        reslevel->log2_prec_width) -

                (reslevel->coord[0][0] >> reslevel->log2_prec_width);



        if (reslevel->coord[1][1] == reslevel->coord[1][0])

            reslevel->num_precincts_y = 0;

        else

            reslevel->num_precincts_y =

                ff_jpeg2000_ceildivpow2(reslevel->coord[1][1],

                                        reslevel->log2_prec_height) -

                (reslevel->coord[1][0] >> reslevel->log2_prec_height);



        reslevel->band = av_malloc_array(reslevel->nbands, sizeof(*reslevel->band));

        if (!reslevel->band)

            return AVERROR(ENOMEM);



        for (VAR_8 = 0; VAR_8 < reslevel->nbands; VAR_8++, VAR_9++) {

            Jpeg2000Band *band = reslevel->band + VAR_8;

            int cblkno, precno;

            int nb_precincts;



            



            switch (VAR_2->quantsty) {

                uint8_t gain;

                int numbps;

            case JPEG2000_QSTY_NONE:

                

                numbps = VAR_3 +

                         lut_gain[VAR_1->transform][VAR_8 + VAR_7 > 0];

                band->stepsize = (float)SHL(2048 + VAR_2->mant[VAR_9],

                                            2 + numbps - VAR_2->expn[VAR_9]);

                break;

            case JPEG2000_QSTY_SI:

                

                band->stepsize = (float) (1 << 13);

                break;

            case JPEG2000_QSTY_SE:

                









                





                gain            = VAR_3;

                band->stepsize  = pow(2.0, gain - VAR_2->expn[VAR_9]);

                band->stepsize *= (float)VAR_2->mant[VAR_9] / 2048.0 + 1.0;

                



                band->stepsize *= 0.5;

                break;

            default:

                band->stepsize = 0;

                av_log(VAR_6, AV_LOG_ERROR, "Unknown quantization format\n");

                break;

            }

            

            if (VAR_6->flags & CODEC_FLAG_BITEXACT)

                band->stepsize = (int32_t)(band->stepsize * (1 << 16));



            







            if (VAR_7 == 0) {

                

                for (VAR_11 = 0; VAR_11 < 2; VAR_11++)

                    for (VAR_12 = 0; VAR_12 < 2; VAR_12++)

                        band->coord[VAR_11][VAR_12] =

                            ff_jpeg2000_ceildivpow2(VAR_0->coord_o[VAR_11][VAR_12],

                                                    declvl - 1);



                log2_band_prec_width  = reslevel->log2_prec_width;

                log2_band_prec_height = reslevel->log2_prec_height;

                

                band->log2_cblk_width  = FFMIN(VAR_1->log2_cblk_width,

                                               reslevel->log2_prec_width);

                band->log2_cblk_height = FFMIN(VAR_1->log2_cblk_height,

                                               reslevel->log2_prec_height);

            } else {

                

                

                for (VAR_11 = 0; VAR_11 < 2; VAR_11++)

                    for (VAR_12 = 0; VAR_12 < 2; VAR_12++)

                        

                        band->coord[VAR_11][VAR_12] =

                            ff_jpeg2000_ceildivpow2(VAR_0->coord_o[VAR_11][VAR_12] -

                                                    (((VAR_8 + 1 >> VAR_11) & 1) << declvl - 1),

                                                    declvl);

                





                

                band->log2_cblk_width  = FFMIN(VAR_1->log2_cblk_width,

                                               reslevel->log2_prec_width - 1);

                band->log2_cblk_height = FFMIN(VAR_1->log2_cblk_height,

                                               reslevel->log2_prec_height - 1);



                log2_band_prec_width  = reslevel->log2_prec_width  - 1;

                log2_band_prec_height = reslevel->log2_prec_height - 1;

            }



            band->prec = av_malloc_array(reslevel->num_precincts_x *

                                         reslevel->num_precincts_y,

                                         sizeof(*band->prec));

            if (!band->prec)

                return AVERROR(ENOMEM);



            nb_precincts = reslevel->num_precincts_x * reslevel->num_precincts_y;



            for (precno = 0; precno < nb_precincts; precno++) {

                Jpeg2000Prec *prec = band->prec + precno;



                

                



                

                prec->coord[0][0] = (precno % reslevel->num_precincts_x) *

                                    (1 << log2_band_prec_width);

                prec->coord[0][0] = FFMAX(prec->coord[0][0], band->coord[0][0]);



                

                prec->coord[1][0] = (precno / reslevel->num_precincts_x) *

                                    (1 << log2_band_prec_height);

                prec->coord[1][0] = FFMAX(prec->coord[1][0], band->coord[1][0]);



                

                prec->coord[0][1] = prec->coord[0][0] +

                                    (1 << log2_band_prec_width);

                prec->coord[0][1] = FFMIN(prec->coord[0][1], band->coord[0][1]);



                

                prec->coord[1][1] = prec->coord[1][0] +

                                    (1 << log2_band_prec_height);

                prec->coord[1][1] = FFMIN(prec->coord[1][1], band->coord[1][1]);



                prec->nb_codeblocks_width =

                    ff_jpeg2000_ceildivpow2(prec->coord[0][1] -

                                            prec->coord[0][0],

                                            band->log2_cblk_width);

                prec->nb_codeblocks_height =

                    ff_jpeg2000_ceildivpow2(prec->coord[1][1] -

                                            prec->coord[1][0],

                                            band->log2_cblk_height);



                

                prec->cblkincl =

                    ff_jpeg2000_tag_tree_init(prec->nb_codeblocks_width,

                                              prec->nb_codeblocks_height);

                if (!prec->cblkincl)

                    return AVERROR(ENOMEM);



                prec->zerobits =

                    ff_jpeg2000_tag_tree_init(prec->nb_codeblocks_width,

                                              prec->nb_codeblocks_height);

                if (!prec->zerobits)

                    return AVERROR(ENOMEM);



                prec->cblk = av_malloc_array(prec->nb_codeblocks_width *

                                             prec->nb_codeblocks_height,

                                             sizeof(*prec->cblk));

                if (!prec->cblk)

                    return AVERROR(ENOMEM);

                for (cblkno = 0; cblkno < prec->nb_codeblocks_width * prec->nb_codeblocks_height; cblkno++) {

                    Jpeg2000Cblk *cblk = prec->cblk + cblkno;

                    uint16_t Cx0, Cy0;



                    

                    

                    Cx0 = (prec->coord[0][0] >> band->log2_cblk_width) << band->log2_cblk_width;

                    Cx0 = Cx0 + ((cblkno % prec->nb_codeblocks_width)  << band->log2_cblk_width);

                    cblk->coord[0][0] = FFMAX(Cx0, prec->coord[0][0]);



                    

                    Cy0 = (prec->coord[1][0] >> band->log2_cblk_height) << band->log2_cblk_height;

                    Cy0 = Cy0 + ((cblkno / prec->nb_codeblocks_width)   << band->log2_cblk_height);

                    cblk->coord[1][0] = FFMAX(Cy0, prec->coord[1][0]);



                    

                    cblk->coord[0][1] = FFMIN(Cx0 + (1 << band->log2_cblk_width),

                                              prec->coord[0][1]);



                    

                    cblk->coord[1][1] = FFMIN(Cy0 + (1 << band->log2_cblk_height),

                                              prec->coord[1][1]);

                    cblk->zero      = 0;

                    cblk->lblock    = 3;

                    cblk->length    = 0;

                    cblk->lengthinc = 0;

                    cblk->npasses   = 0;

                }

            }

        }

    }

    return 0;

}

int FUNC_0(Jpeg2000Component *VAR_0,
Jpeg2000CodingStyle *VAR_1,
Jpeg2000QuantStyle *VAR_2,
int VAR_3, int VAR_4, int VAR_5,
AVCodecContext *VAR_6)
{ uint8_t log2_band_prec_width, log2_band_prec_height; int VAR_7, VAR_8, VAR_9 = 0, VAR_10, VAR_11, VAR_12; uint32_t csize = 1; if (!VAR_1->nreslevels2decode) { av_log(VAR_6, AV_LOG_ERROR, "nreslevels2decode uninitialized\n"); return AVERROR_INVALIDDATA; } if (VAR_10 = ff_jpeg2000_dwt_init(&VAR_0->dwt, VAR_0->coord,
VAR_1->nreslevels2decode - 1,
VAR_1->transform))
return VAR_10; csize = (VAR_0->coord[0][1] - VAR_0->coord[0][0]) *
(VAR_0->coord[1][1] - VAR_0->coord[1][0]); VAR_0->data = av_malloc_array(csize, sizeof(*VAR_0->data)); if (!VAR_0->data)
return AVERROR(ENOMEM); VAR_0->reslevel = av_malloc_array(VAR_1->nreslevels, sizeof(*VAR_0->reslevel)); if (!VAR_0->reslevel)
return AVERROR(ENOMEM); for (VAR_7 = 0; VAR_7 < VAR_1->nreslevels; VAR_7++) { int declvl = VAR_1->nreslevels - VAR_7; Jpeg2000ResLevel *reslevel = VAR_0->reslevel + VAR_7; for (VAR_11 = 0; VAR_11 < 2; VAR_11++) for (VAR_12 = 0; VAR_12 < 2; VAR_12++) reslevel->coord[VAR_11][VAR_12] =
ff_jpeg2000_ceildivpow2(VAR_0->coord_o[VAR_11][VAR_12], declvl - 1); reslevel->log2_prec_width  = VAR_1->log2_prec_widths[VAR_7]; reslevel->log2_prec_height = VAR_1->log2_prec_heights[VAR_7]; if (VAR_7 == 0)
reslevel->nbands = 1; else
reslevel->nbands = 3; if (reslevel->coord[0][1] == reslevel->coord[0][0])
reslevel->num_precincts_x = 0; else
reslevel->num_precincts_x =
ff_jpeg2000_ceildivpow2(reslevel->coord[0][1],
reslevel->log2_prec_width) -
(reslevel->coord[0][0] >> reslevel->log2_prec_width); if (reslevel->coord[1][1] == reslevel->coord[1][0])
reslevel->num_precincts_y = 0; else
reslevel->num_precincts_y =
ff_jpeg2000_ceildivpow2(reslevel->coord[1][1],
reslevel->log2_prec_height) -
(reslevel->coord[1][0] >> reslevel->log2_prec_height); reslevel->band = av_malloc_array(reslevel->nbands, sizeof(*reslevel->band)); if (!reslevel->band)
return AVERROR(ENOMEM); for (VAR_8 = 0; VAR_8 < reslevel->nbands; VAR_8++, VAR_9++) { Jpeg2000Band *band = reslevel->band + VAR_8; int cblkno, precno; int nb_precincts; switch (VAR_2->quantsty) { uint8_t gain; int numbps; case JPEG2000_QSTY_NONE:
numbps = VAR_3 +
lut_gain[VAR_1->transform][VAR_8 + VAR_7 > 0]; band->stepsize = (float)SHL(2048 + VAR_2->mant[VAR_9],
2 + numbps - VAR_2->expn[VAR_9]); break; case JPEG2000_QSTY_SI:
band->stepsize = (float) (1 << 13); break; case JPEG2000_QSTY_SE:
gain            = VAR_3; band->stepsize  = pow(2.0, gain - VAR_2->expn[VAR_9]); band->stepsize *= (float)VAR_2->mant[VAR_9] / 2048.0 + 1.0; band->stepsize *= 0.5; break; default:
band->stepsize = 0; av_log(VAR_6, AV_LOG_ERROR, "Unknown quantization format\n"); break; } if (VAR_6->flags & CODEC_FLAG_BITEXACT)
band->stepsize = (int32_t)(band->stepsize * (1 << 16)); if (VAR_7 == 0) { for (VAR_11 = 0; VAR_11 < 2; VAR_11++) for (VAR_12 = 0; VAR_12 < 2; VAR_12++) band->coord[VAR_11][VAR_12] =
ff_jpeg2000_ceildivpow2(VAR_0->coord_o[VAR_11][VAR_12],
declvl - 1); log2_band_prec_width  = reslevel->log2_prec_width; log2_band_prec_height = reslevel->log2_prec_height; band->log2_cblk_width  = FFMIN(VAR_1->log2_cblk_width,
reslevel->log2_prec_width); band->log2_cblk_height = FFMIN(VAR_1->log2_cblk_height,
reslevel->log2_prec_height); } else { for (VAR_11 = 0; VAR_11 < 2; VAR_11++) for (VAR_12 = 0; VAR_12 < 2; VAR_12++) band->coord[VAR_11][VAR_12] =
ff_jpeg2000_ceildivpow2(VAR_0->coord_o[VAR_11][VAR_12] -
(((VAR_8 + 1 >> VAR_11) & 1) << declvl - 1),
declvl); band->log2_cblk_width  = FFMIN(VAR_1->log2_cblk_width,
reslevel->log2_prec_width - 1); band->log2_cblk_height = FFMIN(VAR_1->log2_cblk_height,
reslevel->log2_prec_height - 1); log2_band_prec_width  = reslevel->log2_prec_width  - 1; log2_band_prec_height = reslevel->log2_prec_height - 1; } band->prec = av_malloc_array(reslevel->num_precincts_x *
reslevel->num_precincts_y,
sizeof(*band->prec)); if (!band->prec)
return AVERROR(ENOMEM); nb_precincts = reslevel->num_precincts_x * reslevel->num_precincts_y; for (precno = 0; precno < nb_precincts; precno++) { Jpeg2000Prec *prec = band->prec + precno; prec->coord[0][0] = (precno % reslevel->num_precincts_x) *
(1 << log2_band_prec_width); prec->coord[0][0] = FFMAX(prec->coord[0][0], band->coord[0][0]); prec->coord[1][0] = (precno / reslevel->num_precincts_x) *
(1 << log2_band_prec_height); prec->coord[1][0] = FFMAX(prec->coord[1][0], band->coord[1][0]); prec->coord[0][1] = prec->coord[0][0] +
(1 << log2_band_prec_width); prec->coord[0][1] = FFMIN(prec->coord[0][1], band->coord[0][1]); prec->coord[1][1] = prec->coord[1][0] +
(1 << log2_band_prec_height); prec->coord[1][1] = FFMIN(prec->coord[1][1], band->coord[1][1]); prec->nb_codeblocks_width =
ff_jpeg2000_ceildivpow2(prec->coord[0][1] -
prec->coord[0][0],
band->log2_cblk_width); prec->nb_codeblocks_height =
ff_jpeg2000_ceildivpow2(prec->coord[1][1] -
prec->coord[1][0],
band->log2_cblk_height); prec->cblkincl =
ff_jpeg2000_tag_tree_init(prec->nb_codeblocks_width,
prec->nb_codeblocks_height); if (!prec->cblkincl)
return AVERROR(ENOMEM); prec->zerobits =
ff_jpeg2000_tag_tree_init(prec->nb_codeblocks_width,
prec->nb_codeblocks_height); if (!prec->zerobits)
return AVERROR(ENOMEM); prec->cblk = av_malloc_array(prec->nb_codeblocks_width *
prec->nb_codeblocks_height,
sizeof(*prec->cblk)); if (!prec->cblk)
return AVERROR(ENOMEM); for (cblkno = 0; cblkno < prec->nb_codeblocks_width * prec->nb_codeblocks_height; cblkno++) { Jpeg2000Cblk *cblk = prec->cblk + cblkno; uint16_t Cx0, Cy0; Cx0 = (prec->coord[0][0] >> band->log2_cblk_width) << band->log2_cblk_width; Cx0 = Cx0 + ((cblkno % prec->nb_codeblocks_width)  << band->log2_cblk_width); cblk->coord[0][0] = FFMAX(Cx0, prec->coord[0][0]); Cy0 = (prec->coord[1][0] >> band->log2_cblk_height) << band->log2_cblk_height; Cy0 = Cy0 + ((cblkno / prec->nb_codeblocks_width)   << band->log2_cblk_height); cblk->coord[1][0] = FFMAX(Cy0, prec->coord[1][0]); cblk->coord[0][1] = FFMIN(Cx0 + (1 << band->log2_cblk_width),
prec->coord[0][1]); cblk->coord[1][1] = FFMIN(Cy0 + (1 << band->log2_cblk_height),
prec->coord[1][1]); cblk->zero      = 0; cblk->lblock    = 3; cblk->length    = 0; cblk->lengthinc = 0; cblk->npasses   = 0; } } } } return 0; }