@gmod/cram 3.0.6 → 3.0.7
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/cram-bundle.js +1 -1
- package/dist/cramFile/codecs/byteArrayStop.d.ts +2 -2
- package/dist/cramFile/file.js +8 -10
- package/dist/cramFile/file.js.map +1 -1
- package/dist/cramFile/util.js.map +1 -1
- package/dist/htscodecs/arith_gen.d.ts +26 -0
- package/dist/htscodecs/arith_gen.js +558 -0
- package/dist/htscodecs/arith_gen.js.map +1 -0
- package/dist/htscodecs/arith_sh.d.ts +16 -0
- package/dist/htscodecs/arith_sh.js +128 -0
- package/dist/htscodecs/arith_sh.js.map +1 -0
- package/dist/htscodecs/byte_model.d.ts +11 -0
- package/dist/htscodecs/byte_model.js +113 -0
- package/dist/htscodecs/byte_model.js.map +1 -0
- package/dist/htscodecs/fqzcomp.d.ts +2 -0
- package/dist/htscodecs/fqzcomp.js +744 -0
- package/dist/htscodecs/fqzcomp.js.map +1 -0
- package/dist/htscodecs/index.d.ts +5 -0
- package/dist/htscodecs/index.js +70 -0
- package/dist/htscodecs/index.js.map +1 -0
- package/dist/htscodecs/iostream.d.ts +26 -0
- package/dist/htscodecs/iostream.js +242 -0
- package/dist/htscodecs/iostream.js.map +1 -0
- package/dist/htscodecs/main_arith_gen.d.ts +1 -0
- package/dist/htscodecs/main_arith_gen.js +86 -0
- package/dist/htscodecs/main_arith_gen.js.map +1 -0
- package/dist/htscodecs/main_fqzcomp.d.ts +1 -0
- package/dist/htscodecs/main_fqzcomp.js +112 -0
- package/dist/htscodecs/main_fqzcomp.js.map +1 -0
- package/dist/htscodecs/main_rans.d.ts +1 -0
- package/dist/htscodecs/main_rans.js +83 -0
- package/dist/htscodecs/main_rans.js.map +1 -0
- package/dist/htscodecs/main_rans4x16.d.ts +1 -0
- package/dist/htscodecs/main_rans4x16.js +82 -0
- package/dist/htscodecs/main_rans4x16.js.map +1 -0
- package/dist/htscodecs/main_tok3.d.ts +1 -0
- package/dist/htscodecs/main_tok3.js +84 -0
- package/dist/htscodecs/main_tok3.js.map +1 -0
- package/dist/htscodecs/rans.d.ts +2 -0
- package/dist/htscodecs/rans.js +480 -0
- package/dist/htscodecs/rans.js.map +1 -0
- package/dist/htscodecs/rans4x16.d.ts +2 -0
- package/dist/htscodecs/rans4x16.js +896 -0
- package/dist/htscodecs/rans4x16.js.map +1 -0
- package/dist/htscodecs/tok3.d.ts +2 -0
- package/dist/htscodecs/tok3.js +347 -0
- package/dist/htscodecs/tok3.js.map +1 -0
- package/esm/cramFile/codecs/byteArrayStop.d.ts +2 -2
- package/esm/cramFile/file.js +8 -10
- package/esm/cramFile/file.js.map +1 -1
- package/esm/cramFile/util.js.map +1 -1
- package/esm/htscodecs/arith_gen.d.ts +26 -0
- package/esm/htscodecs/arith_gen.js +558 -0
- package/esm/htscodecs/arith_gen.js.map +1 -0
- package/esm/htscodecs/arith_sh.d.ts +16 -0
- package/esm/htscodecs/arith_sh.js +128 -0
- package/esm/htscodecs/arith_sh.js.map +1 -0
- package/esm/htscodecs/byte_model.d.ts +11 -0
- package/esm/htscodecs/byte_model.js +113 -0
- package/esm/htscodecs/byte_model.js.map +1 -0
- package/esm/htscodecs/fqzcomp.d.ts +2 -0
- package/esm/htscodecs/fqzcomp.js +744 -0
- package/esm/htscodecs/fqzcomp.js.map +1 -0
- package/esm/htscodecs/index.d.ts +5 -0
- package/esm/htscodecs/index.js +70 -0
- package/esm/htscodecs/index.js.map +1 -0
- package/esm/htscodecs/iostream.d.ts +26 -0
- package/esm/htscodecs/iostream.js +242 -0
- package/esm/htscodecs/iostream.js.map +1 -0
- package/esm/htscodecs/main_arith_gen.d.ts +1 -0
- package/esm/htscodecs/main_arith_gen.js +86 -0
- package/esm/htscodecs/main_arith_gen.js.map +1 -0
- package/esm/htscodecs/main_fqzcomp.d.ts +1 -0
- package/esm/htscodecs/main_fqzcomp.js +112 -0
- package/esm/htscodecs/main_fqzcomp.js.map +1 -0
- package/esm/htscodecs/main_rans.d.ts +1 -0
- package/esm/htscodecs/main_rans.js +83 -0
- package/esm/htscodecs/main_rans.js.map +1 -0
- package/esm/htscodecs/main_rans4x16.d.ts +1 -0
- package/esm/htscodecs/main_rans4x16.js +82 -0
- package/esm/htscodecs/main_rans4x16.js.map +1 -0
- package/esm/htscodecs/main_tok3.d.ts +1 -0
- package/esm/htscodecs/main_tok3.js +84 -0
- package/esm/htscodecs/main_tok3.js.map +1 -0
- package/esm/htscodecs/rans.d.ts +2 -0
- package/esm/htscodecs/rans.js +480 -0
- package/esm/htscodecs/rans.js.map +1 -0
- package/esm/htscodecs/rans4x16.d.ts +2 -0
- package/esm/htscodecs/rans4x16.js +896 -0
- package/esm/htscodecs/rans4x16.js.map +1 -0
- package/esm/htscodecs/tok3.d.ts +2 -0
- package/esm/htscodecs/tok3.js +347 -0
- package/esm/htscodecs/tok3.js.map +1 -0
- package/package.json +1 -2
- package/src/cramFile/file.ts +8 -10
- package/src/cramFile/util.ts +0 -1
- package/src/htscodecs/Makefile +142 -0
- package/src/htscodecs/README.md +64 -0
- package/src/htscodecs/arith_gen.js +607 -0
- package/src/htscodecs/arith_sh.js +138 -0
- package/src/htscodecs/byte_model.js +126 -0
- package/src/htscodecs/fqzcomp.js +834 -0
- package/src/htscodecs/index.js +79 -0
- package/src/htscodecs/iostream.js +256 -0
- package/src/htscodecs/main_arith_gen.js +96 -0
- package/src/htscodecs/main_fqzcomp.js +113 -0
- package/src/htscodecs/main_rans.js +88 -0
- package/src/htscodecs/main_rans4x16.js +87 -0
- package/src/htscodecs/main_tok3.js +86 -0
- package/src/htscodecs/rans.js +545 -0
- package/src/htscodecs/rans4x16.js +1003 -0
- package/src/htscodecs/tok3.js +396 -0
- package/src/typings/htscodecs.d.ts +0 -6
|
@@ -0,0 +1,128 @@
|
|
|
1
|
+
"use strict";
|
|
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
|
+
module.exports = 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 (var i = 0; i < 5; i++)
|
|
54
|
+
this.code = (this.code << 8) + src.ReadByte();
|
|
55
|
+
this.code &= 0xffffffff;
|
|
56
|
+
this.code >>>= 0; // force to be +ve int
|
|
57
|
+
}
|
|
58
|
+
RangeGetFrequency(tot_freq) {
|
|
59
|
+
this.range = Math.floor(this.range / tot_freq);
|
|
60
|
+
//return this.code / this.range;
|
|
61
|
+
return Math.floor(this.code / this.range);
|
|
62
|
+
// Conceptual scenario; return freq only and don't modify range yet
|
|
63
|
+
//return Math.floor(this.code / (Math.floor(this.range / tot_freq)));
|
|
64
|
+
}
|
|
65
|
+
RangeDecode(src, sym_low, sym_freq, tot_freq) {
|
|
66
|
+
// Conceptually we divide range here, but in practice we cached it earlier
|
|
67
|
+
//this.range = Math.floor(this.range / tot_freq);
|
|
68
|
+
this.code -= sym_low * this.range;
|
|
69
|
+
this.range *= sym_freq;
|
|
70
|
+
while (this.range < 1 << 24) {
|
|
71
|
+
this.range *= 256;
|
|
72
|
+
this.code = this.code * 256 + src.ReadByte();
|
|
73
|
+
}
|
|
74
|
+
}
|
|
75
|
+
RangeShiftLow(dst) {
|
|
76
|
+
// We know range is < (1<<24) as we got here. We already have a
|
|
77
|
+
// cached copy of 8 bits from low. Is this correct, or does it need
|
|
78
|
+
// fixing? Possible scenarios.
|
|
79
|
+
// 1. Low < 0xff000000 thus low+range < 0xffffffff and cache
|
|
80
|
+
// cannot possibly change. Output cache and as many ffs as needed.
|
|
81
|
+
// 2. We already detected an overflow in RangeEncode, setting carry.
|
|
82
|
+
// In this case output cached byte + 1 and any 00s needed.
|
|
83
|
+
// 3. Neither case - range is low but we haven't yet detected if we're
|
|
84
|
+
// XXffffff or XY000000 scenario. Increase counter for ff/00s.
|
|
85
|
+
if ((this.low < 0xff000000) | this.carry) {
|
|
86
|
+
// cached byte if no overflow, byte+1 otherwise
|
|
87
|
+
dst.WriteByte(this.cache + this.carry);
|
|
88
|
+
// Flush any tracked FFs (no carry) or 00s (carry).
|
|
89
|
+
while (this.FFnum) {
|
|
90
|
+
dst.WriteByte(this.carry - 1);
|
|
91
|
+
this.FFnum--;
|
|
92
|
+
}
|
|
93
|
+
// Take a copy of top byte ready for next flush
|
|
94
|
+
this.cache = this.low >>> 24;
|
|
95
|
+
this.carry = 0;
|
|
96
|
+
}
|
|
97
|
+
else {
|
|
98
|
+
this.FFnum++; // keep track of number of trailing ff/00 bytes to write
|
|
99
|
+
}
|
|
100
|
+
this.low <<= 8;
|
|
101
|
+
this.low >>>= 0; // force to be +ve int
|
|
102
|
+
}
|
|
103
|
+
RangeEncode(dst, sym_low, sym_freq, tot_freq) {
|
|
104
|
+
var old_low = this.low;
|
|
105
|
+
this.range = Math.floor(this.range / tot_freq);
|
|
106
|
+
this.low += sym_low * this.range;
|
|
107
|
+
this.low >>>= 0; // Truncate to +ve int so we can spot overflow
|
|
108
|
+
this.range *= sym_freq;
|
|
109
|
+
// "low + sym*range < old_low" means we overflow; set carry.
|
|
110
|
+
// NB: can this.low < old_low occur twice before range < (1<<24)?
|
|
111
|
+
// We claim not, but prove it!
|
|
112
|
+
if (this.low < old_low) {
|
|
113
|
+
if (this.carry != 0)
|
|
114
|
+
console.log('ERROR: Multiple carry');
|
|
115
|
+
this.carry = 1;
|
|
116
|
+
}
|
|
117
|
+
// Renormalise if range gets too small
|
|
118
|
+
while (this.range < 1 << 24) {
|
|
119
|
+
this.range *= 256;
|
|
120
|
+
this.RangeShiftLow(dst);
|
|
121
|
+
}
|
|
122
|
+
}
|
|
123
|
+
RangeFinishEncode(dst) {
|
|
124
|
+
for (var i = 0; i < 5; i++)
|
|
125
|
+
this.RangeShiftLow(dst);
|
|
126
|
+
}
|
|
127
|
+
};
|
|
128
|
+
//# sourceMappingURL=arith_sh.js.map
|
|
@@ -0,0 +1 @@
|
|
|
1
|
+
{"version":3,"file":"arith_sh.js","sourceRoot":"","sources":["../../src/htscodecs/arith_sh.js"],"names":[],"mappings":";AAAA;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;GA+BG;AAEH,sEAAsE;AACtE,iCAAiC;AACjC,EAAE;AACF,6CAA6C;AAC7C,kEAAkE;AAClE,qDAAqD;AACrD,wEAAwE;AAExE,wEAAwE;AACxE,2BAA2B;AAC3B,MAAM,CAAC,OAAO,GAAG,MAAM,UAAU;IAC/B,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;YAAE,IAAI,CAAC,IAAI,GAAG,CAAC,IAAI,CAAC,IAAI,IAAI,CAAC,CAAC,GAAG,GAAG,CAAC,QAAQ,EAAE,CAAA;QACzE,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,gCAAgC;QAChC,OAAO,IAAI,CAAC,KAAK,CAAC,IAAI,CAAC,IAAI,GAAG,IAAI,CAAC,KAAK,CAAC,CAAA;QAEzC,mEAAmE;QACnE,qEAAqE;IACvE,CAAC;IAED,WAAW,CAAC,GAAG,EAAE,OAAO,EAAE,QAAQ,EAAE,QAAQ;QAC1C,0EAA0E;QAC1E,iDAAiD;QAEjD,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,IAAI,OAAO,GAAG,IAAI,CAAC,GAAG,CAAA;QACtB,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;gBAAE,OAAO,CAAC,GAAG,CAAC,uBAAuB,CAAC,CAAA;YACzD,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;YAAE,IAAI,CAAC,aAAa,CAAC,GAAG,CAAC,CAAA;IACrD,CAAC;CACF,CAAA"}
|
|
@@ -0,0 +1,11 @@
|
|
|
1
|
+
export = ByteModel;
|
|
2
|
+
declare class ByteModel {
|
|
3
|
+
constructor(max_sym?: number);
|
|
4
|
+
total_freq: number;
|
|
5
|
+
max_sym: number;
|
|
6
|
+
S: any[];
|
|
7
|
+
F: any[];
|
|
8
|
+
ModelDecode(src: any, rc: any): any;
|
|
9
|
+
ModelRenormalise(): void;
|
|
10
|
+
ModelEncode(dst: any, rc: any, sym: any): void;
|
|
11
|
+
}
|
|
@@ -0,0 +1,113 @@
|
|
|
1
|
+
"use strict";
|
|
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 adaptive probability model for encoding and decoding of symbols
|
|
35
|
+
// within a given alphabet, using the range coder to get/put the
|
|
36
|
+
// compressed data.
|
|
37
|
+
const MAX_FREQ = (1 << 16) - 17;
|
|
38
|
+
const STEP = 16;
|
|
39
|
+
module.exports = class ByteModel {
|
|
40
|
+
constructor(max_sym = 256) {
|
|
41
|
+
this.total_freq = max_sym;
|
|
42
|
+
this.max_sym = max_sym - 1;
|
|
43
|
+
this.S = new Array();
|
|
44
|
+
this.F = new Array();
|
|
45
|
+
for (var i = 0; i <= this.max_sym; i++) {
|
|
46
|
+
this.S[i] = i;
|
|
47
|
+
this.F[i] = 1;
|
|
48
|
+
}
|
|
49
|
+
}
|
|
50
|
+
ModelDecode(src, rc) {
|
|
51
|
+
// Find symbol
|
|
52
|
+
var freq = rc.RangeGetFrequency(this.total_freq);
|
|
53
|
+
// Linear scan to find cumulative frequency 'freq'
|
|
54
|
+
var acc = 0;
|
|
55
|
+
var x = 0;
|
|
56
|
+
while (acc + this.F[x] <= freq)
|
|
57
|
+
acc += this.F[x++];
|
|
58
|
+
// for (var acc = 0; (acc += this.F[x]) <= freq; x++)
|
|
59
|
+
// ;
|
|
60
|
+
// acc -= this.F[x];
|
|
61
|
+
// Update range coder
|
|
62
|
+
rc.RangeDecode(src, acc, this.F[x], this.total_freq);
|
|
63
|
+
// Update model
|
|
64
|
+
this.F[x] += STEP;
|
|
65
|
+
this.total_freq += STEP;
|
|
66
|
+
if (this.total_freq > MAX_FREQ)
|
|
67
|
+
this.ModelRenormalise();
|
|
68
|
+
// Keep symbols approximately frequency sorted
|
|
69
|
+
var sym = this.S[x];
|
|
70
|
+
if (x > 0 && this.F[x] > this.F[x - 1]) {
|
|
71
|
+
var tmp = this.F[x];
|
|
72
|
+
this.F[x] = this.F[x - 1];
|
|
73
|
+
this.F[x - 1] = tmp;
|
|
74
|
+
tmp = this.S[x];
|
|
75
|
+
this.S[x] = this.S[x - 1];
|
|
76
|
+
this.S[x - 1] = tmp;
|
|
77
|
+
}
|
|
78
|
+
return sym;
|
|
79
|
+
}
|
|
80
|
+
ModelRenormalise() {
|
|
81
|
+
// Halve all the frequencies, being careful not to hit zero
|
|
82
|
+
this.total_freq = 0;
|
|
83
|
+
for (var i = 0; i <= this.max_sym; i++) {
|
|
84
|
+
this.F[i] -= Math.floor(this.F[i] / 2);
|
|
85
|
+
this.total_freq += this.F[i];
|
|
86
|
+
}
|
|
87
|
+
}
|
|
88
|
+
ModelEncode(dst, rc, sym) {
|
|
89
|
+
// Find cumulative frequency
|
|
90
|
+
var acc = 0;
|
|
91
|
+
for (var x = 0; this.S[x] != sym; x++)
|
|
92
|
+
acc += this.F[x];
|
|
93
|
+
// Encode
|
|
94
|
+
rc.RangeEncode(dst, acc, this.F[x], this.total_freq);
|
|
95
|
+
// Update model
|
|
96
|
+
this.F[x] += STEP;
|
|
97
|
+
this.total_freq += STEP;
|
|
98
|
+
if (this.total_freq > MAX_FREQ)
|
|
99
|
+
// FIXME x2
|
|
100
|
+
this.ModelRenormalise();
|
|
101
|
+
// Keep symbols approximately frequency sorted
|
|
102
|
+
var sym = this.S[x];
|
|
103
|
+
if (x > 0 && this.F[x] > this.F[x - 1]) {
|
|
104
|
+
var tmp = this.F[x];
|
|
105
|
+
this.F[x] = this.F[x - 1];
|
|
106
|
+
this.F[x - 1] = tmp;
|
|
107
|
+
tmp = this.S[x];
|
|
108
|
+
this.S[x] = this.S[x - 1];
|
|
109
|
+
this.S[x - 1] = tmp;
|
|
110
|
+
}
|
|
111
|
+
}
|
|
112
|
+
};
|
|
113
|
+
//# sourceMappingURL=byte_model.js.map
|
|
@@ -0,0 +1 @@
|
|
|
1
|
+
{"version":3,"file":"byte_model.js","sourceRoot":"","sources":["../../src/htscodecs/byte_model.js"],"names":[],"mappings":";AAAA;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;GA+BG;AAEH,qEAAqE;AACrE,gEAAgE;AAChE,mBAAmB;AAEnB,MAAM,QAAQ,GAAG,CAAC,CAAC,IAAI,EAAE,CAAC,GAAG,EAAE,CAAA;AAC/B,MAAM,IAAI,GAAG,EAAE,CAAA;AAEf,MAAM,CAAC,OAAO,GAAG,MAAM,SAAS;IAC9B,YAAY,OAAO,GAAG,GAAG;QACvB,IAAI,CAAC,UAAU,GAAG,OAAO,CAAA;QACzB,IAAI,CAAC,OAAO,GAAG,OAAO,GAAG,CAAC,CAAA;QAC1B,IAAI,CAAC,CAAC,GAAG,IAAI,KAAK,EAAE,CAAA;QACpB,IAAI,CAAC,CAAC,GAAG,IAAI,KAAK,EAAE,CAAA;QAEpB,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,IAAI,IAAI,CAAC,OAAO,EAAE,CAAC,EAAE,EAAE,CAAC;YACvC,IAAI,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAA;YACb,IAAI,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAA;QACf,CAAC;IACH,CAAC;IAED,WAAW,CAAC,GAAG,EAAE,EAAE;QACjB,cAAc;QACd,IAAI,IAAI,GAAG,EAAE,CAAC,iBAAiB,CAAC,IAAI,CAAC,UAAU,CAAC,CAAA;QAEhD,kDAAkD;QAClD,IAAI,GAAG,GAAG,CAAC,CAAA;QACX,IAAI,CAAC,GAAG,CAAC,CAAA;QACT,OAAO,GAAG,GAAG,IAAI,CAAC,CAAC,CAAC,CAAC,CAAC,IAAI,IAAI;YAAE,GAAG,IAAI,IAAI,CAAC,CAAC,CAAC,CAAC,EAAE,CAAC,CAAA;QAElD,qDAAqD;QACrD,QAAQ;QACR,oBAAoB;QAEpB,qBAAqB;QACrB,EAAE,CAAC,WAAW,CAAC,GAAG,EAAE,GAAG,EAAE,IAAI,CAAC,CAAC,CAAC,CAAC,CAAC,EAAE,IAAI,CAAC,UAAU,CAAC,CAAA;QAEpD,eAAe;QACf,IAAI,CAAC,CAAC,CAAC,CAAC,CAAC,IAAI,IAAI,CAAA;QACjB,IAAI,CAAC,UAAU,IAAI,IAAI,CAAA;QACvB,IAAI,IAAI,CAAC,UAAU,GAAG,QAAQ;YAAE,IAAI,CAAC,gBAAgB,EAAE,CAAA;QAEvD,8CAA8C;QAC9C,IAAI,GAAG,GAAG,IAAI,CAAC,CAAC,CAAC,CAAC,CAAC,CAAA;QACnB,IAAI,CAAC,GAAG,CAAC,IAAI,IAAI,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,EAAE,CAAC;YACvC,IAAI,GAAG,GAAG,IAAI,CAAC,CAAC,CAAC,CAAC,CAAC,CAAA;YACnB,IAAI,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAA;YACzB,IAAI,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,GAAG,CAAA;YAEnB,GAAG,GAAG,IAAI,CAAC,CAAC,CAAC,CAAC,CAAC,CAAA;YACf,IAAI,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAA;YACzB,IAAI,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,GAAG,CAAA;QACrB,CAAC;QAED,OAAO,GAAG,CAAA;IACZ,CAAC;IAED,gBAAgB;QACd,2DAA2D;QAC3D,IAAI,CAAC,UAAU,GAAG,CAAC,CAAA;QACnB,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,IAAI,IAAI,CAAC,OAAO,EAAE,CAAC,EAAE,EAAE,CAAC;YACvC,IAAI,CAAC,CAAC,CAAC,CAAC,CAAC,IAAI,IAAI,CAAC,KAAK,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAA;YACtC,IAAI,CAAC,UAAU,IAAI,IAAI,CAAC,CAAC,CAAC,CAAC,CAAC,CAAA;QAC9B,CAAC;IACH,CAAC;IAED,WAAW,CAAC,GAAG,EAAE,EAAE,EAAE,GAAG;QACtB,4BAA4B;QAC5B,IAAI,GAAG,GAAG,CAAC,CAAA;QACX,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,IAAI,CAAC,CAAC,CAAC,CAAC,CAAC,IAAI,GAAG,EAAE,CAAC,EAAE;YAAE,GAAG,IAAI,IAAI,CAAC,CAAC,CAAC,CAAC,CAAC,CAAA;QAEvD,SAAS;QACT,EAAE,CAAC,WAAW,CAAC,GAAG,EAAE,GAAG,EAAE,IAAI,CAAC,CAAC,CAAC,CAAC,CAAC,EAAE,IAAI,CAAC,UAAU,CAAC,CAAA;QAEpD,eAAe;QACf,IAAI,CAAC,CAAC,CAAC,CAAC,CAAC,IAAI,IAAI,CAAA;QACjB,IAAI,CAAC,UAAU,IAAI,IAAI,CAAA;QACvB,IAAI,IAAI,CAAC,UAAU,GAAG,QAAQ;YAC5B,WAAW;YACX,IAAI,CAAC,gBAAgB,EAAE,CAAA;QAEzB,8CAA8C;QAC9C,IAAI,GAAG,GAAG,IAAI,CAAC,CAAC,CAAC,CAAC,CAAC,CAAA;QACnB,IAAI,CAAC,GAAG,CAAC,IAAI,IAAI,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,EAAE,CAAC;YACvC,IAAI,GAAG,GAAG,IAAI,CAAC,CAAC,CAAC,CAAC,CAAC,CAAA;YACnB,IAAI,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAA;YACzB,IAAI,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,GAAG,CAAA;YAEnB,GAAG,GAAG,IAAI,CAAC,CAAC,CAAC,CAAC,CAAC,CAAA;YACf,IAAI,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAA;YACzB,IAAI,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,GAAG,CAAA;QACrB,CAAC;IACH,CAAC;CACF,CAAA"}
|