@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.
Files changed (113) hide show
  1. package/dist/cram-bundle.js +1 -1
  2. package/dist/cramFile/codecs/byteArrayStop.d.ts +2 -2
  3. package/dist/cramFile/file.js +8 -10
  4. package/dist/cramFile/file.js.map +1 -1
  5. package/dist/cramFile/util.js.map +1 -1
  6. package/dist/htscodecs/arith_gen.d.ts +26 -0
  7. package/dist/htscodecs/arith_gen.js +558 -0
  8. package/dist/htscodecs/arith_gen.js.map +1 -0
  9. package/dist/htscodecs/arith_sh.d.ts +16 -0
  10. package/dist/htscodecs/arith_sh.js +128 -0
  11. package/dist/htscodecs/arith_sh.js.map +1 -0
  12. package/dist/htscodecs/byte_model.d.ts +11 -0
  13. package/dist/htscodecs/byte_model.js +113 -0
  14. package/dist/htscodecs/byte_model.js.map +1 -0
  15. package/dist/htscodecs/fqzcomp.d.ts +2 -0
  16. package/dist/htscodecs/fqzcomp.js +744 -0
  17. package/dist/htscodecs/fqzcomp.js.map +1 -0
  18. package/dist/htscodecs/index.d.ts +5 -0
  19. package/dist/htscodecs/index.js +70 -0
  20. package/dist/htscodecs/index.js.map +1 -0
  21. package/dist/htscodecs/iostream.d.ts +26 -0
  22. package/dist/htscodecs/iostream.js +242 -0
  23. package/dist/htscodecs/iostream.js.map +1 -0
  24. package/dist/htscodecs/main_arith_gen.d.ts +1 -0
  25. package/dist/htscodecs/main_arith_gen.js +86 -0
  26. package/dist/htscodecs/main_arith_gen.js.map +1 -0
  27. package/dist/htscodecs/main_fqzcomp.d.ts +1 -0
  28. package/dist/htscodecs/main_fqzcomp.js +112 -0
  29. package/dist/htscodecs/main_fqzcomp.js.map +1 -0
  30. package/dist/htscodecs/main_rans.d.ts +1 -0
  31. package/dist/htscodecs/main_rans.js +83 -0
  32. package/dist/htscodecs/main_rans.js.map +1 -0
  33. package/dist/htscodecs/main_rans4x16.d.ts +1 -0
  34. package/dist/htscodecs/main_rans4x16.js +82 -0
  35. package/dist/htscodecs/main_rans4x16.js.map +1 -0
  36. package/dist/htscodecs/main_tok3.d.ts +1 -0
  37. package/dist/htscodecs/main_tok3.js +84 -0
  38. package/dist/htscodecs/main_tok3.js.map +1 -0
  39. package/dist/htscodecs/rans.d.ts +2 -0
  40. package/dist/htscodecs/rans.js +480 -0
  41. package/dist/htscodecs/rans.js.map +1 -0
  42. package/dist/htscodecs/rans4x16.d.ts +2 -0
  43. package/dist/htscodecs/rans4x16.js +896 -0
  44. package/dist/htscodecs/rans4x16.js.map +1 -0
  45. package/dist/htscodecs/tok3.d.ts +2 -0
  46. package/dist/htscodecs/tok3.js +347 -0
  47. package/dist/htscodecs/tok3.js.map +1 -0
  48. package/esm/cramFile/codecs/byteArrayStop.d.ts +2 -2
  49. package/esm/cramFile/file.js +8 -10
  50. package/esm/cramFile/file.js.map +1 -1
  51. package/esm/cramFile/util.js.map +1 -1
  52. package/esm/htscodecs/arith_gen.d.ts +26 -0
  53. package/esm/htscodecs/arith_gen.js +558 -0
  54. package/esm/htscodecs/arith_gen.js.map +1 -0
  55. package/esm/htscodecs/arith_sh.d.ts +16 -0
  56. package/esm/htscodecs/arith_sh.js +128 -0
  57. package/esm/htscodecs/arith_sh.js.map +1 -0
  58. package/esm/htscodecs/byte_model.d.ts +11 -0
  59. package/esm/htscodecs/byte_model.js +113 -0
  60. package/esm/htscodecs/byte_model.js.map +1 -0
  61. package/esm/htscodecs/fqzcomp.d.ts +2 -0
  62. package/esm/htscodecs/fqzcomp.js +744 -0
  63. package/esm/htscodecs/fqzcomp.js.map +1 -0
  64. package/esm/htscodecs/index.d.ts +5 -0
  65. package/esm/htscodecs/index.js +70 -0
  66. package/esm/htscodecs/index.js.map +1 -0
  67. package/esm/htscodecs/iostream.d.ts +26 -0
  68. package/esm/htscodecs/iostream.js +242 -0
  69. package/esm/htscodecs/iostream.js.map +1 -0
  70. package/esm/htscodecs/main_arith_gen.d.ts +1 -0
  71. package/esm/htscodecs/main_arith_gen.js +86 -0
  72. package/esm/htscodecs/main_arith_gen.js.map +1 -0
  73. package/esm/htscodecs/main_fqzcomp.d.ts +1 -0
  74. package/esm/htscodecs/main_fqzcomp.js +112 -0
  75. package/esm/htscodecs/main_fqzcomp.js.map +1 -0
  76. package/esm/htscodecs/main_rans.d.ts +1 -0
  77. package/esm/htscodecs/main_rans.js +83 -0
  78. package/esm/htscodecs/main_rans.js.map +1 -0
  79. package/esm/htscodecs/main_rans4x16.d.ts +1 -0
  80. package/esm/htscodecs/main_rans4x16.js +82 -0
  81. package/esm/htscodecs/main_rans4x16.js.map +1 -0
  82. package/esm/htscodecs/main_tok3.d.ts +1 -0
  83. package/esm/htscodecs/main_tok3.js +84 -0
  84. package/esm/htscodecs/main_tok3.js.map +1 -0
  85. package/esm/htscodecs/rans.d.ts +2 -0
  86. package/esm/htscodecs/rans.js +480 -0
  87. package/esm/htscodecs/rans.js.map +1 -0
  88. package/esm/htscodecs/rans4x16.d.ts +2 -0
  89. package/esm/htscodecs/rans4x16.js +896 -0
  90. package/esm/htscodecs/rans4x16.js.map +1 -0
  91. package/esm/htscodecs/tok3.d.ts +2 -0
  92. package/esm/htscodecs/tok3.js +347 -0
  93. package/esm/htscodecs/tok3.js.map +1 -0
  94. package/package.json +1 -2
  95. package/src/cramFile/file.ts +8 -10
  96. package/src/cramFile/util.ts +0 -1
  97. package/src/htscodecs/Makefile +142 -0
  98. package/src/htscodecs/README.md +64 -0
  99. package/src/htscodecs/arith_gen.js +607 -0
  100. package/src/htscodecs/arith_sh.js +138 -0
  101. package/src/htscodecs/byte_model.js +126 -0
  102. package/src/htscodecs/fqzcomp.js +834 -0
  103. package/src/htscodecs/index.js +79 -0
  104. package/src/htscodecs/iostream.js +256 -0
  105. package/src/htscodecs/main_arith_gen.js +96 -0
  106. package/src/htscodecs/main_fqzcomp.js +113 -0
  107. package/src/htscodecs/main_rans.js +88 -0
  108. package/src/htscodecs/main_rans4x16.js +87 -0
  109. package/src/htscodecs/main_tok3.js +86 -0
  110. package/src/htscodecs/rans.js +545 -0
  111. package/src/htscodecs/rans4x16.js +1003 -0
  112. package/src/htscodecs/tok3.js +396 -0
  113. 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"}
@@ -0,0 +1,2 @@
1
+ export function decode(src: any, q_lens: any): any;
2
+ export function encode(src: any, q_lens: any, q_dirs: any): any;