+#define WRITE_FRAME(D, PIXEL) \
+ static inline void write_frame_ ## D(Jpeg2000DecoderContext * s, Jpeg2000Tile * tile, \
+ AVFrame * picture, int precision) \
+ { \
+ const AVPixFmtDescriptor *pixdesc = av_pix_fmt_desc_get(s->avctx->pix_fmt); \
+ int planar = !!(pixdesc->flags & AV_PIX_FMT_FLAG_PLANAR); \
+ int pixelsize = planar ? 1 : pixdesc->nb_components; \
+ \
+ int compno; \
+ int x, y; \
+ \
+ for (compno = 0; compno < s->ncomponents; compno++) { \
+ Jpeg2000Component *comp = tile->comp + compno; \
+ Jpeg2000CodingStyle *codsty = tile->codsty + compno; \
+ PIXEL *line; \
+ float *datap = comp->f_data; \
+ int32_t *i_datap = comp->i_data; \
+ int cbps = s->cbps[compno]; \
+ int w = tile->comp[compno].coord[0][1] - s->image_offset_x; \
+ int plane = 0; \
+ \
+ if (planar) \
+ plane = s->cdef[compno] ? s->cdef[compno]-1 : (s->ncomponents-1); \
+ \
+ y = tile->comp[compno].coord[1][0] - s->image_offset_y / s->cdy[compno]; \
+ line = (PIXEL *)picture->data[plane] + y * (picture->linesize[plane] / sizeof(PIXEL));\
+ for (; y < tile->comp[compno].coord[1][1] - s->image_offset_y; y++) { \
+ PIXEL *dst; \
+ \
+ x = tile->comp[compno].coord[0][0] - s->image_offset_x / s->cdx[compno]; \
+ dst = line + x * pixelsize + compno*!planar; \
+ \
+ if (codsty->transform == FF_DWT97) { \
+ for (; x < w; x++) { \
+ int val = lrintf(*datap) + (1 << (cbps - 1)); \
+ /* DC level shift and clip see ISO 15444-1:2002 G.1.2 */ \
+ val = av_clip(val, 0, (1 << cbps) - 1); \
+ *dst = val << (precision - cbps); \
+ datap++; \
+ dst += pixelsize; \
+ } \
+ } else { \
+ for (; x < w; x++) { \
+ int val = *i_datap + (1 << (cbps - 1)); \
+ /* DC level shift and clip see ISO 15444-1:2002 G.1.2 */ \
+ val = av_clip(val, 0, (1 << cbps) - 1); \
+ *dst = val << (precision - cbps); \
+ i_datap++; \
+ dst += pixelsize; \
+ } \
+ } \
+ line += picture->linesize[plane] / sizeof(PIXEL); \
+ } \
+ } \
+ \
+ }
+
+WRITE_FRAME(8, uint8_t)
+WRITE_FRAME(16, uint16_t)
+
+#undef WRITE_FRAME
+