rapier-jxl 2.0.0 → 2.1.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/AGENTS.md +5 -1
- package/ENCODER-COMPARISON.md +4 -4
- package/README.md +45 -14
- package/admit.mjs +5 -1
- package/ans.mjs +157 -0
- package/coefficient-ans.mjs +26 -0
- package/coefficient-effort.mjs +62 -0
- package/effort-level.mjs +10 -0
- package/effort.d.mts +3 -2
- package/effort.min.mjs +2 -2
- package/effort.mjs +31 -8
- package/frame.mjs +4 -1
- package/index.mjs +2 -2
- package/jfif.mjs +11 -0
- package/jpeg-ans.d.mts +3 -0
- package/jpeg-ans.min.mjs +2 -0
- package/jpeg-ans.mjs +9 -0
- package/jpeg-job.mjs +29 -0
- package/jpeg.d.mts +6 -2
- package/jpeg.min.mjs +2 -2
- package/jpeg.mjs +12 -36
- package/local.mjs +197 -0
- package/package.json +9 -1
- package/photo-ans.d.mts +3 -0
- package/photo-ans.min.mjs +2 -0
- package/photo-ans.mjs +9 -0
- package/photo-dct.mjs +83 -0
- package/photo-job.mjs +27 -0
- package/photo-quant.mjs +150 -0
- package/photo.d.mts +5 -1
- package/photo.min.mjs +2 -2
- package/photo.mjs +7 -101
- package/rapier-jxl.min.mjs +2 -2
- package/sizes.json +162 -53
- package/vardct.mjs +45 -8
- package/weighted.mjs +20 -10
package/effort.mjs
CHANGED
|
@@ -5,6 +5,8 @@
|
|
|
5
5
|
// 2. the specification's self-correcting (weighted) predictor, one context per channel;
|
|
6
6
|
// 3. the weighted predictor with its channel's tokens split by the predictor's own error, neighbouring intervals of
|
|
7
7
|
// libjxl's cut points merged wherever a shared prefix code costs less.
|
|
8
|
+
// 4. palette indices with local trees and histograms, the predictor and error intervals learned per group;
|
|
9
|
+
// 6. direct and palette planes with splits on the unclamped gradient and west-minus-northwest difference.
|
|
8
10
|
// A lossy request is effort 1's. Above effort 1 an exact picture's job spends its first half (fractions up to 0.5)
|
|
9
11
|
// writing effort 1's stream and its second searching, so a caller's clock can tell the search's own pace; `hurry`
|
|
10
12
|
// ends the search at its next step with the smallest stream written so far. Every choice is integer arithmetic: the
|
|
@@ -17,6 +19,7 @@ import {inspectPixels, losslessSteps} from './lossless.mjs';
|
|
|
17
19
|
import {lossySteps} from './lossy.mjs';
|
|
18
20
|
import {WEIGHTED_PREDICTOR, WEIGHTED_PROPERTY, WEIGHTED_CUTS, codeWeighted} from './weighted.mjs';
|
|
19
21
|
import {fault, admitOptions, admitPixels, job} from './admit.mjs';
|
|
22
|
+
import {localSteps} from './local.mjs';
|
|
20
23
|
|
|
21
24
|
export {LIMITS};
|
|
22
25
|
|
|
@@ -32,7 +35,7 @@ export function encodeSteps(data, width, height, options) {
|
|
|
32
35
|
}
|
|
33
36
|
|
|
34
37
|
// Effort 1 is the core's work (index.mjs), step for step; the rungs run after it and keep a smaller stream if they
|
|
35
|
-
// find one. A search that runs out of memory leaves
|
|
38
|
+
// find one. A search that runs out of memory leaves the smallest completed stream standing.
|
|
36
39
|
function* effortSteps(data, width, height, quality, colorSpace, effort) {
|
|
37
40
|
const shape = inspectPixels(data, width, height);
|
|
38
41
|
if (quality < 100) {
|
|
@@ -40,15 +43,32 @@ function* effortSteps(data, width, height, quality, colorSpace, effort) {
|
|
|
40
43
|
const exact = yield* part(losslessSteps(data, width, height, {shape, colorSpace}), 0, 2);
|
|
41
44
|
let bytes;
|
|
42
45
|
try { bytes = yield* part(lossySteps(data, width, height, {quality, shape, colorSpace}), 1, 2); }
|
|
43
|
-
catch (error) { if (!(error instanceof RangeError) || error.code) throw error; bytes = exact; }
|
|
46
|
+
catch (error) { if (error.code !== 'JXL_SIZE' && (!(error instanceof RangeError) || error.code)) throw error; bytes = exact; }
|
|
44
47
|
return exact.length <= bytes.length ? exact : bytes;
|
|
45
48
|
}
|
|
46
49
|
if (effort < 2) return yield* losslessSteps(data, width, height, {shape, colorSpace});
|
|
47
50
|
let best = yield* part(losslessSteps(data, width, height, {shape, colorSpace}), 0, 2);
|
|
51
|
+
if (effort >= 4) {
|
|
52
|
+
const searches = [searchSteps(data, width, height, shape.palette ? {...shape, palette: null} : shape, colorSpace, 3)];
|
|
53
|
+
if (shape.palette) searches.push(localSteps(data, width, height, shape, colorSpace, 4, true));
|
|
54
|
+
if (effort >= 6) {
|
|
55
|
+
searches.push(localSteps(data, width, height, shape, colorSpace, 6));
|
|
56
|
+
if (shape.palette) searches.push(localSteps(data, width, height, shape, colorSpace, 6, true));
|
|
57
|
+
}
|
|
58
|
+
for (let i = 0; i < searches.length; i++) {
|
|
59
|
+
try {
|
|
60
|
+
let step, hurried = false;
|
|
61
|
+
while (!(step = searches[i].next(hurried)).done) hurried = yield 0.5 + (i + step.value) / (2 * searches.length);
|
|
62
|
+
if (step.value && step.value.length < best.length) best = step.value;
|
|
63
|
+
if (hurried) return best;
|
|
64
|
+
} catch (error) { if (error.code !== 'JXL_SIZE' && (!(error instanceof RangeError) || error.code)) throw error; }
|
|
65
|
+
}
|
|
66
|
+
return best;
|
|
67
|
+
}
|
|
48
68
|
try {
|
|
49
69
|
const bytes = yield* part(searchSteps(data, width, height, shape.palette ? {...shape, palette: null} : shape, colorSpace, effort), 1, 2);
|
|
50
70
|
if (bytes && bytes.length < best.length) best = bytes;
|
|
51
|
-
} catch (error) { if (!(error instanceof RangeError) || error.code) throw error; }
|
|
71
|
+
} catch (error) { if (error.code !== 'JXL_SIZE' && (!(error instanceof RangeError) || error.code)) throw error; }
|
|
52
72
|
return best;
|
|
53
73
|
}
|
|
54
74
|
|
|
@@ -165,18 +185,21 @@ function* searchSteps(rgba, width, height, shape, colorSpace, effort) {
|
|
|
165
185
|
const [a, b] = candidates.map(price), margin = 27 * (layout.single ? 1 : groups + 1) + 8;
|
|
166
186
|
if (b + margin <= a) chosen = [candidates[1]];
|
|
167
187
|
else if (a + margin <= b) chosen = [candidates[0]];
|
|
188
|
+
else if (b < a) chosen = [candidates[1], candidates[0]];
|
|
168
189
|
}
|
|
169
|
-
let written = 0;
|
|
190
|
+
let written = 0, hurried = false;
|
|
170
191
|
const next = () => 0.5 + ++written / (2 * chosen.length * groups);
|
|
171
192
|
let smallest = null;
|
|
172
193
|
for (const plans of chosen) {
|
|
173
194
|
const bytes = yield* write(plans);
|
|
174
195
|
if (!bytes) return smallest;
|
|
175
|
-
|
|
196
|
+
// Price the likely winner first, but retain rung 2's byte choice on an equal-length completed pair.
|
|
197
|
+
if (!smallest || bytes.length < smallest.length || (bytes.length === smallest.length && plans === candidates[0])) smallest = bytes;
|
|
198
|
+
if (hurried) return smallest;
|
|
176
199
|
}
|
|
177
200
|
return smallest;
|
|
178
201
|
|
|
179
|
-
// A plan's stream, a group per step;
|
|
202
|
+
// A plan's stream, a group per step; a hurry at its final group keeps the already-completed candidate.
|
|
180
203
|
function* write(plans) {
|
|
181
204
|
const {tree, leaves, freqs} = assemble(plans);
|
|
182
205
|
const transforms = channels >= 3 ? [{type: 'rct', beginC: 0, rctType: 6}] : [];
|
|
@@ -197,7 +220,7 @@ function* searchSteps(rgba, width, height, shape, colorSpace, effort) {
|
|
|
197
220
|
const sections = [];
|
|
198
221
|
if (layout.single) {
|
|
199
222
|
group(0, (w, h) => { for (let c = 0; c < channels; c++) code(global, c, w, h); });
|
|
200
|
-
|
|
223
|
+
hurried = yield next();
|
|
201
224
|
sections.push(global.finish());
|
|
202
225
|
} else {
|
|
203
226
|
sections.push(global.finish());
|
|
@@ -209,7 +232,7 @@ function* searchSteps(rgba, width, height, shape, colorSpace, effort) {
|
|
|
209
232
|
for (let c = 0; c < channels; c++) code(section, c, w, h);
|
|
210
233
|
sections.push(section.finish());
|
|
211
234
|
});
|
|
212
|
-
if (yield next()) return null;
|
|
235
|
+
if ((hurried = yield next()) && g + 1 < groups) return null;
|
|
213
236
|
}
|
|
214
237
|
}
|
|
215
238
|
return assembleCodestream(header, sections);
|
package/frame.mjs
CHANGED
|
@@ -1,6 +1,7 @@
|
|
|
1
1
|
// Rapier's JPEG XL encoder: the codestream around a frame. MIT (LICENSE).
|
|
2
2
|
// A bare codestream (no container): signature, size header, image metadata, one frame with its table of contents.
|
|
3
3
|
import {BitWriter} from './bits.mjs';
|
|
4
|
+
import {admitOutputSize} from './admit.mjs';
|
|
4
5
|
|
|
5
6
|
export const GROUP_DIM = 256, DC_GROUP_DIM = 2048;
|
|
6
7
|
|
|
@@ -84,7 +85,9 @@ export function assembleCodestream(header, sections) {
|
|
|
84
85
|
const sizes = sections.map(section => section.length);
|
|
85
86
|
writeTOC(header, sizes);
|
|
86
87
|
const head = header.finish();
|
|
87
|
-
const
|
|
88
|
+
const length = head.length + sizes.reduce((a, b) => a + b, 0);
|
|
89
|
+
admitOutputSize(length);
|
|
90
|
+
const out = new Uint8Array(length);
|
|
88
91
|
out.set(head);
|
|
89
92
|
let at = head.length;
|
|
90
93
|
for (const section of sections) { out.set(section, at); at += section.length; }
|
package/index.mjs
CHANGED
|
@@ -32,10 +32,10 @@ function* pixelSteps(data, width, height, quality, colorSpace) {
|
|
|
32
32
|
if (quality >= 100) return yield* losslessSteps(data, width, height, {shape, colorSpace});
|
|
33
33
|
if (!shape.palette) return yield* lossySteps(data, width, height, {quality, shape, colorSpace});
|
|
34
34
|
// The exact stream of a few-colour picture is cheap to make (256-pixel sections), so it is made first: a lossy
|
|
35
|
-
// attempt that runs out of memory
|
|
35
|
+
// attempt that runs out of memory or stream space still leaves that completed answer.
|
|
36
36
|
const exact = yield* part(losslessSteps(data, width, height, {shape, colorSpace}), 0, 2);
|
|
37
37
|
let bytes;
|
|
38
38
|
try { bytes = yield* part(lossySteps(data, width, height, {quality, shape, colorSpace}), 1, 2); }
|
|
39
|
-
catch (error) { if (!(error instanceof RangeError) || error.code) throw error; bytes = exact; }
|
|
39
|
+
catch (error) { if (error.code !== 'JXL_SIZE' && (!(error instanceof RangeError) || error.code)) throw error; bytes = exact; }
|
|
40
40
|
return exact.length <= bytes.length ? exact : bytes;
|
|
41
41
|
}
|
package/jfif.mjs
CHANGED
|
@@ -42,6 +42,14 @@ const P3 = [0.5151, 0.2412, -0.001, 0.292, 0.6922, 0.0419, 0.1571, 0.0666, 0.784
|
|
|
42
42
|
const CURVE = [0, 0.00516, 0.01435, 0.02934, 0.05088, 0.07958, 0.11602, 0.16068, 0.21404, 0.27652, 0.34851, 0.43039,
|
|
43
43
|
0.52252, 0.62523, 0.73884, 0.86367, 1];
|
|
44
44
|
const PARAMETERS = [2.4, 1 / 1.055, 0.055 / 1.055, 1 / 12.92, 0.04045];
|
|
45
|
+
function srgbLinear(value) {
|
|
46
|
+
if (value <= 0.04045) return value / 12.92;
|
|
47
|
+
const base = (value + 0.055) / 1.055;
|
|
48
|
+
// base ** 2.4 as fixed arithmetic: the fifth root squared, times base squared. All engines round alike.
|
|
49
|
+
let root = 1;
|
|
50
|
+
for (let i = 0; i < 24; i++) root = (4 * root + base / (root * root * root * root)) / 5;
|
|
51
|
+
return base * base * root * root;
|
|
52
|
+
}
|
|
45
53
|
export function profileSpace(profile) {
|
|
46
54
|
const size = profile.length;
|
|
47
55
|
const u16 = at => (profile[at] << 8) | profile[at + 1];
|
|
@@ -64,6 +72,9 @@ export function profileSpace(profile) {
|
|
|
64
72
|
}
|
|
65
73
|
const n = sig(at) === 'curv' ? u32(at + 8) : 0;
|
|
66
74
|
if (n < 2 || at + 12 + 2 * n > size) return false;
|
|
75
|
+
// Matching a few knots cannot identify a transfer curve: a monotonic table can differ between them.
|
|
76
|
+
// Check every stored value as well as the interpolated knots below, which also reject coarse linear tables.
|
|
77
|
+
for (let i = 0; i < n; i++) if (Math.abs(u16(at + 12 + 2 * i) / 65535 - srgbLinear(i / (n - 1))) > 0.001) return false;
|
|
67
78
|
return CURVE.every((want, k) => {
|
|
68
79
|
const p = k * (n - 1) / 16, i = Math.floor(p), j = Math.min(n - 1, i + 1), f = p - i;
|
|
69
80
|
return Math.abs((u16(at + 12 + 2 * i) * (1 - f) + u16(at + 12 + 2 * j) * f) / 65535 - want) <= 0.001;
|
package/jpeg-ans.d.mts
ADDED
package/jpeg-ans.min.mjs
ADDED
|
@@ -0,0 +1,2 @@
|
|
|
1
|
+
/*! Rapier JXL 2.1.0 | MIT License Copyright (c) 2026 rapier.website Permission is hereby granted, free of charge, to any person obtaining a copy of this software and associated documentation files (the "Software"), to deal in the Software without restriction, including without limitation the rights to use, copy, modify, merge, publish, distribute, sublicense, and/or sell copies of the Software, and to permit persons to whom the Software is furnished to do so, subject to the following conditions: The above copyright notice and this permission notice shall be included in all copies or substantial portions of the Software. THE SOFTWARE IS PROVIDED "AS IS", WITHOUT WARRANTY OF ANY KIND, EXPRESS OR IMPLIED, INCLUDING BUT NOT LIMITED TO THE WARRANTIES OF MERCHANTABILITY, FITNESS FOR A PARTICULAR PURPOSE AND NONINFRINGEMENT. IN NO EVENT SHALL THE AUTHORS OR COPYRIGHT HOLDERS BE LIABLE FOR ANY CLAIM, DAMAGES OR OTHER LIABILITY, WHETHER IN AN ACTION OF CONTRACT, TORT OR OTHERWISE, ARISING FROM, OUT OF OR IN CONNECTION WITH THE SOFTWARE OR THE USE OR OTHER DEALINGS IN THE SOFTWARE. */
|
|
2
|
+
const t=Object.freeze({bytes:16777216,pixels:24e6,edge:16384}),o=Object.freeze({bytes:16777216,pixels:64e6,edge:16384}),n=[1];for(let t=1;t<=32;t++)n.push(2*n[t-1]);class r{constructor(t=4096){this.bytes=new Uint8Array(t),this.at=0,this.t=0,this.pending=0}write(t,o){if(!(o>=0&&o<n[t]))throw new Error("bit write out of range: "+t+" bits, "+o);let r=this.t+o*(1<<this.pending),e=this.pending+t;this.at+5>=this.bytes.length&&this.grow();const i=this.bytes;for(;e>=8;)i[this.at++]=255&r,r=Math.floor(r/256),e-=8;this.t=r,this.pending=e}o(t,o){for(let r=0;r<4;r++){const[e,i]=t[r];if(o>=i&&o-i<n[e])return this.write(2,r),void(e&&this.write(e,o-i))}throw new Error("U32 value out of range: "+o)}i(){this.pending&&(this.at+1>=this.bytes.length&&this.grow(),this.bytes[this.at++]=this.t,this.t=0,this.pending=0)}get u(){return 8*this.at+this.pending}grow(o=0){const n=t.bytes+64;if(this.at+o+16>n)throw Object.assign(new Error("The encoded picture exceeds 16 MiB."),{code:"JXL_SIZE"});const r=new Uint8Array(Math.min(n,Math.max(2*this.bytes.length,this.at+o+16)));r.set(this.bytes.subarray(0,this.at)),this.bytes=r}append(t){if(this.at+t.at+8>=this.bytes.length&&this.grow(t.at+8),this.pending)for(let o=0;o<t.at;o++)this.write(8,t.bytes[o]);else this.bytes.set(t.bytes.subarray(0,t.at),this.at),this.at+=t.at;t.pending&&this.write(t.pending,t.t)}finish(){const t=this.at+(this.pending?1:0),o=this.bytes.slice(0,t);return this.pending&&(o[this.at]=this.t),o}}function e(t){return t>=0?2*t:2*-t-1}function i(t){return 31-Math.clz32(t)}function s(t){return t<=1?0:32-Math.clz32(t-1)}function f(t){if(0===t)return 0;const o=t<0?32768:0;if(!((t=Math.abs(t))<65520))throw new Error("half float out of range: "+t);let n=0,r=t;for(;r>=2;)r/=2,n++;for(;r<1;)r*=2,n--;n<-14?(r=16384*t,n=-15):r-=1;let e=Math.round(1024*r);return 1024===e&&(e=0,n++),o|n+15<<10|e}const c=[0,1,8,16,9,2,3,10,17,24,32,25,18,11,4,5,12,19,26,33,40,48,41,34,27,20,13,6,7,14,21,28,35,42,49,56,57,50,43,36,29,22,15,23,30,37,44,51,58,59,52,45,38,31,39,46,53,60,61,54,47,55,62,63];function a(t,o="JXL_JPEG"){return Object.assign(new Error("JPEG: "+t),{code:o})}function l(t,o){const n=new Uint16Array(65536);let r=0,e=0;for(let i=1;i<=16;i++){for(let s=0;s<t[i-1];s++){const t=o[e++],s=r<<16-i,f=1<<16-i;if(s+f>=65536)throw a("bad Huffman table");n.fill(i<<8|t,s,s+f),r++}r<<=1}return n}const h=[.436,.2225,.0139,.3851,.7169,.0971,.143,.0606,.7139],u=[.5151,.2412,-.001,.292,.6922,.0419,.1571,.0666,.7841],w=[0,.00516,.01435,.02934,.05088,.07958,.11602,.16068,.21404,.27652,.34851,.43039,.52252,.62523,.73884,.86367,1],d=[2.4,1/1.055,.055/1.055,1/12.92,.04045];function y(t){if(t<=.04045)return t/12.92;const o=(t+.055)/1.055;let n=1;for(let t=0;t<24;t++)n=(4*n+o/(n*n*n*n))/5;return o*o*n*n}function p(t,o,n){const r=77===t[o]&&77===t[o+1];if(!r&&(73!==t[o]||73!==t[o+1]))return 1;const e=o=>o+2>n?0:r?t[o]<<8|t[o+1]:t[o]|t[o+1]<<8;if(42!==e(o+2))return 1;const i=o+(o=>o+4>n?0:r?(t[o]<<24|t[o+1]<<16|t[o+2]<<8|t[o+3])>>>0:(t[o]|t[o+1]<<8|t[o+2]<<16|t[o+3]<<24)>>>0)(o+4),s=e(i);for(let t=0,o=i+2;t<s&&o+12<=n;t++,o+=12){if(274!==e(o))continue;const t=e(o+8);return t>=1&&t<=8?t:1}return 1}function b(t,o,n,r,e,i,s,f,l){!function(t,o,n,r,e,i){if(!t.progressive){if(0!==n||63!==r||0!==e||0!==i)throw a("a bad spectral selection");for(const t of o){if(t.component.coverage.some(t=>-1!==t))throw a("a component coded twice");t.component.coverage.fill(0)}return}if(i>13||e&&e!==i+1)throw a("a bad successive approximation");if(0===n){if(0!==r)throw a("a bad spectral selection");for(const{component:t}of o){if(t.coverage[0]!==(e||-1))throw a(e?"a refinement scan out of order":"a DC scan twice");t.coverage[0]=i}return}if(1!==o.length)throw a("an interleaved AC scan");if(r>63||n>r)throw a("a bad spectral selection");const{coverage:s}=o[0].component;if(-1===s[0])throw a("an AC scan before the DC scan");for(let t=n;t<=r;t++){if(s[t]!==(e||-1))throw a(e?"a refinement scan out of order":"a band coded twice");s[t]=i}}(n,r,e,i,s,f);let h=o,u=0,w=0,d=0,y=!1;const p=()=>{for(;w<=24;){let o=0;if(y||h>=t.length)y=!0,d+=8;else if(o=t[h],255===o){const n=t[h+1];if(0===n)h+=2;else{if(255===n){h++;continue}y=!0,o=0,d+=8}}else h++;u=(u<<8|o)>>>0,w+=8}},b=()=>a("a scan cut short"),m=t=>{if(!t)return 0;if(p(),t>w-d)throw b();const o=u>>>w-t&(1<<t)-1;return w-=t,o},M=(t,o)=>t<1<<o-1?t-(1<<o)+1:t,g=t=>{if(!t)throw a("a scan without its Huffman table");p();const o=t[u>>>w-16&65535];if(!o)throw a("bad Huffman code");if(o>>8>w-d)throw b();return w-=o>>8,255&o},A=()=>{const o=w-d;if(o>=8)throw a("bytes left over in a scan");if(v)throw a("an end-of-band run past its interval");if(o&&(u>>>d&(1<<o)-1)!=(1<<o)-1)throw a("bad scan padding");for(;h+1<t.length&&255===t[h]&&255===t[h+1];)h++},U=o=>{if(A(),255!==t[h]||t[h+1]!==208+(7&o))throw a("a restart marker missing or out of order");h+=2,u=0,w=0,d=0,y=!1;for(const t of r)t.component.pred=0;v=0},I=(t,o)=>{if(o>11)throw a("a DC value out of range");const n=o?M(m(o),o):0;if(t.component.pred+=n,t.component.pred>2047||t.component.pred<-2048)throw a("a DC value out of range");return t.component.pred},x=(t,o)=>{if(t>10)throw a("an AC value out of range");const n=M(m(t),t)*(1<<o);if(n>1023||n<-1023)throw a("an AC value out of range");return n};let v=0;const E=1===r.length,S=n.progressive,C=()=>a("coefficients past the "+(S?"band":"block")),J=(t,o,n)=>{if(v>0)v--;else for(let r=e||1;r<=i;){const e=g(t.ac),s=e>>4,l=15&e;if(l){if(r+=s,r>i)throw C();o[n+c[r]]=x(l,f),r++}else{if(s<15){if(s&&!S)throw a("an end-of-band run in a sequential scan");v=(1<<s)-1+(s?m(s):0);break}if(r+=16,r>i+1)throw C()}}},L=S?0===e?0===s?(t,o,n)=>{const r=I(t,g(t.dc))*(1<<f);if(r>2047||r<-2048)throw a("a DC value out of range");o[n]=r}:(t,o,n)=>{m(1)&&(o[n]|=1<<f)}:0===s?J:(t,o,n)=>{const r=1<<f,s=-1<<f;let l=e;if(v<=0)for(;l<=i;l++){const e=g(t.ac);let f=e>>4,h=0;const u=15&e;if(u){if(1!==u)throw a("a bad refinement code");h=m(1)?r:s}else if(15!==f){v=1<<f,f&&(v+=m(f));break}for(;l<=i;l++){const t=n+c[l];if(0!==o[t])m(1)&&0===(o[t]&r)&&(o[t]+=o[t]>=0?r:s);else if(--f<0)break}if(f>=0||h&&l>i)throw a("coefficients past the band");h&&(o[n+c[l]]=h)}if(v>0){for(;l<=i;l++){const t=n+c[l];0!==o[t]&&m(1)&&0===(o[t]&r)&&(o[t]+=o[t]>=0?r:s)}v--}}:(t,o,n)=>{o[n]=I(t,g(t.dc)),J(t,o,n)};for(const t of r)t.component.pred=0;const X=r[0].component,k=E?X.blocksW*X.blocksH:n.mcusX*n.mcusY;for(let t=0;t<k;t++)if(l&&t>0&&t%l===0&&U(t/l-1),E){const o=t/X.blocksW|0,n=t-o*X.blocksW;L(r[0],X.coeffs,64*(o*X.stride+n))}else{const o=t/n.mcusX|0,e=t-o*n.mcusX;for(const t of r){const n=t.component;for(let r=0;r<n.v;r++)for(let i=0;i<n.h;i++)L(t,n.coeffs,64*((o*n.v+r)*n.stride+e*n.h+i))}}if(A(),255===t[h]&&t[h+1]>=208&&t[h+1]<=215){if(!l||k%l||t[h+1]!==208+(k/l-1&7))throw a("a restart marker out of order");for(h+=2;h+1<t.length&&255===t[h]&&255===t[h+1];)h++}if(h<t.length&&!(255===t[h]&&h+1<t.length&&0!==t[h+1]&&255!==t[h+1]))throw a("bytes left over after a scan");return h}const m=(t,o)=>Object.assign(new Error(o),{code:t});function M(o){if(o>t.bytes)throw m("JXL_SIZE","The encoded picture exceeds 16 MiB.")}function g(t){try{return t()}catch(t){if(t instanceof RangeError&&!t.code)throw m("JXL_MEMORY","Not enough memory for this picture.");throw t}}function A(n,r,e){const i=function(t={}){if(null===t||"object"!=typeof t)throw m("JXL_INPUT","Options are an object.");const o=void 0===t.effort?1:t.effort;if(!Number.isInteger(o)||o<1||o>9)throw m("JXL_INPUT","Effort is a whole number from 1 to 9.");return o}(r);if(!(n instanceof Uint8Array&&n.length))throw m("JXL_INPUT","A JPEG is a non-empty Uint8Array.");if(n.length>o.bytes)throw m("JXL_SIZE","A JPEG is at most 16 MiB.");const s=function(t){const o=function*(){let n,r;for(;!(n=g(()=>t.next(o.hurry))).done;)yield r=n.value;return r<1&&(yield 1),o.bytes=function(t){return M(t.length),t}(n.value)}();return o.bytes=null,o.hurry=!1,o}(function*(n,r,e,i){const f=function(t){if(!(t.length>4)||255!==t[0]||216!==t[1])throw a("not a JPEG");const n=[],r=[[],[]],e=new Set,i={count:0,p:[]};let s=null,f=0,m=-1,M=!1,g=1,A=2;const U=o=>t[o]<<8|t[o+1],I=(o,n)=>{for(let r=0;r<n.length;r++)if(t[o+r]!==n.charCodeAt(r))return!1;return!0};for(;;){if(A>=t.length){if(s&&s.scans)break;throw a("truncated")}if(255!==t[A])throw a("marker expected");for(;255===t[A];)A++;const h=t[A++];if(217===h)break;if(216===h)throw a("a second start");if(h>=208&&h<=215)throw a("a restart marker outside a scan");if(1===h)continue;if(A+2>t.length)throw a("truncated");const u=U(A),w=A+2,d=A+u;if(u<2||d>t.length)throw a("bad segment");if(A=d,219===h)for(let o=w;o<d;){const r=t[o]>>4,e=15&t[o];if(r>1||e>3||o+1+(r?128:64)>d)throw a("bad quantisation table");o++;const i=new Int32Array(64);for(let n=0;n<64;n++)if(i[c[n]]=r?U(o):t[o],o+=r?2:1,!i[c[n]])throw a("a zero quantisation step");if(s?.progressive&&s.components.some(t=>t.tq===e&&t.quant&&t.coverage.some(t=>0!==t)&&t.quant.some((t,o)=>t!==i[o])))throw a("a quantisation table changed between progressive scans");n[e]=i}else if(196===h)for(let o=w;o<d;){const n=t[o]>>4,i=15&t[o];if(n>1||i>3||o+17>d)throw a("bad Huffman table");o++;const s=t.subarray(o,o+16);let f=0;for(const t of s)f+=t;if(f>256||o+16+f>d)throw a("bad Huffman table");const c=t.subarray(o+16,o+16+f);r[n][i]=l(s,c),0===n&&(c.some(t=>t>15)?e.add(i):e.delete(i)),o+=16+f}else if(192===h||193===h||194===h){if(s)throw a("two frames");if(d-w<6)throw a("bad frame header");if(8!==t[w])throw a("only 8-bit samples");const n=U(w+1),r=U(w+3),e=t[w+5];if(!r||!n)throw a("a height left to a DNL marker");if(1!==e&&3!==e)throw a("only grey or three components");if(d-w<6+3*e)throw a("bad frame header");if(r>o.edge||n>o.edge)throw a("a picture past "+o.edge+" pixels on a side","JXL_DIMENSIONS");if(r*n>o.pixels)throw a("a picture past "+o.pixels+" pixels","JXL_DIMENSIONS");const i=[];let f=0;for(let o=0,n=w+6;o<e;o++,n+=3){const o=t[n+1]>>4,r=15&t[n+1],e=t[n+2];if(o<1||o>4||r<1||r>4)throw a("bad sampling");if(e>3)throw a("bad quantisation table");if(i.some(o=>o.id===t[n]))throw a("two components with one id");f+=o*r,i.push({id:t[n],h:o,v:r,tq:e,pred:0,coverage:new Int8Array(64).fill(-1)})}if(f>10)throw a("bad sampling");const c=Math.max(...i.map(t=>t.h)),l=Math.max(...i.map(t=>t.v)),u=Math.ceil(r/(8*c)),y=Math.ceil(n/(8*l));for(const t of i)t.blocksW=Math.ceil(Math.ceil(r*t.h/c)/8),t.blocksH=Math.ceil(Math.ceil(n*t.v/l)/8),t.stride=u*t.h,t.rows=y*t.v,t.coeffs=new Int16Array(t.stride*t.rows*64);s={progressive:194===h,width:r,height:n,components:i,hmax:c,vmax:l,mcusX:u,mcusY:y,scans:0}}else{if(195===h||h>=197&&h<=199||h>=201&&h<=203||h>=205&&h<=207)throw a("lossless, hierarchical or arithmetic coding");if(221===h){if(u<4)throw a("bad segment");f=U(w)}else{if(220===h)throw a("a DNL marker");if(218===h){if(!s)throw a("a scan before the frame");const o=t[w],i=[];if(o<1||o>4||w+1+2*o+3>d)throw a("bad scan header");let c=w+1;for(let e=0;e<o;e++,c+=2){const o=s.components.find(o=>o.id===t[c]);if(!o)throw a("a scan of an unknown component");if(i.some(t=>t.component===o))throw a("a component twice in one scan");if(o.quant||(o.quant=n[o.tq]),!o.quant)throw a("a scan without its quantisation table");i.push({component:o,dc:r[0][t[c+1]>>4],ac:r[1][15&t[c+1]],dcSlot:t[c+1]>>4})}const l=t[c],h=t[c+1],u=t[c+2]>>4,y=15&t[c+2];if((!s.progressive||0===l&&0===u)&&i.some(t=>e.has(t.dcSlot)))throw a("bad Huffman table");A=b(t,d,s,i,l,h,u,y,f),s.scans++}else if(224===h)I(w,"JFIF")&&(M=!0);else if(225===h)1===g&&I(w,"Exif")&&(g=p(t,w+6,d));else if(226===h){if(u>=16&&I(w,"ICC_PROFILE\0")){const o=t[w+12],n=t[w+13];if(!o||!n||o>n||i.count&&i.count!==n||i.p[o])throw a("a colour profile cut short");i.count=n,i.p[o]=t.subarray(w+14,d)}}else 238===h&&u>=14&&I(w,"Adobe")&&(m=t[w+11])}}}if(!s||!s.scans)throw a("no picture");for(const t of s.components)for(let o=0;o<64;o++)if(0!==t.coverage[o])throw a(s.progressive?"a progressive JPEG whose scans do not finish every coefficient":"a component without its scan");let x="srgb";if(i.count){const t=[];for(let o=1;o<=i.count;o++){if(!i.p[o])throw a("a colour profile cut short");t.push(i.p[o])}const o=new Uint8Array(t.reduce((t,o)=>t+o.length,0));for(let n=0,r=0;n<t.length;r+=t[n].length,n++)o.set(t[n],r);if(x=function(t){const o=t.length,n=o=>t[o]<<8|t[o+1],r=n=>n+4<=o?(t[n]<<24|t[n+1]<<16|t[n+2]<<8|t[n+3])>>>0:0,e=n=>n+4<=o?String.fromCharCode(t[n],t[n+1],t[n+2],t[n+3]):"",i=t=>(0|r(t))/65536,s=r(128);if(o<132||132+12*s>o||"acsp"!==e(36)||"XYZ "!==e(20))return null;const f=t=>{for(let n=0;n<s;n++){const i=132+12*n,s=r(i+4);if(e(i)===t)return s+r(i+8)<=o?s:-1}return-1},c=t=>{const s=f(t);if(s<0)return!1;if("para"===e(s)){const t=n(s+8),o=[1,3,4,5,7][t];if(3!==t&&4!==t)return!1;for(let t=0;t<o;t++)if(Math.abs(i(s+12+4*t)-(t<5?d[t]:0))>(t?.002:.01))return!1;return!0}const c="curv"===e(s)?r(s+8):0;if(c<2||s+12+2*c>o)return!1;for(let t=0;t<c;t++)if(Math.abs(n(s+12+2*t)/65535-y(t/(c-1)))>.001)return!1;return w.every((t,o)=>{const r=o*(c-1)/16,e=Math.floor(r),i=Math.min(c-1,e+1),f=r-e;return Math.abs((n(s+12+2*e)*(1-f)+n(s+12+2*i)*f)/65535-t)<=.001})};if(["A2B0","A2B1","A2B2"].some(t=>f(t)>=0))return null;if("GRAY"===e(16))return c("kTRC")?"srgb":null;if("RGB "!==e(16)||!["rTRC","gTRC","bTRC"].every(c))return null;const a=["rXYZ","gXYZ","bXYZ"].flatMap(t=>{const o=f(t);return o>=0&&"XYZ "===e(o)?[i(o+8),i(o+12),i(o+16)]:[9,9,9]}),l=t=>t.every((t,o)=>Math.abs(a[o]-t)<=.004);return l(h)?"srgb":l(u)?"display-p3":null}(o),!x)throw a("a colour profile other than sRGB or Display P3 (decode the picture and encode its pixels)")}let v=!1;return 3!==s.components.length||M||(v=m>=0?0===m:82===s.components[0].id&&71===s.components[1].id&&66===s.components[2].id),{width:s.width,height:s.height,components:s.components,progressive:s.progressive,ycbcr:!v,orientation:g,colorSpace:x}}(n);!function(o,n,r=t){if(!Number.isInteger(o)||!Number.isInteger(n)||o<1||n<1)throw m("JXL_INPUT","A picture is at least one pixel wide and high, in whole pixels.");if(o>r.edge||n>r.edge)throw m("JXL_DIMENSIONS","A picture is at most "+r.edge+" pixels on a side.");if(o*n>r.pixels)throw m("JXL_DIMENSIONS","A picture is at most "+r.pixels.toLocaleString("en-US")+" pixels.")}(f.width,f.height,o);const M=f.orientation>=5;return(t=>{Object.assign(s,t)})({width:M?f.height:f.width,height:M?f.width:f.height,orientation:f.orientation}),yield*i(f,r)}(n,i,0,e));return s}const U=256;function I(t,o){t.o([[9,1],[13,1],[18,1],[30,1]],o)}const x=[1,2,3,4,0,5,17,6,16,7,8,9,10,11,12,13,14,15],v=[2,4,3,2,2,4],E=[0,7,3,2,1,15];function S(t,o){const n=new Uint8Array(t.length),r=[];for(let o=0;o<t.length;o++)t[o]>0&&r.push(o);if(!r.length)return n;if(1===r.length)return n[r[0]]=1,n;for(let e=1;;e*=2){const i=r.map(o=>({w:Math.max(t[o],e),m:o,l:null,r:null})).sort((t,o)=>t.w-o.w||t.m-o.m),s=[];let f=0,c=0;const a=()=>f<i.length&&(c>=s.length||i[f].w<=s[c].w)?i[f++]:s[c++];for(;i.length-f+(s.length-c)>1;){const t=a(),o=a();s.push({w:t.w+o.w,m:-1,l:t,r:o})}const l=[[a(),0]];let h=0;for(;l.length;){const[t,o]=l.pop();t.m>=0?(n[t.m]=o,o>h&&(h=o)):l.push([t.l,o+1],[t.r,o+1])}if(h<=o)return n}}function C(t,o,n,r,e,i){const s=[];for(;s.push(r&e),r>>=i;)r--;for(let r=s.length-1;r>=0;r--)t.push(n),o.push(s[r])}function J(t,o){let n=0;for(let r=0;r<o;r++)n=n<<1|1&t,t>>=1;return n}function L(t){const o=new Uint32Array(17),n=new Uint32Array(17),r=new Uint16Array(t.length);for(const n of t)o[n]++;o[0]=0;for(let t=1,r=0;t<=16;t++)r=r+o[t-1]<<1,n[t]=r;for(let o=0;o<t.length;o++)t[o]&&(r[o]=J(n[t[o]]++,t[o]));return r}function X(t){const o=[];for(let n=0;n<t.length;n++)t[n]>0&&o.push(n);const n=o.length,r=o[n-1]+1||1,e=n<=4?new Uint8Array(r):S(t,15).subarray(0,r);let i=0;if(n<=4){o.sort((o,n)=>t[n]-t[o]||o-n),i=4===n&&t[o[0]]>t[o[2]]+t[o[3]]?1:0;const r=n<2?[0]:2===n?[1,1]:3===n?[1,2,2]:i?[1,2,3,3]:[2,2,2,2];o.forEach((t,o)=>{e[t]=r[o]}),o.sort((t,o)=>e[t]-e[o]||t-o)}return{M:r,A:e,U:L(e),I:n&&n<=4?o:null,S:i}}function k(t,o){if(o.M<=1)return;if(o.I){const n=i(o.M-1)+1;t.write(2,1),t.write(2,o.I.length-1);for(const r of o.I)t.write(n,r);return void(4===o.I.length&&t.write(1,o.S))}const n=o.A,r=[],e=[];let s=n.length-1;for(;s>0&&!n[s];)s--;let f=8;for(let t=0;t<=s;){const o=n[t];let i=1;for(;t+i<=s&&n[t+i]===o;)i++;let c=i;if(o&&o!==f&&(r.push(o),e.push(0),c--,f=o),c<3)for(let t=0;t<c;t++)r.push(o),e.push(0);else C(r,e,o?16:17,c-3,o?3:7,o?2:3);t+=i}const c=new Uint32Array(18);for(const t of r)c[t]++;let a=0;for(const t of c)t&&a++;1===a&&(c[r[0]?0:1]=c[r[0]]);const l=S(c,5),h=L(l);let u=18;for(;u>0&&!l[x[u-1]];)u--;t.write(2,0);for(let o=0;o<u;o++){const n=l[x[o]];t.write(v[n],E[n])}for(let o=0;o<r.length;o++){const n=r[o];t.write(l[n],h[n]),16===n?t.write(2,e[o]):17===n&&t.write(3,e[o])}}function N(t,o=0,n=0){return{split:t,C:o,J:n,L:1<<t}}function O(t,o,n=15){t.write(s(n+1),o.split),o.split!==n&&(t.write(s(o.split+1),o.C),t.write(s(o.split-o.C+1),o.J))}function _(t,o,n){if(o<t.L)return n[0]=o,n[1]=0,void(n[2]=0);const r=i(o),e=o-(1<<r),s=r-t.C-t.J;n[0]=t.L+(r-t.split<<t.C+t.J|e>>r-t.C<<t.J|e&(1<<t.J)-1),n[1]=s,n[2]=o>>t.J&(1<<s)-1}function T(t,o,n,r=0){_(t,o,D),n[r+D[0]]++}const D=[0,0,0];function P(t,o,n,r,e=0){_(n,r,D);const i=e+D[0];t.write(o.A[i],o.U[i]),D[1]&&t.write(D[1],D[2])}function G(t,o){let n=0;for(const t of o)n=Math.max(n,t);const r=n?s(n+1):0;if(r<=3){if(t.write(1,1),t.write(2,r),r)for(const n of o)t.write(r,n);return}t.write(1,0),t.write(1,0);const e=N(8),i=N(4),f=o.length>=64,c=new Uint32Array(f?257:n+1),a=[];for(let t=0;t<o.length;){let n=1;for(;f&&t+n<o.length&&o[t+n]===o[t];)n++;if(n>=8)a.push([o[t],n-1]),c[o[t]]++,T(i,n-1-7,c,224),t+=n;else{for(let r=0;r<n;r++)a.push([o[t],0]),c[o[t]]++;t+=n}}const l=X(c);if(!f){R(t,{X:new Uint8Array(1),k:[{config:e,code:l}]});for(const[o]of a)P(t,l,e,o);return}const h=new Uint32Array(1);h[0]=1,R(t,{N:{O:224,minLength:7,_:i},X:new Uint8Array([1,0]),k:[{config:N(0),code:X(h)},{config:e,code:l}]});for(const[o,n]of a)P(t,l,e,o),n&&P(t,l,i,n-7,224)}function R(t,{N:o=null,X:n,k:r}){t.write(1,o?1:0),o&&(t.o([[0,224],[0,512],[0,4096],[15,8]],o.O),t.o([[0,3],[0,4],[2,5],[8,9]],o.minLength),O(t,o._,8)),n.length>1&&G(t,n),t.write(1,1);for(const o of r)O(t,o.config,15);for(const o of r){const n=o.code.M-1;if(n){const o=i(n);t.write(1,1),t.write(4,o),t.write(o,n-(1<<o))}else t.write(1,0)}for(const o of r)k(t,o.code)}const Z=N(4),j=Object.freeze({O:224,minLength:7,_:Z}),q=N(0);function B(t,o){const n=[[3,0],[6,8],[10,72],[13,1096]];if("rct"===o.type)t.write(2,0),t.o(n,o.T),t.o([[0,6],[2,0],[4,2],[6,10]],o.D??6);else if("palette"===o.type)t.write(2,1),t.o(n,o.T),t.o([[0,1],[0,3],[0,4],[13,1]],o.P),t.o([[8,0],[10,256],[12,1280],[16,5376]],o.G),t.o([[0,0],[8,1],[10,257],[16,1281]],o.R??0),t.write(4,o.Z??0);else{if("squeeze"!==o.type)throw new Error("unknown transform "+o.type);{t.write(2,2);const r=o.j??[];t.o([[0,0],[4,1],[6,9],[8,41]],r.length);for(const o of r)t.write(1,o.q?1:0),t.write(1,o.B?1:0),t.o(n,o.T),t.o([[0,1],[0,2],[0,3],[4,4]],o.P)}}}function z(t,o,n,r,i,s){const f=s.Z,c=s.offset,a=s.Y,l=q,h=Z;let u=0;const w=(n,r=l,e=0)=>{t?P(t,o,r,n,e):T(r,n,o,e)},d=()=>{if(u){if(u>=8)w(0),w(u-7-1,h,224);else for(let t=0;t<u;t++)w(0);u=0}},y=(t,o)=>{const r=n[t];let i;1===a?i=r-o-c:(i=Math.round((r-o-c)/a),n[t]=o+c+i*a),0!==i?(d(),w(e(i))):u++};for(let t=0,o=0;t<i;t++)for(let e=0;e<r;e++,o++){let i;if(0===f)i=0;else{const s=e?n[o-1]:t?n[o-r]:0,c=t?n[o-r]:s;if(1===f)i=s;else if(2===f)i=c;else if(3===f)i=(s+c)/2|0;else{const a=s+c-(e&&t?n[o-r-1]:s);if(5===f){const t=s<c?s:c,o=s<c?c:s;i=a<t?t:a>o?o:a}else{if(4!==f)throw new Error("predictor "+f+" is not coded here");i=Math.abs(a-s)<Math.abs(a-c)?s:c}}}y(o,i)}d()}const Y=[N(0),N(1),N(2),N(3),N(4),N(4,1,1),N(5,1,1),N(2,0,1),N(3,0,1)],H=N(4,1,1);class F{constructor(t){this.H=t,this.small=new Uint32Array(64*t),this.F=Array.from({length:t},()=>null),this.K=new Float64Array(t)}add(t,o){if(this.K[t]++,o<64)return void this.small[64*t+o]++;const n=this.F[t]||(this.F[t]=new Map);n.set(o,(n.get(o)||0)+1)}V(t,o,n=64){const r=new Uint32Array(n);for(const n of t){const t=64*n;for(let n=0;n<64;n++)this.small[t+n]&&(r[Q(o,n)]+=this.small[t+n]);const e=this.F[n];if(e)for(const[t,n]of e)r[Q(o,t)]+=n}return r}W(t,o){let n=$(this.V(t,o));for(const r of t){const t=64*r;for(let r=o.L;r<64;r++)this.small[t+r]&&(n+=this.small[t+r]*(i(r)-o.C-o.J));const e=this.F[r];if(e)for(const[t,r]of e)n+=r*(i(t)-o.C-o.J)}return n}}const K=[0,0,0];function Q(t,o){return _(t,o,K),K[0]}const V=[1];for(let t=1;t<32;t++)V.push(V[t-1]/2);function W(t){let o=0;for(t>=2&&t<4294967296&&(o=31-Math.clz32(t),t*=V[o]);t>=1.4142135623730951;)t/=2,o++;for(;t<.7071067811865476;)t*=2,o--;const n=(t-1)/(t+1),r=n*n;let e=n,i=n;for(let t=3;t<30;t+=2)i+=(e*=r)/t;return o+2.885390081777927*i}function $(t){let o=0,n=0;for(const n of t)o+=n;for(const r of t)r&&(n+=r*W(o/r));return n}function tt(t,{maxClusters:o=48,newClusterCost:n=320}={}){const r=[];for(let o=0;o<t.H;o++)t.K[o]>0&&r.push(o);r.sort((o,n)=>t.K[n]-t.K[o]);const e=[];for(const i of r){const r=t.V([i],H),s=$(r);let f=-1,c=1/0;for(let t=0;t<e.length;t++){const o=e[t],n=new Uint32Array(64);for(let t=0;t<64;t++)n[t]=o.$[t]+r[t];const i=$(n)-o.tt-s;i<c&&(c=i,f=t)}if(f<0||e.length<o&&c>n)e.push({H:[i],$:r,tt:s});else{const t=e[f];t.H.push(i);for(let o=0;o<64;o++)t.$[o]+=r[o];t.tt=$(t.$)}}e.length||e.push({H:[],$:new Uint32Array(64),tt:0});const i=new Uint8Array(t.H);let s=2*t.H;const f=e.map((o,n)=>{for(const t of o.H)i[t]=n;let r=Y[0],e=1/0;for(const n of Y){const i=t.W(o.H,n);i<e&&(e=i,r=n)}const f=X(t.V(o.H,r,128));return s+=e+3*f.M+40,{config:r,code:f}});return{X:i,k:f,write:(t,o,n)=>{const r=f[i[o]];P(t,r.code,r.config,n)},tt:s}}const ot=[0,0,1,2,3,4,5,6,7,8,9,10,11,12,13,14,15,15,16,16,17,17,18,18,19,19,20,20,21,21,22,22,23,23,23,23,24,24,24,24,25,25,25,25,26,26,26,26,27,27,27,27,28,28,28,28,29,29,29,29,30,30,30,30],nt=[0,0,31,62,62,93,93,93,93,123,123,123,123,152,152,152,152,152,152,152,152,180,180,180,180,180,180,180,180,180,180,180,180,206,206,206,206,206,206,206,206,206,206,206,206,206,206,206,206,206,206,206,206,206,206,206,206,206,206,206,206,206,206,206],rt=c.map(t=>(7&t)<<3|t>>3),et=(t,o,n)=>2*(nt[t]+ot[o])+n;function it(t,o,n){const r=[];if(o>0){const n=new Uint32Array(2048);let e=0;for(const o of t)n[Math.min(2047,Math.max(0,o+1024))]++,e++;let i=0,s=e/(o+1);for(let t=0;t<2048&&r.length<o;t++)i+=n[t],i>s&&(r.push(t-1025),s=e*(r.length+1)/(o+1))}const e=r.length+1,i=new Uint8Array(39*e);for(let t=0;t<e;t++)for(let o=0;o<13;o++)i[(0+o)*e+t]=t,i[(13+o)*e+t]=n?e:e+(t>>1),i[(26+o)*e+t]=n?e:e+(e-1>>1)+1+(t>>1);let s=0;for(const t of i)s=Math.max(s,t+1);return{thresholds:r,ot:e,map:i,nt:s,rt:t=>{let o=0;for(const n of r)t>n&&o++;return o},et:(t,o)=>i[13*(t<2?1^t:2)*e+o],it:(t,o)=>(Math.min(t,64)<8?Math.min(t,64):4+(Math.min(t,64)>>1))*s+o,st:t=>37*s+458*t,H:495*s}}function st(t,o,n){const r=o.map(()=>function(t,o=0,n=1){return{Z:t,offset:o,Y:n,context:-1}}(n)),i=r.map(()=>new Uint32Array(257));o.forEach((t,o)=>z(null,i[o],t.data,t.w,t.h,r[o])),function(t,{ft:o=!0,ct:n=[]}={}){t.write(1,o?1:0),t.write(1,1),t.o([[0,0],[0,1],[4,2],[8,18]],n.length);for(const o of n)B(t,o)}(t,{ft:!1,ct:[]});const s=function(t,o){const n=[o],r=[],i=[];for(const t of n){if(t.left){r.push(t.lt+1,e(t.ht)),n.push(t.left,t.right);continue}t.context=i.length,i.push(t);const o=1===t.Y?0:31-Math.clz32(t.Y&-t.Y);r.push(0,t.Z,e(t.offset),o,(t.Y>>o)-1)}const s=new Uint32Array(64);for(const t of r)T(q,t,s);const f=X(s);R(t,{X:new Uint8Array(6),k:[{config:q,code:f}]});for(const o of r)P(t,f,q,o);return i}(t,function(t){return function(t,o){const n=(r,e)=>{if(r===e)return o?t[r].ut:t[r];const i=r+e>>1;return function(t,o,n,r){return{lt:t,ht:o,left:n,right:r}}(o,o?t[i].wt:i,n(i+1,e),n(r,i))};return n(0,t.length-1)}(t,0)}(r)),f=function(t,o,n,r=t=>t.context){const e=new Uint8Array(o.length+1);for(const t of o)e[t.context]=r(t)+1;const i=new Uint32Array(2);i[1]=1;const s=[{config:q,code:X(i)}],f=Math.max(...e);for(let t=0;t<f;t++)s.push({config:q,code:X(n[t])});return R(t,{N:j,X:e,k:s}),s}(t,s,i,t=>r.indexOf(t));o.forEach((o,n)=>z(t,f[n+1].code,o.data,o.w,o.h,r[n]))}function*ft(t,o={}){const{width:n,height:i,components:c,alpha:l}=t,h=1===c.length,u=t.quantScale||1;let w=t.quantFieldBase||1;for(;w/(2040*u)<1/16384;)w*=2;const d=t.ycbcr||h,y=h?[0,0,0]:d?[1,0,2]:[0,1,2],p=y.map(t=>c[t]),b=o.orders?y.map(t=>o.orders[t].map(t=>rt[t])):[rt,rt,rt],m=t=>1===t?0:2===t?1:-1,g=p.map(t=>m(t.h)),A=p.map(t=>m(t.v));if(g.includes(-1)||A.includes(-1))throw a("sampling factors other than one or two");const x=Math.max(...g),v=Math.max(...A),E=t=>x-g[t],S=t=>v-A[t];if(!d&&(x||v))throw a("an RGB JPEG with subsampling");const C=Math.ceil(n/(8<<x))<<x,J=Math.ceil(i/(8<<v))<<v,L=Math.ceil(n/U),k=Math.ceil(i/U),O=L*k,_=Math.ceil(C/256),D=_*Math.ceil(J/256),Z=1===O,j=(t,o,n)=>{const r=p[t];return h&&1!==t?-1:64*(n*r.stride+o)},q=t=>{let o=1024;for(;t%o;)o>>=1;return o},B=p.map(t=>d?1:t.quant[0]/q(t.quant[0])),z=p.map(t=>d?0:1024/q(t.quant[0])),Y=new Int32Array(C*J);for(let t=0;t<J;t++)for(let o=0;o<C;o++)Y[t*C+o]=p[1].coeffs[j(1,o,t)]*B[1]+z[1];const H=p[1].quant,K=(H[8]+H[16]+H[24]+H[32]+H[40])/u,Q=C*J,V=Q<256?0:Math.max(1,Math.min(7,s(Q)-s(K)-7));let W=it(Y,V,h);const $=new r(256);if(function(t,o,n,r,e,{xyb:i=!1,orientation:s=1,colorSpace:f}={}){t.write(16,2815),t.write(1,0),I(t,n),t.write(3,0),I(t,o),t.write(1,0);const c=1!==s;t.write(1,c?1:0),c&&t.write(6,s-1),t.write(4,8),e?t.write(3,5):t.write(2,0),t.write(1,i?1:0),3===r&&"display-p3"===f?t.write(21,903568):3===r?t.write(1,1):t.write(15,14100),c&&t.write(1,1),t.write(2,0),t.write(1,1),t.i()}($,n,i,h?1:3,!!l,{orientation:t.orientation,colorSpace:t.colorSpace}),$.write(1,0),$.write(2,0),$.write(1,0),$.write(2,2),$.write(8,111),$.write(1,d?1:0),d)for(let t=0;t<3;t++)$.write(2,0===g[t]&&0===A[t]?0:1===g[t]&&1===A[t]?1:1===g[t]?2:3);$.write(2,0),l&&$.write(2,0),function(t,o){t.write(o?7:5,0),t.write(11,1)}($,l);const ot=[0,1,2].map(()=>new Int32Array(1024)),nt=(t,o)=>{const n=t%L*32,r=32*(t/L|0),i=Math.min(32,C-n),s=Math.min(32,J-r);for(let t=0;t<s;t++)for(let s=0;s<i;s++)for(const i of[1,0,2]){const f=E(i),c=S(i),a=s>>f,l=t>>c;if(a<<f!==s||l<<c!==t)continue;const h=j(i,(n>>f)+a,(r>>c)+l),u=p[i].coeffs;let w=0;if(h>=0)for(let t=1;t<64;t++)u[h+t]&&w++;const d=ot[i],y=0===a?0===l?32:d[32*(l-1)]:0===l?d[32*l+a-1]:d[32*(l-1)+a]+d[32*l+a-1]+1>>1;d[32*l+a]=w;const m=W.et(i,W.rt(Y[(r+t)*C+n+s]));o(W.it(y,m),w);const M=W.st(m);let g=w>4?0:1,A=w;for(let t=1;t<64&&A;t++){const n=u[h+b[i][t]];o(M+et(A,t,g),e(n)),g=n?1:0,A-=g}}},ft=(V>0?3:2)*O+D;let ct=0;const at=function*(t){W=t;const n=new F(t.H);for(let t=0;t<O;t++)nt(t,(t,o)=>n.add(t,o)),yield++ct/ft;return(o.coding||tt)(n,o.clusters)};let lt=yield*at(W);if(V>0){const t=W,o=it(Y,0,h),n=yield*at(o);n.tt<=lt.tt?lt=n:W=t}const ht=Array.from({length:Z?1:2+D+O},()=>null),ut=t=>ht[Z?0:t]||(ht[Z?0:t]=new r(4096)),wt=new Uint32Array(ht.length);let dt=0;const yt=t=>{t=Z?0:t;const o=Math.ceil(ht[t].u/8);dt+=o-wt[t],wt[t]=o,M(dt)},pt=ut(0);pt.write(1,0);const bt=(t,o)=>{const r=t*U,e=o*U,s=Math.min(U,n-r),f=Math.min(U,i-e),c=new Int32Array(s*f);for(let t=0;t<f;t++)for(let o=0;o<s;o++)c[t*s+o]=l[4*((e+t)*n+r+o)+3];return{w:s,h:f,data:c}};for(let t=0;t<3;t++)pt.write(16,f(p[t].quant[0]/B[t]/(2040*u)*128));pt.write(2,3),pt.write(16,57343),pt.write(2,1),pt.write(5,0),pt.write(1,0),pt.write(4,0),pt.write(4,W.thresholds.length);for(const t of W.thresholds)pt.o([[4,0],[8,16],[16,272],[32,65808]],e(t));pt.write(4,0),pt.write(4,0),G(pt,W.map),pt.write(1,0),pt.write(2,0),pt.write(16,0),pt.write(16,0),pt.write(8,128),pt.write(8,128),pt.write(1,0),l&&st(pt,[Z?bt(0,0):{w:0,h:0,data:new Int32Array(0)}],5),yt(0);for(let o=0;o<D;o++){const n=ut(1+o),r=o%_*256,e=256*(o/_|0),i=Math.min(256,C-r),f=Math.min(256,J-e);n.write(2,0),st(n,[1,0,2].map(t=>{const o=i>>E(t),n=f>>S(t),s=new Int32Array(o*n),c=B[t],a=z[t];for(let i=0;i<n;i++)for(let n=0;n<o;n++){const f=j(t,(r>>E(t))+n,(e>>S(t))+i);s[i*o+n]=f<0?0:p[t].coeffs[f]*c+a}return{w:o,h:n,data:s}}),5);const c=i*f,a=s(c);a&&n.write(a,c-1);const l=i+7>>3,h=f+7>>3,u=new Int32Array(2*c).fill(w-1,c);if(t.quantFields)for(let o=0;o<f;o++)for(let n=0;n<i;n++)u[c+o*i+n]=t.quantFields[(e+o)*C+r+n]-1;st(n,[{w:l,h:h,data:new Int32Array(l*h)},{w:l,h:h,data:new Int32Array(l*h)},{w:c,h:2,data:u},{w:i,h:f,data:new Int32Array(c)}],w>1?5:0),yt(1+o),yield++ct/ft}const mt=ut(1+D);mt.write(1,0),mt.write(3,7),mt.write(16,f(w/(2040*u))),st(mt,[0,1,2].map(t=>{const o=new Int32Array(64),n=p[t].quant;for(let t=0;t<8;t++)for(let r=0;r<8;r++)o[8*r+t]=n[8*t+r];return{w:8,h:8,data:o}}),5);for(let t=1;t<17;t++)mt.write(3,0);if(O>1&&mt.write(s(O),0),o.orders){mt.o([[0,95],[0,19],[0,0],[13,0]],1);const t=[];for(const n of y){const r=Array.from({length:64},(t,o)=>o),e=o.orders[n].map(t=>{const o=r.indexOf(t);return r.splice(o,1),o});let i=64;for(;i>1&&!e[i-1];)i--;t.push(i-1,...e.slice(1,i))}const n=N(0),r=new Uint32Array(64);for(const o of t)T(n,o,r);const e=X(r);R(mt,{X:new Uint8Array(8),k:[{config:n,code:e}]});for(const o of t)P(mt,e,n,o)}else mt.write(2,2);lt.writeHistograms?lt.writeHistograms(mt):R(mt,{X:lt.X,k:lt.k}),yt(1+D);for(let t=0;t<O;t++){const o=ut(2+D+t);nt(t,(t,n)=>lt.write(o,t,n)),lt.flush&<.flush(o),l&&!Z&&st(o,[bt(t%L,Math.floor(t/L))],5),yt(2+D+t),yield++ct/ft}return function(t,o){const n=o.map(t=>t.length);!function(t,o){t.write(1,0),t.i();for(const n of o)t.o([[10,0],[14,1024],[22,17408],[30,4211712]],n);t.i()}(t,n);const r=t.finish(),e=r.length+n.reduce((t,o)=>t+o,0);M(e);const i=new Uint8Array(e);i.set(r);let s=r.length;for(const t of o)i.set(t,s),s+=t.length;return i}($,function(t){return t.map(t=>t?t.finish():new Uint8Array(0))}(ht))}function*ct(t,o,n){if(o<3&&!n)return yield*ft(t);let r=n;try{const n=ft(t);let e,i;for(;!(e=n.next()).done;)if(i=yield e.value/2,i&&r&&e.value<1)return n.return(),r;if((!r||e.value.length<r.length)&&(r=e.value),i||o<3)return r;const s=o<4?1:3,f=ft(t,{clusters:{maxClusters:32,newClusterCost:160}});for(;!(e=f.next()).done;)if(i=yield.5+e.value/(2*s),i&&e.value<1)return f.return(),r;if(e.value.length<r.length&&(r=e.value),i||o<4)return r;const a=function*(t){const o=c.map(t=>(7&t)<<3|t>>3),n=t.map(()=>new Uint32Array(64)),r=t.reduce((t,o)=>t+o.coeffs.length,0);let e=0;for(let i=0;i<t.length;i++){const{coeffs:s}=t[i],f=n[i];for(let t=0;t<s.length;t+=65536){const n=Math.min(t+65536,s.length);for(let r=t;r<n;r+=64)for(let t=1;t<64;t++)s[r+o[t]]&&f[t]++;if(e+=n-t,yield e/r)return null}}return n.map(t=>[0,...Array.from({length:63},(t,o)=>o+1).sort((o,n)=>Math.floor(t[n]/8)-Math.floor(t[o]/8)||o-n)])}(t.components);for(;!(e=a.next()).done;)if(yield.5+(1+e.value)/6)return a.return(),r;const l=ft(t,{orders:e.value});for(;!(e=l.next()).done;)if(i=yield.5+(2+e.value)/6,i&&e.value<1)return l.return(),r;return e.value.length<r.length?e.value:r}catch(t){if(!r||(!(t instanceof RangeError)||t.code)&&"JXL_SIZE"!==t.code)throw t;return r}}const at=[5,4,4,4,4,4,3,3,3,3,3,6,7],lt=[17,11,15,3,9,7,4,2,5,6,0,33,1],ht=196608;function ut(t,o){if(!o)return void t.write(1,0);const n=i(o);t.write(1,1),t.write(3,n),t.write(n,o-(1<<n))}function wt(t){let o=0,n=0;for(let r=0;r<t.length;r++)o+=t[r],t[r]&&(n=r+1);const r=new Uint16Array(n||1);if(!o)return r[0]=4096,r;let e=0;for(let i=0;i<n;i++)t[i]&&(e+=r[i]=Math.max(1,Math.round(4096*t[i]/o)));for(;4096!==e;){const i=e<4096?1:-1;let s=-1,f=-1/0;for(let e=0;e<n;e++)if(t[e]&&(i>0||r[e]>1)){const n=i*(4096*t[e]-r[e]*o);n>f&&(s=e,f=n)}r[s]+=i,e+=i}return r}function dt(t,o){const n=[];let r=0;for(let t=0;t<o.length;t++)o[t]&&n.push(t),o[t]>o[r]&&(r=t);if(n.length<=2){t.write(1,1),t.write(1,n.length-1);for(const o of n)ut(t,o);return void(2===n.length&&t.write(12,o[n[0]]))}t.write(1,0),t.write(1,0),t.write(3,7),t.write(3,5),ut(t,o.length-3);const e=new Uint8Array(o.length);let s=0;for(let t=0;t<o.length;t++)t!==r&&o[t]&&(e[t]=i(o[t])+1,s=Math.max(s,e[t]+(t<r?1:0)));e[r]=s;for(const o of e)t.write(at[o],lt[o]);for(let n=0;n<o.length;n++)if(n!==r&&e[n]>1){const r=e[n]-1;t.write(r,o[n]-(1<<r))}}function yt(t,o){const n=1<<o,r=4096>>o,e=new Uint16Array(t.length+1),i=new Uint16Array(4096);let s=-1;for(let o=0;o<t.length;o++)e[o+1]=e[o]+t[o],4096===t[o]&&(s=o);if(s>=0){for(let t=0;t<4096;t++)i[t]=t;return{starts:e,reverse:i}}const f=new Int32Array(n),c=new Uint16Array(n),a=new Int32Array(n),l=[],h=[];f.set(t);for(let t=0;t<n;t++)f[t]<r?l.push(t):f[t]>r&&h.push(t);for(;h.length;){const t=h.pop(),o=l.pop(),n=r-f[o];f[t]-=n,c[o]=t,a[o]=f[t]-f[o],f[t]<r?l.push(t):f[t]>r&&h.push(t)}for(let t=0;t<n;t++)for(let o=0;o<r;o++){const n=o<f[t],s=n?t:c[t],l=n?o:o+a[t];i[e[s]+l]=t*r+o}return{starts:e,reverse:i}}function pt(t,o){const n=tt(t,o),e=n.k.map(()=>[]);for(let o=0;o<t.H;o++)t.K[o]&&e[n.X[o]].push(o);const i=n.k.map(({config:o},n)=>({config:o,counts:wt(t.V(e[n],o,256))})),f=Math.max(5,s(Math.max(...i.map(t=>t.counts.length))));for(const t of i)Object.assign(t,yt(t.counts,f));const c=new r;c.write(1,0),n.X.length>1&&G(c,n.X),c.write(1,0),c.write(2,f-5);for(const t of i)O(c,t.config,f);for(const t of i)dt(c,t.counts);const a=new Uint32Array(ht),l=new Uint8Array(ht),h=new Uint16Array(ht),u=[0,0,0];let w=0;return{tt:n.tt,write:(t,o,r)=>{a[w]=r,l[w++]=n.X[o]},flush:t=>{let o=1245184;for(let t=w-1;t>=0;t--){const n=i[l[t]];_(n.config,a[t],u);const r=u[0],e=n.counts[r];a[t]=u[2],l[t]=u[1],o>>>20>=e&&(h[t]=65535&o,l[t]|=64,o>>>=16);const s=Math.floor(o/e);o=4096*s+n.reverse[n.starts[r]+o-s*e]}t.write(32,o);for(let o=0;o<w;o++){64&l[o]&&t.write(16,h[o]);const n=63&l[o];n&&t.write(n,a[o])}w=0},writeHistograms:t=>t.append(c)}}function*bt(t,o,n){if(o<2)return yield*ct(t,o,n);let r=yield*function*(t){let o,n;for(;!(o=t.next(n)).done;)n=yield(0+o.value)/2;return o.value}(ct(t,o,n));if(yield.5)return r;try{const o=ft(t,{coding:pt});let n;for(;!(n=o.next()).done;)if((yield.5+n.value/2)&&n.value<1)return o.return(),r;n.value.length<r.length&&(r=n.value)}catch(t){if((!(t instanceof RangeError)||t.code)&&"JXL_SIZE"!==t.code)throw t}return r}function mt(t,o){return function(t){return function(t){let o;for(;!(o=t.next()).done;);o.value}(t),{bytes:t.bytes,width:t.width,height:t.height,orientation:t.orientation}}(Mt(t,o))}function Mt(t,o){return A(t,o,bt)}export{o as LIMITS,mt as transcode,Mt as transcodeSteps};
|
package/jpeg-ans.mjs
ADDED
|
@@ -0,0 +1,9 @@
|
|
|
1
|
+
// SPDX-License-Identifier: MIT
|
|
2
|
+
// JPEG coefficients carried without reconstruction boxes; one checked job owner.
|
|
3
|
+
import {JPEG_LIMITS} from './bits.mjs';
|
|
4
|
+
import {jpegAnswer, jpegJob} from './jpeg-job.mjs';
|
|
5
|
+
import {coefficientAnsSteps} from './coefficient-ans.mjs';
|
|
6
|
+
|
|
7
|
+
export {JPEG_LIMITS as LIMITS};
|
|
8
|
+
export function transcode(jpeg, options) { return jpegAnswer(transcodeSteps(jpeg, options)); }
|
|
9
|
+
export function transcodeSteps(jpeg, options) { return jpegJob(jpeg, options, coefficientAnsSteps); }
|
package/jpeg-job.mjs
ADDED
|
@@ -0,0 +1,29 @@
|
|
|
1
|
+
// SPDX-License-Identifier: MIT
|
|
2
|
+
// One checked JPEG job shared by the prefix and optional ANS doors.
|
|
3
|
+
import {parseJPEG} from './jfif.mjs';
|
|
4
|
+
import {admitEffort} from './effort-level.mjs';
|
|
5
|
+
import {JPEG_LIMITS, complete} from './bits.mjs';
|
|
6
|
+
import {fault, admitSize, job} from './admit.mjs';
|
|
7
|
+
|
|
8
|
+
export function jpegAnswer(it) {
|
|
9
|
+
complete(it);
|
|
10
|
+
return {bytes: it.bytes, width: it.width, height: it.height, orientation: it.orientation};
|
|
11
|
+
}
|
|
12
|
+
|
|
13
|
+
export function jpegJob(jpeg, options, coefficients) {
|
|
14
|
+
const effort = admitEffort(options);
|
|
15
|
+
if (!(jpeg instanceof Uint8Array) || !jpeg.length) throw fault('JXL_INPUT', 'A JPEG is a non-empty Uint8Array.');
|
|
16
|
+
if (jpeg.length > JPEG_LIMITS.bytes) throw fault('JXL_SIZE', 'A JPEG is at most 16 MiB.');
|
|
17
|
+
const it = job(carry(jpeg, effort, shown => Object.assign(it, shown), coefficients));
|
|
18
|
+
return it;
|
|
19
|
+
}
|
|
20
|
+
|
|
21
|
+
// The first step reads the scans and sets the shown dimensions and orientation;
|
|
22
|
+
// the selected coefficient writer then owns the remaining steps. Admission is shared.
|
|
23
|
+
function* carry(jpeg, effort, shown, coefficients) {
|
|
24
|
+
const parsed = parseJPEG(jpeg);
|
|
25
|
+
admitSize(parsed.width, parsed.height, JPEG_LIMITS);
|
|
26
|
+
const swapped = parsed.orientation >= 5;
|
|
27
|
+
shown({width: swapped ? parsed.height : parsed.width, height: swapped ? parsed.width : parsed.height, orientation: parsed.orientation});
|
|
28
|
+
return yield* coefficients(parsed, effort);
|
|
29
|
+
}
|
package/jpeg.d.mts
CHANGED
|
@@ -9,7 +9,11 @@ export interface Transcoded {
|
|
|
9
9
|
/** The Exif orientation, 1 to 8, kept in the JPEG XL header. */ orientation: number;
|
|
10
10
|
}
|
|
11
11
|
/** A JPEG's bytes to JPEG XL, its coefficients carried; a JPEG this path does not take is refused as JXL_JPEG. */
|
|
12
|
-
export
|
|
12
|
+
export interface TranscodeOptions {
|
|
13
|
+
/** Whole number 1 to 9; default 1. From 3, try a 32-cluster budget; from 4, also a learned coefficient order. */
|
|
14
|
+
effort?: number;
|
|
15
|
+
}
|
|
16
|
+
export function transcode(jpeg: Uint8Array, options?: TranscodeOptions): Transcoded;
|
|
13
17
|
/** transcode's work as steps; the first step reads the scans and sets the picture's size as shown and orientation. */
|
|
14
18
|
export interface TranscodeJob extends Job { width?: number; height?: number; orientation?: number; }
|
|
15
|
-
export function transcodeSteps(jpeg: Uint8Array): TranscodeJob;
|
|
19
|
+
export function transcodeSteps(jpeg: Uint8Array, options?: TranscodeOptions): TranscodeJob;
|
package/jpeg.min.mjs
CHANGED
|
@@ -1,2 +1,2 @@
|
|
|
1
|
-
/*! Rapier JXL 2.
|
|
2
|
-
const t=Object.freeze({bytes:16777216,pixels:24e6,edge:16384}),o=Object.freeze({bytes:16777216,pixels:64e6,edge:16384}),n=[1];for(let t=1;t<=32;t++)n.push(2*n[t-1]);class r{constructor(t=4096){this.bytes=new Uint8Array(t),this.at=0,this.t=0,this.pending=0}write(t,o){if(!(o>=0&&o<n[t]))throw new Error("bit write out of range: "+t+" bits, "+o);let r=this.t+o*(1<<this.pending),e=this.pending+t;this.at+5>=this.bytes.length&&this.grow();const s=this.bytes;for(;e>=8;)s[this.at++]=255&r,r=Math.floor(r/256),e-=8;this.t=r,this.pending=e}o(t,o){for(let r=0;r<4;r++){const[e,s]=t[r];if(o>=s&&o-s<n[e])return this.write(2,r),void(e&&this.write(e,o-s))}throw new Error("U32 value out of range: "+o)}i(){this.pending&&(this.at+1>=this.bytes.length&&this.grow(),this.bytes[this.at++]=this.t,this.t=0,this.pending=0)}get u(){return 8*this.at+this.pending}grow(o=0){const n=t.bytes+64;if(this.at+o+16>n)throw Object.assign(new Error("The encoded picture exceeds 16 MiB."),{code:"JXL_SIZE"});const r=new Uint8Array(Math.min(n,Math.max(2*this.bytes.length,this.at+o+16)));r.set(this.bytes.subarray(0,this.at)),this.bytes=r}append(t){if(this.at+t.at+8>=this.bytes.length&&this.grow(t.at+8),this.pending)for(let o=0;o<t.at;o++)this.write(8,t.bytes[o]);else this.bytes.set(t.bytes.subarray(0,t.at),this.at),this.at+=t.at;t.pending&&this.write(t.pending,t.t)}finish(){const t=this.at+(this.pending?1:0),o=this.bytes.slice(0,t);return this.pending&&(o[this.at]=this.t),o}}function e(t){return t>=0?2*t:2*-t-1}function s(t){return 31-Math.clz32(t)}function i(t){return t<=1?0:32-Math.clz32(t-1)}function f(t){if(0===t)return 0;const o=t<0?32768:0;if(!((t=Math.abs(t))<65520))throw new Error("half float out of range: "+t);let n=0,r=t;for(;r>=2;)r/=2,n++;for(;r<1;)r*=2,n--;n<-14?(r=16384*t,n=-15):r-=1;let e=Math.round(1024*r);return 1024===e&&(e=0,n++),o|n+15<<10|e}const c=[0,1,8,16,9,2,3,10,17,24,32,25,18,11,4,5,12,19,26,33,40,48,41,34,27,20,13,6,7,14,21,28,35,42,49,56,57,50,43,36,29,22,15,23,30,37,44,51,58,59,52,45,38,31,39,46,53,60,61,54,47,55,62,63];function a(t,o="JXL_JPEG"){return Object.assign(new Error("JPEG: "+t),{code:o})}function h(t,o){const n=new Uint16Array(65536);let r=0,e=0;for(let s=1;s<=16;s++){for(let i=0;i<t[s-1];i++){const t=o[e++],i=r<<16-s,f=1<<16-s;if(i+f>=65536)throw a("bad Huffman table");n.fill(s<<8|t,i,i+f),r++}r<<=1}return n}const l=[.436,.2225,.0139,.3851,.7169,.0971,.143,.0606,.7139],u=[.5151,.2412,-.001,.292,.6922,.0419,.1571,.0666,.7841],w=[0,.00516,.01435,.02934,.05088,.07958,.11602,.16068,.21404,.27652,.34851,.43039,.52252,.62523,.73884,.86367,1],d=[2.4,1/1.055,.055/1.055,1/12.92,.04045];function p(t,o,n){const r=77===t[o]&&77===t[o+1];if(!r&&(73!==t[o]||73!==t[o+1]))return 1;const e=o=>o+2>n?0:r?t[o]<<8|t[o+1]:t[o]|t[o+1]<<8;if(42!==e(o+2))return 1;const s=o+(o=>o+4>n?0:r?(t[o]<<24|t[o+1]<<16|t[o+2]<<8|t[o+3])>>>0:(t[o]|t[o+1]<<8|t[o+2]<<16|t[o+3]<<24)>>>0)(o+4),i=e(s);for(let t=0,o=s+2;t<i&&o+12<=n;t++,o+=12){if(274!==e(o))continue;const t=e(o+8);return t>=1&&t<=8?t:1}return 1}function b(t,o,n,r,e,s,i,f,h){!function(t,o,n,r,e,s){if(!t.progressive){if(0!==n||63!==r||0!==e||0!==s)throw a("a bad spectral selection");for(const t of o){if(t.component.coverage.some(t=>-1!==t))throw a("a component coded twice");t.component.coverage.fill(0)}return}if(s>13||e&&e!==s+1)throw a("a bad successive approximation");if(0===n){if(0!==r)throw a("a bad spectral selection");for(const{component:t}of o){if(t.coverage[0]!==(e||-1))throw a(e?"a refinement scan out of order":"a DC scan twice");t.coverage[0]=s}return}if(1!==o.length)throw a("an interleaved AC scan");if(r>63||n>r)throw a("a bad spectral selection");const{coverage:i}=o[0].component;if(-1===i[0])throw a("an AC scan before the DC scan");for(let t=n;t<=r;t++){if(i[t]!==(e||-1))throw a(e?"a refinement scan out of order":"a band coded twice");i[t]=s}}(n,r,e,s,i,f);let l=o,u=0,w=0,d=0,p=!1;const b=()=>{for(;w<=24;){let o=0;if(p||l>=t.length)p=!0,d+=8;else if(o=t[l],255===o){const n=t[l+1];if(0===n)l+=2;else{if(255===n){l++;continue}p=!0,o=0,d+=8}}else l++;u=(u<<8|o)>>>0,w+=8}},y=()=>a("a scan cut short"),m=t=>{if(!t)return 0;if(b(),t>w-d)throw y();const o=u>>>w-t&(1<<t)-1;return w-=t,o},g=(t,o)=>t<1<<o-1?t-(1<<o)+1:t,M=t=>{if(!t)throw a("a scan without its Huffman table");b();const o=t[u>>>w-16&65535];if(!o)throw a("bad Huffman code");if(o>>8>w-d)throw y();return w-=o>>8,255&o},A=()=>{const o=w-d;if(o>=8)throw a("bytes left over in a scan");if(v)throw a("an end-of-band run past its interval");if(o&&(u>>>d&(1<<o)-1)!=(1<<o)-1)throw a("bad scan padding");for(;l+1<t.length&&255===t[l]&&255===t[l+1];)l++},U=o=>{if(A(),255!==t[l]||t[l+1]!==208+(7&o))throw a("a restart marker missing or out of order");l+=2,u=0,w=0,d=0,p=!1;for(const t of r)t.component.pred=0;v=0},x=(t,o)=>{if(o>11)throw a("a DC value out of range");const n=o?g(m(o),o):0;if(t.component.pred+=n,t.component.pred>2047||t.component.pred<-2048)throw a("a DC value out of range");return t.component.pred},I=(t,o)=>{if(t>10)throw a("an AC value out of range");const n=g(m(t),t)*(1<<o);if(n>1023||n<-1023)throw a("an AC value out of range");return n};let v=0;const E=1===r.length,S=n.progressive,C=()=>a("coefficients past the "+(S?"band":"block")),J=(t,o,n)=>{if(v>0)v--;else for(let r=e||1;r<=s;){const e=M(t.ac),i=e>>4,h=15&e;if(h){if(r+=i,r>s)throw C();o[n+c[r]]=I(h,f),r++}else{if(i<15){if(i&&!S)throw a("an end-of-band run in a sequential scan");v=(1<<i)-1+(i?m(i):0);break}if(r+=16,r>s+1)throw C()}}},k=S?0===e?0===i?(t,o,n)=>{const r=x(t,M(t.dc))*(1<<f);if(r>2047||r<-2048)throw a("a DC value out of range");o[n]=r}:(t,o,n)=>{m(1)&&(o[n]|=1<<f)}:0===i?J:(t,o,n)=>{const r=1<<f,i=-1<<f;let h=e;if(v<=0)for(;h<=s;h++){const e=M(t.ac);let f=e>>4,l=0;const u=15&e;if(u){if(1!==u)throw a("a bad refinement code");l=m(1)?r:i}else if(15!==f){v=1<<f,f&&(v+=m(f));break}for(;h<=s;h++){const t=n+c[h];if(0!==o[t])m(1)&&0===(o[t]&r)&&(o[t]+=o[t]>=0?r:i);else if(--f<0)break}if(f>=0||l&&h>s)throw a("coefficients past the band");l&&(o[n+c[h]]=l)}if(v>0){for(;h<=s;h++){const t=n+c[h];0!==o[t]&&m(1)&&0===(o[t]&r)&&(o[t]+=o[t]>=0?r:i)}v--}}:(t,o,n)=>{o[n]=x(t,M(t.dc)),J(t,o,n)};for(const t of r)t.component.pred=0;const L=r[0].component,X=E?L.blocksW*L.blocksH:n.mcusX*n.mcusY;for(let t=0;t<X;t++)if(h&&t>0&&t%h===0&&U(t/h-1),E){const o=t/L.blocksW|0,n=t-o*L.blocksW;k(r[0],L.coeffs,64*(o*L.stride+n))}else{const o=t/n.mcusX|0,e=t-o*n.mcusX;for(const t of r){const n=t.component;for(let r=0;r<n.v;r++)for(let s=0;s<n.h;s++)k(t,n.coeffs,64*((o*n.v+r)*n.stride+e*n.h+s))}}if(A(),255===t[l]&&t[l+1]>=208&&t[l+1]<=215){if(!h||X%h||t[l+1]!==208+(X/h-1&7))throw a("a restart marker out of order");for(l+=2;l+1<t.length&&255===t[l]&&255===t[l+1];)l++}if(l<t.length&&!(255===t[l]&&l+1<t.length&&0!==t[l+1]&&255!==t[l+1]))throw a("bytes left over after a scan");return l}const y=256;function m(t,o){t.o([[9,1],[13,1],[18,1],[30,1]],o)}const g=[1,2,3,4,0,5,17,6,16,7,8,9,10,11,12,13,14,15],M=[2,4,3,2,2,4],A=[0,7,3,2,1,15];function U(t,o){const n=new Uint8Array(t.length),r=[];for(let o=0;o<t.length;o++)t[o]>0&&r.push(o);if(!r.length)return n;if(1===r.length)return n[r[0]]=1,n;for(let e=1;;e*=2){const s=r.map(o=>({w:Math.max(t[o],e),p:o,l:null,r:null})).sort((t,o)=>t.w-o.w||t.p-o.p),i=[];let f=0,c=0;const a=()=>f<s.length&&(c>=i.length||s[f].w<=i[c].w)?s[f++]:i[c++];for(;s.length-f+(i.length-c)>1;){const t=a(),o=a();i.push({w:t.w+o.w,p:-1,l:t,r:o})}const h=[[a(),0]];let l=0;for(;h.length;){const[t,o]=h.pop();t.p>=0?(n[t.p]=o,o>l&&(l=o)):h.push([t.l,o+1],[t.r,o+1])}if(l<=o)return n}}function x(t,o,n,r,e,s){const i=[];for(;i.push(r&e),r>>=s;)r--;for(let r=i.length-1;r>=0;r--)t.push(n),o.push(i[r])}function I(t,o){let n=0;for(let r=0;r<o;r++)n=n<<1|1&t,t>>=1;return n}function v(t){const o=new Uint32Array(17),n=new Uint32Array(17),r=new Uint16Array(t.length);for(const n of t)o[n]++;o[0]=0;for(let t=1,r=0;t<=16;t++)r=r+o[t-1]<<1,n[t]=r;for(let o=0;o<t.length;o++)t[o]&&(r[o]=I(n[t[o]]++,t[o]));return r}function E(t){const o=[];for(let n=0;n<t.length;n++)t[n]>0&&o.push(n);const n=o.length,r=o[n-1]+1||1,e=n<=4?new Uint8Array(r):U(t,15).subarray(0,r);let s=0;if(n<=4){o.sort((o,n)=>t[n]-t[o]||o-n),s=4===n&&t[o[0]]>t[o[2]]+t[o[3]]?1:0;const r=n<2?[0]:2===n?[1,1]:3===n?[1,2,2]:s?[1,2,3,3]:[2,2,2,2];o.forEach((t,o)=>{e[t]=r[o]}),o.sort((t,o)=>e[t]-e[o]||t-o)}return{m:r,M:e,A:v(e),U:n&&n<=4?o:null,I:s}}function S(t,o){if(o.m<=1)return;if(o.U){const n=s(o.m-1)+1;t.write(2,1),t.write(2,o.U.length-1);for(const r of o.U)t.write(n,r);return void(4===o.U.length&&t.write(1,o.I))}const n=o.M,r=[],e=[];let i=n.length-1;for(;i>0&&!n[i];)i--;let f=8;for(let t=0;t<=i;){const o=n[t];let s=1;for(;t+s<=i&&n[t+s]===o;)s++;let c=s;if(o&&o!==f&&(r.push(o),e.push(0),c--,f=o),c<3)for(let t=0;t<c;t++)r.push(o),e.push(0);else x(r,e,o?16:17,c-3,o?3:7,o?2:3);t+=s}const c=new Uint32Array(18);for(const t of r)c[t]++;let a=0;for(const t of c)t&&a++;1===a&&(c[r[0]?0:1]=c[r[0]]);const h=U(c,5),l=v(h);let u=18;for(;u>0&&!h[g[u-1]];)u--;t.write(2,0);for(let o=0;o<u;o++){const n=h[g[o]];t.write(M[n],A[n])}for(let o=0;o<r.length;o++){const n=r[o];t.write(h[n],l[n]),16===n?t.write(2,e[o]):17===n&&t.write(3,e[o])}}function C(t,o=0,n=0){return{split:t,S:o,C:n,J:1<<t}}function J(t,o,n=15){t.write(i(n+1),o.split),o.split!==n&&(t.write(i(o.split+1),o.S),t.write(i(o.split-o.S+1),o.C))}function k(t,o,n){if(o<t.J)return n[0]=o,n[1]=0,void(n[2]=0);const r=s(o),e=o-(1<<r),i=r-t.S-t.C;n[0]=t.J+(r-t.split<<t.S+t.C|e>>r-t.S<<t.C|e&(1<<t.C)-1),n[1]=i,n[2]=o>>t.C&(1<<i)-1}function L(t,o,n,r=0){k(t,o,X),n[r+X[0]]++}const X=[0,0,0];function O(t,o,n,r,e=0){k(n,r,X);const s=e+X[0];t.write(o.M[s],o.A[s]),X[1]&&t.write(X[1],X[2])}function N(t,o){let n=0;for(const t of o)n=Math.max(n,t);const r=n?i(n+1):0;if(r<=3){if(t.write(1,1),t.write(2,r),r)for(const n of o)t.write(r,n);return}t.write(1,0),t.write(1,0);const e=C(8),s=C(4),f=o.length>=64,c=new Uint32Array(f?257:n+1),a=[];for(let t=0;t<o.length;){let n=1;for(;f&&t+n<o.length&&o[t+n]===o[t];)n++;if(n>=8)a.push([o[t],n-1]),c[o[t]]++,L(s,n-1-7,c,224),t+=n;else{for(let r=0;r<n;r++)a.push([o[t],0]),c[o[t]]++;t+=n}}const h=E(c);if(!f){D(t,{k:new Uint8Array(1),L:[{config:e,code:h}]});for(const[o]of a)O(t,h,e,o);return}const l=new Uint32Array(1);l[0]=1,D(t,{X:{O:224,minLength:7,N:s},k:new Uint8Array([1,0]),L:[{config:C(0),code:E(l)},{config:e,code:h}]});for(const[o,n]of a)O(t,h,e,o),n&&O(t,h,s,n-7,224)}function D(t,{X:o=null,k:n,L:r}){t.write(1,o?1:0),o&&(t.o([[0,224],[0,512],[0,4096],[15,8]],o.O),t.o([[0,3],[0,4],[2,5],[8,9]],o.minLength),J(t,o.N,8)),n.length>1&&N(t,n),t.write(1,1);for(const o of r)J(t,o.config,15);for(const o of r){const n=o.code.m-1;if(n){const o=s(n);t.write(1,1),t.write(4,o),t.write(o,n-(1<<o))}else t.write(1,0)}for(const o of r)S(t,o.code)}const G=C(4),T=Object.freeze({O:224,minLength:7,N:G}),P=C(0);function _(t,o){const n=[[3,0],[6,8],[10,72],[13,1096]];if("rct"===o.type)t.write(2,0),t.o(n,o.D),t.o([[0,6],[2,0],[4,2],[6,10]],o.G??6);else if("palette"===o.type)t.write(2,1),t.o(n,o.D),t.o([[0,1],[0,3],[0,4],[13,1]],o.T),t.o([[8,0],[10,256],[12,1280],[16,5376]],o.P),t.o([[0,0],[8,1],[10,257],[16,1281]],o._??0),t.write(4,o.R??0);else{if("squeeze"!==o.type)throw new Error("unknown transform "+o.type);{t.write(2,2);const r=o.q??[];t.o([[0,0],[4,1],[6,9],[8,41]],r.length);for(const o of r)t.write(1,o.B?1:0),t.write(1,o.Z?1:0),t.o(n,o.D),t.o([[0,1],[0,2],[0,3],[4,4]],o.T)}}}function R(t,o,n,r,s,i){const f=i.R,c=i.offset,a=i.Y,h=P,l=G;let u=0;const w=(n,r=h,e=0)=>{t?O(t,o,r,n,e):L(r,n,o,e)},d=()=>{if(u){if(u>=8)w(0),w(u-7-1,l,224);else for(let t=0;t<u;t++)w(0);u=0}},p=(t,o)=>{const r=n[t];let s;1===a?s=r-o-c:(s=Math.round((r-o-c)/a),n[t]=o+c+s*a),0!==s?(d(),w(e(s))):u++};for(let t=0,o=0;t<s;t++)for(let e=0;e<r;e++,o++){let s;if(0===f)s=0;else{const i=e?n[o-1]:t?n[o-r]:0,c=t?n[o-r]:i;if(1===f)s=i;else if(2===f)s=c;else if(3===f)s=(i+c)/2|0;else{const a=i+c-(e&&t?n[o-r-1]:i);if(5===f){const t=i<c?i:c,o=i<c?c:i;s=a<t?t:a>o?o:a}else{if(4!==f)throw new Error("predictor "+f+" is not coded here");s=Math.abs(a-i)<Math.abs(a-c)?i:c}}}p(o,s)}d()}const q=[C(0),C(1),C(2),C(3),C(4),C(4,1,1),C(5,1,1),C(2,0,1),C(3,0,1)],B=C(4,1,1);class Z{constructor(t){this.j=t,this.small=new Uint32Array(64*t),this.H=Array.from({length:t},()=>null),this.F=new Float64Array(t)}add(t,o){if(this.F[t]++,o<64)return void this.small[64*t+o]++;const n=this.H[t]||(this.H[t]=new Map);n.set(o,(n.get(o)||0)+1)}K(t,o,n=64){const r=new Uint32Array(n);for(const n of t){const t=64*n;for(let n=0;n<64;n++)this.small[t+n]&&(r[Y(o,n)]+=this.small[t+n]);const e=this.H[n];if(e)for(const[t,n]of e)r[Y(o,t)]+=n}return r}V(t,o){let n=F(this.K(t,o));for(const r of t){const t=64*r;for(let r=o.J;r<64;r++)this.small[t+r]&&(n+=this.small[t+r]*(s(r)-o.S-o.C));const e=this.H[r];if(e)for(const[t,r]of e)n+=r*(s(t)-o.S-o.C)}return n}}const z=[0,0,0];function Y(t,o){return k(t,o,z),z[0]}const j=[1];for(let t=1;t<32;t++)j.push(j[t-1]/2);function H(t){let o=0;for(t>=2&&t<4294967296&&(o=31-Math.clz32(t),t*=j[o]);t>=1.4142135623730951;)t/=2,o++;for(;t<.7071067811865476;)t*=2,o--;const n=(t-1)/(t+1),r=n*n;let e=n,s=n;for(let t=3;t<30;t+=2)s+=(e*=r)/t;return o+2.885390081777927*s}function F(t){let o=0,n=0;for(const n of t)o+=n;for(const r of t)r&&(n+=r*H(o/r));return n}const K=[0,0,1,2,3,4,5,6,7,8,9,10,11,12,13,14,15,15,16,16,17,17,18,18,19,19,20,20,21,21,22,22,23,23,23,23,24,24,24,24,25,25,25,25,26,26,26,26,27,27,27,27,28,28,28,28,29,29,29,29,30,30,30,30],Q=[0,0,31,62,62,93,93,93,93,123,123,123,123,152,152,152,152,152,152,152,152,180,180,180,180,180,180,180,180,180,180,180,180,206,206,206,206,206,206,206,206,206,206,206,206,206,206,206,206,206,206,206,206,206,206,206,206,206,206,206,206,206,206,206],V=c.map(t=>(7&t)<<3|t>>3),W=(t,o,n)=>2*(Q[t]+K[o])+n;function $(t,o,n){const r=[];if(o>0){const n=new Uint32Array(2048);let e=0;for(const o of t)n[Math.min(2047,Math.max(0,o+1024))]++,e++;let s=0,i=e/(o+1);for(let t=0;t<2048&&r.length<o;t++)s+=n[t],s>i&&(r.push(t-1025),i=e*(r.length+1)/(o+1))}const e=r.length+1,s=new Uint8Array(39*e);for(let t=0;t<e;t++)for(let o=0;o<13;o++)s[(0+o)*e+t]=t,s[(13+o)*e+t]=n?e:e+(t>>1),s[(26+o)*e+t]=n?e:e+(e-1>>1)+1+(t>>1);let i=0;for(const t of s)i=Math.max(i,t+1);return{thresholds:r,W:e,map:s,$:i,tt:t=>{let o=0;for(const n of r)t>n&&o++;return o},ot:(t,o)=>s[13*(t<2?1^t:2)*e+o],nt:(t,o)=>(Math.min(t,64)<8?Math.min(t,64):4+(Math.min(t,64)>>1))*i+o,rt:t=>37*i+458*t,j:495*i}}function tt(t,o,n){const r=o.map(()=>function(t,o=0,n=1){return{R:t,offset:o,Y:n,context:-1}}(n)),s=r.map(()=>new Uint32Array(257));o.forEach((t,o)=>R(null,s[o],t.data,t.w,t.h,r[o])),function(t,{et:o=!0,st:n=[]}={}){t.write(1,o?1:0),t.write(1,1),t.o([[0,0],[0,1],[4,2],[8,18]],n.length);for(const o of n)_(t,o)}(t,{et:!1,st:[]});const i=function(t,o){const n=[o],r=[],s=[];for(const t of n){if(t.left){r.push(t.it+1,e(t.ft)),n.push(t.left,t.right);continue}t.context=s.length,s.push(t);const o=1===t.Y?0:31-Math.clz32(t.Y&-t.Y);r.push(0,t.R,e(t.offset),o,(t.Y>>o)-1)}const i=new Uint32Array(64);for(const t of r)L(P,t,i);const f=E(i);D(t,{k:new Uint8Array(6),L:[{config:P,code:f}]});for(const o of r)O(t,f,P,o);return s}(t,function(t){return function(t,o){const n=(r,e)=>{if(r===e)return o?t[r].ct:t[r];const s=r+e>>1;return function(t,o,n,r){return{it:t,ft:o,left:n,right:r}}(o,o?t[s].ht:s,n(s+1,e),n(r,s))};return n(0,t.length-1)}(t,0)}(r)),f=function(t,o,n,r=t=>t.context){const e=new Uint8Array(o.length+1);for(const t of o)e[t.context]=r(t)+1;const s=new Uint32Array(2);s[1]=1;const i=[{config:P,code:E(s)}],f=Math.max(...e);for(let t=0;t<f;t++)i.push({config:P,code:E(n[t])});return D(t,{X:T,k:e,L:i}),i}(t,i,s,t=>r.indexOf(t));o.forEach((o,n)=>R(t,f[n+1].code,o.data,o.w,o.h,r[n]))}const ot=(t,o)=>Object.assign(new Error(o),{code:t});function nt(t){try{return t()}catch(t){if(t instanceof RangeError&&!t.code)throw ot("JXL_MEMORY","Not enough memory for this picture.");throw t}}function rt(t){const o=et(t);return function(t){let o;for(;!(o=t.next()).done;);o.value}(o),{bytes:o.bytes,width:o.width,height:o.height,orientation:o.orientation}}function et(n){if(!(n instanceof Uint8Array&&n.length))throw ot("JXL_INPUT","A JPEG is a non-empty Uint8Array.");if(n.length>o.bytes)throw ot("JXL_SIZE","A JPEG is at most 16 MiB.");const s=function(o){const n=function*(){let r,e;for(;!(r=nt(()=>o.next(n.hurry))).done;)yield e=r.value;return e<1&&(yield 1),n.bytes=function(o){if(o.length>t.bytes)throw ot("JXL_SIZE","The encoded picture exceeds 16 MiB.");return o}(r.value)}();return n.bytes=null,n.hurry=!1,n}(function*(n){const g=function(t){if(!(t.length>4)||255!==t[0]||216!==t[1])throw a("not a JPEG");const n=[],r=[[],[]],e=new Set,s={count:0,lt:[]};let i=null,f=0,y=-1,m=!1,g=1,M=2;const A=o=>t[o]<<8|t[o+1],U=(o,n)=>{for(let r=0;r<n.length;r++)if(t[o+r]!==n.charCodeAt(r))return!1;return!0};for(;;){if(M>=t.length){if(i&&i.scans)break;throw a("truncated")}if(255!==t[M])throw a("marker expected");for(;255===t[M];)M++;const l=t[M++];if(217===l)break;if(216===l)throw a("a second start");if(l>=208&&l<=215)throw a("a restart marker outside a scan");if(1===l)continue;if(M+2>t.length)throw a("truncated");const u=A(M),w=M+2,d=M+u;if(u<2||d>t.length)throw a("bad segment");if(M=d,219===l)for(let o=w;o<d;){const r=t[o]>>4,e=15&t[o];if(r>1||e>3||o+1+(r?128:64)>d)throw a("bad quantisation table");o++;const s=new Int32Array(64);for(let n=0;n<64;n++)if(s[c[n]]=r?A(o):t[o],o+=r?2:1,!s[c[n]])throw a("a zero quantisation step");if(i?.progressive&&i.components.some(t=>t.tq===e&&t.quant&&t.coverage.some(t=>0!==t)&&t.quant.some((t,o)=>t!==s[o])))throw a("a quantisation table changed between progressive scans");n[e]=s}else if(196===l)for(let o=w;o<d;){const n=t[o]>>4,s=15&t[o];if(n>1||s>3||o+17>d)throw a("bad Huffman table");o++;const i=t.subarray(o,o+16);let f=0;for(const t of i)f+=t;if(f>256||o+16+f>d)throw a("bad Huffman table");const c=t.subarray(o+16,o+16+f);r[n][s]=h(i,c),0===n&&(c.some(t=>t>15)?e.add(s):e.delete(s)),o+=16+f}else if(192===l||193===l||194===l){if(i)throw a("two frames");if(d-w<6)throw a("bad frame header");if(8!==t[w])throw a("only 8-bit samples");const n=A(w+1),r=A(w+3),e=t[w+5];if(!r||!n)throw a("a height left to a DNL marker");if(1!==e&&3!==e)throw a("only grey or three components");if(d-w<6+3*e)throw a("bad frame header");if(r>o.edge||n>o.edge)throw a("a picture past "+o.edge+" pixels on a side","JXL_DIMENSIONS");if(r*n>o.pixels)throw a("a picture past "+o.pixels+" pixels","JXL_DIMENSIONS");const s=[];let f=0;for(let o=0,n=w+6;o<e;o++,n+=3){const o=t[n+1]>>4,r=15&t[n+1],e=t[n+2];if(o<1||o>4||r<1||r>4)throw a("bad sampling");if(e>3)throw a("bad quantisation table");if(s.some(o=>o.id===t[n]))throw a("two components with one id");f+=o*r,s.push({id:t[n],h:o,v:r,tq:e,pred:0,coverage:new Int8Array(64).fill(-1)})}if(f>10)throw a("bad sampling");const c=Math.max(...s.map(t=>t.h)),h=Math.max(...s.map(t=>t.v)),u=Math.ceil(r/(8*c)),p=Math.ceil(n/(8*h));for(const t of s)t.blocksW=Math.ceil(Math.ceil(r*t.h/c)/8),t.blocksH=Math.ceil(Math.ceil(n*t.v/h)/8),t.stride=u*t.h,t.rows=p*t.v,t.coeffs=new Int16Array(t.stride*t.rows*64);i={progressive:194===l,width:r,height:n,components:s,hmax:c,vmax:h,mcusX:u,mcusY:p,scans:0}}else{if(195===l||l>=197&&l<=199||l>=201&&l<=203||l>=205&&l<=207)throw a("lossless, hierarchical or arithmetic coding");if(221===l){if(u<4)throw a("bad segment");f=A(w)}else{if(220===l)throw a("a DNL marker");if(218===l){if(!i)throw a("a scan before the frame");const o=t[w],s=[];if(o<1||o>4||w+1+2*o+3>d)throw a("bad scan header");let c=w+1;for(let e=0;e<o;e++,c+=2){const o=i.components.find(o=>o.id===t[c]);if(!o)throw a("a scan of an unknown component");if(s.some(t=>t.component===o))throw a("a component twice in one scan");if(o.quant||(o.quant=n[o.tq]),!o.quant)throw a("a scan without its quantisation table");s.push({component:o,dc:r[0][t[c+1]>>4],ac:r[1][15&t[c+1]],dcSlot:t[c+1]>>4})}const h=t[c],l=t[c+1],u=t[c+2]>>4,p=15&t[c+2];if((!i.progressive||0===h&&0===u)&&s.some(t=>e.has(t.dcSlot)))throw a("bad Huffman table");M=b(t,d,i,s,h,l,u,p,f),i.scans++}else if(224===l)U(w,"JFIF")&&(m=!0);else if(225===l)1===g&&U(w,"Exif")&&(g=p(t,w+6,d));else if(226===l){if(u>=16&&U(w,"ICC_PROFILE\0")){const o=t[w+12],n=t[w+13];if(!o||!n||o>n||s.count&&s.count!==n||s.lt[o])throw a("a colour profile cut short");s.count=n,s.lt[o]=t.subarray(w+14,d)}}else 238===l&&u>=14&&U(w,"Adobe")&&(y=t[w+11])}}}if(!i||!i.scans)throw a("no picture");for(const t of i.components)for(let o=0;o<64;o++)if(0!==t.coverage[o])throw a(i.progressive?"a progressive JPEG whose scans do not finish every coefficient":"a component without its scan");let x="srgb";if(s.count){const t=[];for(let o=1;o<=s.count;o++){if(!s.lt[o])throw a("a colour profile cut short");t.push(s.lt[o])}const o=new Uint8Array(t.reduce((t,o)=>t+o.length,0));for(let n=0,r=0;n<t.length;r+=t[n].length,n++)o.set(t[n],r);if(x=function(t){const o=t.length,n=o=>t[o]<<8|t[o+1],r=n=>n+4<=o?(t[n]<<24|t[n+1]<<16|t[n+2]<<8|t[n+3])>>>0:0,e=n=>n+4<=o?String.fromCharCode(t[n],t[n+1],t[n+2],t[n+3]):"",s=t=>(0|r(t))/65536,i=r(128);if(o<132||132+12*i>o||"acsp"!==e(36)||"XYZ "!==e(20))return null;const f=t=>{for(let n=0;n<i;n++){const s=132+12*n,i=r(s+4);if(e(s)===t)return i+r(s+8)<=o?i:-1}return-1},c=t=>{const i=f(t);if(i<0)return!1;if("para"===e(i)){const t=n(i+8),o=[1,3,4,5,7][t];if(3!==t&&4!==t)return!1;for(let t=0;t<o;t++)if(Math.abs(s(i+12+4*t)-(t<5?d[t]:0))>(t?.002:.01))return!1;return!0}const c="curv"===e(i)?r(i+8):0;return!(c<2||i+12+2*c>o)&&w.every((t,o)=>{const r=o*(c-1)/16,e=Math.floor(r),s=Math.min(c-1,e+1),f=r-e;return Math.abs((n(i+12+2*e)*(1-f)+n(i+12+2*s)*f)/65535-t)<=.001})};if(["A2B0","A2B1","A2B2"].some(t=>f(t)>=0))return null;if("GRAY"===e(16))return c("kTRC")?"srgb":null;if("RGB "!==e(16)||!["rTRC","gTRC","bTRC"].every(c))return null;const a=["rXYZ","gXYZ","bXYZ"].flatMap(t=>{const o=f(t);return o>=0&&"XYZ "===e(o)?[s(o+8),s(o+12),s(o+16)]:[9,9,9]}),h=t=>t.every((t,o)=>Math.abs(a[o]-t)<=.004);return h(l)?"srgb":h(u)?"display-p3":null}(o),!x)throw a("a colour profile other than sRGB or Display P3 (decode the picture and encode its pixels)")}let I=!1;return 3!==i.components.length||m||(I=y>=0?0===y:82===i.components[0].id&&71===i.components[1].id&&66===i.components[2].id),{width:i.width,height:i.height,components:i.components,progressive:i.progressive,ycbcr:!I,orientation:g,colorSpace:x}}(n);!function(o,n,r=t){if(!Number.isInteger(o)||!Number.isInteger(n)||o<1||n<1)throw ot("JXL_INPUT","A picture is at least one pixel wide and high, in whole pixels.");if(o>r.edge||n>r.edge)throw ot("JXL_DIMENSIONS","A picture is at most "+r.edge+" pixels on a side.");if(o*n>r.pixels)throw ot("JXL_DIMENSIONS","A picture is at most "+r.pixels.toLocaleString("en-US")+" pixels.")}(g.width,g.height,o);const M=g.orientation>=5;return(t=>{Object.assign(s,t)})({width:M?g.height:g.width,height:M?g.width:g.height,orientation:g.orientation}),yield*function*(t){const{width:o,height:n,components:s,alpha:c}=t,h=1===s.length,l=t.quantScale||1;let u=1;for(;u/(2040*l)<1/16384;)u*=2;const w=t.ycbcr||h,d=(h?[0,0,0]:w?[1,0,2]:[0,1,2]).map(t=>s[t]),p=t=>1===t?0:2===t?1:-1,b=d.map(t=>p(t.h)),g=d.map(t=>p(t.v));if(b.includes(-1)||g.includes(-1))throw a("sampling factors other than one or two");const M=Math.max(...b),A=Math.max(...g),U=t=>M-b[t],x=t=>A-g[t];if(!w&&(M||A))throw a("an RGB JPEG with subsampling");const I=Math.ceil(o/(8<<M))<<M,v=Math.ceil(n/(8<<A))<<A,S=Math.ceil(o/y),C=Math.ceil(n/y),J=S*C,k=Math.ceil(I/256),L=k*Math.ceil(v/256),X=1===J,G=(t,o,n)=>{const r=d[t];return h&&1!==t?-1:64*(n*r.stride+o)},T=t=>{let o=1024;for(;t%o;)o>>=1;return o},P=d.map(t=>w?1:t.quant[0]/T(t.quant[0])),_=d.map(t=>w?0:1024/T(t.quant[0])),R=new Int32Array(I*v);for(let t=0;t<v;t++)for(let o=0;o<I;o++)R[t*I+o]=d[1].coeffs[G(1,o,t)]*P[1]+_[1];const z=d[1].quant,Y=(z[8]+z[16]+z[24]+z[32]+z[40])/l,j=I*v,H=j<256?0:Math.max(1,Math.min(7,i(j)-i(Y)-7));let K=$(R,H,h);const Q=new r(256);if(function(t,o,n,r,e,{xyb:s=!1,orientation:i=1,colorSpace:f}={}){t.write(16,2815),t.write(1,0),m(t,n),t.write(3,0),m(t,o),t.write(1,0);const c=1!==i;t.write(1,c?1:0),c&&t.write(6,i-1),t.write(4,8),e?t.write(3,5):t.write(2,0),t.write(1,s?1:0),3===r&&"display-p3"===f?t.write(21,903568):3===r?t.write(1,1):t.write(15,14100),c&&t.write(1,1),t.write(2,0),t.write(1,1),t.i()}(Q,o,n,h?1:3,!!c,{orientation:t.orientation,colorSpace:t.colorSpace}),Q.write(1,0),Q.write(2,0),Q.write(1,0),Q.write(2,2),Q.write(8,111),Q.write(1,w?1:0),w)for(let t=0;t<3;t++)Q.write(2,0===b[t]&&0===g[t]?0:1===b[t]&&1===g[t]?1:1===b[t]?2:3);Q.write(2,0),c&&Q.write(2,0),function(t,o){t.write(o?7:5,0),t.write(11,1)}(Q,c);const ot=[0,1,2].map(()=>new Int32Array(1024)),nt=(t,o)=>{const n=t%S*32,r=32*(t/S|0),s=Math.min(32,I-n),i=Math.min(32,v-r);for(let t=0;t<i;t++)for(let i=0;i<s;i++)for(const s of[1,0,2]){const f=U(s),c=x(s),a=i>>f,h=t>>c;if(a<<f!==i||h<<c!==t)continue;const l=G(s,(n>>f)+a,(r>>c)+h),u=d[s].coeffs;let w=0;if(l>=0)for(let t=1;t<64;t++)u[l+t]&&w++;const p=ot[s],b=0===a?0===h?32:p[32*(h-1)]:0===h?p[32*h+a-1]:p[32*(h-1)+a]+p[32*h+a-1]+1>>1;p[32*h+a]=w;const y=K.ot(s,K.tt(R[(r+t)*I+n+i]));o(K.nt(b,y),w);const m=K.rt(y);let g=w>4?0:1,M=w;for(let t=1;t<64&&M;t++){const n=u[l+V[t]];o(m+W(M,t,g),e(n)),g=n?1:0,M-=g}}},rt=(H>0?3:2)*J+L;let et=0;const st=function*(t){K=t;const o=new Z(t.j);for(let t=0;t<J;t++)nt(t,(t,n)=>o.add(t,n)),yield++et/rt;return function(t,{maxClusters:o=48,newClusterCost:n=320}={}){const r=[];for(let o=0;o<t.j;o++)t.F[o]>0&&r.push(o);r.sort((o,n)=>t.F[n]-t.F[o]);const e=[];for(const s of r){const r=t.K([s],B),i=F(r);let f=-1,c=1/0;for(let t=0;t<e.length;t++){const o=e[t],n=new Uint32Array(64);for(let t=0;t<64;t++)n[t]=o.ut[t]+r[t];const s=F(n)-o.wt-i;s<c&&(c=s,f=t)}if(f<0||e.length<o&&c>n)e.push({j:[s],ut:r,wt:i});else{const t=e[f];t.j.push(s);for(let o=0;o<64;o++)t.ut[o]+=r[o];t.wt=F(t.ut)}}e.length||e.push({j:[],ut:new Uint32Array(64),wt:0});const s=new Uint8Array(t.j);let i=2*t.j;const f=e.map((o,n)=>{for(const t of o.j)s[t]=n;let r=q[0],e=1/0;for(const n of q){const s=t.V(o.j,n);s<e&&(e=s,r=n)}const f=E(t.K(o.j,r,128));return i+=e+3*f.m+40,{config:r,code:f}});return{k:s,L:f,write:(t,o,n)=>{const r=f[s[o]];O(t,r.code,r.config,n)},wt:i}}(o)};let it=yield*st(K);if(H>0){const t=K,o=$(R,0,h),n=yield*st(o);n.wt<=it.wt?it=n:K=t}const ft=Array.from({length:X?1:2+L+J},()=>null),ct=t=>ft[X?0:t]||(ft[X?0:t]=new r(4096)),at=ct(0);at.write(1,0);const ht=(t,r)=>{const e=t*y,s=r*y,i=Math.min(y,o-e),f=Math.min(y,n-s),a=new Int32Array(i*f);for(let t=0;t<f;t++)for(let n=0;n<i;n++)a[t*i+n]=c[4*((s+t)*o+e+n)+3];return{w:i,h:f,data:a}};for(let t=0;t<3;t++)at.write(16,f(d[t].quant[0]/P[t]/(2040*l)*128));at.write(2,3),at.write(16,57343),at.write(2,1),at.write(5,0),at.write(1,0),at.write(4,0),at.write(4,K.thresholds.length);for(const t of K.thresholds)at.o([[4,0],[8,16],[16,272],[32,65808]],e(t));at.write(4,0),at.write(4,0),N(at,K.map),at.write(1,0),at.write(2,0),at.write(16,0),at.write(16,0),at.write(8,128),at.write(8,128),at.write(1,0),c&&tt(at,[X?ht(0,0):{w:0,h:0,data:new Int32Array(0)}],5);for(let t=0;t<L;t++){const o=ct(1+t),n=t%k*256,r=256*(t/k|0),e=Math.min(256,I-n),s=Math.min(256,v-r);o.write(2,0),tt(o,[1,0,2].map(t=>{const o=e>>U(t),i=s>>x(t),f=new Int32Array(o*i),c=P[t],a=_[t];for(let e=0;e<i;e++)for(let s=0;s<o;s++){const i=G(t,(n>>U(t))+s,(r>>x(t))+e);f[e*o+s]=i<0?0:d[t].coeffs[i]*c+a}return{w:o,h:i,data:f}}),5);const f=e*s,c=i(f);c&&o.write(c,f-1);const a=e+7>>3,h=s+7>>3,l=new Int32Array(2*f).fill(u-1,f);tt(o,[{w:a,h:h,data:new Int32Array(a*h)},{w:a,h:h,data:new Int32Array(a*h)},{w:f,h:2,data:l},{w:e,h:s,data:new Int32Array(f)}],u>1?5:0),yield++et/rt}const lt=ct(1+L);lt.write(1,0),lt.write(3,7),lt.write(16,f(u/(2040*l))),tt(lt,[0,1,2].map(t=>{const o=new Int32Array(64),n=d[t].quant;for(let t=0;t<8;t++)for(let r=0;r<8;r++)o[8*r+t]=n[8*t+r];return{w:8,h:8,data:o}}),5);for(let t=1;t<17;t++)lt.write(3,0);J>1&<.write(i(J),0),lt.write(2,2),D(lt,{k:it.k,L:it.L});for(let t=0;t<J;t++){const o=ct(2+L+t);nt(t,(t,n)=>it.write(o,t,n)),c&&!X&&tt(o,[ht(t%S,Math.floor(t/S))],5),yield++et/rt}return function(t,o){const n=o.map(t=>t.length);!function(t,o){t.write(1,0),t.i();for(const n of o)t.o([[10,0],[14,1024],[22,17408],[30,4211712]],n);t.i()}(t,n);const r=t.finish(),e=new Uint8Array(r.length+n.reduce((t,o)=>t+o,0));e.set(r);let s=r.length;for(const t of o)e.set(t,s),s+=t.length;return e}(Q,function(t){return t.map(t=>t?t.finish():new Uint8Array(0))}(ft))}(g)}(n));return s}export{o as LIMITS,rt as transcode,et as transcodeSteps};
|
|
1
|
+
/*! Rapier JXL 2.1.0 | MIT License Copyright (c) 2026 rapier.website Permission is hereby granted, free of charge, to any person obtaining a copy of this software and associated documentation files (the "Software"), to deal in the Software without restriction, including without limitation the rights to use, copy, modify, merge, publish, distribute, sublicense, and/or sell copies of the Software, and to permit persons to whom the Software is furnished to do so, subject to the following conditions: The above copyright notice and this permission notice shall be included in all copies or substantial portions of the Software. THE SOFTWARE IS PROVIDED "AS IS", WITHOUT WARRANTY OF ANY KIND, EXPRESS OR IMPLIED, INCLUDING BUT NOT LIMITED TO THE WARRANTIES OF MERCHANTABILITY, FITNESS FOR A PARTICULAR PURPOSE AND NONINFRINGEMENT. IN NO EVENT SHALL THE AUTHORS OR COPYRIGHT HOLDERS BE LIABLE FOR ANY CLAIM, DAMAGES OR OTHER LIABILITY, WHETHER IN AN ACTION OF CONTRACT, TORT OR OTHERWISE, ARISING FROM, OUT OF OR IN CONNECTION WITH THE SOFTWARE OR THE USE OR OTHER DEALINGS IN THE SOFTWARE. */
|
|
2
|
+
const t=Object.freeze({bytes:16777216,pixels:24e6,edge:16384}),o=Object.freeze({bytes:16777216,pixels:64e6,edge:16384}),n=[1];for(let t=1;t<=32;t++)n.push(2*n[t-1]);class r{constructor(t=4096){this.bytes=new Uint8Array(t),this.at=0,this.t=0,this.pending=0}write(t,o){if(!(o>=0&&o<n[t]))throw new Error("bit write out of range: "+t+" bits, "+o);let r=this.t+o*(1<<this.pending),e=this.pending+t;this.at+5>=this.bytes.length&&this.grow();const i=this.bytes;for(;e>=8;)i[this.at++]=255&r,r=Math.floor(r/256),e-=8;this.t=r,this.pending=e}o(t,o){for(let r=0;r<4;r++){const[e,i]=t[r];if(o>=i&&o-i<n[e])return this.write(2,r),void(e&&this.write(e,o-i))}throw new Error("U32 value out of range: "+o)}i(){this.pending&&(this.at+1>=this.bytes.length&&this.grow(),this.bytes[this.at++]=this.t,this.t=0,this.pending=0)}get u(){return 8*this.at+this.pending}grow(o=0){const n=t.bytes+64;if(this.at+o+16>n)throw Object.assign(new Error("The encoded picture exceeds 16 MiB."),{code:"JXL_SIZE"});const r=new Uint8Array(Math.min(n,Math.max(2*this.bytes.length,this.at+o+16)));r.set(this.bytes.subarray(0,this.at)),this.bytes=r}append(t){if(this.at+t.at+8>=this.bytes.length&&this.grow(t.at+8),this.pending)for(let o=0;o<t.at;o++)this.write(8,t.bytes[o]);else this.bytes.set(t.bytes.subarray(0,t.at),this.at),this.at+=t.at;t.pending&&this.write(t.pending,t.t)}finish(){const t=this.at+(this.pending?1:0),o=this.bytes.slice(0,t);return this.pending&&(o[this.at]=this.t),o}}function e(t){return t>=0?2*t:2*-t-1}function i(t){return 31-Math.clz32(t)}function s(t){return t<=1?0:32-Math.clz32(t-1)}function f(t){if(0===t)return 0;const o=t<0?32768:0;if(!((t=Math.abs(t))<65520))throw new Error("half float out of range: "+t);let n=0,r=t;for(;r>=2;)r/=2,n++;for(;r<1;)r*=2,n--;n<-14?(r=16384*t,n=-15):r-=1;let e=Math.round(1024*r);return 1024===e&&(e=0,n++),o|n+15<<10|e}const c=[0,1,8,16,9,2,3,10,17,24,32,25,18,11,4,5,12,19,26,33,40,48,41,34,27,20,13,6,7,14,21,28,35,42,49,56,57,50,43,36,29,22,15,23,30,37,44,51,58,59,52,45,38,31,39,46,53,60,61,54,47,55,62,63];function a(t,o="JXL_JPEG"){return Object.assign(new Error("JPEG: "+t),{code:o})}function l(t,o){const n=new Uint16Array(65536);let r=0,e=0;for(let i=1;i<=16;i++){for(let s=0;s<t[i-1];s++){const t=o[e++],s=r<<16-i,f=1<<16-i;if(s+f>=65536)throw a("bad Huffman table");n.fill(i<<8|t,s,s+f),r++}r<<=1}return n}const h=[.436,.2225,.0139,.3851,.7169,.0971,.143,.0606,.7139],u=[.5151,.2412,-.001,.292,.6922,.0419,.1571,.0666,.7841],w=[0,.00516,.01435,.02934,.05088,.07958,.11602,.16068,.21404,.27652,.34851,.43039,.52252,.62523,.73884,.86367,1],d=[2.4,1/1.055,.055/1.055,1/12.92,.04045];function p(t){if(t<=.04045)return t/12.92;const o=(t+.055)/1.055;let n=1;for(let t=0;t<24;t++)n=(4*n+o/(n*n*n*n))/5;return o*o*n*n}function y(t,o,n){const r=77===t[o]&&77===t[o+1];if(!r&&(73!==t[o]||73!==t[o+1]))return 1;const e=o=>o+2>n?0:r?t[o]<<8|t[o+1]:t[o]|t[o+1]<<8;if(42!==e(o+2))return 1;const i=o+(o=>o+4>n?0:r?(t[o]<<24|t[o+1]<<16|t[o+2]<<8|t[o+3])>>>0:(t[o]|t[o+1]<<8|t[o+2]<<16|t[o+3]<<24)>>>0)(o+4),s=e(i);for(let t=0,o=i+2;t<s&&o+12<=n;t++,o+=12){if(274!==e(o))continue;const t=e(o+8);return t>=1&&t<=8?t:1}return 1}function b(t,o,n,r,e,i,s,f,l){!function(t,o,n,r,e,i){if(!t.progressive){if(0!==n||63!==r||0!==e||0!==i)throw a("a bad spectral selection");for(const t of o){if(t.component.coverage.some(t=>-1!==t))throw a("a component coded twice");t.component.coverage.fill(0)}return}if(i>13||e&&e!==i+1)throw a("a bad successive approximation");if(0===n){if(0!==r)throw a("a bad spectral selection");for(const{component:t}of o){if(t.coverage[0]!==(e||-1))throw a(e?"a refinement scan out of order":"a DC scan twice");t.coverage[0]=i}return}if(1!==o.length)throw a("an interleaved AC scan");if(r>63||n>r)throw a("a bad spectral selection");const{coverage:s}=o[0].component;if(-1===s[0])throw a("an AC scan before the DC scan");for(let t=n;t<=r;t++){if(s[t]!==(e||-1))throw a(e?"a refinement scan out of order":"a band coded twice");s[t]=i}}(n,r,e,i,s,f);let h=o,u=0,w=0,d=0,p=!1;const y=()=>{for(;w<=24;){let o=0;if(p||h>=t.length)p=!0,d+=8;else if(o=t[h],255===o){const n=t[h+1];if(0===n)h+=2;else{if(255===n){h++;continue}p=!0,o=0,d+=8}}else h++;u=(u<<8|o)>>>0,w+=8}},b=()=>a("a scan cut short"),m=t=>{if(!t)return 0;if(y(),t>w-d)throw b();const o=u>>>w-t&(1<<t)-1;return w-=t,o},M=(t,o)=>t<1<<o-1?t-(1<<o)+1:t,g=t=>{if(!t)throw a("a scan without its Huffman table");y();const o=t[u>>>w-16&65535];if(!o)throw a("bad Huffman code");if(o>>8>w-d)throw b();return w-=o>>8,255&o},A=()=>{const o=w-d;if(o>=8)throw a("bytes left over in a scan");if(E)throw a("an end-of-band run past its interval");if(o&&(u>>>d&(1<<o)-1)!=(1<<o)-1)throw a("bad scan padding");for(;h+1<t.length&&255===t[h]&&255===t[h+1];)h++},U=o=>{if(A(),255!==t[h]||t[h+1]!==208+(7&o))throw a("a restart marker missing or out of order");h+=2,u=0,w=0,d=0,p=!1;for(const t of r)t.component.pred=0;E=0},x=(t,o)=>{if(o>11)throw a("a DC value out of range");const n=o?M(m(o),o):0;if(t.component.pred+=n,t.component.pred>2047||t.component.pred<-2048)throw a("a DC value out of range");return t.component.pred},I=(t,o)=>{if(t>10)throw a("an AC value out of range");const n=M(m(t),t)*(1<<o);if(n>1023||n<-1023)throw a("an AC value out of range");return n};let E=0;const v=1===r.length,C=n.progressive,S=()=>a("coefficients past the "+(C?"band":"block")),J=(t,o,n)=>{if(E>0)E--;else for(let r=e||1;r<=i;){const e=g(t.ac),s=e>>4,l=15&e;if(l){if(r+=s,r>i)throw S();o[n+c[r]]=I(l,f),r++}else{if(s<15){if(s&&!C)throw a("an end-of-band run in a sequential scan");E=(1<<s)-1+(s?m(s):0);break}if(r+=16,r>i+1)throw S()}}},L=C?0===e?0===s?(t,o,n)=>{const r=x(t,g(t.dc))*(1<<f);if(r>2047||r<-2048)throw a("a DC value out of range");o[n]=r}:(t,o,n)=>{m(1)&&(o[n]|=1<<f)}:0===s?J:(t,o,n)=>{const r=1<<f,s=-1<<f;let l=e;if(E<=0)for(;l<=i;l++){const e=g(t.ac);let f=e>>4,h=0;const u=15&e;if(u){if(1!==u)throw a("a bad refinement code");h=m(1)?r:s}else if(15!==f){E=1<<f,f&&(E+=m(f));break}for(;l<=i;l++){const t=n+c[l];if(0!==o[t])m(1)&&0===(o[t]&r)&&(o[t]+=o[t]>=0?r:s);else if(--f<0)break}if(f>=0||h&&l>i)throw a("coefficients past the band");h&&(o[n+c[l]]=h)}if(E>0){for(;l<=i;l++){const t=n+c[l];0!==o[t]&&m(1)&&0===(o[t]&r)&&(o[t]+=o[t]>=0?r:s)}E--}}:(t,o,n)=>{o[n]=x(t,g(t.dc)),J(t,o,n)};for(const t of r)t.component.pred=0;const X=r[0].component,k=v?X.blocksW*X.blocksH:n.mcusX*n.mcusY;for(let t=0;t<k;t++)if(l&&t>0&&t%l===0&&U(t/l-1),v){const o=t/X.blocksW|0,n=t-o*X.blocksW;L(r[0],X.coeffs,64*(o*X.stride+n))}else{const o=t/n.mcusX|0,e=t-o*n.mcusX;for(const t of r){const n=t.component;for(let r=0;r<n.v;r++)for(let i=0;i<n.h;i++)L(t,n.coeffs,64*((o*n.v+r)*n.stride+e*n.h+i))}}if(A(),255===t[h]&&t[h+1]>=208&&t[h+1]<=215){if(!l||k%l||t[h+1]!==208+(k/l-1&7))throw a("a restart marker out of order");for(h+=2;h+1<t.length&&255===t[h]&&255===t[h+1];)h++}if(h<t.length&&!(255===t[h]&&h+1<t.length&&0!==t[h+1]&&255!==t[h+1]))throw a("bytes left over after a scan");return h}const m=(t,o)=>Object.assign(new Error(o),{code:t});function M(o){if(o>t.bytes)throw m("JXL_SIZE","The encoded picture exceeds 16 MiB.")}function g(t){try{return t()}catch(t){if(t instanceof RangeError&&!t.code)throw m("JXL_MEMORY","Not enough memory for this picture.");throw t}}function A(n,r,e){const i=function(t={}){if(null===t||"object"!=typeof t)throw m("JXL_INPUT","Options are an object.");const o=void 0===t.effort?1:t.effort;if(!Number.isInteger(o)||o<1||o>9)throw m("JXL_INPUT","Effort is a whole number from 1 to 9.");return o}(r);if(!(n instanceof Uint8Array&&n.length))throw m("JXL_INPUT","A JPEG is a non-empty Uint8Array.");if(n.length>o.bytes)throw m("JXL_SIZE","A JPEG is at most 16 MiB.");const s=function(t){const o=function*(){let n,r;for(;!(n=g(()=>t.next(o.hurry))).done;)yield r=n.value;return r<1&&(yield 1),o.bytes=function(t){return M(t.length),t}(n.value)}();return o.bytes=null,o.hurry=!1,o}(function*(n,r,e,i){const f=function(t){if(!(t.length>4)||255!==t[0]||216!==t[1])throw a("not a JPEG");const n=[],r=[[],[]],e=new Set,i={count:0,p:[]};let s=null,f=0,m=-1,M=!1,g=1,A=2;const U=o=>t[o]<<8|t[o+1],x=(o,n)=>{for(let r=0;r<n.length;r++)if(t[o+r]!==n.charCodeAt(r))return!1;return!0};for(;;){if(A>=t.length){if(s&&s.scans)break;throw a("truncated")}if(255!==t[A])throw a("marker expected");for(;255===t[A];)A++;const h=t[A++];if(217===h)break;if(216===h)throw a("a second start");if(h>=208&&h<=215)throw a("a restart marker outside a scan");if(1===h)continue;if(A+2>t.length)throw a("truncated");const u=U(A),w=A+2,d=A+u;if(u<2||d>t.length)throw a("bad segment");if(A=d,219===h)for(let o=w;o<d;){const r=t[o]>>4,e=15&t[o];if(r>1||e>3||o+1+(r?128:64)>d)throw a("bad quantisation table");o++;const i=new Int32Array(64);for(let n=0;n<64;n++)if(i[c[n]]=r?U(o):t[o],o+=r?2:1,!i[c[n]])throw a("a zero quantisation step");if(s?.progressive&&s.components.some(t=>t.tq===e&&t.quant&&t.coverage.some(t=>0!==t)&&t.quant.some((t,o)=>t!==i[o])))throw a("a quantisation table changed between progressive scans");n[e]=i}else if(196===h)for(let o=w;o<d;){const n=t[o]>>4,i=15&t[o];if(n>1||i>3||o+17>d)throw a("bad Huffman table");o++;const s=t.subarray(o,o+16);let f=0;for(const t of s)f+=t;if(f>256||o+16+f>d)throw a("bad Huffman table");const c=t.subarray(o+16,o+16+f);r[n][i]=l(s,c),0===n&&(c.some(t=>t>15)?e.add(i):e.delete(i)),o+=16+f}else if(192===h||193===h||194===h){if(s)throw a("two frames");if(d-w<6)throw a("bad frame header");if(8!==t[w])throw a("only 8-bit samples");const n=U(w+1),r=U(w+3),e=t[w+5];if(!r||!n)throw a("a height left to a DNL marker");if(1!==e&&3!==e)throw a("only grey or three components");if(d-w<6+3*e)throw a("bad frame header");if(r>o.edge||n>o.edge)throw a("a picture past "+o.edge+" pixels on a side","JXL_DIMENSIONS");if(r*n>o.pixels)throw a("a picture past "+o.pixels+" pixels","JXL_DIMENSIONS");const i=[];let f=0;for(let o=0,n=w+6;o<e;o++,n+=3){const o=t[n+1]>>4,r=15&t[n+1],e=t[n+2];if(o<1||o>4||r<1||r>4)throw a("bad sampling");if(e>3)throw a("bad quantisation table");if(i.some(o=>o.id===t[n]))throw a("two components with one id");f+=o*r,i.push({id:t[n],h:o,v:r,tq:e,pred:0,coverage:new Int8Array(64).fill(-1)})}if(f>10)throw a("bad sampling");const c=Math.max(...i.map(t=>t.h)),l=Math.max(...i.map(t=>t.v)),u=Math.ceil(r/(8*c)),p=Math.ceil(n/(8*l));for(const t of i)t.blocksW=Math.ceil(Math.ceil(r*t.h/c)/8),t.blocksH=Math.ceil(Math.ceil(n*t.v/l)/8),t.stride=u*t.h,t.rows=p*t.v,t.coeffs=new Int16Array(t.stride*t.rows*64);s={progressive:194===h,width:r,height:n,components:i,hmax:c,vmax:l,mcusX:u,mcusY:p,scans:0}}else{if(195===h||h>=197&&h<=199||h>=201&&h<=203||h>=205&&h<=207)throw a("lossless, hierarchical or arithmetic coding");if(221===h){if(u<4)throw a("bad segment");f=U(w)}else{if(220===h)throw a("a DNL marker");if(218===h){if(!s)throw a("a scan before the frame");const o=t[w],i=[];if(o<1||o>4||w+1+2*o+3>d)throw a("bad scan header");let c=w+1;for(let e=0;e<o;e++,c+=2){const o=s.components.find(o=>o.id===t[c]);if(!o)throw a("a scan of an unknown component");if(i.some(t=>t.component===o))throw a("a component twice in one scan");if(o.quant||(o.quant=n[o.tq]),!o.quant)throw a("a scan without its quantisation table");i.push({component:o,dc:r[0][t[c+1]>>4],ac:r[1][15&t[c+1]],dcSlot:t[c+1]>>4})}const l=t[c],h=t[c+1],u=t[c+2]>>4,p=15&t[c+2];if((!s.progressive||0===l&&0===u)&&i.some(t=>e.has(t.dcSlot)))throw a("bad Huffman table");A=b(t,d,s,i,l,h,u,p,f),s.scans++}else if(224===h)x(w,"JFIF")&&(M=!0);else if(225===h)1===g&&x(w,"Exif")&&(g=y(t,w+6,d));else if(226===h){if(u>=16&&x(w,"ICC_PROFILE\0")){const o=t[w+12],n=t[w+13];if(!o||!n||o>n||i.count&&i.count!==n||i.p[o])throw a("a colour profile cut short");i.count=n,i.p[o]=t.subarray(w+14,d)}}else 238===h&&u>=14&&x(w,"Adobe")&&(m=t[w+11])}}}if(!s||!s.scans)throw a("no picture");for(const t of s.components)for(let o=0;o<64;o++)if(0!==t.coverage[o])throw a(s.progressive?"a progressive JPEG whose scans do not finish every coefficient":"a component without its scan");let I="srgb";if(i.count){const t=[];for(let o=1;o<=i.count;o++){if(!i.p[o])throw a("a colour profile cut short");t.push(i.p[o])}const o=new Uint8Array(t.reduce((t,o)=>t+o.length,0));for(let n=0,r=0;n<t.length;r+=t[n].length,n++)o.set(t[n],r);if(I=function(t){const o=t.length,n=o=>t[o]<<8|t[o+1],r=n=>n+4<=o?(t[n]<<24|t[n+1]<<16|t[n+2]<<8|t[n+3])>>>0:0,e=n=>n+4<=o?String.fromCharCode(t[n],t[n+1],t[n+2],t[n+3]):"",i=t=>(0|r(t))/65536,s=r(128);if(o<132||132+12*s>o||"acsp"!==e(36)||"XYZ "!==e(20))return null;const f=t=>{for(let n=0;n<s;n++){const i=132+12*n,s=r(i+4);if(e(i)===t)return s+r(i+8)<=o?s:-1}return-1},c=t=>{const s=f(t);if(s<0)return!1;if("para"===e(s)){const t=n(s+8),o=[1,3,4,5,7][t];if(3!==t&&4!==t)return!1;for(let t=0;t<o;t++)if(Math.abs(i(s+12+4*t)-(t<5?d[t]:0))>(t?.002:.01))return!1;return!0}const c="curv"===e(s)?r(s+8):0;if(c<2||s+12+2*c>o)return!1;for(let t=0;t<c;t++)if(Math.abs(n(s+12+2*t)/65535-p(t/(c-1)))>.001)return!1;return w.every((t,o)=>{const r=o*(c-1)/16,e=Math.floor(r),i=Math.min(c-1,e+1),f=r-e;return Math.abs((n(s+12+2*e)*(1-f)+n(s+12+2*i)*f)/65535-t)<=.001})};if(["A2B0","A2B1","A2B2"].some(t=>f(t)>=0))return null;if("GRAY"===e(16))return c("kTRC")?"srgb":null;if("RGB "!==e(16)||!["rTRC","gTRC","bTRC"].every(c))return null;const a=["rXYZ","gXYZ","bXYZ"].flatMap(t=>{const o=f(t);return o>=0&&"XYZ "===e(o)?[i(o+8),i(o+12),i(o+16)]:[9,9,9]}),l=t=>t.every((t,o)=>Math.abs(a[o]-t)<=.004);return l(h)?"srgb":l(u)?"display-p3":null}(o),!I)throw a("a colour profile other than sRGB or Display P3 (decode the picture and encode its pixels)")}let E=!1;return 3!==s.components.length||M||(E=m>=0?0===m:82===s.components[0].id&&71===s.components[1].id&&66===s.components[2].id),{width:s.width,height:s.height,components:s.components,progressive:s.progressive,ycbcr:!E,orientation:g,colorSpace:I}}(n);!function(o,n,r=t){if(!Number.isInteger(o)||!Number.isInteger(n)||o<1||n<1)throw m("JXL_INPUT","A picture is at least one pixel wide and high, in whole pixels.");if(o>r.edge||n>r.edge)throw m("JXL_DIMENSIONS","A picture is at most "+r.edge+" pixels on a side.");if(o*n>r.pixels)throw m("JXL_DIMENSIONS","A picture is at most "+r.pixels.toLocaleString("en-US")+" pixels.")}(f.width,f.height,o);const M=f.orientation>=5;return(t=>{Object.assign(s,t)})({width:M?f.height:f.width,height:M?f.width:f.height,orientation:f.orientation}),yield*i(f,r)}(n,i,0,e));return s}const U=256;function x(t,o){t.o([[9,1],[13,1],[18,1],[30,1]],o)}const I=[1,2,3,4,0,5,17,6,16,7,8,9,10,11,12,13,14,15],E=[2,4,3,2,2,4],v=[0,7,3,2,1,15];function C(t,o){const n=new Uint8Array(t.length),r=[];for(let o=0;o<t.length;o++)t[o]>0&&r.push(o);if(!r.length)return n;if(1===r.length)return n[r[0]]=1,n;for(let e=1;;e*=2){const i=r.map(o=>({w:Math.max(t[o],e),m:o,l:null,r:null})).sort((t,o)=>t.w-o.w||t.m-o.m),s=[];let f=0,c=0;const a=()=>f<i.length&&(c>=s.length||i[f].w<=s[c].w)?i[f++]:s[c++];for(;i.length-f+(s.length-c)>1;){const t=a(),o=a();s.push({w:t.w+o.w,m:-1,l:t,r:o})}const l=[[a(),0]];let h=0;for(;l.length;){const[t,o]=l.pop();t.m>=0?(n[t.m]=o,o>h&&(h=o)):l.push([t.l,o+1],[t.r,o+1])}if(h<=o)return n}}function S(t,o,n,r,e,i){const s=[];for(;s.push(r&e),r>>=i;)r--;for(let r=s.length-1;r>=0;r--)t.push(n),o.push(s[r])}function J(t,o){let n=0;for(let r=0;r<o;r++)n=n<<1|1&t,t>>=1;return n}function L(t){const o=new Uint32Array(17),n=new Uint32Array(17),r=new Uint16Array(t.length);for(const n of t)o[n]++;o[0]=0;for(let t=1,r=0;t<=16;t++)r=r+o[t-1]<<1,n[t]=r;for(let o=0;o<t.length;o++)t[o]&&(r[o]=J(n[t[o]]++,t[o]));return r}function X(t){const o=[];for(let n=0;n<t.length;n++)t[n]>0&&o.push(n);const n=o.length,r=o[n-1]+1||1,e=n<=4?new Uint8Array(r):C(t,15).subarray(0,r);let i=0;if(n<=4){o.sort((o,n)=>t[n]-t[o]||o-n),i=4===n&&t[o[0]]>t[o[2]]+t[o[3]]?1:0;const r=n<2?[0]:2===n?[1,1]:3===n?[1,2,2]:i?[1,2,3,3]:[2,2,2,2];o.forEach((t,o)=>{e[t]=r[o]}),o.sort((t,o)=>e[t]-e[o]||t-o)}return{M:r,A:e,U:L(e),I:n&&n<=4?o:null,C:i}}function k(t,o){if(o.M<=1)return;if(o.I){const n=i(o.M-1)+1;t.write(2,1),t.write(2,o.I.length-1);for(const r of o.I)t.write(n,r);return void(4===o.I.length&&t.write(1,o.C))}const n=o.A,r=[],e=[];let s=n.length-1;for(;s>0&&!n[s];)s--;let f=8;for(let t=0;t<=s;){const o=n[t];let i=1;for(;t+i<=s&&n[t+i]===o;)i++;let c=i;if(o&&o!==f&&(r.push(o),e.push(0),c--,f=o),c<3)for(let t=0;t<c;t++)r.push(o),e.push(0);else S(r,e,o?16:17,c-3,o?3:7,o?2:3);t+=i}const c=new Uint32Array(18);for(const t of r)c[t]++;let a=0;for(const t of c)t&&a++;1===a&&(c[r[0]?0:1]=c[r[0]]);const l=C(c,5),h=L(l);let u=18;for(;u>0&&!l[I[u-1]];)u--;t.write(2,0);for(let o=0;o<u;o++){const n=l[I[o]];t.write(E[n],v[n])}for(let o=0;o<r.length;o++){const n=r[o];t.write(l[n],h[n]),16===n?t.write(2,e[o]):17===n&&t.write(3,e[o])}}function N(t,o=0,n=0){return{split:t,S:o,J:n,L:1<<t}}function O(t,o,n=15){t.write(s(n+1),o.split),o.split!==n&&(t.write(s(o.split+1),o.S),t.write(s(o.split-o.S+1),o.J))}function T(t,o,n){if(o<t.L)return n[0]=o,n[1]=0,void(n[2]=0);const r=i(o),e=o-(1<<r),s=r-t.S-t.J;n[0]=t.L+(r-t.split<<t.S+t.J|e>>r-t.S<<t.J|e&(1<<t.J)-1),n[1]=s,n[2]=o>>t.J&(1<<s)-1}function _(t,o,n,r=0){T(t,o,D),n[r+D[0]]++}const D=[0,0,0];function P(t,o,n,r,e=0){T(n,r,D);const i=e+D[0];t.write(o.A[i],o.U[i]),D[1]&&t.write(D[1],D[2])}function G(t,o){let n=0;for(const t of o)n=Math.max(n,t);const r=n?s(n+1):0;if(r<=3){if(t.write(1,1),t.write(2,r),r)for(const n of o)t.write(r,n);return}t.write(1,0),t.write(1,0);const e=N(8),i=N(4),f=o.length>=64,c=new Uint32Array(f?257:n+1),a=[];for(let t=0;t<o.length;){let n=1;for(;f&&t+n<o.length&&o[t+n]===o[t];)n++;if(n>=8)a.push([o[t],n-1]),c[o[t]]++,_(i,n-1-7,c,224),t+=n;else{for(let r=0;r<n;r++)a.push([o[t],0]),c[o[t]]++;t+=n}}const l=X(c);if(!f){R(t,{X:new Uint8Array(1),k:[{config:e,code:l}]});for(const[o]of a)P(t,l,e,o);return}const h=new Uint32Array(1);h[0]=1,R(t,{N:{O:224,minLength:7,T:i},X:new Uint8Array([1,0]),k:[{config:N(0),code:X(h)},{config:e,code:l}]});for(const[o,n]of a)P(t,l,e,o),n&&P(t,l,i,n-7,224)}function R(t,{N:o=null,X:n,k:r}){t.write(1,o?1:0),o&&(t.o([[0,224],[0,512],[0,4096],[15,8]],o.O),t.o([[0,3],[0,4],[2,5],[8,9]],o.minLength),O(t,o.T,8)),n.length>1&&G(t,n),t.write(1,1);for(const o of r)O(t,o.config,15);for(const o of r){const n=o.code.M-1;if(n){const o=i(n);t.write(1,1),t.write(4,o),t.write(o,n-(1<<o))}else t.write(1,0)}for(const o of r)k(t,o.code)}const q=N(4),B=Object.freeze({O:224,minLength:7,T:q}),Z=N(0);function j(t,o){const n=[[3,0],[6,8],[10,72],[13,1096]];if("rct"===o.type)t.write(2,0),t.o(n,o._),t.o([[0,6],[2,0],[4,2],[6,10]],o.D??6);else if("palette"===o.type)t.write(2,1),t.o(n,o._),t.o([[0,1],[0,3],[0,4],[13,1]],o.P),t.o([[8,0],[10,256],[12,1280],[16,5376]],o.G),t.o([[0,0],[8,1],[10,257],[16,1281]],o.R??0),t.write(4,o.q??0);else{if("squeeze"!==o.type)throw new Error("unknown transform "+o.type);{t.write(2,2);const r=o.B??[];t.o([[0,0],[4,1],[6,9],[8,41]],r.length);for(const o of r)t.write(1,o.Z?1:0),t.write(1,o.j?1:0),t.o(n,o._),t.o([[0,1],[0,2],[0,3],[4,4]],o.P)}}}function z(t,o,n,r,i,s){const f=s.q,c=s.offset,a=s.Y,l=Z,h=q;let u=0;const w=(n,r=l,e=0)=>{t?P(t,o,r,n,e):_(r,n,o,e)},d=()=>{if(u){if(u>=8)w(0),w(u-7-1,h,224);else for(let t=0;t<u;t++)w(0);u=0}},p=(t,o)=>{const r=n[t];let i;1===a?i=r-o-c:(i=Math.round((r-o-c)/a),n[t]=o+c+i*a),0!==i?(d(),w(e(i))):u++};for(let t=0,o=0;t<i;t++)for(let e=0;e<r;e++,o++){let i;if(0===f)i=0;else{const s=e?n[o-1]:t?n[o-r]:0,c=t?n[o-r]:s;if(1===f)i=s;else if(2===f)i=c;else if(3===f)i=(s+c)/2|0;else{const a=s+c-(e&&t?n[o-r-1]:s);if(5===f){const t=s<c?s:c,o=s<c?c:s;i=a<t?t:a>o?o:a}else{if(4!==f)throw new Error("predictor "+f+" is not coded here");i=Math.abs(a-s)<Math.abs(a-c)?s:c}}}p(o,i)}d()}const Y=[N(0),N(1),N(2),N(3),N(4),N(4,1,1),N(5,1,1),N(2,0,1),N(3,0,1)],H=N(4,1,1);class F{constructor(t){this.H=t,this.small=new Uint32Array(64*t),this.F=Array.from({length:t},()=>null),this.K=new Float64Array(t)}add(t,o){if(this.K[t]++,o<64)return void this.small[64*t+o]++;const n=this.F[t]||(this.F[t]=new Map);n.set(o,(n.get(o)||0)+1)}V(t,o,n=64){const r=new Uint32Array(n);for(const n of t){const t=64*n;for(let n=0;n<64;n++)this.small[t+n]&&(r[Q(o,n)]+=this.small[t+n]);const e=this.F[n];if(e)for(const[t,n]of e)r[Q(o,t)]+=n}return r}W(t,o){let n=$(this.V(t,o));for(const r of t){const t=64*r;for(let r=o.L;r<64;r++)this.small[t+r]&&(n+=this.small[t+r]*(i(r)-o.S-o.J));const e=this.F[r];if(e)for(const[t,r]of e)n+=r*(i(t)-o.S-o.J)}return n}}const K=[0,0,0];function Q(t,o){return T(t,o,K),K[0]}const V=[1];for(let t=1;t<32;t++)V.push(V[t-1]/2);function W(t){let o=0;for(t>=2&&t<4294967296&&(o=31-Math.clz32(t),t*=V[o]);t>=1.4142135623730951;)t/=2,o++;for(;t<.7071067811865476;)t*=2,o--;const n=(t-1)/(t+1),r=n*n;let e=n,i=n;for(let t=3;t<30;t+=2)i+=(e*=r)/t;return o+2.885390081777927*i}function $(t){let o=0,n=0;for(const n of t)o+=n;for(const r of t)r&&(n+=r*W(o/r));return n}function tt(t,{maxClusters:o=48,newClusterCost:n=320}={}){const r=[];for(let o=0;o<t.H;o++)t.K[o]>0&&r.push(o);r.sort((o,n)=>t.K[n]-t.K[o]);const e=[];for(const i of r){const r=t.V([i],H),s=$(r);let f=-1,c=1/0;for(let t=0;t<e.length;t++){const o=e[t],n=new Uint32Array(64);for(let t=0;t<64;t++)n[t]=o.$[t]+r[t];const i=$(n)-o.tt-s;i<c&&(c=i,f=t)}if(f<0||e.length<o&&c>n)e.push({H:[i],$:r,tt:s});else{const t=e[f];t.H.push(i);for(let o=0;o<64;o++)t.$[o]+=r[o];t.tt=$(t.$)}}e.length||e.push({H:[],$:new Uint32Array(64),tt:0});const i=new Uint8Array(t.H);let s=2*t.H;const f=e.map((o,n)=>{for(const t of o.H)i[t]=n;let r=Y[0],e=1/0;for(const n of Y){const i=t.W(o.H,n);i<e&&(e=i,r=n)}const f=X(t.V(o.H,r,128));return s+=e+3*f.M+40,{config:r,code:f}});return{X:i,k:f,write:(t,o,n)=>{const r=f[i[o]];P(t,r.code,r.config,n)},tt:s}}const ot=[0,0,1,2,3,4,5,6,7,8,9,10,11,12,13,14,15,15,16,16,17,17,18,18,19,19,20,20,21,21,22,22,23,23,23,23,24,24,24,24,25,25,25,25,26,26,26,26,27,27,27,27,28,28,28,28,29,29,29,29,30,30,30,30],nt=[0,0,31,62,62,93,93,93,93,123,123,123,123,152,152,152,152,152,152,152,152,180,180,180,180,180,180,180,180,180,180,180,180,206,206,206,206,206,206,206,206,206,206,206,206,206,206,206,206,206,206,206,206,206,206,206,206,206,206,206,206,206,206,206],rt=c.map(t=>(7&t)<<3|t>>3),et=(t,o,n)=>2*(nt[t]+ot[o])+n;function it(t,o,n){const r=[];if(o>0){const n=new Uint32Array(2048);let e=0;for(const o of t)n[Math.min(2047,Math.max(0,o+1024))]++,e++;let i=0,s=e/(o+1);for(let t=0;t<2048&&r.length<o;t++)i+=n[t],i>s&&(r.push(t-1025),s=e*(r.length+1)/(o+1))}const e=r.length+1,i=new Uint8Array(39*e);for(let t=0;t<e;t++)for(let o=0;o<13;o++)i[(0+o)*e+t]=t,i[(13+o)*e+t]=n?e:e+(t>>1),i[(26+o)*e+t]=n?e:e+(e-1>>1)+1+(t>>1);let s=0;for(const t of i)s=Math.max(s,t+1);return{thresholds:r,ot:e,map:i,nt:s,rt:t=>{let o=0;for(const n of r)t>n&&o++;return o},et:(t,o)=>i[13*(t<2?1^t:2)*e+o],it:(t,o)=>(Math.min(t,64)<8?Math.min(t,64):4+(Math.min(t,64)>>1))*s+o,st:t=>37*s+458*t,H:495*s}}function st(t,o,n){const r=o.map(()=>function(t,o=0,n=1){return{q:t,offset:o,Y:n,context:-1}}(n)),i=r.map(()=>new Uint32Array(257));o.forEach((t,o)=>z(null,i[o],t.data,t.w,t.h,r[o])),function(t,{ft:o=!0,ct:n=[]}={}){t.write(1,o?1:0),t.write(1,1),t.o([[0,0],[0,1],[4,2],[8,18]],n.length);for(const o of n)j(t,o)}(t,{ft:!1,ct:[]});const s=function(t,o){const n=[o],r=[],i=[];for(const t of n){if(t.left){r.push(t.lt+1,e(t.ht)),n.push(t.left,t.right);continue}t.context=i.length,i.push(t);const o=1===t.Y?0:31-Math.clz32(t.Y&-t.Y);r.push(0,t.q,e(t.offset),o,(t.Y>>o)-1)}const s=new Uint32Array(64);for(const t of r)_(Z,t,s);const f=X(s);R(t,{X:new Uint8Array(6),k:[{config:Z,code:f}]});for(const o of r)P(t,f,Z,o);return i}(t,function(t){return function(t,o){const n=(r,e)=>{if(r===e)return o?t[r].ut:t[r];const i=r+e>>1;return function(t,o,n,r){return{lt:t,ht:o,left:n,right:r}}(o,o?t[i].wt:i,n(i+1,e),n(r,i))};return n(0,t.length-1)}(t,0)}(r)),f=function(t,o,n,r=t=>t.context){const e=new Uint8Array(o.length+1);for(const t of o)e[t.context]=r(t)+1;const i=new Uint32Array(2);i[1]=1;const s=[{config:Z,code:X(i)}],f=Math.max(...e);for(let t=0;t<f;t++)s.push({config:Z,code:X(n[t])});return R(t,{N:B,X:e,k:s}),s}(t,s,i,t=>r.indexOf(t));o.forEach((o,n)=>z(t,f[n+1].code,o.data,o.w,o.h,r[n]))}function*ft(t,o={}){const{width:n,height:i,components:c,alpha:l}=t,h=1===c.length,u=t.quantScale||1;let w=t.quantFieldBase||1;for(;w/(2040*u)<1/16384;)w*=2;const d=t.ycbcr||h,p=h?[0,0,0]:d?[1,0,2]:[0,1,2],y=p.map(t=>c[t]),b=o.orders?p.map(t=>o.orders[t].map(t=>rt[t])):[rt,rt,rt],m=t=>1===t?0:2===t?1:-1,g=y.map(t=>m(t.h)),A=y.map(t=>m(t.v));if(g.includes(-1)||A.includes(-1))throw a("sampling factors other than one or two");const I=Math.max(...g),E=Math.max(...A),v=t=>I-g[t],C=t=>E-A[t];if(!d&&(I||E))throw a("an RGB JPEG with subsampling");const S=Math.ceil(n/(8<<I))<<I,J=Math.ceil(i/(8<<E))<<E,L=Math.ceil(n/U),k=Math.ceil(i/U),O=L*k,T=Math.ceil(S/256),D=T*Math.ceil(J/256),q=1===O,B=(t,o,n)=>{const r=y[t];return h&&1!==t?-1:64*(n*r.stride+o)},Z=t=>{let o=1024;for(;t%o;)o>>=1;return o},j=y.map(t=>d?1:t.quant[0]/Z(t.quant[0])),z=y.map(t=>d?0:1024/Z(t.quant[0])),Y=new Int32Array(S*J);for(let t=0;t<J;t++)for(let o=0;o<S;o++)Y[t*S+o]=y[1].coeffs[B(1,o,t)]*j[1]+z[1];const H=y[1].quant,K=(H[8]+H[16]+H[24]+H[32]+H[40])/u,Q=S*J,V=Q<256?0:Math.max(1,Math.min(7,s(Q)-s(K)-7));let W=it(Y,V,h);const $=new r(256);if(function(t,o,n,r,e,{xyb:i=!1,orientation:s=1,colorSpace:f}={}){t.write(16,2815),t.write(1,0),x(t,n),t.write(3,0),x(t,o),t.write(1,0);const c=1!==s;t.write(1,c?1:0),c&&t.write(6,s-1),t.write(4,8),e?t.write(3,5):t.write(2,0),t.write(1,i?1:0),3===r&&"display-p3"===f?t.write(21,903568):3===r?t.write(1,1):t.write(15,14100),c&&t.write(1,1),t.write(2,0),t.write(1,1),t.i()}($,n,i,h?1:3,!!l,{orientation:t.orientation,colorSpace:t.colorSpace}),$.write(1,0),$.write(2,0),$.write(1,0),$.write(2,2),$.write(8,111),$.write(1,d?1:0),d)for(let t=0;t<3;t++)$.write(2,0===g[t]&&0===A[t]?0:1===g[t]&&1===A[t]?1:1===g[t]?2:3);$.write(2,0),l&&$.write(2,0),function(t,o){t.write(o?7:5,0),t.write(11,1)}($,l);const ot=[0,1,2].map(()=>new Int32Array(1024)),nt=(t,o)=>{const n=t%L*32,r=32*(t/L|0),i=Math.min(32,S-n),s=Math.min(32,J-r);for(let t=0;t<s;t++)for(let s=0;s<i;s++)for(const i of[1,0,2]){const f=v(i),c=C(i),a=s>>f,l=t>>c;if(a<<f!==s||l<<c!==t)continue;const h=B(i,(n>>f)+a,(r>>c)+l),u=y[i].coeffs;let w=0;if(h>=0)for(let t=1;t<64;t++)u[h+t]&&w++;const d=ot[i],p=0===a?0===l?32:d[32*(l-1)]:0===l?d[32*l+a-1]:d[32*(l-1)+a]+d[32*l+a-1]+1>>1;d[32*l+a]=w;const m=W.et(i,W.rt(Y[(r+t)*S+n+s]));o(W.it(p,m),w);const M=W.st(m);let g=w>4?0:1,A=w;for(let t=1;t<64&&A;t++){const n=u[h+b[i][t]];o(M+et(A,t,g),e(n)),g=n?1:0,A-=g}}},ft=(V>0?3:2)*O+D;let ct=0;const at=function*(t){W=t;const n=new F(t.H);for(let t=0;t<O;t++)nt(t,(t,o)=>n.add(t,o)),yield++ct/ft;return(o.coding||tt)(n,o.clusters)};let lt=yield*at(W);if(V>0){const t=W,o=it(Y,0,h),n=yield*at(o);n.tt<=lt.tt?lt=n:W=t}const ht=Array.from({length:q?1:2+D+O},()=>null),ut=t=>ht[q?0:t]||(ht[q?0:t]=new r(4096)),wt=new Uint32Array(ht.length);let dt=0;const pt=t=>{t=q?0:t;const o=Math.ceil(ht[t].u/8);dt+=o-wt[t],wt[t]=o,M(dt)},yt=ut(0);yt.write(1,0);const bt=(t,o)=>{const r=t*U,e=o*U,s=Math.min(U,n-r),f=Math.min(U,i-e),c=new Int32Array(s*f);for(let t=0;t<f;t++)for(let o=0;o<s;o++)c[t*s+o]=l[4*((e+t)*n+r+o)+3];return{w:s,h:f,data:c}};for(let t=0;t<3;t++)yt.write(16,f(y[t].quant[0]/j[t]/(2040*u)*128));yt.write(2,3),yt.write(16,57343),yt.write(2,1),yt.write(5,0),yt.write(1,0),yt.write(4,0),yt.write(4,W.thresholds.length);for(const t of W.thresholds)yt.o([[4,0],[8,16],[16,272],[32,65808]],e(t));yt.write(4,0),yt.write(4,0),G(yt,W.map),yt.write(1,0),yt.write(2,0),yt.write(16,0),yt.write(16,0),yt.write(8,128),yt.write(8,128),yt.write(1,0),l&&st(yt,[q?bt(0,0):{w:0,h:0,data:new Int32Array(0)}],5),pt(0);for(let o=0;o<D;o++){const n=ut(1+o),r=o%T*256,e=256*(o/T|0),i=Math.min(256,S-r),f=Math.min(256,J-e);n.write(2,0),st(n,[1,0,2].map(t=>{const o=i>>v(t),n=f>>C(t),s=new Int32Array(o*n),c=j[t],a=z[t];for(let i=0;i<n;i++)for(let n=0;n<o;n++){const f=B(t,(r>>v(t))+n,(e>>C(t))+i);s[i*o+n]=f<0?0:y[t].coeffs[f]*c+a}return{w:o,h:n,data:s}}),5);const c=i*f,a=s(c);a&&n.write(a,c-1);const l=i+7>>3,h=f+7>>3,u=new Int32Array(2*c).fill(w-1,c);if(t.quantFields)for(let o=0;o<f;o++)for(let n=0;n<i;n++)u[c+o*i+n]=t.quantFields[(e+o)*S+r+n]-1;st(n,[{w:l,h:h,data:new Int32Array(l*h)},{w:l,h:h,data:new Int32Array(l*h)},{w:c,h:2,data:u},{w:i,h:f,data:new Int32Array(c)}],w>1?5:0),pt(1+o),yield++ct/ft}const mt=ut(1+D);mt.write(1,0),mt.write(3,7),mt.write(16,f(w/(2040*u))),st(mt,[0,1,2].map(t=>{const o=new Int32Array(64),n=y[t].quant;for(let t=0;t<8;t++)for(let r=0;r<8;r++)o[8*r+t]=n[8*t+r];return{w:8,h:8,data:o}}),5);for(let t=1;t<17;t++)mt.write(3,0);if(O>1&&mt.write(s(O),0),o.orders){mt.o([[0,95],[0,19],[0,0],[13,0]],1);const t=[];for(const n of p){const r=Array.from({length:64},(t,o)=>o),e=o.orders[n].map(t=>{const o=r.indexOf(t);return r.splice(o,1),o});let i=64;for(;i>1&&!e[i-1];)i--;t.push(i-1,...e.slice(1,i))}const n=N(0),r=new Uint32Array(64);for(const o of t)_(n,o,r);const e=X(r);R(mt,{X:new Uint8Array(8),k:[{config:n,code:e}]});for(const o of t)P(mt,e,n,o)}else mt.write(2,2);lt.writeHistograms?lt.writeHistograms(mt):R(mt,{X:lt.X,k:lt.k}),pt(1+D);for(let t=0;t<O;t++){const o=ut(2+D+t);nt(t,(t,n)=>lt.write(o,t,n)),lt.flush&<.flush(o),l&&!q&&st(o,[bt(t%L,Math.floor(t/L))],5),pt(2+D+t),yield++ct/ft}return function(t,o){const n=o.map(t=>t.length);!function(t,o){t.write(1,0),t.i();for(const n of o)t.o([[10,0],[14,1024],[22,17408],[30,4211712]],n);t.i()}(t,n);const r=t.finish(),e=r.length+n.reduce((t,o)=>t+o,0);M(e);const i=new Uint8Array(e);i.set(r);let s=r.length;for(const t of o)i.set(t,s),s+=t.length;return i}($,function(t){return t.map(t=>t?t.finish():new Uint8Array(0))}(ht))}function*ct(t,o,n){if(o<3&&!n)return yield*ft(t);let r=n;try{const n=ft(t);let e,i;for(;!(e=n.next()).done;)if(i=yield e.value/2,i&&r&&e.value<1)return n.return(),r;if((!r||e.value.length<r.length)&&(r=e.value),i||o<3)return r;const s=o<4?1:3,f=ft(t,{clusters:{maxClusters:32,newClusterCost:160}});for(;!(e=f.next()).done;)if(i=yield.5+e.value/(2*s),i&&e.value<1)return f.return(),r;if(e.value.length<r.length&&(r=e.value),i||o<4)return r;const a=function*(t){const o=c.map(t=>(7&t)<<3|t>>3),n=t.map(()=>new Uint32Array(64)),r=t.reduce((t,o)=>t+o.coeffs.length,0);let e=0;for(let i=0;i<t.length;i++){const{coeffs:s}=t[i],f=n[i];for(let t=0;t<s.length;t+=65536){const n=Math.min(t+65536,s.length);for(let r=t;r<n;r+=64)for(let t=1;t<64;t++)s[r+o[t]]&&f[t]++;if(e+=n-t,yield e/r)return null}}return n.map(t=>[0,...Array.from({length:63},(t,o)=>o+1).sort((o,n)=>Math.floor(t[n]/8)-Math.floor(t[o]/8)||o-n)])}(t.components);for(;!(e=a.next()).done;)if(yield.5+(1+e.value)/6)return a.return(),r;const l=ft(t,{orders:e.value});for(;!(e=l.next()).done;)if(i=yield.5+(2+e.value)/6,i&&e.value<1)return l.return(),r;return e.value.length<r.length?e.value:r}catch(t){if(!r||(!(t instanceof RangeError)||t.code)&&"JXL_SIZE"!==t.code)throw t;return r}}function at(t,o){return function(t){return function(t){let o;for(;!(o=t.next()).done;);o.value}(t),{bytes:t.bytes,width:t.width,height:t.height,orientation:t.orientation}}(lt(t,o))}function lt(t,o){return A(t,o,ct)}export{o as LIMITS,at as transcode,lt as transcodeSteps};
|
package/jpeg.mjs
CHANGED
|
@@ -1,39 +1,15 @@
|
|
|
1
1
|
// SPDX-License-Identifier: MIT
|
|
2
|
-
// Rapier JXL's JPEG door: a JPEG carried as its
|
|
3
|
-
// quantisation tables, colour and subsampling
|
|
4
|
-
//
|
|
5
|
-
//
|
|
6
|
-
// Not carried:
|
|
7
|
-
//
|
|
8
|
-
//
|
|
9
|
-
import {
|
|
10
|
-
import {
|
|
11
|
-
import {
|
|
12
|
-
import {fault, admitSize, job} from './admit.mjs';
|
|
2
|
+
// Rapier JXL's JPEG door: a JPEG carried as its quantised DCT coefficients,
|
|
3
|
+
// quantisation tables, colour and subsampling, with only the entropy coding changed.
|
|
4
|
+
// Baseline, extended sequential and progressive scans, restarts; 8-bit grey, YCbCr or RGB;
|
|
5
|
+
// Exif orientation and an admitted sRGB or Display P3 declaration are kept.
|
|
6
|
+
// Not carried: JPEG reconstruction data, ICC bytes, Exif beyond orientation or XMP.
|
|
7
|
+
// Arithmetic, 12-bit, lossless, CMYK, DNL, another colour profile, truncation or a table
|
|
8
|
+
// libjpeg refuses is JXL_JPEG: decode it and encode its pixels instead.
|
|
9
|
+
import {JPEG_LIMITS} from './bits.mjs';
|
|
10
|
+
import {jpegAnswer, jpegJob} from './jpeg-job.mjs';
|
|
11
|
+
import {coefficientEffortSteps} from './coefficient-effort.mjs';
|
|
13
12
|
|
|
14
|
-
// This door's limits (bits.mjs): 64 million pixels, its memory 6.4 bytes a pixel at 4:4:4 and 3.3 at 4:2:0 besides the
|
|
15
|
-
// JPEG at that size, so a phone's 24 and 48 megapixel JPEGs are carried.
|
|
16
13
|
export {JPEG_LIMITS as LIMITS};
|
|
17
|
-
|
|
18
|
-
export function
|
|
19
|
-
const it = transcodeSteps(jpeg);
|
|
20
|
-
complete(it);
|
|
21
|
-
return {bytes: it.bytes, width: it.width, height: it.height, orientation: it.orientation};
|
|
22
|
-
}
|
|
23
|
-
|
|
24
|
-
// The same work as a job (admit.mjs): the scans read in the first step, which also sets the job's width, height and
|
|
25
|
-
// orientation (the picture as shown), then the VarDCT writer's steps.
|
|
26
|
-
export function transcodeSteps(jpeg) {
|
|
27
|
-
if (!(jpeg instanceof Uint8Array) || !jpeg.length) throw fault('JXL_INPUT', 'A JPEG is a non-empty Uint8Array.');
|
|
28
|
-
if (jpeg.length > JPEG_LIMITS.bytes) throw fault('JXL_SIZE', 'A JPEG is at most 16 MiB.');
|
|
29
|
-
const it = job(carry(jpeg, shown => Object.assign(it, shown)));
|
|
30
|
-
return it;
|
|
31
|
-
}
|
|
32
|
-
|
|
33
|
-
function* carry(jpeg, shown) {
|
|
34
|
-
const parsed = parseJPEG(jpeg);
|
|
35
|
-
admitSize(parsed.width, parsed.height, JPEG_LIMITS);
|
|
36
|
-
const swapped = parsed.orientation >= 5;
|
|
37
|
-
shown({width: swapped ? parsed.height : parsed.width, height: swapped ? parsed.width : parsed.height, orientation: parsed.orientation});
|
|
38
|
-
return yield* varDCTSteps(parsed);
|
|
39
|
-
}
|
|
14
|
+
export function transcode(jpeg, options) { return jpegAnswer(transcodeSteps(jpeg, options)); }
|
|
15
|
+
export function transcodeSteps(jpeg, options) { return jpegJob(jpeg, options, coefficientEffortSteps); }
|