2 * Copyright (c) 1999-2000 Image Power, Inc. and the University of
4 * Copyright (c) 2001-2003 Michael David Adams.
8 /* __START_OF_JASPER_LICENSE__
10 * JasPer License Version 2.0
12 * Copyright (c) 1999-2000 Image Power, Inc.
13 * Copyright (c) 1999-2000 The University of British Columbia
14 * Copyright (c) 2001-2003 Michael David Adams
16 * All rights reserved.
18 * Permission is hereby granted, free of charge, to any person (the
19 * "User") obtaining a copy of this software and associated documentation
20 * files (the "Software"), to deal in the Software without restriction,
21 * including without limitation the rights to use, copy, modify, merge,
22 * publish, distribute, and/or sell copies of the Software, and to permit
23 * persons to whom the Software is furnished to do so, subject to the
24 * following conditions:
26 * 1. The above copyright notices and this permission notice (which
27 * includes the disclaimer below) shall be included in all copies or
28 * substantial portions of the Software.
30 * 2. The name of a copyright holder shall not be used to endorse or
31 * promote products derived from the Software without specific prior
34 * THIS DISCLAIMER OF WARRANTY CONSTITUTES AN ESSENTIAL PART OF THIS
35 * LICENSE. NO USE OF THE SOFTWARE IS AUTHORIZED HEREUNDER EXCEPT UNDER
36 * THIS DISCLAIMER. THE SOFTWARE IS PROVIDED BY THE COPYRIGHT HOLDERS
37 * "AS IS", WITHOUT WARRANTY OF ANY KIND, EXPRESS OR IMPLIED, INCLUDING
38 * BUT NOT LIMITED TO THE WARRANTIES OF MERCHANTABILITY, FITNESS FOR A
39 * PARTICULAR PURPOSE AND NONINFRINGEMENT OF THIRD PARTY RIGHTS. IN NO
40 * EVENT SHALL THE COPYRIGHT HOLDERS BE LIABLE FOR ANY CLAIM, OR ANY SPECIAL
41 * INDIRECT OR CONSEQUENTIAL DAMAGES, OR ANY DAMAGES WHATSOEVER RESULTING
42 * FROM LOSS OF USE, DATA OR PROFITS, WHETHER IN AN ACTION OF CONTRACT,
43 * NEGLIGENCE OR OTHER TORTIOUS ACTION, ARISING OUT OF OR IN CONNECTION
44 * WITH THE USE OR PERFORMANCE OF THIS SOFTWARE. NO ASSURANCES ARE
45 * PROVIDED BY THE COPYRIGHT HOLDERS THAT THE SOFTWARE DOES NOT INFRINGE
46 * THE PATENT OR OTHER INTELLECTUAL PROPERTY RIGHTS OF ANY OTHER ENTITY.
47 * EACH COPYRIGHT HOLDER DISCLAIMS ANY LIABILITY TO THE USER FOR CLAIMS
48 * BROUGHT BY ANY OTHER ENTITY BASED ON INFRINGEMENT OF INTELLECTUAL
49 * PROPERTY RIGHTS OR OTHERWISE. AS A CONDITION TO EXERCISING THE RIGHTS
50 * GRANTED HEREUNDER, EACH USER HEREBY ASSUMES SOLE RESPONSIBILITY TO SECURE
51 * ANY OTHER INTELLECTUAL PROPERTY RIGHTS NEEDED, IF ANY. THE SOFTWARE
52 * IS NOT FAULT-TOLERANT AND IS NOT INTENDED FOR USE IN MISSION-CRITICAL
53 * SYSTEMS, SUCH AS THOSE USED IN THE OPERATION OF NUCLEAR FACILITIES,
54 * AIRCRAFT NAVIGATION OR COMMUNICATION SYSTEMS, AIR TRAFFIC CONTROL
55 * SYSTEMS, DIRECT LIFE SUPPORT MACHINES, OR WEAPONS SYSTEMS, IN WHICH
56 * THE FAILURE OF THE SOFTWARE OR SYSTEM COULD LEAD DIRECTLY TO DEATH,
57 * PERSONAL INJURY, OR SEVERE PHYSICAL OR ENVIRONMENTAL DAMAGE ("HIGH
58 * RISK ACTIVITIES"). THE COPYRIGHT HOLDERS SPECIFICALLY DISCLAIM ANY
59 * EXPRESS OR IMPLIED WARRANTY OF FITNESS FOR HIGH RISK ACTIVITIES.
61 * __END_OF_JASPER_LICENSE__
65 * $Id: jpc_t1cod.c,v 1.1 2005/05/22 18:33:05 malaterre Exp $
68 /******************************************************************************\
70 \******************************************************************************/
77 #include "jasper/jas_types.h"
78 #include "jasper/jas_math.h"
83 #include "jpc_mqcod.h"
84 #include "jpc_t1cod.h"
87 double jpc_pow2i(int n);
89 /******************************************************************************\
91 \******************************************************************************/
93 int jpc_zcctxnolut[4 * 256];
95 int jpc_scctxnolut[256];
96 int jpc_magctxnolut[4096];
98 jpc_fix_t jpc_signmsedec[1 << JPC_NMSEDEC_BITS];
99 jpc_fix_t jpc_refnmsedec[1 << JPC_NMSEDEC_BITS];
100 jpc_fix_t jpc_signmsedec0[1 << JPC_NMSEDEC_BITS];
101 jpc_fix_t jpc_refnmsedec0[1 << JPC_NMSEDEC_BITS];
103 jpc_mqctx_t jpc_mqctxs[JPC_NUMCTXS];
105 /******************************************************************************\
107 \******************************************************************************/
109 void jpc_initmqctxs(void);
111 /******************************************************************************\
113 \******************************************************************************/
115 int JPC_PASSTYPE(int passno)
118 switch (passno % 3) {
120 passtype = JPC_CLNPASS;
123 passtype = JPC_SIGPASS;
126 passtype = JPC_REFPASS;
136 int JPC_NOMINALGAIN(int qmfbid, int numlvls, int lvlno, int orient)
138 /* Avoid compiler warnings about unused parameters. */
141 if (qmfbid == JPC_COX_INS) {
144 assert(qmfbid == JPC_COX_RFT);
146 assert(orient == JPC_TSFB_LL);
162 /******************************************************************************\
163 * Coding pass related functions.
164 \******************************************************************************/
166 int JPC_SEGTYPE(int passno, int firstpassno, int bypass)
170 passtype = JPC_PASSTYPE(passno);
171 if (passtype == JPC_CLNPASS) {
174 return ((passno < firstpassno + 10) ? JPC_SEG_MQ : JPC_SEG_RAW);
180 int JPC_SEGPASSCNT(int passno, int firstpassno, int numpasses, int bypass, int termall)
188 if (passno < firstpassno + 10) {
189 ret = 10 - (passno - firstpassno);
191 passtype = JPC_PASSTYPE(passno);
209 ret = JPC_PREC * 3 - 2;
211 ret = JAS_MIN(ret, numpasses - passno);
215 int JPC_ISTERMINATED(int passno, int firstpassno, int numpasses, int termall,
220 if (passno - firstpassno == numpasses - 1) {
223 n = JPC_SEGPASSCNT(passno, firstpassno, numpasses, lazy, termall);
224 ret = (n <= 1) ? 1 : 0;
230 /******************************************************************************\
232 \******************************************************************************/
246 for (orient = 0; orient < 4; ++orient) {
247 for (i = 0; i < 256; ++i) {
248 jpc_zcctxnolut[(orient << 8) | i] = jpc_getzcctxno(i, orient);
252 for (i = 0; i < 256; ++i) {
253 jpc_spblut[i] = jpc_getspb(i << 4);
256 for (i = 0; i < 256; ++i) {
257 jpc_scctxnolut[i] = jpc_getscctxno(i << 4);
260 for (refine = 0; refine < 2; ++refine) {
261 for (i = 0; i < 2048; ++i) {
262 jpc_magctxnolut[(refine << 11) + i] = jpc_getmagctxno((refine ? JPC_REFINE : 0) | i);
266 for (i = 0; i < (1 << JPC_NMSEDEC_BITS); ++i) {
267 t = i * jpc_pow2i(-JPC_NMSEDEC_FRACBITS);
270 jpc_signmsedec[i] = jpc_dbltofix(floor((u * u - v * v) * jpc_pow2i(JPC_NMSEDEC_FRACBITS) + 0.5) / jpc_pow2i(JPC_NMSEDEC_FRACBITS));
271 /* XXX - this calc is not correct */
272 jpc_signmsedec0[i] = jpc_dbltofix(floor((u * u) * jpc_pow2i(JPC_NMSEDEC_FRACBITS) + 0.5) / jpc_pow2i(JPC_NMSEDEC_FRACBITS));
274 if (i & (1 << (JPC_NMSEDEC_BITS - 1))) {
279 jpc_refnmsedec[i] = jpc_dbltofix(floor((u * u - v * v) * jpc_pow2i(JPC_NMSEDEC_FRACBITS) + 0.5) / jpc_pow2i(JPC_NMSEDEC_FRACBITS));
280 /* XXX - this calc is not correct */
281 jpc_refnmsedec0[i] = jpc_dbltofix(floor((u * u) * jpc_pow2i(JPC_NMSEDEC_FRACBITS) + 0.5) / jpc_pow2i(JPC_NMSEDEC_FRACBITS));
285 jpc_fix_t jpc_getsignmsedec_func(jpc_fix_t x, int bitpos)
288 assert(!(x & (~JAS_ONES(bitpos + 1))));
289 y = jpc_getsignmsedec_macro(x, bitpos);
293 int jpc_getzcctxno(int f, int orient)
302 /* Avoid compiler warning. */
305 h = ((f & JPC_WSIG) != 0) + ((f & JPC_ESIG) != 0);
306 v = ((f & JPC_NSIG) != 0) + ((f & JPC_SSIG) != 0);
307 d = ((f & JPC_NWSIG) != 0) + ((f & JPC_NESIG) != 0) + ((f & JPC_SESIG) != 0) + ((f & JPC_SWSIG) != 0);
348 } else if (hv == 1) {
356 } else if (hv == 1) {
372 assert(n < JPC_NUMZCCTXS);
373 return JPC_ZCCTXNO + n;
376 int jpc_getspb(int f)
382 hc = JAS_MIN(((f & (JPC_ESIG | JPC_ESGN)) == JPC_ESIG) + ((f & (JPC_WSIG | JPC_WSGN)) == JPC_WSIG), 1) -
383 JAS_MIN(((f & (JPC_ESIG | JPC_ESGN)) == (JPC_ESIG | JPC_ESGN)) + ((f & (JPC_WSIG | JPC_WSGN)) == (JPC_WSIG | JPC_WSGN)), 1);
384 vc = JAS_MIN(((f & (JPC_NSIG | JPC_NSGN)) == JPC_NSIG) + ((f & (JPC_SSIG | JPC_SSGN)) == JPC_SSIG), 1) -
385 JAS_MIN(((f & (JPC_NSIG | JPC_NSGN)) == (JPC_NSIG | JPC_NSGN)) + ((f & (JPC_SSIG | JPC_SSGN)) == (JPC_SSIG | JPC_SSGN)), 1);
389 n = (!(hc > 0 || (!hc && vc > 0)));
394 int jpc_getscctxno(int f)
400 /* Avoid compiler warning. */
403 hc = JAS_MIN(((f & (JPC_ESIG | JPC_ESGN)) == JPC_ESIG) + ((f & (JPC_WSIG | JPC_WSGN)) == JPC_WSIG),
404 1) - JAS_MIN(((f & (JPC_ESIG | JPC_ESGN)) == (JPC_ESIG | JPC_ESGN)) +
405 ((f & (JPC_WSIG | JPC_WSGN)) == (JPC_WSIG | JPC_WSGN)), 1);
406 vc = JAS_MIN(((f & (JPC_NSIG | JPC_NSGN)) == JPC_NSIG) + ((f & (JPC_SSIG | JPC_SSGN)) == JPC_SSIG),
407 1) - JAS_MIN(((f & (JPC_NSIG | JPC_NSGN)) == (JPC_NSIG | JPC_NSGN)) +
408 ((f & (JPC_SSIG | JPC_SSGN)) == (JPC_SSIG | JPC_SSGN)), 1);
409 assert(hc >= -1 && hc <= 1 && vc >= -1 && vc <= 1);
422 } else if (hc == 1) {
431 assert(n < JPC_NUMSCCTXS);
432 return JPC_SCCTXNO + n;
435 int jpc_getmagctxno(int f)
439 if (!(f & JPC_REFINE)) {
440 n = (f & (JPC_OTHSIGMSK)) ? 1 : 0;
445 assert(n < JPC_NUMMAGCTXS);
446 return JPC_MAGCTXNO + n;
449 void jpc_initctxs(jpc_mqctx_t *ctxs)
455 for (i = 0; i < JPC_NUMCTXS; ++i) {
475 void jpc_initmqctxs()
477 jpc_initctxs(jpc_mqctxs);
480 /* Calculate the real quantity exp2(n), where x is an integer. */
481 double jpc_pow2i(int n)