}	dr = jas_matrix_getref(data, 0, 0);	drs = jas_matrix_rowstep(data);	for (i = 0; i < height; ++i, dr += drs) {		d = dr;		if (jas_stream_seek(cmpt->stream_, (cmpt->width_ * (y + i) + x)		  * cmpt->cps_, SEEK_SET) < 0) {			return -1;		}		for (j = width; j > 0; --j, ++d) {			v = inttobits(*d, cmpt->prec_, cmpt->sgnd_);			for (k = cmpt->cps_; k > 0; --k) {				c = (v >> (8 * (cmpt->cps_ - 1))) & 0xff;				if (jas_stream_putc(cmpt->stream_,				  (unsigned char) c) == EOF) {					return -1;				}				v <<= 8;			}		}	}	return 0;}/******************************************************************************\* File format operations.\******************************************************************************/void jas_image_clearfmts(){	int i;	jas_image_fmtinfo_t *fmtinfo;	for (i = 0; i < jas_image_numfmts; ++i) {		fmtinfo = &jas_image_fmtinfos[i];		if (fmtinfo->name) {			jas_free(fmtinfo->name);			fmtinfo->name = 0;		}		if (fmtinfo->ext) {			jas_free(fmtinfo->ext);			fmtinfo->ext = 0;		}		if (fmtinfo->desc) {			jas_free(fmtinfo->desc);			fmtinfo->desc = 0;		}	}	jas_image_numfmts = 0;}int jas_image_addfmt(int id, char *name, char *ext, char *desc,  jas_image_fmtops_t *ops){	jas_image_fmtinfo_t *fmtinfo;	assert(id >= 0 && name && ext && ops);	if (jas_image_numfmts >= JAS_IMAGE_MAXFMTS) {		return -1;	}	fmtinfo = &jas_image_fmtinfos[jas_image_numfmts];	fmtinfo->id = id;	if (!(fmtinfo->name = jas_strdup(name))) {		return -1;	}	if (!(fmtinfo->ext = jas_strdup(ext))) {		jas_free(fmtinfo->name);		return -1;	}	if (!(fmtinfo->desc = jas_strdup(desc))) {		jas_free(fmtinfo->name);		jas_free(fmtinfo->ext);		return -1;	}	fmtinfo->ops = *ops;	++jas_image_numfmts;	return 0;}int jas_image_strtofmt(char *name){	jas_image_fmtinfo_t *fmtinfo;	if (!(fmtinfo = jas_image_lookupfmtbyname(name))) {		return -1;	}	return fmtinfo->id;}char *jas_image_fmttostr(int fmt){	jas_image_fmtinfo_t *fmtinfo;	if (!(fmtinfo = jas_image_lookupfmtbyid(fmt))) {		return 0;	}	return fmtinfo->name;}int jas_image_getfmt(jas_stream_t *in){	jas_image_fmtinfo_t *fmtinfo;	int found;	int i;	/* Check for data in each of the supported formats. */	found = 0;	for (i = 0, fmtinfo = jas_image_fmtinfos; i < jas_image_numfmts; ++i,	  ++fmtinfo) {		if (fmtinfo->ops.validate) {			/* Is the input data valid for this format? */			if (!(*fmtinfo->ops.validate)(in)) {				found = 1;				break;			}		}	}	return found ? fmtinfo->id : (-1);}int jas_image_fmtfromname(char *name){	int i;	char *ext;	jas_image_fmtinfo_t *fmtinfo;	/* Get the file name extension. */	if (!(ext = strrchr(name, '.'))) {		return -1;	}	++ext;	/* Try to find a format that uses this extension. */		for (i = 0, fmtinfo = jas_image_fmtinfos; i < jas_image_numfmts; ++i,	  ++fmtinfo) {		/* Do we have a match? */		if (!strcmp(ext, fmtinfo->ext)) {			return fmtinfo->id;		}	}	return -1;}/******************************************************************************\* Miscellaneous operations.\******************************************************************************/uint_fast32_t jas_image_rawsize(jas_image_t *image){	uint_fast32_t rawsize;	int cmptno;	jas_image_cmpt_t *cmpt;	rawsize = 0;	for (cmptno = 0; cmptno < image->numcmpts_; ++cmptno) {		cmpt = image->cmpts_[cmptno];		rawsize += (cmpt->width_ * cmpt->height_ * cmpt->prec_ +		  7) / 8;	}	return rawsize;}void jas_image_delcmpt(jas_image_t *image, int cmptno){	if (cmptno >= image->numcmpts_) {		return;	}	jas_image_cmpt_destroy(image->cmpts_[cmptno]);	if (cmptno < image->numcmpts_) {		memmove(&image->cmpts_[cmptno], &image->cmpts_[cmptno + 1],		  (image->numcmpts_ - 1 - cmptno) * sizeof(jas_image_cmpt_t *));	}	--image->numcmpts_;	jas_image_setbbox(image);}int jas_image_addcmpt(jas_image_t *image, int cmptno,  jas_image_cmptparm_t *cmptparm){	jas_image_cmpt_t *newcmpt;	if (cmptno < 0)		cmptno = image->numcmpts_;	assert(cmptno >= 0 && cmptno <= image->numcmpts_);	if (image->numcmpts_ >= image->maxcmpts_) {		if (jas_image_growcmpts(image, image->maxcmpts_ + 128)) {			return -1;		}	}	if (!(newcmpt = jas_image_cmpt_create(cmptparm->tlx,	  cmptparm->tly, cmptparm->hstep, cmptparm->vstep,	  cmptparm->width, cmptparm->height, cmptparm->prec,	  cmptparm->sgnd, 1))) {		return -1;	}	if (cmptno < image->numcmpts_) {		memmove(&image->cmpts_[cmptno + 1], &image->cmpts_[cmptno],		  (image->numcmpts_ - cmptno) * sizeof(jas_image_cmpt_t *));	}	image->cmpts_[cmptno] = newcmpt;	++image->numcmpts_;	jas_image_setbbox(image);	return 0;}jas_image_fmtinfo_t *jas_image_lookupfmtbyid(int id){	int i;	jas_image_fmtinfo_t *fmtinfo;	for (i = 0, fmtinfo = jas_image_fmtinfos; i < jas_image_numfmts; ++i, ++fmtinfo) {		if (fmtinfo->id == id) {			return fmtinfo;		}	}	return 0;}jas_image_fmtinfo_t *jas_image_lookupfmtbyname(const char *name){	int i;	jas_image_fmtinfo_t *fmtinfo;	for (i = 0, fmtinfo = jas_image_fmtinfos; i < jas_image_numfmts; ++i, ++fmtinfo) {		if (!strcmp(fmtinfo->name, name)) {			return fmtinfo;		}	}	return 0;}static uint_fast32_t inttobits(jas_seqent_t v, int prec, bool sgnd){	uint_fast32_t ret;	ret = ((sgnd && v < 0) ? ((1 << prec) + v) : v) & JAS_ONES(prec);	return ret;}static jas_seqent_t bitstoint(uint_fast32_t v, int prec, bool sgnd){	jas_seqent_t ret;	v &= JAS_ONES(prec);	ret = (sgnd && (v & (1 << (prec - 1)))) ? (v - (1 << prec)) : v;	return ret;}static void jas_image_setbbox(jas_image_t *image){	jas_image_cmpt_t *cmpt;	int cmptno;	int_fast32_t x;	int_fast32_t y;	if (image->numcmpts_ > 0) {		/* Determine the bounding box for all of the components on the		  reference grid (i.e., the image area) */		cmpt = image->cmpts_[0];		image->tlx_ = cmpt->tlx_;		image->tly_ = cmpt->tly_;		image->brx_ = cmpt->tlx_ + cmpt->hstep_ * (cmpt->width_ - 1) + 1;		image->bry_ = cmpt->tly_ + cmpt->vstep_ * (cmpt->height_ - 1) + 1;		for (cmptno = 1; cmptno < image->numcmpts_; ++cmptno) {			cmpt = image->cmpts_[cmptno];			if (image->tlx_ > cmpt->tlx_) {				image->tlx_ = cmpt->tlx_;			}			if (image->tly_ > cmpt->tly_) {				image->tly_ = cmpt->tly_;			}			x = cmpt->tlx_ + cmpt->hstep_ * (cmpt->width_ - 1) + 1;			if (image->brx_ < x) {				image->brx_ = x;			}			y = cmpt->tly_ + cmpt->vstep_ * (cmpt->height_ - 1) + 1;			if (image->bry_ < y) {				image->bry_ = y;			}		}	} else {		image->tlx_ = 0;		image->tly_ = 0;		image->brx_ = 0;		image->bry_ = 0;	}}static int jas_image_growcmpts(jas_image_t *image, int maxcmpts){	jas_image_cmpt_t **newcmpts;	int cmptno;	newcmpts = (!image->cmpts_) ? jas_malloc(maxcmpts * sizeof(jas_image_cmpt_t *)) :	  jas_realloc(image->cmpts_, maxcmpts * sizeof(jas_image_cmpt_t *));	if (!newcmpts) {		return -1;	}	image->cmpts_ = newcmpts;	image->maxcmpts_ = maxcmpts;	for (cmptno = image->numcmpts_; cmptno < image->maxcmpts_; ++cmptno) {		image->cmpts_[cmptno] = 0;	}	return 0;}int jas_image_copycmpt(jas_image_t *dstimage, int dstcmptno, jas_image_t *srcimage,  int srccmptno){	jas_image_cmpt_t *newcmpt;	if (dstimage->numcmpts_ >= dstimage->maxcmpts_) {		if (jas_image_growcmpts(dstimage, dstimage->maxcmpts_ + 128)) {			return -1;		}	}	if (!(newcmpt = jas_image_cmpt_copy(srcimage->cmpts_[srccmptno]))) {		return -1;	}	if (dstcmptno < dstimage->numcmpts_) {		memmove(&dstimage->cmpts_[dstcmptno + 1], &dstimage->cmpts_[dstcmptno],		  (dstimage->numcmpts_ - dstcmptno) * sizeof(jas_image_cmpt_t *));	}	dstimage->cmpts_[dstcmptno] = newcmpt;	++dstimage->numcmpts_;	jas_image_setbbox(dstimage);	return 0;}void jas_image_dump(jas_image_t *image, FILE *out){	long buf[1024];	int cmptno;	int n;	int i;	int width;	int height;	jas_image_cmpt_t *cmpt;	for (cmptno = 0; cmptno < image->numcmpts_; ++cmptno) {		cmpt = image->cmpts_[cmptno];		fprintf(out, "prec=%d, sgnd=%d, cmpttype=%d\n", cmpt->prec_,		  cmpt->sgnd_, cmpt->type_);		width = jas_image_cmptwidth(image, cmptno);		height = jas_image_cmptheight(image, cmptno);		n = JAS_MIN(16, width);		if (jas_image_readcmpt2(image, cmptno, 0, 0, n, 1, buf)) {			abort();		}		for (i = 0; i < n; ++i) {			fprintf(out, " f(%d,%d)=%ld", i, 0, buf[i]);		}		fprintf(out, "\n");		if (jas_image_readcmpt2(image, cmptno, width - n, height - 1, n, 1, buf)) {			abort();		}		for (i = 0; i < n; ++i) {			fprintf(out, " f(%d,%d)=%ld", width - n + i, height - 1, buf[i]);		}		fprintf(out, "\n");	}}int jas_image_depalettize(jas_image_t *image, int cmptno, int numlutents,  int_fast32_t *lutents, int dtype, int newcmptno){	jas_image_cmptparm_t cmptparms;	int_fast32_t v;	int i;	int j;	jas_image_cmpt_t *cmpt;	cmpt = image->cmpts_[cmptno];	cmptparms.tlx = cmpt->tlx_;	cmptparms.tly = cmpt->tly_;	cmptparms.hstep = cmpt->hstep_;	cmptparms.vstep = cmpt->vstep_;	cmptparms.width = cmpt->width_;	cmptparms.height = cmpt->height_;	cmptparms.prec = JAS_IMAGE_CDT_GETPREC(dtype);	cmptparms.sgnd = JAS_IMAGE_CDT_GETSGND(dtype);	if (jas_image_addcmpt(image, newcmptno, &cmptparms)) {		return -1;	}	if (newcmptno <= cmptno) {		++cmptno;		cmpt = image->cmpts_[cmptno];	}	for (j = 0; j < cmpt->height_; ++j) {		for (i = 0; i < cmpt->width_; ++i) {			v = jas_image_readcmptsample(image, cmptno, i, j);			if (v < 0) {				v = 0;			} else if (v >= numlutents) {				v = numlutents - 1;			}			jas_image_writecmptsample(image, newcmptno, i, j,			  lutents[v]);		}	}	return 0;}int jas_image_readcmptsample(jas_image_t *image, int cmptno, int x, int y){	jas_image_cmpt_t *cmpt;	uint_fast32_t v;	int k;	int c;	cmpt = image->cmpts_[cmptno];	if (jas_stream_seek(cmpt->stream_, (cmpt->width_ * y + x) * cmpt->cps_,	  SEEK_SET) < 0) {		return -1;	}	v = 0;	for (k = cmpt->cps_; k > 0; --k) {		if ((c = jas_stream_getc(cmpt->stream_)) == EOF) {			return -1;		}		v = (v << 8) | (c & 0xff);	}	return bitstoint(v, cmpt->prec_, cmpt->sgnd_);}void jas_image_writecmptsample(jas_image_t *image, int cmptno, int x, int y,  int_fast32_t v){	jas_image_cmpt_t *cmpt;	uint_fast32_t t;	int k;	int c;	cmpt = image->cmpts_[cmptno];	if (jas_stream_seek(cmpt->stream_, (cmpt->width_ * y + x) * cmpt->cps_,	  SEEK_SET) < 0) {		return;	}	t = inttobits(v, cmpt->prec_, cmpt->sgnd_);	for (k = cmpt->cps_; k > 0; --k) {		c = (t >> (8 * (cmpt->cps_ - 1))) & 0xff;		if (jas_stream_putc(cmpt->stream_, (unsigned char) c) == EOF) {			return;		}		t <<= 8;	}}int jas_image_getcmptbytype(jas_image_t *image, int ctype){	int cmptno;	for (cmptno = 0; cmptno < image->numcmpts_; ++cmptno) {		if (image->cmpts_[cmptno]->type_ == ctype) {			return cmptno;		}	}	return -1;}/***********************************************//***********************************************//***********************************************//***********************************************/