]> Creatis software - gdcm.git/blob - src/gdcmjpegls/Encoder/bitio.c
Reverse order sorting now works, even with user supplied function.
[gdcm.git] / src / gdcmjpegls / Encoder / bitio.c
1 /* SPMG/JPEG-LS IMPLEMENTATION V.2.1
2    =====================================
3    These programs are Copyright (c) University of British Columbia. All rights reserved.
4    They may be freely redistributed in their entirety provided that this copyright
5    notice is not removed. THEY MAY NOT BE SOLD FOR PROFIT OR INCORPORATED IN
6    COMMERCIAL PROGRAMS WITHOUT THE WRITTEN PERMISSION OF THE COPYRIGHT HOLDER.
7    Each program is provided as is, without any express or implied warranty,
8    without even the warranty of fitness for a particular purpose.
9
10    =========================================================
11    THIS SOFTWARE IS BASED ON HP's implementation of jpeg-ls:
12    =========================================================
13
14    LOCO-I/JPEG-LS IMPLEMENTATION V.0.90
15    -------------------------------------------------------------------------------
16    (c) COPYRIGHT HEWLETT-PACKARD COMPANY, 1995-1999.
17         HEWLETT-PACKARD COMPANY ("HP") DOES NOT WARRANT THE ACCURACY OR
18    COMPLETENESS OF THE INFORMATION GIVEN HERE.  ANY USE MADE OF, OR
19    RELIANCE ON, SUCH INFORMATION IS ENTIRELY AT USER'S OWN RISK.
20         BY DOWNLOADING THE LOCO-I/JPEG-LS COMPRESSORS/DECOMPRESSORS
21    ("THE SOFTWARE") YOU AGREE TO BE BOUND BY THE TERMS AND CONDITIONS
22    OF THIS LICENSING AGREEMENT.
23         YOU MAY DOWNLOAD AND USE THE SOFTWARE FOR NON-COMMERCIAL PURPOSES
24    FREE OF CHARGE OR FURTHER OBLIGATION.  YOU MAY NOT, DIRECTLY OR
25    INDIRECTLY, DISTRIBUTE THE SOFTWARE FOR A FEE, INCORPORATE THIS
26    SOFTWARE INTO ANY PRODUCT OFFERED FOR SALE, OR USE THE SOFTWARE
27    TO PROVIDE A SERVICE FOR WHICH A FEE IS CHARGED.
28         YOU MAY MAKE COPIES OF THE SOFTWARE AND DISTRIBUTE SUCH COPIES TO
29    OTHER PERSONS PROVIDED THAT SUCH COPIES ARE ACCOMPANIED BY
30    HEWLETT-PACKARD'S COPYRIGHT NOTICE AND THIS AGREEMENT AND THAT
31    SUCH OTHER PERSONS AGREE TO BE BOUND BY THE TERMS OF THIS AGREEMENT.
32         THE SOFTWARE IS NOT OF PRODUCT QUALITY AND MAY HAVE ERRORS OR DEFECTS.
33    THE JPEG-LS STANDARD IS STILL UNDER DEVELOPMENT. THE SOFTWARE IS NOT A
34    FINAL OR FULL IMPLEMENTATION OF THE STANDARD.  HP GIVES NO EXPRESS OR
35    IMPLIED WARRANTY OF ANY KIND AND ANY IMPLIED WARRANTIES OF
36    MERCHANTABILITY AND FITNESS FOR PURPOSE ARE DISCLAIMED.
37         HP SHALL NOT BE LIABLE FOR ANY DIRECT, INDIRECT, SPECIAL, INCIDENTAL,
38    OR CONSEQUENTIAL DAMAGES ARISING OUT OF ANY USE OF THE SOFTWARE.
39    -------------------------------------------------------------------------------
40 */
41
42 /*
43  * bitio.c --- for I/O routines
44  *
45  * Initial code by Alex Jakulin,  Aug. 1995
46  *
47  * Modified and optimized: Gadiel Seroussi, October 1995
48  *
49  * Modified and added Restart marker and input tables by:
50  * David Cheng-Hsiu Chu, and Ismail R. Ismail march 1999
51  */
52
53 #include "global.h"
54 #include "bitio.h"
55
56 extern FILE *in, *out;
57
58 byte negbuff[BUFSIZE+4];    /* byte I/O buffer, allowing for 4 "negative"
59              locations  */
60
61  /*
62  'buff' is defined as 'rawbuff+4' in bitio.h, so that buff[-4]..buff[-1]
63  are well defined. Those locations are used to "return" data to
64  the byte buffer when flushing the input bit buffer .
65  */
66
67 int fp;     /* index into byte buffer */
68 int truebufsize;        /* true size of byte buffer ( <= BUFSIZE) */
69 int foundeof;
70
71
72 /* BIT I/O variables */
73 dword reg;     /* BIT buffer for input/output */
74 int bits;          /* number of bits free in bit buffer (on output) */
75        /* (number of bits free)-8 in bit buffer (on input)*/
76
77
78
79 /****************************************************************************
80  *  OUTPUT ROUTINES
81  *  note: some routines are implemented as preprocessor macros. See bitio.h.
82  ****************************************************************************/
83
84 void flushbuff(FILE *fil) {
85   /* fwrite must work correctly, even if fp is equal to 0 */
86     fwrite(buff, 1, fp, fil);
87     fp = 0;
88 }
89
90
91
92 /* Flushes the bit output buffer and the byte output buffer */
93 void bitoflush() {
94   register unsigned int outbyte;
95     
96     while (bits < 32) {
97     outbyte = reg >> 24;
98         myputc(outbyte, out);
99     if ( outbyte == 0xff ) {
100       bits += 7;
101       reg <<= 7;
102       reg &= ~(1<<(8*sizeof(reg)-1)); /* stuff a 0 at MSB */
103     } else {
104         bits += 8;
105         reg <<= 8;
106     }
107   }
108   flushbuff(out);
109   bitoinit();
110 }
111
112
113
114 /* Initializes the bit output routines */
115 void bitoinit() {
116   bits = 32;
117   reg = 0;
118   fp = 0;
119 }