@gmod/cram 6.0.0 → 7.0.0

This diff represents the content of publicly available package versions that have been released to one of the supported registries. The information contained in this diff is provided for informational purposes only and reflects changes between package versions as they appear in their respective public registries.
Files changed (178) hide show
  1. package/dist/craiIndex.js +2 -2
  2. package/dist/craiIndex.js.map +1 -1
  3. package/dist/cram-bundle.js +1 -1
  4. package/dist/cramFile/file.d.ts +2 -2
  5. package/dist/cramFile/file.js +9 -13
  6. package/dist/cramFile/file.js.map +1 -1
  7. package/dist/htscodecs/index.d.ts +1 -5
  8. package/dist/htscodecs/index.js +11 -72
  9. package/dist/htscodecs/index.js.map +1 -1
  10. package/dist/htscodecs-wasm.d.ts +7 -0
  11. package/dist/htscodecs-wasm.js +183 -0
  12. package/dist/htscodecs-wasm.js.map +1 -0
  13. package/dist/unzip.d.ts +1 -1
  14. package/dist/unzip.js +3 -3
  15. package/dist/unzip.js.map +1 -1
  16. package/dist/wasm/htscodecs.d.ts +5 -0
  17. package/dist/wasm/htscodecs.js +0 -0
  18. package/dist/wasm/htscodecs.js.map +1 -0
  19. package/esm/craiIndex.js +2 -2
  20. package/esm/craiIndex.js.map +1 -1
  21. package/esm/cramFile/file.d.ts +2 -2
  22. package/esm/cramFile/file.js +7 -11
  23. package/esm/cramFile/file.js.map +1 -1
  24. package/esm/htscodecs/index.d.ts +1 -5
  25. package/esm/htscodecs/index.js +3 -31
  26. package/esm/htscodecs/index.js.map +1 -1
  27. package/esm/htscodecs-wasm.d.ts +7 -0
  28. package/esm/htscodecs-wasm.js +171 -0
  29. package/esm/htscodecs-wasm.js.map +1 -0
  30. package/esm/unzip.d.ts +1 -1
  31. package/esm/unzip.js +3 -3
  32. package/esm/unzip.js.map +1 -1
  33. package/esm/wasm/htscodecs.d.ts +5 -0
  34. package/esm/wasm/htscodecs.js +0 -0
  35. package/esm/wasm/htscodecs.js.map +1 -0
  36. package/package.json +9 -9
  37. package/src/craiIndex.ts +4 -2
  38. package/src/cramFile/file.ts +7 -11
  39. package/src/htscodecs/index.ts +12 -38
  40. package/src/htscodecs-wasm.ts +210 -0
  41. package/src/unzip.ts +3 -3
  42. package/src/wasm/htscodecs.d.ts +58 -0
  43. package/src/wasm/htscodecs.js +0 -0
  44. package/dist/htscodecs/arith_gen.d.ts +0 -16
  45. package/dist/htscodecs/arith_gen.js +0 -323
  46. package/dist/htscodecs/arith_gen.js.map +0 -1
  47. package/dist/htscodecs/arith_sh.d.ts +0 -9
  48. package/dist/htscodecs/arith_sh.js +0 -134
  49. package/dist/htscodecs/arith_sh.js.map +0 -1
  50. package/dist/htscodecs/byte_model.d.ts +0 -6
  51. package/dist/htscodecs/byte_model.js +0 -121
  52. package/dist/htscodecs/byte_model.js.map +0 -1
  53. package/dist/htscodecs/fqzcomp.d.ts +0 -1
  54. package/dist/htscodecs/fqzcomp.js +0 -346
  55. package/dist/htscodecs/fqzcomp.js.map +0 -1
  56. package/dist/htscodecs/iostream.d.ts +0 -16
  57. package/dist/htscodecs/iostream.js +0 -148
  58. package/dist/htscodecs/iostream.js.map +0 -1
  59. package/dist/htscodecs/rans.d.ts +0 -1
  60. package/dist/htscodecs/rans.js +0 -224
  61. package/dist/htscodecs/rans.js.map +0 -1
  62. package/dist/htscodecs/rans4x16.d.ts +0 -1
  63. package/dist/htscodecs/rans4x16.js +0 -445
  64. package/dist/htscodecs/rans4x16.js.map +0 -1
  65. package/dist/htscodecs/tok3.d.ts +0 -1
  66. package/dist/htscodecs/tok3.js +0 -206
  67. package/dist/htscodecs/tok3.js.map +0 -1
  68. package/dist/rans/constants.d.ts +0 -4
  69. package/dist/rans/constants.js +0 -10
  70. package/dist/rans/constants.js.map +0 -1
  71. package/dist/rans/d04.d.ts +0 -1
  72. package/dist/rans/d04.js +0 -100
  73. package/dist/rans/d04.js.map +0 -1
  74. package/dist/rans/d14.d.ts +0 -1
  75. package/dist/rans/d14.js +0 -96
  76. package/dist/rans/d14.js.map +0 -1
  77. package/dist/rans/decoding.d.ts +0 -27
  78. package/dist/rans/decoding.js +0 -119
  79. package/dist/rans/decoding.js.map +0 -1
  80. package/dist/rans/frequencies.d.ts +0 -2
  81. package/dist/rans/frequencies.js +0 -112
  82. package/dist/rans/frequencies.js.map +0 -1
  83. package/dist/rans/index.d.ts +0 -1
  84. package/dist/rans/index.js +0 -198
  85. package/dist/rans/index.js.map +0 -1
  86. package/dist/seek-bzip/bitreader.d.ts +0 -27
  87. package/dist/seek-bzip/bitreader.js +0 -101
  88. package/dist/seek-bzip/bitreader.js.map +0 -1
  89. package/dist/seek-bzip/crc32.d.ts +0 -26
  90. package/dist/seek-bzip/crc32.js +0 -115
  91. package/dist/seek-bzip/crc32.js.map +0 -1
  92. package/dist/seek-bzip/index.d.ts +0 -9
  93. package/dist/seek-bzip/index.js +0 -585
  94. package/dist/seek-bzip/index.js.map +0 -1
  95. package/dist/seek-bzip/stream.d.ts +0 -8
  96. package/dist/seek-bzip/stream.js +0 -37
  97. package/dist/seek-bzip/stream.js.map +0 -1
  98. package/dist/seek-bzip/toHex.d.ts +0 -1
  99. package/dist/seek-bzip/toHex.js +0 -36
  100. package/dist/seek-bzip/toHex.js.map +0 -1
  101. package/esm/htscodecs/arith_gen.d.ts +0 -16
  102. package/esm/htscodecs/arith_gen.js +0 -317
  103. package/esm/htscodecs/arith_gen.js.map +0 -1
  104. package/esm/htscodecs/arith_sh.d.ts +0 -9
  105. package/esm/htscodecs/arith_sh.js +0 -131
  106. package/esm/htscodecs/arith_sh.js.map +0 -1
  107. package/esm/htscodecs/byte_model.d.ts +0 -6
  108. package/esm/htscodecs/byte_model.js +0 -118
  109. package/esm/htscodecs/byte_model.js.map +0 -1
  110. package/esm/htscodecs/fqzcomp.d.ts +0 -1
  111. package/esm/htscodecs/fqzcomp.js +0 -340
  112. package/esm/htscodecs/fqzcomp.js.map +0 -1
  113. package/esm/htscodecs/iostream.d.ts +0 -16
  114. package/esm/htscodecs/iostream.js +0 -145
  115. package/esm/htscodecs/iostream.js.map +0 -1
  116. package/esm/htscodecs/rans.d.ts +0 -1
  117. package/esm/htscodecs/rans.js +0 -218
  118. package/esm/htscodecs/rans.js.map +0 -1
  119. package/esm/htscodecs/rans4x16.d.ts +0 -1
  120. package/esm/htscodecs/rans4x16.js +0 -439
  121. package/esm/htscodecs/rans4x16.js.map +0 -1
  122. package/esm/htscodecs/tok3.d.ts +0 -1
  123. package/esm/htscodecs/tok3.js +0 -167
  124. package/esm/htscodecs/tok3.js.map +0 -1
  125. package/esm/rans/constants.d.ts +0 -4
  126. package/esm/rans/constants.js +0 -5
  127. package/esm/rans/constants.js.map +0 -1
  128. package/esm/rans/d04.d.ts +0 -1
  129. package/esm/rans/d04.js +0 -94
  130. package/esm/rans/d04.js.map +0 -1
  131. package/esm/rans/d14.d.ts +0 -1
  132. package/esm/rans/d14.js +0 -93
  133. package/esm/rans/d14.js.map +0 -1
  134. package/esm/rans/decoding.d.ts +0 -27
  135. package/esm/rans/decoding.js +0 -117
  136. package/esm/rans/decoding.js.map +0 -1
  137. package/esm/rans/frequencies.d.ts +0 -2
  138. package/esm/rans/frequencies.js +0 -105
  139. package/esm/rans/frequencies.js.map +0 -1
  140. package/esm/rans/index.d.ts +0 -1
  141. package/esm/rans/index.js +0 -192
  142. package/esm/rans/index.js.map +0 -1
  143. package/esm/seek-bzip/bitreader.d.ts +0 -27
  144. package/esm/seek-bzip/bitreader.js +0 -98
  145. package/esm/seek-bzip/bitreader.js.map +0 -1
  146. package/esm/seek-bzip/crc32.d.ts +0 -26
  147. package/esm/seek-bzip/crc32.js +0 -112
  148. package/esm/seek-bzip/crc32.js.map +0 -1
  149. package/esm/seek-bzip/index.d.ts +0 -9
  150. package/esm/seek-bzip/index.js +0 -579
  151. package/esm/seek-bzip/index.js.map +0 -1
  152. package/esm/seek-bzip/stream.d.ts +0 -8
  153. package/esm/seek-bzip/stream.js +0 -34
  154. package/esm/seek-bzip/stream.js.map +0 -1
  155. package/esm/seek-bzip/toHex.d.ts +0 -1
  156. package/esm/seek-bzip/toHex.js +0 -33
  157. package/esm/seek-bzip/toHex.js.map +0 -1
  158. package/src/htscodecs/arith_gen.ts +0 -371
  159. package/src/htscodecs/arith_sh.ts +0 -146
  160. package/src/htscodecs/byte_model.ts +0 -136
  161. package/src/htscodecs/fqzcomp.ts +0 -394
  162. package/src/htscodecs/iostream.ts +0 -159
  163. package/src/htscodecs/rans.ts +0 -250
  164. package/src/htscodecs/rans4x16.ts +0 -521
  165. package/src/htscodecs/tok3.ts +0 -197
  166. package/src/rans/README.md +0 -1
  167. package/src/rans/constants.ts +0 -5
  168. package/src/rans/d04.ts +0 -112
  169. package/src/rans/d14.ts +0 -106
  170. package/src/rans/decoding.ts +0 -138
  171. package/src/rans/frequencies.ts +0 -115
  172. package/src/rans/index.ts +0 -245
  173. package/src/seek-bzip/README.md +0 -6
  174. package/src/seek-bzip/bitreader.ts +0 -113
  175. package/src/seek-bzip/crc32.ts +0 -121
  176. package/src/seek-bzip/index.ts +0 -615
  177. package/src/seek-bzip/stream.ts +0 -33
  178. package/src/seek-bzip/toHex.ts +0 -33
@@ -1,36 +0,0 @@
1
- "use strict";
2
- // from https://www.xaymar.com/articles/2020/12/08/fastest-uint8array-to-hex-string-conversion-in-javascript/
3
- Object.defineProperty(exports, "__esModule", { value: true });
4
- exports.toHex = toHex;
5
- // Pre-Init
6
- const LUT_HEX_4b = [
7
- '0',
8
- '1',
9
- '2',
10
- '3',
11
- '4',
12
- '5',
13
- '6',
14
- '7',
15
- '8',
16
- '9',
17
- 'A',
18
- 'B',
19
- 'C',
20
- 'D',
21
- 'E',
22
- 'F',
23
- ];
24
- const LUT_HEX_8b = new Array(0x100);
25
- for (let n = 0; n < 0x100; n++) {
26
- LUT_HEX_8b[n] = `${LUT_HEX_4b[(n >>> 4) & 0xf]}${LUT_HEX_4b[n & 0xf]}`;
27
- }
28
- // End Pre-Init
29
- function toHex(buffer) {
30
- let out = '';
31
- for (let idx = 0, edx = buffer.length; idx < edx; idx++) {
32
- out += LUT_HEX_8b[buffer[idx]];
33
- }
34
- return out;
35
- }
36
- //# sourceMappingURL=toHex.js.map
@@ -1 +0,0 @@
1
- {"version":3,"file":"toHex.js","sourceRoot":"","sources":["../../src/seek-bzip/toHex.ts"],"names":[],"mappings":";AAAA,6GAA6G;;AA0B7G,sBAMC;AA9BD,WAAW;AACX,MAAM,UAAU,GAAG;IACjB,GAAG;IACH,GAAG;IACH,GAAG;IACH,GAAG;IACH,GAAG;IACH,GAAG;IACH,GAAG;IACH,GAAG;IACH,GAAG;IACH,GAAG;IACH,GAAG;IACH,GAAG;IACH,GAAG;IACH,GAAG;IACH,GAAG;IACH,GAAG;CACJ,CAAA;AACD,MAAM,UAAU,GAAG,IAAI,KAAK,CAAC,KAAK,CAAC,CAAA;AACnC,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,KAAK,EAAE,CAAC,EAAE,EAAE,CAAC;IAC/B,UAAU,CAAC,CAAC,CAAC,GAAG,GAAG,UAAU,CAAC,CAAC,CAAC,KAAK,CAAC,CAAC,GAAG,GAAG,CAAC,GAAG,UAAU,CAAC,CAAC,GAAG,GAAG,CAAC,EAAE,CAAA;AACxE,CAAC;AACD,eAAe;AACf,SAAgB,KAAK,CAAC,MAAkB;IACtC,IAAI,GAAG,GAAG,EAAE,CAAA;IACZ,KAAK,IAAI,GAAG,GAAG,CAAC,EAAE,GAAG,GAAG,MAAM,CAAC,MAAM,EAAE,GAAG,GAAG,GAAG,EAAE,GAAG,EAAE,EAAE,CAAC;QACxD,GAAG,IAAI,UAAU,CAAC,MAAM,CAAC,GAAG,CAAE,CAAC,CAAA;IACjC,CAAC;IACD,OAAO,GAAG,CAAA;AACZ,CAAC"}
@@ -1,16 +0,0 @@
1
- import IOStream from './iostream.ts';
2
- export default class RangeCoderGen {
3
- stream: IOStream;
4
- decode(src: Uint8Array): Uint8Array<any>;
5
- decodeStream(stream: IOStream, n_out?: number): Uint8Array<any>;
6
- decode0(stream: any, n_out: any): Uint8Array<any>;
7
- decode1(stream: any, n_out: any): Uint8Array<any>;
8
- decodeExt(stream: any, n_out: any): any;
9
- decodeRLE0(stream: any, n_out: any): Uint8Array<any>;
10
- decodeRLE1(stream: any, n_out: any): Uint8Array<any>;
11
- decodePackMeta(stream: any): any[];
12
- decodePack(data: any, M: any, len: any): any;
13
- packMeta(src: any): (number | any[] | IOStream)[];
14
- decodeStripe(stream: any, len: any): Uint8Array<any>;
15
- decodeCat(stream: any, len: any): Uint8Array<any>;
16
- }
@@ -1,317 +0,0 @@
1
- /* eslint-disable no-var */
2
- // @ts-nocheck
3
- /*
4
- * Copyright (c) 2019,2020 Genome Research Ltd.
5
- * Author(s): James Bonfield
6
- *
7
- * Redistribution and use in source and binary forms, with or without
8
- * modification, are permitted provided that the following conditions are met:
9
- *
10
- * 1. Redistributions of source code must retain the above copyright notice,
11
- * this list of conditions and the following disclaimer.
12
- *
13
- * 2. Redistributions in binary form must reproduce the above
14
- * copyright notice, this list of conditions and the following
15
- * disclaimer in the documentation and/or other materials provided
16
- * with the distribution.
17
- *
18
- * 3. Neither the names Genome Research Ltd and Wellcome Trust Sanger
19
- * Institute nor the names of its contributors may be used to endorse
20
- * or promote products derived from this software without specific
21
- * prior written permission.
22
- *
23
- * THIS SOFTWARE IS PROVIDED BY GENOME RESEARCH LTD AND CONTRIBUTORS "AS
24
- * IS" AND ANY EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT LIMITED
25
- * TO, THE IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS FOR A
26
- * PARTICULAR PURPOSE ARE DISCLAIMED. IN NO EVENT SHALL GENOME RESEARCH
27
- * LTD OR CONTRIBUTORS BE LIABLE FOR ANY DIRECT, INDIRECT, INCIDENTAL,
28
- * SPECIAL, EXEMPLARY, OR CONSEQUENTIAL DAMAGES (INCLUDING, BUT NOT
29
- * LIMITED TO, PROCUREMENT OF SUBSTITUTE GOODS OR SERVICES; LOSS OF USE,
30
- * DATA, OR PROFITS; OR BUSINESS INTERRUPTION) HOWEVER CAUSED AND ON ANY
31
- * THEORY OF LIABILITY, WHETHER IN CONTRACT, STRICT LIABILITY, OR TORT
32
- * (INCLUDING NEGLIGENCE OR OTHERWISE) ARISING IN ANY WAY OUT OF THE USE
33
- * OF THIS SOFTWARE, EVEN IF ADVISED OF THE POSSIBILITY OF SUCH DAMAGE.
34
- */
35
- import { decode } from "../seek-bzip/index.js";
36
- import RangeCoder from "./arith_sh.js";
37
- import ByteModel from "./byte_model.js";
38
- import IOStream from "./iostream.js";
39
- const ARITH_ORDER = 1;
40
- const ARITH_EXT = 4;
41
- const ARITH_STRIPE = 8;
42
- const ARITH_NOSIZE = 16;
43
- const ARITH_CAT = 32;
44
- const ARITH_RLE = 64;
45
- const ARITH_PACK = 128;
46
- export default class RangeCoderGen {
47
- decode(src) {
48
- this.stream = new IOStream(src);
49
- return this.decodeStream(this.stream);
50
- }
51
- decodeStream(stream, n_out = 0) {
52
- const flags = this.stream.ReadByte();
53
- if (!(flags & ARITH_NOSIZE)) {
54
- n_out = this.stream.ReadUint7();
55
- }
56
- let e_len = n_out;
57
- const order = flags & ARITH_ORDER;
58
- // 4-way recursion
59
- if (flags & ARITH_STRIPE) {
60
- return this.decodeStripe(this.stream, n_out);
61
- }
62
- // Meta data
63
- if (flags & ARITH_PACK) {
64
- var P;
65
- [P, e_len] = this.decodePackMeta(this.stream);
66
- }
67
- // NOP, useful for tiny blocks
68
- if (flags & ARITH_CAT) {
69
- var data = this.decodeCat(this.stream, e_len);
70
- }
71
- // Entropy decode
72
- else if (flags & ARITH_EXT) {
73
- var data = this.decodeExt(this.stream, e_len);
74
- }
75
- else if (flags & ARITH_RLE) {
76
- var data = order
77
- ? this.decodeRLE1(this.stream, e_len)
78
- : this.decodeRLE0(this.stream, e_len);
79
- }
80
- else {
81
- var data = order
82
- ? this.decode1(this.stream, e_len)
83
- : this.decode0(this.stream, e_len);
84
- }
85
- // Transforms
86
- if (flags & ARITH_PACK) {
87
- data = this.decodePack(data, P, n_out);
88
- }
89
- return data;
90
- }
91
- // ----------------------------------------------------------------------
92
- // Order-0 codec
93
- decode0(stream, n_out) {
94
- const output = new Uint8Array(n_out);
95
- let max_sym = stream.ReadByte();
96
- if (max_sym == 0) {
97
- max_sym = 256;
98
- }
99
- const byte_model = new ByteModel(max_sym);
100
- const rc = new RangeCoder(stream);
101
- rc.RangeStartDecode(stream);
102
- for (let i = 0; i < n_out; i++) {
103
- output[i] = byte_model.ModelDecode(stream, rc);
104
- }
105
- return output;
106
- }
107
- // ----------------------------------------------------------------------
108
- // Order-1 codec
109
- decode1(stream, n_out) {
110
- const output = new Uint8Array(n_out);
111
- let max_sym = stream.ReadByte();
112
- if (max_sym == 0) {
113
- max_sym = 256;
114
- }
115
- const byte_model = new Array(max_sym);
116
- for (var i = 0; i < max_sym; i++) {
117
- byte_model[i] = new ByteModel(max_sym);
118
- }
119
- const rc = new RangeCoder(stream);
120
- rc.RangeStartDecode(stream);
121
- let last = 0;
122
- for (var i = 0; i < n_out; i++) {
123
- output[i] = byte_model[last].ModelDecode(stream, rc);
124
- last = output[i];
125
- }
126
- return output;
127
- }
128
- // ----------------------------------------------------------------------
129
- // External codec
130
- decodeExt(stream, n_out) {
131
- return decode(stream.buf.slice(stream.pos));
132
- }
133
- // ----------------------------------------------------------------------
134
- // Order-0 RLE codec
135
- decodeRLE0(stream, n_out) {
136
- const output = new Uint8Array(n_out);
137
- let max_sym = stream.ReadByte();
138
- if (max_sym == 0) {
139
- max_sym = 256;
140
- }
141
- const model_lit = new ByteModel(max_sym);
142
- const model_run = new Array(258);
143
- for (var i = 0; i <= 257; i++) {
144
- model_run[i] = new ByteModel(4);
145
- }
146
- const rc = new RangeCoder(stream);
147
- rc.RangeStartDecode(stream);
148
- var i = 0;
149
- while (i < n_out) {
150
- output[i] = model_lit.ModelDecode(stream, rc);
151
- let part = model_run[output[i]].ModelDecode(stream, rc);
152
- let run = part;
153
- let rctx = 256;
154
- while (part == 3) {
155
- part = model_run[rctx].ModelDecode(stream, rc);
156
- rctx = 257;
157
- run += part;
158
- }
159
- for (let j = 1; j <= run; j++) {
160
- output[i + j] = output[i];
161
- }
162
- i += run + 1;
163
- }
164
- return output;
165
- }
166
- // ----------------------------------------------------------------------
167
- // Order-1 RLE codec
168
- decodeRLE1(stream, n_out) {
169
- const output = new Uint8Array(n_out);
170
- let max_sym = stream.ReadByte();
171
- if (max_sym == 0) {
172
- max_sym = 256;
173
- }
174
- const model_lit = new Array(max_sym);
175
- for (var i = 0; i < max_sym; i++) {
176
- model_lit[i] = new ByteModel(max_sym);
177
- }
178
- const model_run = new Array(258);
179
- for (var i = 0; i <= 257; i++) {
180
- model_run[i] = new ByteModel(4);
181
- }
182
- const rc = new RangeCoder(stream);
183
- rc.RangeStartDecode(stream);
184
- let last = 0;
185
- var i = 0;
186
- while (i < n_out) {
187
- output[i] = model_lit[last].ModelDecode(stream, rc);
188
- last = output[i];
189
- let part = model_run[output[i]].ModelDecode(stream, rc);
190
- let run = part;
191
- let rctx = 256;
192
- while (part == 3) {
193
- part = model_run[rctx].ModelDecode(stream, rc);
194
- rctx = 257;
195
- run += part;
196
- }
197
- for (let j = 1; j <= run; j++) {
198
- output[i + j] = output[i];
199
- }
200
- i += run + 1;
201
- }
202
- return output;
203
- }
204
- // ----------------------------------------------------------------------
205
- // Pack method
206
- decodePackMeta(stream) {
207
- this.nsym = stream.ReadByte();
208
- const M = new Array(this.nsym);
209
- for (let i = 0; i < this.nsym; i++) {
210
- M[i] = stream.ReadByte();
211
- }
212
- const e_len = stream.ReadUint7(); // Could be derived data from nsym and n_out
213
- return [M, e_len];
214
- }
215
- decodePack(data, M, len) {
216
- const out = new Uint8Array(len);
217
- if (this.nsym <= 1) {
218
- // Constant value
219
- for (var i = 0; i < len; i++) {
220
- out[i] = M[0];
221
- }
222
- }
223
- else if (this.nsym <= 2) {
224
- // 1 bit per value
225
- for (var i = 0, j = 0; i < len; i++) {
226
- if (i % 8 == 0) {
227
- var v = data[j++];
228
- }
229
- out[i] = M[v & 1];
230
- v >>= 1;
231
- }
232
- }
233
- else if (this.nsym <= 4) {
234
- // 2 bits per value
235
- for (var i = 0, j = 0; i < len; i++) {
236
- if (i % 4 == 0) {
237
- var v = data[j++];
238
- }
239
- out[i] = M[v & 3];
240
- v >>= 2;
241
- }
242
- }
243
- else if (this.nsym <= 16) {
244
- // 4 bits per value
245
- for (var i = 0, j = 0; i < len; i++) {
246
- if (i % 2 == 0) {
247
- var v = data[j++];
248
- }
249
- out[i] = M[v & 15];
250
- v >>= 4;
251
- }
252
- }
253
- else {
254
- // 8 bits per value: NOP
255
- return data;
256
- }
257
- return out;
258
- }
259
- // Compute M array and return meta-data stream
260
- packMeta(src) {
261
- const stream = new IOStream('', 0, 1024);
262
- // Count symbols
263
- const M = new Array(256);
264
- for (var i = 0; i < src.length; i++) {
265
- M[src[i]] = 1;
266
- }
267
- // Write Map
268
- for (var nsym = 0, i = 0; i < 256; i++) {
269
- if (M[i]) {
270
- M[i] = ++nsym;
271
- }
272
- } // map to 1..N
273
- stream.WriteByte(nsym);
274
- // FIXME: add check for nsym > 16?
275
- // Or just accept it as an inefficient waste of time.
276
- for (var i = 0; i < 256; i++) {
277
- if (M[i]) {
278
- stream.WriteByte(i); // adjust to 0..N-1
279
- M[i]--;
280
- }
281
- }
282
- return [stream, M, nsym];
283
- }
284
- decodeStripe(stream, len) {
285
- const N = stream.ReadByte();
286
- // Retrieve lengths
287
- const clen = new Array(N);
288
- const ulen = new Array(N);
289
- for (var j = 0; j < N; j++) {
290
- clen[j] = stream.ReadUint7();
291
- }
292
- // Decode streams
293
- const T = new Array(N);
294
- for (var j = 0; j < N; j++) {
295
- ulen[j] = Math.floor(len / N) + (len % N > j);
296
- T[j] = this.decodeStream(stream, ulen[j]);
297
- }
298
- // Transpose
299
- const out = new Uint8Array(len);
300
- for (var j = 0; j < N; j++) {
301
- for (let i = 0; i < ulen[j]; i++) {
302
- out[i * N + j] = T[j][i];
303
- }
304
- }
305
- return out;
306
- }
307
- // ----------------------------------------------------------------------
308
- // Cat method
309
- decodeCat(stream, len) {
310
- const out = new Uint8Array(len);
311
- for (let i = 0; i < len; i++) {
312
- out[i] = stream.ReadByte();
313
- }
314
- return out;
315
- }
316
- }
317
- //# sourceMappingURL=arith_gen.js.map
@@ -1 +0,0 @@
1
- {"version":3,"file":"arith_gen.js","sourceRoot":"","sources":["../../src/htscodecs/arith_gen.ts"],"names":[],"mappings":"AAAA,2BAA2B;AAC3B,cAAc;AAEd;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;GA+BG;AAEH,OAAO,EAAE,MAAM,EAAE,MAAM,uBAAuB,CAAA;AAE9C,OAAO,UAAU,MAAM,eAAe,CAAA;AACtC,OAAO,SAAS,MAAM,iBAAiB,CAAA;AACvC,OAAO,QAAQ,MAAM,eAAe,CAAA;AAGpC,MAAM,WAAW,GAAG,CAAC,CAAA;AACrB,MAAM,SAAS,GAAG,CAAC,CAAA;AACnB,MAAM,YAAY,GAAG,CAAC,CAAA;AACtB,MAAM,YAAY,GAAG,EAAE,CAAA;AACvB,MAAM,SAAS,GAAG,EAAE,CAAA;AACpB,MAAM,SAAS,GAAG,EAAE,CAAA;AACpB,MAAM,UAAU,GAAG,GAAG,CAAA;AAEtB,MAAM,CAAC,OAAO,OAAO,aAAa;IAEhC,MAAM,CAAC,GAAe;QACpB,IAAI,CAAC,MAAM,GAAG,IAAI,QAAQ,CAAC,GAAG,CAAC,CAAA;QAC/B,OAAO,IAAI,CAAC,YAAY,CAAC,IAAI,CAAC,MAAM,CAAC,CAAA;IACvC,CAAC;IAED,YAAY,CAAC,MAAgB,EAAE,KAAK,GAAG,CAAC;QACtC,MAAM,KAAK,GAAG,IAAI,CAAC,MAAM,CAAC,QAAQ,EAAE,CAAA;QACpC,IAAI,CAAC,CAAC,KAAK,GAAG,YAAY,CAAC,EAAE,CAAC;YAC5B,KAAK,GAAG,IAAI,CAAC,MAAM,CAAC,SAAS,EAAE,CAAA;QACjC,CAAC;QACD,IAAI,KAAK,GAAG,KAAK,CAAA;QAEjB,MAAM,KAAK,GAAG,KAAK,GAAG,WAAW,CAAA;QAEjC,kBAAkB;QAClB,IAAI,KAAK,GAAG,YAAY,EAAE,CAAC;YACzB,OAAO,IAAI,CAAC,YAAY,CAAC,IAAI,CAAC,MAAM,EAAE,KAAK,CAAC,CAAA;QAC9C,CAAC;QAED,YAAY;QACZ,IAAI,KAAK,GAAG,UAAU,EAAE,CAAC;YACvB,IAAI,CAAC,CACJ;YAAA,CAAC,CAAC,EAAE,KAAK,CAAC,GAAG,IAAI,CAAC,cAAc,CAAC,IAAI,CAAC,MAAM,CAAC,CAAA;QAChD,CAAC;QAED,8BAA8B;QAC9B,IAAI,KAAK,GAAG,SAAS,EAAE,CAAC;YACtB,IAAI,IAAI,GAAG,IAAI,CAAC,SAAS,CAAC,IAAI,CAAC,MAAM,EAAE,KAAK,CAAC,CAAA;QAC/C,CAAC;QACD,iBAAiB;aACZ,IAAI,KAAK,GAAG,SAAS,EAAE,CAAC;YAC3B,IAAI,IAAI,GAAG,IAAI,CAAC,SAAS,CAAC,IAAI,CAAC,MAAM,EAAE,KAAK,CAAC,CAAA;QAC/C,CAAC;aAAM,IAAI,KAAK,GAAG,SAAS,EAAE,CAAC;YAC7B,IAAI,IAAI,GAAG,KAAK;gBACd,CAAC,CAAC,IAAI,CAAC,UAAU,CAAC,IAAI,CAAC,MAAM,EAAE,KAAK,CAAC;gBACrC,CAAC,CAAC,IAAI,CAAC,UAAU,CAAC,IAAI,CAAC,MAAM,EAAE,KAAK,CAAC,CAAA;QACzC,CAAC;aAAM,CAAC;YACN,IAAI,IAAI,GAAG,KAAK;gBACd,CAAC,CAAC,IAAI,CAAC,OAAO,CAAC,IAAI,CAAC,MAAM,EAAE,KAAK,CAAC;gBAClC,CAAC,CAAC,IAAI,CAAC,OAAO,CAAC,IAAI,CAAC,MAAM,EAAE,KAAK,CAAC,CAAA;QACtC,CAAC;QAED,aAAa;QACb,IAAI,KAAK,GAAG,UAAU,EAAE,CAAC;YACvB,IAAI,GAAG,IAAI,CAAC,UAAU,CAAC,IAAI,EAAE,CAAC,EAAE,KAAK,CAAC,CAAA;QACxC,CAAC;QAED,OAAO,IAAI,CAAA;IACb,CAAC;IAED,yEAAyE;IACzE,gBAAgB;IAChB,OAAO,CAAC,MAAM,EAAE,KAAK;QACnB,MAAM,MAAM,GAAG,IAAI,UAAU,CAAC,KAAK,CAAC,CAAA;QAEpC,IAAI,OAAO,GAAG,MAAM,CAAC,QAAQ,EAAE,CAAA;QAC/B,IAAI,OAAO,IAAI,CAAC,EAAE,CAAC;YACjB,OAAO,GAAG,GAAG,CAAA;QACf,CAAC;QAED,MAAM,UAAU,GAAG,IAAI,SAAS,CAAC,OAAO,CAAC,CAAA;QAEzC,MAAM,EAAE,GAAG,IAAI,UAAU,CAAC,MAAM,CAAC,CAAA;QACjC,EAAE,CAAC,gBAAgB,CAAC,MAAM,CAAC,CAAA;QAE3B,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,KAAK,EAAE,CAAC,EAAE,EAAE,CAAC;YAC/B,MAAM,CAAC,CAAC,CAAC,GAAG,UAAU,CAAC,WAAW,CAAC,MAAM,EAAE,EAAE,CAAC,CAAA;QAChD,CAAC;QAED,OAAO,MAAM,CAAA;IACf,CAAC;IAED,yEAAyE;IACzE,gBAAgB;IAEhB,OAAO,CAAC,MAAM,EAAE,KAAK;QACnB,MAAM,MAAM,GAAG,IAAI,UAAU,CAAC,KAAK,CAAC,CAAA;QAEpC,IAAI,OAAO,GAAG,MAAM,CAAC,QAAQ,EAAE,CAAA;QAC/B,IAAI,OAAO,IAAI,CAAC,EAAE,CAAC;YACjB,OAAO,GAAG,GAAG,CAAA;QACf,CAAC;QAED,MAAM,UAAU,GAAG,IAAI,KAAK,CAAC,OAAO,CAAC,CAAA;QACrC,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,OAAO,EAAE,CAAC,EAAE,EAAE,CAAC;YACjC,UAAU,CAAC,CAAC,CAAC,GAAG,IAAI,SAAS,CAAC,OAAO,CAAC,CAAA;QACxC,CAAC;QAED,MAAM,EAAE,GAAG,IAAI,UAAU,CAAC,MAAM,CAAC,CAAA;QACjC,EAAE,CAAC,gBAAgB,CAAC,MAAM,CAAC,CAAA;QAE3B,IAAI,IAAI,GAAG,CAAC,CAAA;QACZ,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,KAAK,EAAE,CAAC,EAAE,EAAE,CAAC;YAC/B,MAAM,CAAC,CAAC,CAAC,GAAG,UAAU,CAAC,IAAI,CAAC,CAAC,WAAW,CAAC,MAAM,EAAE,EAAE,CAAC,CAAA;YACpD,IAAI,GAAG,MAAM,CAAC,CAAC,CAAC,CAAA;QAClB,CAAC;QAED,OAAO,MAAM,CAAA;IACf,CAAC;IAED,yEAAyE;IACzE,iBAAiB;IACjB,SAAS,CAAC,MAAM,EAAE,KAAK;QACrB,OAAO,MAAM,CAAC,MAAM,CAAC,GAAG,CAAC,KAAK,CAAC,MAAM,CAAC,GAAG,CAAC,CAAC,CAAA;IAC7C,CAAC;IAED,yEAAyE;IACzE,oBAAoB;IACpB,UAAU,CAAC,MAAM,EAAE,KAAK;QACtB,MAAM,MAAM,GAAG,IAAI,UAAU,CAAC,KAAK,CAAC,CAAA;QAEpC,IAAI,OAAO,GAAG,MAAM,CAAC,QAAQ,EAAE,CAAA;QAC/B,IAAI,OAAO,IAAI,CAAC,EAAE,CAAC;YACjB,OAAO,GAAG,GAAG,CAAA;QACf,CAAC;QAED,MAAM,SAAS,GAAG,IAAI,SAAS,CAAC,OAAO,CAAC,CAAA;QACxC,MAAM,SAAS,GAAG,IAAI,KAAK,CAAC,GAAG,CAAC,CAAA;QAChC,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,IAAI,GAAG,EAAE,CAAC,EAAE,EAAE,CAAC;YAC9B,SAAS,CAAC,CAAC,CAAC,GAAG,IAAI,SAAS,CAAC,CAAC,CAAC,CAAA;QACjC,CAAC;QAED,MAAM,EAAE,GAAG,IAAI,UAAU,CAAC,MAAM,CAAC,CAAA;QACjC,EAAE,CAAC,gBAAgB,CAAC,MAAM,CAAC,CAAA;QAE3B,IAAI,CAAC,GAAG,CAAC,CAAA;QACT,OAAO,CAAC,GAAG,KAAK,EAAE,CAAC;YACjB,MAAM,CAAC,CAAC,CAAC,GAAG,SAAS,CAAC,WAAW,CAAC,MAAM,EAAE,EAAE,CAAC,CAAA;YAC7C,IAAI,IAAI,GAAG,SAAS,CAAC,MAAM,CAAC,CAAC,CAAC,CAAC,CAAC,WAAW,CAAC,MAAM,EAAE,EAAE,CAAC,CAAA;YACvD,IAAI,GAAG,GAAG,IAAI,CAAA;YACd,IAAI,IAAI,GAAG,GAAG,CAAA;YACd,OAAO,IAAI,IAAI,CAAC,EAAE,CAAC;gBACjB,IAAI,GAAG,SAAS,CAAC,IAAI,CAAC,CAAC,WAAW,CAAC,MAAM,EAAE,EAAE,CAAC,CAAA;gBAC9C,IAAI,GAAG,GAAG,CAAA;gBACV,GAAG,IAAI,IAAI,CAAA;YACb,CAAC;YACD,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,IAAI,GAAG,EAAE,CAAC,EAAE,EAAE,CAAC;gBAC9B,MAAM,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,MAAM,CAAC,CAAC,CAAC,CAAA;YAC3B,CAAC;YACD,CAAC,IAAI,GAAG,GAAG,CAAC,CAAA;QACd,CAAC;QAED,OAAO,MAAM,CAAA;IACf,CAAC;IAED,yEAAyE;IACzE,oBAAoB;IAEpB,UAAU,CAAC,MAAM,EAAE,KAAK;QACtB,MAAM,MAAM,GAAG,IAAI,UAAU,CAAC,KAAK,CAAC,CAAA;QAEpC,IAAI,OAAO,GAAG,MAAM,CAAC,QAAQ,EAAE,CAAA;QAC/B,IAAI,OAAO,IAAI,CAAC,EAAE,CAAC;YACjB,OAAO,GAAG,GAAG,CAAA;QACf,CAAC;QAED,MAAM,SAAS,GAAG,IAAI,KAAK,CAAC,OAAO,CAAC,CAAA;QACpC,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,OAAO,EAAE,CAAC,EAAE,EAAE,CAAC;YACjC,SAAS,CAAC,CAAC,CAAC,GAAG,IAAI,SAAS,CAAC,OAAO,CAAC,CAAA;QACvC,CAAC;QAED,MAAM,SAAS,GAAG,IAAI,KAAK,CAAC,GAAG,CAAC,CAAA;QAChC,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,IAAI,GAAG,EAAE,CAAC,EAAE,EAAE,CAAC;YAC9B,SAAS,CAAC,CAAC,CAAC,GAAG,IAAI,SAAS,CAAC,CAAC,CAAC,CAAA;QACjC,CAAC;QAED,MAAM,EAAE,GAAG,IAAI,UAAU,CAAC,MAAM,CAAC,CAAA;QACjC,EAAE,CAAC,gBAAgB,CAAC,MAAM,CAAC,CAAA;QAE3B,IAAI,IAAI,GAAG,CAAC,CAAA;QACZ,IAAI,CAAC,GAAG,CAAC,CAAA;QACT,OAAO,CAAC,GAAG,KAAK,EAAE,CAAC;YACjB,MAAM,CAAC,CAAC,CAAC,GAAG,SAAS,CAAC,IAAI,CAAC,CAAC,WAAW,CAAC,MAAM,EAAE,EAAE,CAAC,CAAA;YACnD,IAAI,GAAG,MAAM,CAAC,CAAC,CAAC,CAAA;YAChB,IAAI,IAAI,GAAG,SAAS,CAAC,MAAM,CAAC,CAAC,CAAC,CAAC,CAAC,WAAW,CAAC,MAAM,EAAE,EAAE,CAAC,CAAA;YACvD,IAAI,GAAG,GAAG,IAAI,CAAA;YACd,IAAI,IAAI,GAAG,GAAG,CAAA;YACd,OAAO,IAAI,IAAI,CAAC,EAAE,CAAC;gBACjB,IAAI,GAAG,SAAS,CAAC,IAAI,CAAC,CAAC,WAAW,CAAC,MAAM,EAAE,EAAE,CAAC,CAAA;gBAC9C,IAAI,GAAG,GAAG,CAAA;gBACV,GAAG,IAAI,IAAI,CAAA;YACb,CAAC;YACD,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,IAAI,GAAG,EAAE,CAAC,EAAE,EAAE,CAAC;gBAC9B,MAAM,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,MAAM,CAAC,CAAC,CAAC,CAAA;YAC3B,CAAC;YACD,CAAC,IAAI,GAAG,GAAG,CAAC,CAAA;QACd,CAAC;QAED,OAAO,MAAM,CAAA;IACf,CAAC;IAED,yEAAyE;IACzE,cAAc;IACd,cAAc,CAAC,MAAM;QACnB,IAAI,CAAC,IAAI,GAAG,MAAM,CAAC,QAAQ,EAAE,CAAA;QAE7B,MAAM,CAAC,GAAG,IAAI,KAAK,CAAC,IAAI,CAAC,IAAI,CAAC,CAAA;QAC9B,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,IAAI,CAAC,IAAI,EAAE,CAAC,EAAE,EAAE,CAAC;YACnC,CAAC,CAAC,CAAC,CAAC,GAAG,MAAM,CAAC,QAAQ,EAAE,CAAA;QAC1B,CAAC;QAED,MAAM,KAAK,GAAG,MAAM,CAAC,SAAS,EAAE,CAAA,CAAC,4CAA4C;QAE7E,OAAO,CAAC,CAAC,EAAE,KAAK,CAAC,CAAA;IACnB,CAAC;IAED,UAAU,CAAC,IAAI,EAAE,CAAC,EAAE,GAAG;QACrB,MAAM,GAAG,GAAG,IAAI,UAAU,CAAC,GAAG,CAAC,CAAA;QAE/B,IAAI,IAAI,CAAC,IAAI,IAAI,CAAC,EAAE,CAAC;YACnB,iBAAiB;YACjB,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,GAAG,EAAE,CAAC,EAAE,EAAE,CAAC;gBAC7B,GAAG,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAA;YACf,CAAC;QACH,CAAC;aAAM,IAAI,IAAI,CAAC,IAAI,IAAI,CAAC,EAAE,CAAC;YAC1B,kBAAkB;YAClB,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,GAAG,EAAE,CAAC,EAAE,EAAE,CAAC;gBACpC,IAAI,CAAC,GAAG,CAAC,IAAI,CAAC,EAAE,CAAC;oBACf,IAAI,CAAC,GAAG,IAAI,CAAC,CAAC,EAAE,CAAC,CAAA;gBACnB,CAAC;gBACD,GAAG,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAA;gBACjB,CAAC,KAAK,CAAC,CAAA;YACT,CAAC;QACH,CAAC;aAAM,IAAI,IAAI,CAAC,IAAI,IAAI,CAAC,EAAE,CAAC;YAC1B,mBAAmB;YACnB,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,GAAG,EAAE,CAAC,EAAE,EAAE,CAAC;gBACpC,IAAI,CAAC,GAAG,CAAC,IAAI,CAAC,EAAE,CAAC;oBACf,IAAI,CAAC,GAAG,IAAI,CAAC,CAAC,EAAE,CAAC,CAAA;gBACnB,CAAC;gBACD,GAAG,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAA;gBACjB,CAAC,KAAK,CAAC,CAAA;YACT,CAAC;QACH,CAAC;aAAM,IAAI,IAAI,CAAC,IAAI,IAAI,EAAE,EAAE,CAAC;YAC3B,mBAAmB;YACnB,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,GAAG,EAAE,CAAC,EAAE,EAAE,CAAC;gBACpC,IAAI,CAAC,GAAG,CAAC,IAAI,CAAC,EAAE,CAAC;oBACf,IAAI,CAAC,GAAG,IAAI,CAAC,CAAC,EAAE,CAAC,CAAA;gBACnB,CAAC;gBACD,GAAG,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,GAAG,EAAE,CAAC,CAAA;gBAClB,CAAC,KAAK,CAAC,CAAA;YACT,CAAC;QACH,CAAC;aAAM,CAAC;YACN,wBAAwB;YACxB,OAAO,IAAI,CAAA;QACb,CAAC;QAED,OAAO,GAAG,CAAA;IACZ,CAAC;IAED,8CAA8C;IAC9C,QAAQ,CAAC,GAAG;QACV,MAAM,MAAM,GAAG,IAAI,QAAQ,CAAC,EAAE,EAAE,CAAC,EAAE,IAAI,CAAC,CAAA;QAExC,gBAAgB;QAChB,MAAM,CAAC,GAAG,IAAI,KAAK,CAAC,GAAG,CAAC,CAAA;QACxB,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,GAAG,CAAC,MAAM,EAAE,CAAC,EAAE,EAAE,CAAC;YACpC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAA;QACf,CAAC;QAED,YAAY;QACZ,KAAK,IAAI,IAAI,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,GAAG,EAAE,CAAC,EAAE,EAAE,CAAC;YACvC,IAAI,CAAC,CAAC,CAAC,CAAC,EAAE,CAAC;gBACT,CAAC,CAAC,CAAC,CAAC,GAAG,EAAE,IAAI,CAAA;YACf,CAAC;QACH,CAAC,CAAC,cAAc;QAChB,MAAM,CAAC,SAAS,CAAC,IAAI,CAAC,CAAA;QAEtB,kCAAkC;QAClC,qDAAqD;QACrD,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,GAAG,EAAE,CAAC,EAAE,EAAE,CAAC;YAC7B,IAAI,CAAC,CAAC,CAAC,CAAC,EAAE,CAAC;gBACT,MAAM,CAAC,SAAS,CAAC,CAAC,CAAC,CAAA,CAAC,mBAAmB;gBACvC,CAAC,CAAC,CAAC,CAAC,EAAE,CAAA;YACR,CAAC;QACH,CAAC;QAED,OAAO,CAAC,MAAM,EAAE,CAAC,EAAE,IAAI,CAAC,CAAA;IAC1B,CAAC;IAED,YAAY,CAAC,MAAM,EAAE,GAAG;QACtB,MAAM,CAAC,GAAG,MAAM,CAAC,QAAQ,EAAE,CAAA;QAE3B,mBAAmB;QACnB,MAAM,IAAI,GAAG,IAAI,KAAK,CAAC,CAAC,CAAC,CAAA;QACzB,MAAM,IAAI,GAAG,IAAI,KAAK,CAAC,CAAC,CAAC,CAAA;QACzB,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC;YAC3B,IAAI,CAAC,CAAC,CAAC,GAAG,MAAM,CAAC,SAAS,EAAE,CAAA;QAC9B,CAAC;QAED,iBAAiB;QACjB,MAAM,CAAC,GAAG,IAAI,KAAK,CAAC,CAAC,CAAC,CAAA;QACtB,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC;YAC3B,IAAI,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,KAAK,CAAC,GAAG,GAAG,CAAC,CAAC,GAAG,CAAC,GAAG,GAAG,CAAC,GAAG,CAAC,CAAC,CAAA;YAC7C,CAAC,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,YAAY,CAAC,MAAM,EAAE,IAAI,CAAC,CAAC,CAAC,CAAC,CAAA;QAC3C,CAAC;QAED,YAAY;QACZ,MAAM,GAAG,GAAG,IAAI,UAAU,CAAC,GAAG,CAAC,CAAA;QAC/B,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC;YAC3B,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,IAAI,CAAC,CAAC,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC;gBACjC,GAAG,CAAC,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAA;YAC1B,CAAC;QACH,CAAC;QAED,OAAO,GAAG,CAAA;IACZ,CAAC;IAED,yEAAyE;IACzE,aAAa;IACb,SAAS,CAAC,MAAM,EAAE,GAAG;QACnB,MAAM,GAAG,GAAG,IAAI,UAAU,CAAC,GAAG,CAAC,CAAA;QAC/B,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,GAAG,EAAE,CAAC,EAAE,EAAE,CAAC;YAC7B,GAAG,CAAC,CAAC,CAAC,GAAG,MAAM,CAAC,QAAQ,EAAE,CAAA;QAC5B,CAAC;QAED,OAAO,GAAG,CAAA;IACZ,CAAC;CACF"}
@@ -1,9 +0,0 @@
1
- export default class RangeCoder {
2
- constructor(src: any);
3
- RangeStartDecode(src: any): void;
4
- RangeGetFrequency(tot_freq: any): number;
5
- RangeDecode(src: any, sym_low: any, sym_freq: any, tot_freq: any): void;
6
- RangeShiftLow(dst: any): void;
7
- RangeEncode(dst: any, sym_low: any, sym_freq: any, tot_freq: any): void;
8
- RangeFinishEncode(dst: any): void;
9
- }
@@ -1,131 +0,0 @@
1
- // @ts-nocheck
2
- /*
3
- * Copyright (c) 2019 Genome Research Ltd.
4
- * Author(s): James Bonfield
5
- *
6
- * Redistribution and use in source and binary forms, with or without
7
- * modification, are permitted provided that the following conditions are met:
8
- *
9
- * 1. Redistributions of source code must retain the above copyright notice,
10
- * this list of conditions and the following disclaimer.
11
- *
12
- * 2. Redistributions in binary form must reproduce the above
13
- * copyright notice, this list of conditions and the following
14
- * disclaimer in the documentation and/or other materials provided
15
- * with the distribution.
16
- *
17
- * 3. Neither the names Genome Research Ltd and Wellcome Trust Sanger
18
- * Institute nor the names of its contributors may be used to endorse
19
- * or promote products derived from this software without specific
20
- * prior written permission.
21
- *
22
- * THIS SOFTWARE IS PROVIDED BY GENOME RESEARCH LTD AND CONTRIBUTORS "AS
23
- * IS" AND ANY EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT LIMITED
24
- * TO, THE IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS FOR A
25
- * PARTICULAR PURPOSE ARE DISCLAIMED. IN NO EVENT SHALL GENOME RESEARCH
26
- * LTD OR CONTRIBUTORS BE LIABLE FOR ANY DIRECT, INDIRECT, INCIDENTAL,
27
- * SPECIAL, EXEMPLARY, OR CONSEQUENTIAL DAMAGES (INCLUDING, BUT NOT
28
- * LIMITED TO, PROCUREMENT OF SUBSTITUTE GOODS OR SERVICES; LOSS OF USE,
29
- * DATA, OR PROFITS; OR BUSINESS INTERRUPTION) HOWEVER CAUSED AND ON ANY
30
- * THEORY OF LIABILITY, WHETHER IN CONTRACT, STRICT LIABILITY, OR TORT
31
- * (INCLUDING NEGLIGENCE OR OTHERWISE) ARISING IN ANY WAY OUT OF THE USE
32
- * OF THIS SOFTWARE, EVEN IF ADVISED OF THE POSSIBILITY OF SUCH DAMAGE.
33
- */
34
- // An arithmetic coder, based on Eugene Shelwien's reimplementation of
35
- // Michael Schindler range coder.
36
- //
37
- // Order-0 byte stream of ~/scratch/data/q40b
38
- // C: 3.1s decode (approx same vs 32-bit and 64-bit)
39
- // Arith_sh.js 6.7s decode (32-bit with carries)
40
- // Arith.js 317.0s decode (64-bit no carries); int64 crippling it.
41
- // ----------------------------------------------------------------------
42
- // Arithmetic (range) coder
43
- export default class RangeCoder {
44
- constructor(src) {
45
- this.low = 0;
46
- this.range = 0xffffffff;
47
- this.code = 0;
48
- this.FFnum = 0;
49
- this.carry = 0;
50
- this.cache = 0;
51
- }
52
- RangeStartDecode(src) {
53
- for (let i = 0; i < 5; i++) {
54
- this.code = (this.code << 8) + src.ReadByte();
55
- }
56
- this.code &= 0xffffffff;
57
- this.code >>>= 0; // force to be +ve int
58
- }
59
- RangeGetFrequency(tot_freq) {
60
- this.range = Math.floor(this.range / tot_freq);
61
- // return this.code / this.range;
62
- return Math.floor(this.code / this.range);
63
- // Conceptual scenario; return freq only and don't modify range yet
64
- // return Math.floor(this.code / (Math.floor(this.range / tot_freq)));
65
- }
66
- RangeDecode(src, sym_low, sym_freq, tot_freq) {
67
- // Conceptually we divide range here, but in practice we cached it earlier
68
- // this.range = Math.floor(this.range / tot_freq);
69
- this.code -= sym_low * this.range;
70
- this.range *= sym_freq;
71
- while (this.range < 1 << 24) {
72
- this.range *= 256;
73
- this.code = this.code * 256 + src.ReadByte();
74
- }
75
- }
76
- RangeShiftLow(dst) {
77
- // We know range is < (1<<24) as we got here. We already have a
78
- // cached copy of 8 bits from low. Is this correct, or does it need
79
- // fixing? Possible scenarios.
80
- // 1. Low < 0xff000000 thus low+range < 0xffffffff and cache
81
- // cannot possibly change. Output cache and as many ffs as needed.
82
- // 2. We already detected an overflow in RangeEncode, setting carry.
83
- // In this case output cached byte + 1 and any 00s needed.
84
- // 3. Neither case - range is low but we haven't yet detected if we're
85
- // XXffffff or XY000000 scenario. Increase counter for ff/00s.
86
- if ((this.low < 0xff000000) | this.carry) {
87
- // cached byte if no overflow, byte+1 otherwise
88
- dst.WriteByte(this.cache + this.carry);
89
- // Flush any tracked FFs (no carry) or 00s (carry).
90
- while (this.FFnum) {
91
- dst.WriteByte(this.carry - 1);
92
- this.FFnum--;
93
- }
94
- // Take a copy of top byte ready for next flush
95
- this.cache = this.low >>> 24;
96
- this.carry = 0;
97
- }
98
- else {
99
- this.FFnum++; // keep track of number of trailing ff/00 bytes to write
100
- }
101
- this.low <<= 8;
102
- this.low >>>= 0; // force to be +ve int
103
- }
104
- RangeEncode(dst, sym_low, sym_freq, tot_freq) {
105
- const old_low = this.low;
106
- this.range = Math.floor(this.range / tot_freq);
107
- this.low += sym_low * this.range;
108
- this.low >>>= 0; // Truncate to +ve int so we can spot overflow
109
- this.range *= sym_freq;
110
- // "low + sym*range < old_low" means we overflow; set carry.
111
- // NB: can this.low < old_low occur twice before range < (1<<24)?
112
- // We claim not, but prove it!
113
- if (this.low < old_low) {
114
- if (this.carry != 0) {
115
- console.log('ERROR: Multiple carry');
116
- }
117
- this.carry = 1;
118
- }
119
- // Renormalise if range gets too small
120
- while (this.range < 1 << 24) {
121
- this.range *= 256;
122
- this.RangeShiftLow(dst);
123
- }
124
- }
125
- RangeFinishEncode(dst) {
126
- for (let i = 0; i < 5; i++) {
127
- this.RangeShiftLow(dst);
128
- }
129
- }
130
- }
131
- //# sourceMappingURL=arith_sh.js.map
@@ -1 +0,0 @@
1
- {"version":3,"file":"arith_sh.js","sourceRoot":"","sources":["../../src/htscodecs/arith_sh.ts"],"names":[],"mappings":"AAAA,cAAc;AAEd;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;GA+BG;AAEH,sEAAsE;AACtE,iCAAiC;AACjC,EAAE;AACF,6CAA6C;AAC7C,kEAAkE;AAClE,qDAAqD;AACrD,wEAAwE;AAExE,yEAAyE;AACzE,2BAA2B;AAC3B,MAAM,CAAC,OAAO,OAAO,UAAU;IAC7B,YAAY,GAAG;QACb,IAAI,CAAC,GAAG,GAAG,CAAC,CAAA;QACZ,IAAI,CAAC,KAAK,GAAG,UAAU,CAAA;QACvB,IAAI,CAAC,IAAI,GAAG,CAAC,CAAA;QACb,IAAI,CAAC,KAAK,GAAG,CAAC,CAAA;QACd,IAAI,CAAC,KAAK,GAAG,CAAC,CAAA;QACd,IAAI,CAAC,KAAK,GAAG,CAAC,CAAA;IAChB,CAAC;IAED,gBAAgB,CAAC,GAAG;QAClB,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC;YAC3B,IAAI,CAAC,IAAI,GAAG,CAAC,IAAI,CAAC,IAAI,IAAI,CAAC,CAAC,GAAG,GAAG,CAAC,QAAQ,EAAE,CAAA;QAC/C,CAAC;QACD,IAAI,CAAC,IAAI,IAAI,UAAU,CAAA;QACvB,IAAI,CAAC,IAAI,MAAM,CAAC,CAAA,CAAC,sBAAsB;IACzC,CAAC;IAED,iBAAiB,CAAC,QAAQ;QACxB,IAAI,CAAC,KAAK,GAAG,IAAI,CAAC,KAAK,CAAC,IAAI,CAAC,KAAK,GAAG,QAAQ,CAAC,CAAA;QAC9C,iCAAiC;QACjC,OAAO,IAAI,CAAC,KAAK,CAAC,IAAI,CAAC,IAAI,GAAG,IAAI,CAAC,KAAK,CAAC,CAAA;QAEzC,mEAAmE;QACnE,sEAAsE;IACxE,CAAC;IAED,WAAW,CAAC,GAAG,EAAE,OAAO,EAAE,QAAQ,EAAE,QAAQ;QAC1C,0EAA0E;QAC1E,kDAAkD;QAElD,IAAI,CAAC,IAAI,IAAI,OAAO,GAAG,IAAI,CAAC,KAAK,CAAA;QACjC,IAAI,CAAC,KAAK,IAAI,QAAQ,CAAA;QAEtB,OAAO,IAAI,CAAC,KAAK,GAAG,CAAC,IAAI,EAAE,EAAE,CAAC;YAC5B,IAAI,CAAC,KAAK,IAAI,GAAG,CAAA;YACjB,IAAI,CAAC,IAAI,GAAG,IAAI,CAAC,IAAI,GAAG,GAAG,GAAG,GAAG,CAAC,QAAQ,EAAE,CAAA;QAC9C,CAAC;IACH,CAAC;IAED,aAAa,CAAC,GAAG;QACf,gEAAgE;QAChE,oEAAoE;QACpE,+BAA+B;QAC/B,4DAA4D;QAC5D,sEAAsE;QACtE,oEAAoE;QACpE,6DAA6D;QAC7D,sEAAsE;QACtE,kEAAkE;QAElE,IAAI,CAAC,IAAI,CAAC,GAAG,GAAG,UAAU,CAAC,GAAG,IAAI,CAAC,KAAK,EAAE,CAAC;YACzC,+CAA+C;YAC/C,GAAG,CAAC,SAAS,CAAC,IAAI,CAAC,KAAK,GAAG,IAAI,CAAC,KAAK,CAAC,CAAA;YAEtC,mDAAmD;YACnD,OAAO,IAAI,CAAC,KAAK,EAAE,CAAC;gBAClB,GAAG,CAAC,SAAS,CAAC,IAAI,CAAC,KAAK,GAAG,CAAC,CAAC,CAAA;gBAC7B,IAAI,CAAC,KAAK,EAAE,CAAA;YACd,CAAC;YAED,+CAA+C;YAC/C,IAAI,CAAC,KAAK,GAAG,IAAI,CAAC,GAAG,KAAK,EAAE,CAAA;YAC5B,IAAI,CAAC,KAAK,GAAG,CAAC,CAAA;QAChB,CAAC;aAAM,CAAC;YACN,IAAI,CAAC,KAAK,EAAE,CAAA,CAAC,wDAAwD;QACvE,CAAC;QACD,IAAI,CAAC,GAAG,KAAK,CAAC,CAAA;QACd,IAAI,CAAC,GAAG,MAAM,CAAC,CAAA,CAAC,sBAAsB;IACxC,CAAC;IAED,WAAW,CAAC,GAAG,EAAE,OAAO,EAAE,QAAQ,EAAE,QAAQ;QAC1C,MAAM,OAAO,GAAG,IAAI,CAAC,GAAG,CAAA;QACxB,IAAI,CAAC,KAAK,GAAG,IAAI,CAAC,KAAK,CAAC,IAAI,CAAC,KAAK,GAAG,QAAQ,CAAC,CAAA;QAC9C,IAAI,CAAC,GAAG,IAAI,OAAO,GAAG,IAAI,CAAC,KAAK,CAAA;QAChC,IAAI,CAAC,GAAG,MAAM,CAAC,CAAA,CAAC,8CAA8C;QAC9D,IAAI,CAAC,KAAK,IAAI,QAAQ,CAAA;QAEtB,4DAA4D;QAC5D,iEAAiE;QACjE,8BAA8B;QAC9B,IAAI,IAAI,CAAC,GAAG,GAAG,OAAO,EAAE,CAAC;YACvB,IAAI,IAAI,CAAC,KAAK,IAAI,CAAC,EAAE,CAAC;gBACpB,OAAO,CAAC,GAAG,CAAC,uBAAuB,CAAC,CAAA;YACtC,CAAC;YACD,IAAI,CAAC,KAAK,GAAG,CAAC,CAAA;QAChB,CAAC;QAED,sCAAsC;QACtC,OAAO,IAAI,CAAC,KAAK,GAAG,CAAC,IAAI,EAAE,EAAE,CAAC;YAC5B,IAAI,CAAC,KAAK,IAAI,GAAG,CAAA;YACjB,IAAI,CAAC,aAAa,CAAC,GAAG,CAAC,CAAA;QACzB,CAAC;IACH,CAAC;IAED,iBAAiB,CAAC,GAAG;QACnB,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC;YAC3B,IAAI,CAAC,aAAa,CAAC,GAAG,CAAC,CAAA;QACzB,CAAC;IACH,CAAC;CACF"}
@@ -1,6 +0,0 @@
1
- export default class ByteModel {
2
- constructor(max_sym?: number);
3
- ModelDecode(src: any, rc: any): any;
4
- ModelRenormalise(): void;
5
- ModelEncode(dst: any, rc: any, sym: any): void;
6
- }
@@ -1,118 +0,0 @@
1
- /* eslint-disable no-var */
2
- // @ts-nocheck
3
- /*
4
- * Copyright (c) 2019 Genome Research Ltd.
5
- * Author(s): James Bonfield
6
- *
7
- * Redistribution and use in source and binary forms, with or without
8
- * modification, are permitted provided that the following conditions are met:
9
- *
10
- * 1. Redistributions of source code must retain the above copyright notice,
11
- * this list of conditions and the following disclaimer.
12
- *
13
- * 2. Redistributions in binary form must reproduce the above
14
- * copyright notice, this list of conditions and the following
15
- * disclaimer in the documentation and/or other materials provided
16
- * with the distribution.
17
- *
18
- * 3. Neither the names Genome Research Ltd and Wellcome Trust Sanger
19
- * Institute nor the names of its contributors may be used to endorse
20
- * or promote products derived from this software without specific
21
- * prior written permission.
22
- *
23
- * THIS SOFTWARE IS PROVIDED BY GENOME RESEARCH LTD AND CONTRIBUTORS "AS
24
- * IS" AND ANY EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT LIMITED
25
- * TO, THE IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS FOR A
26
- * PARTICULAR PURPOSE ARE DISCLAIMED. IN NO EVENT SHALL GENOME RESEARCH
27
- * LTD OR CONTRIBUTORS BE LIABLE FOR ANY DIRECT, INDIRECT, INCIDENTAL,
28
- * SPECIAL, EXEMPLARY, OR CONSEQUENTIAL DAMAGES (INCLUDING, BUT NOT
29
- * LIMITED TO, PROCUREMENT OF SUBSTITUTE GOODS OR SERVICES; LOSS OF USE,
30
- * DATA, OR PROFITS; OR BUSINESS INTERRUPTION) HOWEVER CAUSED AND ON ANY
31
- * THEORY OF LIABILITY, WHETHER IN CONTRACT, STRICT LIABILITY, OR TORT
32
- * (INCLUDING NEGLIGENCE OR OTHERWISE) ARISING IN ANY WAY OUT OF THE USE
33
- * OF THIS SOFTWARE, EVEN IF ADVISED OF THE POSSIBILITY OF SUCH DAMAGE.
34
- */
35
- // An adaptive probability model for encoding and decoding of symbols
36
- // within a given alphabet, using the range coder to get/put the
37
- // compressed data.
38
- const MAX_FREQ = (1 << 16) - 17;
39
- const STEP = 16;
40
- export default class ByteModel {
41
- constructor(max_sym = 256) {
42
- this.total_freq = max_sym;
43
- this.max_sym = max_sym - 1;
44
- this.S = [];
45
- this.F = [];
46
- for (let i = 0; i <= this.max_sym; i++) {
47
- this.S[i] = i;
48
- this.F[i] = 1;
49
- }
50
- }
51
- ModelDecode(src, rc) {
52
- // Find symbol
53
- const freq = rc.RangeGetFrequency(this.total_freq);
54
- // Linear scan to find cumulative frequency 'freq'
55
- let acc = 0;
56
- let x = 0;
57
- while (acc + this.F[x] <= freq) {
58
- acc += this.F[x++];
59
- }
60
- // for (var acc = 0; (acc += this.F[x]) <= freq; x++)
61
- // ;
62
- // acc -= this.F[x];
63
- // Update range coder
64
- rc.RangeDecode(src, acc, this.F[x], this.total_freq);
65
- // Update model
66
- this.F[x] += STEP;
67
- this.total_freq += STEP;
68
- if (this.total_freq > MAX_FREQ) {
69
- this.ModelRenormalise();
70
- }
71
- // Keep symbols approximately frequency sorted
72
- const sym = this.S[x];
73
- if (x > 0 && this.F[x] > this.F[x - 1]) {
74
- let tmp = this.F[x];
75
- this.F[x] = this.F[x - 1];
76
- this.F[x - 1] = tmp;
77
- tmp = this.S[x];
78
- this.S[x] = this.S[x - 1];
79
- this.S[x - 1] = tmp;
80
- }
81
- return sym;
82
- }
83
- ModelRenormalise() {
84
- // Halve all the frequencies, being careful not to hit zero
85
- this.total_freq = 0;
86
- for (let i = 0; i <= this.max_sym; i++) {
87
- this.F[i] -= Math.floor(this.F[i] / 2);
88
- this.total_freq += this.F[i];
89
- }
90
- }
91
- ModelEncode(dst, rc, sym) {
92
- // Find cumulative frequency
93
- let acc = 0;
94
- for (var x = 0; this.S[x] != sym; x++) {
95
- acc += this.F[x];
96
- }
97
- // Encode
98
- rc.RangeEncode(dst, acc, this.F[x], this.total_freq);
99
- // Update model
100
- this.F[x] += STEP;
101
- this.total_freq += STEP;
102
- if (this.total_freq > MAX_FREQ) {
103
- // FIXME x2
104
- this.ModelRenormalise();
105
- }
106
- // Keep symbols approximately frequency sorted
107
- var sym = this.S[x];
108
- if (x > 0 && this.F[x] > this.F[x - 1]) {
109
- let tmp = this.F[x];
110
- this.F[x] = this.F[x - 1];
111
- this.F[x - 1] = tmp;
112
- tmp = this.S[x];
113
- this.S[x] = this.S[x - 1];
114
- this.S[x - 1] = tmp;
115
- }
116
- }
117
- }
118
- //# sourceMappingURL=byte_model.js.map