@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.
- package/dist/craiIndex.js +2 -2
- package/dist/craiIndex.js.map +1 -1
- package/dist/cram-bundle.js +1 -1
- package/dist/cramFile/file.d.ts +2 -2
- package/dist/cramFile/file.js +9 -13
- package/dist/cramFile/file.js.map +1 -1
- package/dist/htscodecs/index.d.ts +1 -5
- package/dist/htscodecs/index.js +11 -72
- package/dist/htscodecs/index.js.map +1 -1
- package/dist/htscodecs-wasm.d.ts +7 -0
- package/dist/htscodecs-wasm.js +183 -0
- package/dist/htscodecs-wasm.js.map +1 -0
- package/dist/unzip.d.ts +1 -1
- package/dist/unzip.js +3 -3
- package/dist/unzip.js.map +1 -1
- package/dist/wasm/htscodecs.d.ts +5 -0
- package/dist/wasm/htscodecs.js +0 -0
- package/dist/wasm/htscodecs.js.map +1 -0
- package/esm/craiIndex.js +2 -2
- package/esm/craiIndex.js.map +1 -1
- package/esm/cramFile/file.d.ts +2 -2
- package/esm/cramFile/file.js +7 -11
- package/esm/cramFile/file.js.map +1 -1
- package/esm/htscodecs/index.d.ts +1 -5
- package/esm/htscodecs/index.js +3 -31
- package/esm/htscodecs/index.js.map +1 -1
- package/esm/htscodecs-wasm.d.ts +7 -0
- package/esm/htscodecs-wasm.js +171 -0
- package/esm/htscodecs-wasm.js.map +1 -0
- package/esm/unzip.d.ts +1 -1
- package/esm/unzip.js +3 -3
- package/esm/unzip.js.map +1 -1
- package/esm/wasm/htscodecs.d.ts +5 -0
- package/esm/wasm/htscodecs.js +0 -0
- package/esm/wasm/htscodecs.js.map +1 -0
- package/package.json +9 -9
- package/src/craiIndex.ts +4 -2
- package/src/cramFile/file.ts +7 -11
- package/src/htscodecs/index.ts +12 -38
- package/src/htscodecs-wasm.ts +210 -0
- package/src/unzip.ts +3 -3
- package/src/wasm/htscodecs.d.ts +58 -0
- package/src/wasm/htscodecs.js +0 -0
- package/dist/htscodecs/arith_gen.d.ts +0 -16
- package/dist/htscodecs/arith_gen.js +0 -323
- package/dist/htscodecs/arith_gen.js.map +0 -1
- package/dist/htscodecs/arith_sh.d.ts +0 -9
- package/dist/htscodecs/arith_sh.js +0 -134
- package/dist/htscodecs/arith_sh.js.map +0 -1
- package/dist/htscodecs/byte_model.d.ts +0 -6
- package/dist/htscodecs/byte_model.js +0 -121
- package/dist/htscodecs/byte_model.js.map +0 -1
- package/dist/htscodecs/fqzcomp.d.ts +0 -1
- package/dist/htscodecs/fqzcomp.js +0 -346
- package/dist/htscodecs/fqzcomp.js.map +0 -1
- package/dist/htscodecs/iostream.d.ts +0 -16
- package/dist/htscodecs/iostream.js +0 -148
- package/dist/htscodecs/iostream.js.map +0 -1
- package/dist/htscodecs/rans.d.ts +0 -1
- package/dist/htscodecs/rans.js +0 -224
- package/dist/htscodecs/rans.js.map +0 -1
- package/dist/htscodecs/rans4x16.d.ts +0 -1
- package/dist/htscodecs/rans4x16.js +0 -445
- package/dist/htscodecs/rans4x16.js.map +0 -1
- package/dist/htscodecs/tok3.d.ts +0 -1
- package/dist/htscodecs/tok3.js +0 -206
- package/dist/htscodecs/tok3.js.map +0 -1
- package/dist/rans/constants.d.ts +0 -4
- package/dist/rans/constants.js +0 -10
- package/dist/rans/constants.js.map +0 -1
- package/dist/rans/d04.d.ts +0 -1
- package/dist/rans/d04.js +0 -100
- package/dist/rans/d04.js.map +0 -1
- package/dist/rans/d14.d.ts +0 -1
- package/dist/rans/d14.js +0 -96
- package/dist/rans/d14.js.map +0 -1
- package/dist/rans/decoding.d.ts +0 -27
- package/dist/rans/decoding.js +0 -119
- package/dist/rans/decoding.js.map +0 -1
- package/dist/rans/frequencies.d.ts +0 -2
- package/dist/rans/frequencies.js +0 -112
- package/dist/rans/frequencies.js.map +0 -1
- package/dist/rans/index.d.ts +0 -1
- package/dist/rans/index.js +0 -198
- package/dist/rans/index.js.map +0 -1
- package/dist/seek-bzip/bitreader.d.ts +0 -27
- package/dist/seek-bzip/bitreader.js +0 -101
- package/dist/seek-bzip/bitreader.js.map +0 -1
- package/dist/seek-bzip/crc32.d.ts +0 -26
- package/dist/seek-bzip/crc32.js +0 -115
- package/dist/seek-bzip/crc32.js.map +0 -1
- package/dist/seek-bzip/index.d.ts +0 -9
- package/dist/seek-bzip/index.js +0 -585
- package/dist/seek-bzip/index.js.map +0 -1
- package/dist/seek-bzip/stream.d.ts +0 -8
- package/dist/seek-bzip/stream.js +0 -37
- package/dist/seek-bzip/stream.js.map +0 -1
- package/dist/seek-bzip/toHex.d.ts +0 -1
- package/dist/seek-bzip/toHex.js +0 -36
- package/dist/seek-bzip/toHex.js.map +0 -1
- package/esm/htscodecs/arith_gen.d.ts +0 -16
- package/esm/htscodecs/arith_gen.js +0 -317
- package/esm/htscodecs/arith_gen.js.map +0 -1
- package/esm/htscodecs/arith_sh.d.ts +0 -9
- package/esm/htscodecs/arith_sh.js +0 -131
- package/esm/htscodecs/arith_sh.js.map +0 -1
- package/esm/htscodecs/byte_model.d.ts +0 -6
- package/esm/htscodecs/byte_model.js +0 -118
- package/esm/htscodecs/byte_model.js.map +0 -1
- package/esm/htscodecs/fqzcomp.d.ts +0 -1
- package/esm/htscodecs/fqzcomp.js +0 -340
- package/esm/htscodecs/fqzcomp.js.map +0 -1
- package/esm/htscodecs/iostream.d.ts +0 -16
- package/esm/htscodecs/iostream.js +0 -145
- package/esm/htscodecs/iostream.js.map +0 -1
- package/esm/htscodecs/rans.d.ts +0 -1
- package/esm/htscodecs/rans.js +0 -218
- package/esm/htscodecs/rans.js.map +0 -1
- package/esm/htscodecs/rans4x16.d.ts +0 -1
- package/esm/htscodecs/rans4x16.js +0 -439
- package/esm/htscodecs/rans4x16.js.map +0 -1
- package/esm/htscodecs/tok3.d.ts +0 -1
- package/esm/htscodecs/tok3.js +0 -167
- package/esm/htscodecs/tok3.js.map +0 -1
- package/esm/rans/constants.d.ts +0 -4
- package/esm/rans/constants.js +0 -5
- package/esm/rans/constants.js.map +0 -1
- package/esm/rans/d04.d.ts +0 -1
- package/esm/rans/d04.js +0 -94
- package/esm/rans/d04.js.map +0 -1
- package/esm/rans/d14.d.ts +0 -1
- package/esm/rans/d14.js +0 -93
- package/esm/rans/d14.js.map +0 -1
- package/esm/rans/decoding.d.ts +0 -27
- package/esm/rans/decoding.js +0 -117
- package/esm/rans/decoding.js.map +0 -1
- package/esm/rans/frequencies.d.ts +0 -2
- package/esm/rans/frequencies.js +0 -105
- package/esm/rans/frequencies.js.map +0 -1
- package/esm/rans/index.d.ts +0 -1
- package/esm/rans/index.js +0 -192
- package/esm/rans/index.js.map +0 -1
- package/esm/seek-bzip/bitreader.d.ts +0 -27
- package/esm/seek-bzip/bitreader.js +0 -98
- package/esm/seek-bzip/bitreader.js.map +0 -1
- package/esm/seek-bzip/crc32.d.ts +0 -26
- package/esm/seek-bzip/crc32.js +0 -112
- package/esm/seek-bzip/crc32.js.map +0 -1
- package/esm/seek-bzip/index.d.ts +0 -9
- package/esm/seek-bzip/index.js +0 -579
- package/esm/seek-bzip/index.js.map +0 -1
- package/esm/seek-bzip/stream.d.ts +0 -8
- package/esm/seek-bzip/stream.js +0 -34
- package/esm/seek-bzip/stream.js.map +0 -1
- package/esm/seek-bzip/toHex.d.ts +0 -1
- package/esm/seek-bzip/toHex.js +0 -33
- package/esm/seek-bzip/toHex.js.map +0 -1
- package/src/htscodecs/arith_gen.ts +0 -371
- package/src/htscodecs/arith_sh.ts +0 -146
- package/src/htscodecs/byte_model.ts +0 -136
- package/src/htscodecs/fqzcomp.ts +0 -394
- package/src/htscodecs/iostream.ts +0 -159
- package/src/htscodecs/rans.ts +0 -250
- package/src/htscodecs/rans4x16.ts +0 -521
- package/src/htscodecs/tok3.ts +0 -197
- package/src/rans/README.md +0 -1
- package/src/rans/constants.ts +0 -5
- package/src/rans/d04.ts +0 -112
- package/src/rans/d14.ts +0 -106
- package/src/rans/decoding.ts +0 -138
- package/src/rans/frequencies.ts +0 -115
- package/src/rans/index.ts +0 -245
- package/src/seek-bzip/README.md +0 -6
- package/src/seek-bzip/bitreader.ts +0 -113
- package/src/seek-bzip/crc32.ts +0 -121
- package/src/seek-bzip/index.ts +0 -615
- package/src/seek-bzip/stream.ts +0 -33
- package/src/seek-bzip/toHex.ts +0 -33
package/dist/seek-bzip/toHex.js
DELETED
|
@@ -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,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
|