@johnhenry/math 0.1.0 → 0.3.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/README.md +5 -4
- package/dist/Distributions.d.ts +86 -16
- package/dist/Distributions.d.ts.map +1 -1
- package/dist/Distributions.js +340 -89
- package/dist/Distributions.js.map +1 -1
- package/dist/EigenGeneral.d.ts +22 -0
- package/dist/EigenGeneral.d.ts.map +1 -0
- package/dist/EigenGeneral.js +327 -0
- package/dist/EigenGeneral.js.map +1 -0
- package/dist/Frobenius.d.ts +60 -0
- package/dist/Frobenius.d.ts.map +1 -0
- package/dist/Frobenius.js +162 -0
- package/dist/Frobenius.js.map +1 -0
- package/dist/Graph.d.ts +26 -0
- package/dist/Graph.d.ts.map +1 -1
- package/dist/Graph.js +75 -0
- package/dist/Graph.js.map +1 -1
- package/dist/MatrixMath.d.ts +10 -0
- package/dist/MatrixMath.d.ts.map +1 -1
- package/dist/MatrixMath.js +11 -0
- package/dist/MatrixMath.js.map +1 -1
- package/dist/MonteCarlo.d.ts +106 -0
- package/dist/MonteCarlo.d.ts.map +1 -0
- package/dist/MonteCarlo.js +183 -0
- package/dist/MonteCarlo.js.map +1 -0
- package/dist/Random.d.ts +31 -0
- package/dist/Random.d.ts.map +1 -0
- package/dist/Random.js +95 -0
- package/dist/Random.js.map +1 -0
- package/dist/Regression.d.ts +114 -0
- package/dist/Regression.d.ts.map +1 -0
- package/dist/Regression.js +338 -0
- package/dist/Regression.js.map +1 -0
- package/dist/Shader.d.ts +41 -0
- package/dist/Shader.d.ts.map +1 -0
- package/dist/Shader.js +758 -0
- package/dist/Shader.js.map +1 -0
- package/dist/SpecialFunctions.d.ts +3 -3
- package/dist/SpecialFunctions.d.ts.map +1 -1
- package/dist/SpecialFunctions.js +61 -51
- package/dist/SpecialFunctions.js.map +1 -1
- package/dist/Symbolic.d.ts +12 -0
- package/dist/Symbolic.d.ts.map +1 -1
- package/dist/Symbolic.js +16 -0
- package/dist/Symbolic.js.map +1 -1
- package/dist/index.d.ts +7 -1
- package/dist/index.d.ts.map +1 -1
- package/dist/index.js +5 -0
- package/dist/index.js.map +1 -1
- package/package.json +1 -1
- package/src/Distributions.ts +391 -94
- package/src/EigenGeneral.ts +342 -0
- package/src/Frobenius.ts +203 -0
- package/src/Graph.ts +68 -0
- package/src/MatrixMath.ts +13 -0
- package/src/MonteCarlo.ts +257 -0
- package/src/Random.ts +104 -0
- package/src/Regression.ts +423 -0
- package/src/Shader.ts +837 -0
- package/src/SpecialFunctions.ts +52 -48
- package/src/Symbolic.ts +18 -0
- package/src/index.ts +34 -0
package/dist/Distributions.js
CHANGED
|
@@ -1,4 +1,93 @@
|
|
|
1
1
|
import { SpecialFunctions as SF } from "./SpecialFunctions.js";
|
|
2
|
+
/** Build a continuous distribution; `ppf` is derived by bisection unless supplied. */
|
|
3
|
+
function continuous(d) {
|
|
4
|
+
const ppf = d.ppf ??
|
|
5
|
+
((p) => invertContinuous(d.cdf, d.sf, p, d.lower ?? Number.NEGATIVE_INFINITY, d.upper ?? Number.POSITIVE_INFINITY));
|
|
6
|
+
return {
|
|
7
|
+
pdf: d.pdf,
|
|
8
|
+
cdf: d.cdf,
|
|
9
|
+
sf: d.sf,
|
|
10
|
+
ppf,
|
|
11
|
+
quantile: ppf,
|
|
12
|
+
mean: d.mean,
|
|
13
|
+
variance: d.variance,
|
|
14
|
+
std: () => Math.sqrt(d.variance()),
|
|
15
|
+
sample: d.sample,
|
|
16
|
+
rvs: (n) => Array.from({ length: n }, () => d.sample()),
|
|
17
|
+
};
|
|
18
|
+
}
|
|
19
|
+
function discrete(d) {
|
|
20
|
+
const ppf = (p) => {
|
|
21
|
+
if (Number.isNaN(p) || p < 0 || p > 1)
|
|
22
|
+
return Number.NaN;
|
|
23
|
+
if (p === 0)
|
|
24
|
+
return -1;
|
|
25
|
+
const max = d.max ?? Number.POSITIVE_INFINITY;
|
|
26
|
+
if (p === 1)
|
|
27
|
+
return max;
|
|
28
|
+
let hi = Math.max(1, Math.ceil(d.mean() + 4 * Math.sqrt(d.variance())));
|
|
29
|
+
while (hi < max && d.cdf(hi) < p)
|
|
30
|
+
hi *= 2;
|
|
31
|
+
hi = Math.min(hi, max);
|
|
32
|
+
let lo = -1; // cdf(-1) = 0 < p
|
|
33
|
+
while (hi - lo > 1) {
|
|
34
|
+
const mid = Math.floor((lo + hi) / 2);
|
|
35
|
+
if (d.cdf(mid) >= p)
|
|
36
|
+
hi = mid;
|
|
37
|
+
else
|
|
38
|
+
lo = mid;
|
|
39
|
+
}
|
|
40
|
+
return hi;
|
|
41
|
+
};
|
|
42
|
+
return {
|
|
43
|
+
pmf: d.pmf,
|
|
44
|
+
cdf: d.cdf,
|
|
45
|
+
sf: d.sf,
|
|
46
|
+
ppf,
|
|
47
|
+
quantile: ppf,
|
|
48
|
+
mean: d.mean,
|
|
49
|
+
variance: d.variance,
|
|
50
|
+
std: () => Math.sqrt(d.variance()),
|
|
51
|
+
sample: d.sample,
|
|
52
|
+
rvs: (n) => Array.from({ length: n }, () => d.sample()),
|
|
53
|
+
};
|
|
54
|
+
}
|
|
55
|
+
/**
|
|
56
|
+
* Invert a continuous CDF by bracketing + bisection. Uses `cdf` for the lower
|
|
57
|
+
* half and `sf` for the upper half so tail quantiles keep full precision.
|
|
58
|
+
*/
|
|
59
|
+
function invertContinuous(cdf, sf, p, lower, upper) {
|
|
60
|
+
if (Number.isNaN(p) || p < 0 || p > 1)
|
|
61
|
+
return Number.NaN;
|
|
62
|
+
if (p === 0)
|
|
63
|
+
return lower;
|
|
64
|
+
if (p === 1)
|
|
65
|
+
return upper;
|
|
66
|
+
const below = p <= 0.5 ? (x) => cdf(x) < p : (x) => sf(x) > 1 - p;
|
|
67
|
+
let lo = lower;
|
|
68
|
+
if (!Number.isFinite(lo)) {
|
|
69
|
+
lo = -1;
|
|
70
|
+
for (let step = 1; !below(lo); step *= 2)
|
|
71
|
+
lo -= step;
|
|
72
|
+
}
|
|
73
|
+
let hi = Number.isFinite(upper) ? upper : lo + 1;
|
|
74
|
+
if (!Number.isFinite(upper)) {
|
|
75
|
+
for (let step = 1; below(hi); step *= 2) {
|
|
76
|
+
lo = hi;
|
|
77
|
+
hi += step;
|
|
78
|
+
}
|
|
79
|
+
}
|
|
80
|
+
for (let i = 0; i < 400; i++) {
|
|
81
|
+
const mid = lo + (hi - lo) / 2;
|
|
82
|
+
if (mid === lo || mid === hi)
|
|
83
|
+
break;
|
|
84
|
+
if (below(mid))
|
|
85
|
+
lo = mid;
|
|
86
|
+
else
|
|
87
|
+
hi = mid;
|
|
88
|
+
}
|
|
89
|
+
return hi;
|
|
90
|
+
}
|
|
2
91
|
/** The inverse standard-normal CDF (Acklam's rational approximation). */
|
|
3
92
|
function normInv(p) {
|
|
4
93
|
if (p <= 0)
|
|
@@ -34,81 +123,157 @@ function normInv(p) {
|
|
|
34
123
|
return (-(((((c[0] * q + c[1]) * q + c[2]) * q + c[3]) * q + c[4]) * q + c[5]) /
|
|
35
124
|
((((d[0] * q + d[1]) * q + d[2]) * q + d[3]) * q + 1));
|
|
36
125
|
}
|
|
126
|
+
const SQRT_2PI = Math.sqrt(2 * Math.PI);
|
|
37
127
|
export class Distributions {
|
|
38
|
-
/** The inverse standard-normal CDF (quantile / probit function). */
|
|
128
|
+
/** The inverse standard-normal CDF (quantile / probit function), to ~1e-15 after one Halley refinement. */
|
|
39
129
|
static normalQuantile(p, mean = 0, sd = 1) {
|
|
40
|
-
|
|
130
|
+
if (Number.isNaN(p) || p < 0 || p > 1)
|
|
131
|
+
return Number.NaN;
|
|
132
|
+
let x = normInv(p);
|
|
133
|
+
if (Number.isFinite(x)) {
|
|
134
|
+
const e = 0.5 * SF.erfc(-x / Math.SQRT2) - p;
|
|
135
|
+
const u = e * SQRT_2PI * Math.exp((x * x) / 2);
|
|
136
|
+
x -= u / (1 + (x * u) / 2);
|
|
137
|
+
}
|
|
138
|
+
return mean + sd * x;
|
|
41
139
|
}
|
|
140
|
+
/** Normal(`mean`, `sd`) -- `scipy.stats.norm(loc=mean, scale=sd)`. */
|
|
42
141
|
static normal(mean = 0, sd = 1, rng = Math.random) {
|
|
43
|
-
return {
|
|
44
|
-
pdf: (x) => Math.exp(-0.5 * ((x - mean) / sd) ** 2) / (sd *
|
|
142
|
+
return continuous({
|
|
143
|
+
pdf: (x) => Math.exp(-0.5 * ((x - mean) / sd) ** 2) / (sd * SQRT_2PI),
|
|
45
144
|
cdf: (x) => 0.5 * SF.erfc(-(x - mean) / (sd * Math.SQRT2)),
|
|
145
|
+
sf: (x) => 0.5 * SF.erfc((x - mean) / (sd * Math.SQRT2)),
|
|
146
|
+
ppf: (p) => Distributions.normalQuantile(p, mean, sd),
|
|
46
147
|
mean: () => mean,
|
|
47
148
|
variance: () => sd * sd,
|
|
48
149
|
sample: () => mean + sd * Math.sqrt(-2 * Math.log(1 - rng())) * Math.cos(2 * Math.PI * rng()),
|
|
49
|
-
};
|
|
150
|
+
});
|
|
50
151
|
}
|
|
152
|
+
/** Exponential with the given `rate` -- `scipy.stats.expon(scale=1/rate)`. */
|
|
51
153
|
static exponential(rate = 1, rng = Math.random) {
|
|
52
|
-
return {
|
|
154
|
+
return continuous({
|
|
53
155
|
pdf: (x) => (x < 0 ? 0 : rate * Math.exp(-rate * x)),
|
|
54
|
-
cdf: (x) => (x < 0 ? 0 :
|
|
156
|
+
cdf: (x) => (x < 0 ? 0 : -Math.expm1(-rate * x)),
|
|
157
|
+
sf: (x) => (x < 0 ? 1 : Math.exp(-rate * x)),
|
|
158
|
+
ppf: (p) => (p < 0 || p > 1 ? Number.NaN : -Math.log1p(-p) / rate),
|
|
55
159
|
mean: () => 1 / rate,
|
|
56
160
|
variance: () => 1 / (rate * rate),
|
|
57
161
|
sample: () => -Math.log(1 - rng()) / rate,
|
|
58
|
-
};
|
|
162
|
+
});
|
|
59
163
|
}
|
|
164
|
+
/** Uniform on `[a, b]` -- `scipy.stats.uniform(loc=a, scale=b-a)`. */
|
|
60
165
|
static uniform(a = 0, b = 1, rng = Math.random) {
|
|
61
|
-
return {
|
|
166
|
+
return continuous({
|
|
62
167
|
pdf: (x) => (x >= a && x <= b ? 1 / (b - a) : 0),
|
|
63
168
|
cdf: (x) => (x < a ? 0 : x > b ? 1 : (x - a) / (b - a)),
|
|
169
|
+
sf: (x) => (x < a ? 1 : x > b ? 0 : (b - x) / (b - a)),
|
|
170
|
+
ppf: (p) => (p < 0 || p > 1 ? Number.NaN : a + p * (b - a)),
|
|
64
171
|
mean: () => (a + b) / 2,
|
|
65
172
|
variance: () => (b - a) ** 2 / 12,
|
|
66
173
|
sample: () => a + rng() * (b - a),
|
|
67
|
-
};
|
|
174
|
+
});
|
|
68
175
|
}
|
|
176
|
+
/** Gamma with `shape` and `scale` -- `scipy.stats.gamma(a=shape, scale=scale)`. */
|
|
69
177
|
static gamma(shape, scale = 1, rng = Math.random) {
|
|
70
|
-
return {
|
|
71
|
-
pdf: (x) =>
|
|
178
|
+
return continuous({
|
|
179
|
+
pdf: (x) => {
|
|
180
|
+
if (x < 0)
|
|
181
|
+
return 0;
|
|
182
|
+
if (x === 0)
|
|
183
|
+
return shape < 1 ? Number.POSITIVE_INFINITY : shape === 1 ? 1 / scale : 0;
|
|
184
|
+
return Math.exp((shape - 1) * Math.log(x) - x / scale - SF.lnGamma(shape) - shape * Math.log(scale));
|
|
185
|
+
},
|
|
72
186
|
cdf: (x) => (x <= 0 ? 0 : SF.regularizedGammaP(shape, x / scale)),
|
|
187
|
+
sf: (x) => (x <= 0 ? 1 : SF.regularizedGammaQ(shape, x / scale)),
|
|
188
|
+
lower: 0,
|
|
73
189
|
mean: () => shape * scale,
|
|
74
190
|
variance: () => shape * scale * scale,
|
|
75
191
|
sample: () => Distributions.sampleGamma(shape, rng) * scale,
|
|
76
|
-
};
|
|
192
|
+
});
|
|
77
193
|
}
|
|
194
|
+
/** Chi-square with `df` degrees of freedom -- `scipy.stats.chi2(df)`. */
|
|
78
195
|
static chiSquare(df, rng = Math.random) {
|
|
79
196
|
const g = Distributions.gamma(df / 2, 2, rng);
|
|
80
|
-
return { ...g, mean: () => df, variance: () => 2 * df };
|
|
197
|
+
return { ...g, mean: () => df, variance: () => 2 * df, std: () => Math.sqrt(2 * df) };
|
|
81
198
|
}
|
|
199
|
+
/** Student's t with `df` degrees of freedom -- `scipy.stats.t(df)`. */
|
|
82
200
|
static studentT(df, rng = Math.random) {
|
|
83
|
-
const
|
|
84
|
-
|
|
85
|
-
|
|
86
|
-
|
|
87
|
-
|
|
88
|
-
return
|
|
89
|
-
|
|
201
|
+
const lc = SF.lnGamma((df + 1) / 2) - SF.lnGamma(df / 2) - 0.5 * Math.log(df * Math.PI);
|
|
202
|
+
const cdf = (x) => {
|
|
203
|
+
if (x === Number.POSITIVE_INFINITY)
|
|
204
|
+
return 1;
|
|
205
|
+
if (x === Number.NEGATIVE_INFINITY)
|
|
206
|
+
return 0;
|
|
207
|
+
const ib = SF.regularizedIncompleteBeta(df / (df + x * x), df / 2, 0.5);
|
|
208
|
+
return x >= 0 ? 1 - 0.5 * ib : 0.5 * ib;
|
|
209
|
+
};
|
|
210
|
+
return continuous({
|
|
211
|
+
pdf: (x) => Math.exp(lc - ((df + 1) / 2) * Math.log1p((x * x) / df)),
|
|
212
|
+
cdf,
|
|
213
|
+
sf: (x) => cdf(-x),
|
|
90
214
|
mean: () => (df > 1 ? 0 : Number.NaN),
|
|
91
|
-
variance: () => (df > 2 ? df / (df - 2) : Number.POSITIVE_INFINITY),
|
|
215
|
+
variance: () => (df > 2 ? df / (df - 2) : df > 1 ? Number.POSITIVE_INFINITY : Number.NaN),
|
|
92
216
|
sample: () => {
|
|
93
217
|
const z = Distributions.normal(0, 1, rng).sample();
|
|
94
218
|
const chi = Distributions.chiSquare(df, rng).sample();
|
|
95
219
|
return z / Math.sqrt(chi / df);
|
|
96
220
|
},
|
|
97
|
-
};
|
|
221
|
+
});
|
|
222
|
+
}
|
|
223
|
+
/** Fisher-Snedecor F with numerator `d1` and denominator `d2` degrees of freedom -- `scipy.stats.f(d1, d2)`. */
|
|
224
|
+
static fisherF(d1, d2, rng = Math.random) {
|
|
225
|
+
const lbeta = SF.lnGamma(d1 / 2) + SF.lnGamma(d2 / 2) - SF.lnGamma((d1 + d2) / 2);
|
|
226
|
+
return continuous({
|
|
227
|
+
pdf: (x) => {
|
|
228
|
+
if (x < 0)
|
|
229
|
+
return 0;
|
|
230
|
+
if (x === 0)
|
|
231
|
+
return d1 < 2 ? Number.POSITIVE_INFINITY : d1 === 2 ? 1 : 0;
|
|
232
|
+
return Math.exp(0.5 * d1 * Math.log(d1) +
|
|
233
|
+
0.5 * d2 * Math.log(d2) +
|
|
234
|
+
(d1 / 2 - 1) * Math.log(x) -
|
|
235
|
+
((d1 + d2) / 2) * Math.log(d2 + d1 * x) -
|
|
236
|
+
lbeta);
|
|
237
|
+
},
|
|
238
|
+
cdf: (x) => (x <= 0 ? 0 : SF.regularizedIncompleteBeta((d1 * x) / (d1 * x + d2), d1 / 2, d2 / 2)),
|
|
239
|
+
sf: (x) => (x <= 0 ? 1 : SF.regularizedIncompleteBeta(d2 / (d2 + d1 * x), d2 / 2, d1 / 2)),
|
|
240
|
+
lower: 0,
|
|
241
|
+
mean: () => (d2 > 2 ? d2 / (d2 - 2) : Number.NaN),
|
|
242
|
+
variance: () => d2 > 4
|
|
243
|
+
? (2 * d2 * d2 * (d1 + d2 - 2)) / (d1 * (d2 - 2) ** 2 * (d2 - 4))
|
|
244
|
+
: d2 > 2
|
|
245
|
+
? Number.POSITIVE_INFINITY
|
|
246
|
+
: Number.NaN,
|
|
247
|
+
sample: () => (Distributions.sampleGamma(d1 / 2, rng) * 2) / d1 / ((Distributions.sampleGamma(d2 / 2, rng) * 2) / d2),
|
|
248
|
+
});
|
|
98
249
|
}
|
|
250
|
+
/** Binomial(`n`, `p`) -- `scipy.stats.binom(n, p)`. */
|
|
99
251
|
static binomial(n, p, rng = Math.random) {
|
|
100
252
|
const logChoose = (k) => SF.lnGamma(n + 1) - SF.lnGamma(k + 1) - SF.lnGamma(n - k + 1);
|
|
101
|
-
const pmf = (k) =>
|
|
102
|
-
|
|
253
|
+
const pmf = (k) => {
|
|
254
|
+
if (k < 0 || k > n || !Number.isInteger(k))
|
|
255
|
+
return 0;
|
|
256
|
+
if (p === 0)
|
|
257
|
+
return k === 0 ? 1 : 0;
|
|
258
|
+
if (p === 1)
|
|
259
|
+
return k === n ? 1 : 0;
|
|
260
|
+
return Math.exp(logChoose(k) + k * Math.log(p) + (n - k) * Math.log1p(-p));
|
|
261
|
+
};
|
|
262
|
+
return discrete({
|
|
103
263
|
pmf,
|
|
104
264
|
cdf: (k) => {
|
|
105
|
-
|
|
106
|
-
|
|
107
|
-
|
|
108
|
-
|
|
109
|
-
|
|
110
|
-
|
|
265
|
+
if (k < 0)
|
|
266
|
+
return 0;
|
|
267
|
+
const f = Math.floor(k);
|
|
268
|
+
return f >= n ? 1 : SF.regularizedIncompleteBeta(1 - p, n - f, f + 1);
|
|
269
|
+
},
|
|
270
|
+
sf: (k) => {
|
|
271
|
+
if (k < 0)
|
|
272
|
+
return 1;
|
|
273
|
+
const f = Math.floor(k);
|
|
274
|
+
return f >= n ? 0 : SF.regularizedIncompleteBeta(p, f + 1, n - f);
|
|
111
275
|
},
|
|
276
|
+
max: n,
|
|
112
277
|
mean: () => n * p,
|
|
113
278
|
variance: () => n * p * (1 - p),
|
|
114
279
|
sample: () => {
|
|
@@ -118,25 +283,36 @@ export class Distributions {
|
|
|
118
283
|
count++;
|
|
119
284
|
return count;
|
|
120
285
|
},
|
|
121
|
-
};
|
|
286
|
+
});
|
|
122
287
|
}
|
|
288
|
+
/** Poisson(`lambda`) -- `scipy.stats.poisson(mu=lambda)`. */
|
|
123
289
|
static poisson(lambda, rng = Math.random) {
|
|
124
|
-
return {
|
|
290
|
+
return discrete({
|
|
125
291
|
pmf: (k) => (k < 0 || !Number.isInteger(k) ? 0 : Math.exp(k * Math.log(lambda) - lambda - SF.lnGamma(k + 1))),
|
|
126
292
|
cdf: (k) => (k < 0 ? 0 : SF.regularizedGammaQ(Math.floor(k) + 1, lambda)),
|
|
293
|
+
sf: (k) => (k < 0 ? 1 : SF.regularizedGammaP(Math.floor(k) + 1, lambda)),
|
|
127
294
|
mean: () => lambda,
|
|
128
295
|
variance: () => lambda,
|
|
129
296
|
sample: () => {
|
|
130
|
-
|
|
131
|
-
|
|
132
|
-
let
|
|
133
|
-
|
|
134
|
-
|
|
135
|
-
|
|
136
|
-
|
|
137
|
-
|
|
297
|
+
// Knuth's product method underflows exp(-lambda) past ~745, so split
|
|
298
|
+
// large rates into independent Poisson(<=500) pieces and add them.
|
|
299
|
+
let remaining = lambda;
|
|
300
|
+
let total = 0;
|
|
301
|
+
while (remaining > 0) {
|
|
302
|
+
const part = Math.min(remaining, 500);
|
|
303
|
+
remaining -= part;
|
|
304
|
+
const L = Math.exp(-part);
|
|
305
|
+
let k = 0;
|
|
306
|
+
let prod = 1;
|
|
307
|
+
do {
|
|
308
|
+
k++;
|
|
309
|
+
prod *= rng();
|
|
310
|
+
} while (prod > L);
|
|
311
|
+
total += k - 1;
|
|
312
|
+
}
|
|
313
|
+
return total;
|
|
138
314
|
},
|
|
139
|
-
};
|
|
315
|
+
});
|
|
140
316
|
}
|
|
141
317
|
/** Marsaglia-Tsang gamma sampler (shape ≥ 0). */
|
|
142
318
|
static sampleGamma(shape, rng = Math.random) {
|
|
@@ -160,63 +336,138 @@ export class Distributions {
|
|
|
160
336
|
}
|
|
161
337
|
}
|
|
162
338
|
}
|
|
339
|
+
const sumOf = (xs) => xs.reduce((a, b) => a + b, 0);
|
|
340
|
+
const meanOf = (xs) => sumOf(xs) / xs.length;
|
|
341
|
+
const sampleVar = (xs) => {
|
|
342
|
+
const m = meanOf(xs);
|
|
343
|
+
return xs.reduce((acc, x) => acc + (x - m) ** 2, 0) / (xs.length - 1);
|
|
344
|
+
};
|
|
345
|
+
/** p-value of a t statistic under the chosen alternative. */
|
|
346
|
+
function tPValue(t, df, alternative = "two-sided") {
|
|
347
|
+
const dist = Distributions.studentT(df);
|
|
348
|
+
if (alternative === "less")
|
|
349
|
+
return dist.cdf(t);
|
|
350
|
+
if (alternative === "greater")
|
|
351
|
+
return dist.sf(t);
|
|
352
|
+
return Math.min(1, 2 * dist.sf(Math.abs(t)));
|
|
353
|
+
}
|
|
163
354
|
export class HypothesisTests {
|
|
164
|
-
/** One-sample
|
|
165
|
-
static tTestOneSample(sample, mu0 = 0) {
|
|
355
|
+
/** One-sample Student t-test against a hypothesised mean `mu0` -- `scipy.stats.ttest_1samp`. */
|
|
356
|
+
static tTestOneSample(sample, mu0 = 0, options = {}) {
|
|
166
357
|
const n = sample.length;
|
|
167
|
-
const
|
|
168
|
-
const
|
|
169
|
-
const se = Math.sqrt(variance / n);
|
|
170
|
-
const t = (mean - mu0) / se;
|
|
358
|
+
const se = Math.sqrt(sampleVar(sample) / n);
|
|
359
|
+
const t = (meanOf(sample) - mu0) / se;
|
|
171
360
|
const df = n - 1;
|
|
172
|
-
|
|
173
|
-
return { statistic: t, df, pValue: p };
|
|
361
|
+
return { statistic: t, df, pValue: tPValue(t, df, options.alternative) };
|
|
174
362
|
}
|
|
175
|
-
/**
|
|
176
|
-
|
|
363
|
+
/**
|
|
364
|
+
* Welch's two-sample t-test (unequal variances) -- `scipy.stats.ttest_ind(a, b, equal_var=False)`.
|
|
365
|
+
* This is the mallory behaviour; for SciPy's pooled default see {@link tTestIndependent}.
|
|
366
|
+
*/
|
|
367
|
+
static tTestTwoSample(a, b, options = {}) {
|
|
368
|
+
return HypothesisTests.tTestIndependent(a, b, { ...options, equalVar: false });
|
|
369
|
+
}
|
|
370
|
+
/** Alias of Welch's test (`tTestTwoSample`). */
|
|
371
|
+
static welch(a, b, options = {}) {
|
|
372
|
+
return HypothesisTests.tTestTwoSample(a, b, options);
|
|
373
|
+
}
|
|
374
|
+
/** Independent two-sample t-test -- `scipy.stats.ttest_ind`. Pooled-variance (`equalVar: true`) by default, like SciPy. */
|
|
375
|
+
static tTestIndependent(a, b, options = {}) {
|
|
376
|
+
const { equalVar = true, alternative } = options;
|
|
177
377
|
const na = a.length;
|
|
178
378
|
const nb = b.length;
|
|
179
|
-
const
|
|
180
|
-
const
|
|
181
|
-
|
|
182
|
-
|
|
183
|
-
|
|
184
|
-
|
|
185
|
-
|
|
186
|
-
|
|
187
|
-
|
|
188
|
-
|
|
189
|
-
|
|
190
|
-
|
|
191
|
-
|
|
192
|
-
|
|
193
|
-
|
|
194
|
-
return { statistic:
|
|
195
|
-
}
|
|
196
|
-
/**
|
|
379
|
+
const va = sampleVar(a);
|
|
380
|
+
const vb = sampleVar(b);
|
|
381
|
+
let t;
|
|
382
|
+
let df;
|
|
383
|
+
if (equalVar) {
|
|
384
|
+
df = na + nb - 2;
|
|
385
|
+
const pooled = ((na - 1) * va + (nb - 1) * vb) / df;
|
|
386
|
+
t = (meanOf(a) - meanOf(b)) / Math.sqrt(pooled * (1 / na + 1 / nb));
|
|
387
|
+
}
|
|
388
|
+
else {
|
|
389
|
+
const sa = va / na;
|
|
390
|
+
const sb = vb / nb;
|
|
391
|
+
t = (meanOf(a) - meanOf(b)) / Math.sqrt(sa + sb);
|
|
392
|
+
df = (sa + sb) ** 2 / (sa ** 2 / (na - 1) + sb ** 2 / (nb - 1));
|
|
393
|
+
}
|
|
394
|
+
return { statistic: t, df, pValue: tPValue(t, df, alternative) };
|
|
395
|
+
}
|
|
396
|
+
/** Paired t-test on matched samples -- `scipy.stats.ttest_rel`. */
|
|
397
|
+
static tTestPaired(a, b, options = {}) {
|
|
398
|
+
if (a.length !== b.length)
|
|
399
|
+
throw new Error("tTestPaired requires samples of equal length.");
|
|
400
|
+
return HypothesisTests.tTestOneSample(a.map((x, i) => x - b[i]), 0, options);
|
|
401
|
+
}
|
|
402
|
+
/**
|
|
403
|
+
* Pearson chi-square goodness-of-fit test -- `scipy.stats.chisquare`.
|
|
404
|
+
* `expected` defaults to uniform; `ddof` is SciPy's delta degrees of freedom
|
|
405
|
+
* (the p-value uses `k - 1 - ddof`). Throws when the totals disagree.
|
|
406
|
+
*/
|
|
407
|
+
static chiSquareGoodnessOfFit(observed, expected, ddof = 0) {
|
|
408
|
+
const k = observed.length;
|
|
409
|
+
const exp = expected ?? new Array(k).fill(sumOf(observed) / k);
|
|
410
|
+
if (exp.length !== k)
|
|
411
|
+
throw new Error("observed and expected must have the same length.");
|
|
412
|
+
const so = sumOf(observed);
|
|
413
|
+
const se = sumOf(exp);
|
|
414
|
+
if (Math.abs(so - se) > 1e-8 * Math.max(Math.abs(so), Math.abs(se))) {
|
|
415
|
+
throw new Error(`observed (${so}) and expected (${se}) frequencies must sum to the same total.`);
|
|
416
|
+
}
|
|
417
|
+
const chi2 = observed.reduce((s, o, i) => s + (o - exp[i]) ** 2 / exp[i], 0);
|
|
418
|
+
const df = k - 1 - ddof;
|
|
419
|
+
return { statistic: chi2, df, pValue: Distributions.chiSquare(df).sf(chi2) };
|
|
420
|
+
}
|
|
421
|
+
/**
|
|
422
|
+
* Chi-square test of independence on a contingency table --
|
|
423
|
+
* `scipy.stats.chi2_contingency`. Like SciPy, applies Yates' continuity
|
|
424
|
+
* correction by default when there is one degree of freedom.
|
|
425
|
+
*/
|
|
426
|
+
static chiSquareIndependence(table, options = {}) {
|
|
427
|
+
const { correction = true } = options;
|
|
428
|
+
const r = table.length;
|
|
429
|
+
const c = table[0].length;
|
|
430
|
+
const rows = table.map(sumOf);
|
|
431
|
+
const cols = Array.from({ length: c }, (_, j) => sumOf(table.map((row) => row[j])));
|
|
432
|
+
const total = sumOf(rows);
|
|
433
|
+
const expected = table.map((_, i) => cols.map((cj) => (rows[i] * cj) / total));
|
|
434
|
+
const df = (r - 1) * (c - 1);
|
|
435
|
+
let stat = 0;
|
|
436
|
+
for (let i = 0; i < r; i++) {
|
|
437
|
+
for (let j = 0; j < c; j++) {
|
|
438
|
+
let diff = Math.abs(table[i][j] - expected[i][j]);
|
|
439
|
+
if (correction && df === 1)
|
|
440
|
+
diff = Math.max(0, diff - 0.5);
|
|
441
|
+
stat += diff ** 2 / expected[i][j];
|
|
442
|
+
}
|
|
443
|
+
}
|
|
444
|
+
return { statistic: stat, df, pValue: df === 0 ? 1 : Distributions.chiSquare(df).sf(stat), expected };
|
|
445
|
+
}
|
|
446
|
+
/** t-based confidence interval for the mean of a sample. */
|
|
197
447
|
static confidenceIntervalMean(sample, level = 0.95) {
|
|
198
448
|
const n = sample.length;
|
|
199
|
-
const m = sample
|
|
200
|
-
const
|
|
201
|
-
const
|
|
202
|
-
const alpha = 1 - level;
|
|
203
|
-
const tCrit = tQuantile(1 - alpha / 2, n - 1);
|
|
449
|
+
const m = meanOf(sample);
|
|
450
|
+
const se = Math.sqrt(sampleVar(sample) / n);
|
|
451
|
+
const tCrit = Distributions.studentT(n - 1).ppf(1 - (1 - level) / 2);
|
|
204
452
|
return [m - tCrit * se, m + tCrit * se];
|
|
205
453
|
}
|
|
206
|
-
|
|
207
|
-
|
|
208
|
-
|
|
209
|
-
|
|
210
|
-
|
|
211
|
-
|
|
212
|
-
|
|
213
|
-
|
|
214
|
-
|
|
215
|
-
|
|
216
|
-
|
|
217
|
-
|
|
218
|
-
|
|
454
|
+
/** Welch confidence interval for `mean(a) - mean(b)`. */
|
|
455
|
+
static confidenceIntervalDifference(a, b, level = 0.95) {
|
|
456
|
+
const sa = sampleVar(a) / a.length;
|
|
457
|
+
const sb = sampleVar(b) / b.length;
|
|
458
|
+
const df = (sa + sb) ** 2 / (sa ** 2 / (a.length - 1) + sb ** 2 / (b.length - 1));
|
|
459
|
+
const d = meanOf(a) - meanOf(b);
|
|
460
|
+
const half = Distributions.studentT(df).ppf(1 - (1 - level) / 2) * Math.sqrt(sa + sb);
|
|
461
|
+
return [d - half, d + half];
|
|
462
|
+
}
|
|
463
|
+
/** Wilson score interval for a binomial proportion (`scipy.stats.binomtest(...).proportion_ci(method="wilson", ...)` without continuity correction). */
|
|
464
|
+
static confidenceIntervalProportion(successes, n, level = 0.95) {
|
|
465
|
+
const z = Distributions.normalQuantile(1 - (1 - level) / 2);
|
|
466
|
+
const phat = successes / n;
|
|
467
|
+
const denom = 1 + (z * z) / n;
|
|
468
|
+
const centre = (phat + (z * z) / (2 * n)) / denom;
|
|
469
|
+
const half = (z * Math.sqrt((phat * (1 - phat)) / n + (z * z) / (4 * n * n))) / denom;
|
|
470
|
+
return [Math.max(0, centre - half), Math.min(1, centre + half)];
|
|
219
471
|
}
|
|
220
|
-
return (lo + hi) / 2;
|
|
221
472
|
}
|
|
222
473
|
//# sourceMappingURL=Distributions.js.map
|
|
@@ -1 +1 @@
|
|
|
1
|
-
{"version":3,"file":"Distributions.js","sourceRoot":"","sources":["../src/Distributions.ts"],"names":[],"mappings":"AAAA,OAAO,EAAE,gBAAgB,IAAI,EAAE,EAAE,MAAM,uBAAuB,CAAC;AA8B/D,yEAAyE;AACzE,SAAS,OAAO,CAAC,CAAS;IACxB,IAAI,CAAC,IAAI,CAAC;QAAE,OAAO,MAAM,CAAC,iBAAiB,CAAC;IAC5C,IAAI,CAAC,IAAI,CAAC;QAAE,OAAO,MAAM,CAAC,iBAAiB,CAAC;IAC5C,MAAM,CAAC,GAAG;QACR,CAAC,mBAAmB,EAAE,mBAAmB,EAAE,CAAC,mBAAmB,EAAE,kBAAkB,EAAE,CAAC,mBAAmB;QACzG,iBAAiB;KAClB,CAAC;IACF,MAAM,CAAC,GAAG;QACR,CAAC,mBAAmB,EAAE,mBAAmB,EAAE,CAAC,mBAAmB,EAAE,mBAAmB,EAAE,CAAC,mBAAmB;KAC3G,CAAC;IACF,MAAM,CAAC,GAAG;QACR,CAAC,oBAAoB,EAAE,CAAC,oBAAoB,EAAE,CAAC,iBAAiB,EAAE,CAAC,iBAAiB,EAAE,iBAAiB;QACvG,iBAAiB;KAClB,CAAC;IACF,MAAM,CAAC,GAAG,CAAC,oBAAoB,EAAE,oBAAoB,EAAE,iBAAiB,EAAE,iBAAiB,CAAC,CAAC;IAC7F,MAAM,IAAI,GAAG,OAAO,CAAC;IACrB,MAAM,KAAK,GAAG,CAAC,GAAG,IAAI,CAAC;IACvB,IAAI,CAAC,GAAG,IAAI,EAAE,CAAC;QACb,MAAM,CAAC,GAAG,IAAI,CAAC,IAAI,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC;QACtC,OAAO,CACL,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC;YACrE,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,CACtD,CAAC;IACJ,CAAC;IACD,IAAI,CAAC,IAAI,KAAK,EAAE,CAAC;QACf,MAAM,CAAC,GAAG,CAAC,GAAG,GAAG,CAAC;QAClB,MAAM,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC;QAChB,OAAO,CACL,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC;YAC3E,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,CACnE,CAAC;IACJ,CAAC;IACD,MAAM,CAAC,GAAG,IAAI,CAAC,IAAI,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC;IAC1C,OAAO,CACL,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC;QACtE,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,CACtD,CAAC;AACJ,CAAC;AAED,MAAM,OAAO,aAAa;IACxB,oEAAoE;IACpE,MAAM,CAAC,cAAc,CAAC,CAAS,EAAE,IAAI,GAAG,CAAC,EAAE,EAAE,GAAG,CAAC;QAC/C,OAAO,IAAI,GAAG,EAAE,GAAG,OAAO,CAAC,CAAC,CAAC,CAAC;IAChC,CAAC;IAED,MAAM,CAAC,MAAM,CAAC,IAAI,GAAG,CAAC,EAAE,EAAE,GAAG,CAAC,EAAE,MAAoB,IAAI,CAAC,MAAM;QAC7D,OAAO;YACL,GAAG,EAAE,CAAC,CAAC,EAAE,EAAE,CAAC,IAAI,CAAC,GAAG,CAAC,CAAC,GAAG,GAAG,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,GAAG,EAAE,CAAC,IAAI,CAAC,CAAC,GAAG,CAAC,EAAE,GAAG,IAAI,CAAC,IAAI,CAAC,CAAC,GAAG,IAAI,CAAC,EAAE,CAAC,CAAC;YACnF,GAAG,EAAE,CAAC,CAAC,EAAE,EAAE,CAAC,GAAG,GAAG,EAAE,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,EAAE,GAAG,IAAI,CAAC,KAAK,CAAC,CAAC;YAC1D,IAAI,EAAE,GAAG,EAAE,CAAC,IAAI;YAChB,QAAQ,EAAE,GAAG,EAAE,CAAC,EAAE,GAAG,EAAE;YACvB,MAAM,EAAE,GAAG,EAAE,CAAC,IAAI,GAAG,EAAE,GAAG,IAAI,CAAC,IAAI,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,GAAG,GAAG,EAAE,CAAC,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,GAAG,IAAI,CAAC,EAAE,GAAG,GAAG,EAAE,CAAC;SAC9F,CAAC;IACJ,CAAC;IAED,MAAM,CAAC,WAAW,CAAC,IAAI,GAAG,CAAC,EAAE,MAAoB,IAAI,CAAC,MAAM;QAC1D,OAAO;YACL,GAAG,EAAE,CAAC,CAAC,EAAE,EAAE,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,IAAI,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,IAAI,GAAG,CAAC,CAAC,CAAC;YACpD,GAAG,EAAE,CAAC,CAAC,EAAE,EAAE,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,IAAI,GAAG,CAAC,CAAC,CAAC;YACjD,IAAI,EAAE,GAAG,EAAE,CAAC,CAAC,GAAG,IAAI;YACpB,QAAQ,EAAE,GAAG,EAAE,CAAC,CAAC,GAAG,CAAC,IAAI,GAAG,IAAI,CAAC;YACjC,MAAM,EAAE,GAAG,EAAE,CAAC,CAAC,IAAI,CAAC,GAAG,CAAC,CAAC,GAAG,GAAG,EAAE,CAAC,GAAG,IAAI;SAC1C,CAAC;IACJ,CAAC;IAED,MAAM,CAAC,OAAO,CAAC,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,EAAE,MAAoB,IAAI,CAAC,MAAM;QAC1D,OAAO;YACL,GAAG,EAAE,CAAC,CAAC,EAAE,EAAE,CAAC,CAAC,CAAC,IAAI,CAAC,IAAI,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;YAChD,GAAG,EAAE,CAAC,CAAC,EAAE,EAAE,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC;YACvD,IAAI,EAAE,GAAG,EAAE,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC;YACvB,QAAQ,EAAE,GAAG,EAAE,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,IAAI,CAAC,GAAG,EAAE;YACjC,MAAM,EAAE,GAAG,EAAE,CAAC,CAAC,GAAG,GAAG,EAAE,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC;SAClC,CAAC;IACJ,CAAC;IAED,MAAM,CAAC,KAAK,CAAC,KAAa,EAAE,KAAK,GAAG,CAAC,EAAE,MAAoB,IAAI,CAAC,MAAM;QACpE,OAAO;YACL,GAAG,EAAE,CAAC,CAAC,EAAE,EAAE,CAAC,CAAC,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,IAAI,CAAC,KAAK,GAAG,CAAC,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,GAAG,KAAK,CAAC,CAAC,GAAG,CAAC,EAAE,CAAC,KAAK,CAAC,KAAK,CAAC,GAAG,KAAK,IAAI,KAAK,CAAC,CAAC;YACzG,GAAG,EAAE,CAAC,CAAC,EAAE,EAAE,CAAC,CAAC,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,EAAE,CAAC,iBAAiB,CAAC,KAAK,EAAE,CAAC,GAAG,KAAK,CAAC,CAAC;YACjE,IAAI,EAAE,GAAG,EAAE,CAAC,KAAK,GAAG,KAAK;YACzB,QAAQ,EAAE,GAAG,EAAE,CAAC,KAAK,GAAG,KAAK,GAAG,KAAK;YACrC,MAAM,EAAE,GAAG,EAAE,CAAC,aAAa,CAAC,WAAW,CAAC,KAAK,EAAE,GAAG,CAAC,GAAG,KAAK;SAC5D,CAAC;IACJ,CAAC;IAED,MAAM,CAAC,SAAS,CAAC,EAAU,EAAE,MAAoB,IAAI,CAAC,MAAM;QAC1D,MAAM,CAAC,GAAG,aAAa,CAAC,KAAK,CAAC,EAAE,GAAG,CAAC,EAAE,CAAC,EAAE,GAAG,CAAC,CAAC;QAC9C,OAAO,EAAE,GAAG,CAAC,EAAE,IAAI,EAAE,GAAG,EAAE,CAAC,EAAE,EAAE,QAAQ,EAAE,GAAG,EAAE,CAAC,CAAC,GAAG,EAAE,EAAE,CAAC;IAC1D,CAAC;IAED,MAAM,CAAC,QAAQ,CAAC,EAAU,EAAE,MAAoB,IAAI,CAAC,MAAM;QACzD,MAAM,CAAC,GAAG,EAAE,CAAC,KAAK,CAAC,CAAC,EAAE,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,IAAI,CAAC,IAAI,CAAC,EAAE,GAAG,IAAI,CAAC,EAAE,CAAC,GAAG,EAAE,CAAC,KAAK,CAAC,EAAE,GAAG,CAAC,CAAC,CAAC,CAAC;QAChF,OAAO;YACL,GAAG,EAAE,CAAC,CAAC,EAAE,EAAE,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,EAAE,CAAC,IAAI,CAAC,CAAC,CAAC,EAAE,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC;YACrD,GAAG,EAAE,CAAC,CAAC,EAAE,EAAE;gBACT,MAAM,EAAE,GAAG,EAAE,CAAC,yBAAyB,CAAC,EAAE,GAAG,CAAC,EAAE,GAAG,CAAC,GAAG,CAAC,CAAC,EAAE,EAAE,GAAG,CAAC,EAAE,GAAG,CAAC,CAAC;gBACxE,OAAO,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,GAAG,GAAG,EAAE,CAAC,CAAC,CAAC,GAAG,GAAG,EAAE,CAAC;YAC1C,CAAC;YACD,IAAI,EAAE,GAAG,EAAE,CAAC,CAAC,EAAE,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,MAAM,CAAC,GAAG,CAAC;YACrC,QAAQ,EAAE,GAAG,EAAE,CAAC,CAAC,EAAE,GAAG,CAAC,CAAC,CAAC,CAAC,EAAE,GAAG,CAAC,EAAE,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,MAAM,CAAC,iBAAiB,CAAC;YACnE,MAAM,EAAE,GAAG,EAAE;gBACX,MAAM,CAAC,GAAG,aAAa,CAAC,MAAM,CAAC,CAAC,EAAE,CAAC,EAAE,GAAG,CAAC,CAAC,MAAM,EAAE,CAAC;gBACnD,MAAM,GAAG,GAAG,aAAa,CAAC,SAAS,CAAC,EAAE,EAAE,GAAG,CAAC,CAAC,MAAM,EAAE,CAAC;gBACtD,OAAO,CAAC,GAAG,IAAI,CAAC,IAAI,CAAC,GAAG,GAAG,EAAE,CAAC,CAAC;YACjC,CAAC;SACF,CAAC;IACJ,CAAC;IAED,MAAM,CAAC,QAAQ,CAAC,CAAS,EAAE,CAAS,EAAE,MAAoB,IAAI,CAAC,MAAM;QACnE,MAAM,SAAS,GAAG,CAAC,CAAS,EAAE,EAAE,CAAC,EAAE,CAAC,OAAO,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,EAAE,CAAC,OAAO,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,EAAE,CAAC,OAAO,CAAC,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,CAAC;QAC/F,MAAM,GAAG,GAAG,CAAC,CAAS,EAAU,EAAE,CAChC,CAAC,GAAG,CAAC,IAAI,CAAC,GAAG,CAAC,IAAI,CAAC,MAAM,CAAC,SAAS,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,IAAI,CAAC,GAAG,CAAC,SAAS,CAAC,CAAC,CAAC,GAAG,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC;QACpH,OAAO;YACL,GAAG;YACH,GAAG,EAAE,CAAC,CAAC,EAAE,EAAE;gBACT,uEAAuE;gBACvE,mEAAmE;gBACnE,IAAI,CAAC,GAAG,CAAC,CAAC;gBACV,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,IAAI,IAAI,CAAC,KAAK,CAAC,CAAC,CAAC,EAAE,CAAC,EAAE;oBAAE,CAAC,IAAI,GAAG,CAAC,CAAC,CAAC,CAAC;gBACrD,OAAO,IAAI,CAAC,GAAG,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC;YACxB,CAAC;YACD,IAAI,EAAE,GAAG,EAAE,CAAC,CAAC,GAAG,CAAC;YACjB,QAAQ,EAAE,GAAG,EAAE,CAAC,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC;YAC/B,MAAM,EAAE,GAAG,EAAE;gBACX,IAAI,KAAK,GAAG,CAAC,CAAC;gBACd,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,EAAE,CAAC,EAAE;oBAAE,IAAI,GAAG,EAAE,GAAG,CAAC;wBAAE,KAAK,EAAE,CAAC;gBACnD,OAAO,KAAK,CAAC;YACf,CAAC;SACF,CAAC;IACJ,CAAC;IAED,MAAM,CAAC,OAAO,CAAC,MAAc,EAAE,MAAoB,IAAI,CAAC,MAAM;QAC5D,OAAO;YACL,GAAG,EAAE,CAAC,CAAC,EAAE,EAAE,CAAC,CAAC,CAAC,GAAG,CAAC,IAAI,CAAC,MAAM,CAAC,SAAS,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,IAAI,CAAC,GAAG,CAAC,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,MAAM,CAAC,GAAG,MAAM,GAAG,EAAE,CAAC,OAAO,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC;YAC7G,GAAG,EAAE,CAAC,CAAC,EAAE,EAAE,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,EAAE,CAAC,iBAAiB,CAAC,IAAI,CAAC,KAAK,CAAC,CAAC,CAAC,GAAG,CAAC,EAAE,MAAM,CAAC,CAAC;YACzE,IAAI,EAAE,GAAG,EAAE,CAAC,MAAM;YAClB,QAAQ,EAAE,GAAG,EAAE,CAAC,MAAM;YACtB,MAAM,EAAE,GAAG,EAAE;gBACX,MAAM,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,MAAM,CAAC,CAAC;gBAC5B,IAAI,CAAC,GAAG,CAAC,CAAC;gBACV,IAAI,IAAI,GAAG,CAAC,CAAC;gBACb,GAAG,CAAC;oBACF,CAAC,EAAE,CAAC;oBACJ,IAAI,IAAI,GAAG,EAAE,CAAC;gBAChB,CAAC,QAAQ,IAAI,GAAG,CAAC,EAAE;gBACnB,OAAO,CAAC,GAAG,CAAC,CAAC;YACf,CAAC;SACF,CAAC;IACJ,CAAC;IAED,iDAAiD;IACzC,MAAM,CAAC,WAAW,CAAC,KAAa,EAAE,MAAoB,IAAI,CAAC,MAAM;QACvE,IAAI,KAAK,GAAG,CAAC;YAAE,OAAO,aAAa,CAAC,WAAW,CAAC,KAAK,GAAG,CAAC,EAAE,GAAG,CAAC,GAAG,GAAG,EAAE,IAAI,CAAC,CAAC,GAAG,KAAK,CAAC,CAAC;QACvF,MAAM,CAAC,GAAG,KAAK,GAAG,CAAC,GAAG,CAAC,CAAC;QACxB,MAAM,CAAC,GAAG,CAAC,GAAG,IAAI,CAAC,IAAI,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC;QAC/B,SAAS,CAAC;YACR,IAAI,CAAS,CAAC;YACd,IAAI,CAAS,CAAC;YACd,GAAG,CAAC;gBACF,CAAC,GAAG,aAAa,CAAC,MAAM,CAAC,CAAC,EAAE,CAAC,EAAE,GAAG,CAAC,CAAC,MAAM,EAAE,CAAC;gBAC7C,CAAC,GAAG,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC;YAChB,CAAC,QAAQ,CAAC,IAAI,CAAC,EAAE;YACjB,CAAC,GAAG,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC;YACd,MAAM,CAAC,GAAG,GAAG,EAAE,CAAC;YAChB,IAAI,CAAC,GAAG,CAAC,GAAG,MAAM,GAAG,CAAC,IAAI,CAAC;gBAAE,OAAO,CAAC,GAAG,CAAC,CAAC;YAC1C,IAAI,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,GAAG,GAAG,GAAG,CAAC,GAAG,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC;gBAAE,OAAO,CAAC,GAAG,CAAC,CAAC;QAC1E,CAAC;IACH,CAAC;CACF;AAQD,MAAM,OAAO,eAAe;IAC1B,8EAA8E;IAC9E,MAAM,CAAC,cAAc,CAAC,MAAgB,EAAE,GAAG,GAAG,CAAC;QAC7C,MAAM,CAAC,GAAG,MAAM,CAAC,MAAM,CAAC;QACxB,MAAM,IAAI,GAAG,MAAM,CAAC,MAAM,CAAC,CAAC,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC,CAAC,GAAG,CAAC,EAAE,CAAC,CAAC,GAAG,CAAC,CAAC;QACnD,MAAM,QAAQ,GAAG,MAAM,CAAC,MAAM,CAAC,CAAC,GAAG,EAAE,CAAC,EAAE,EAAE,CAAC,GAAG,GAAG,CAAC,CAAC,GAAG,IAAI,CAAC,IAAI,CAAC,EAAE,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC;QAC/E,MAAM,EAAE,GAAG,IAAI,CAAC,IAAI,CAAC,QAAQ,GAAG,CAAC,CAAC,CAAC;QACnC,MAAM,CAAC,GAAG,CAAC,IAAI,GAAG,GAAG,CAAC,GAAG,EAAE,CAAC;QAC5B,MAAM,EAAE,GAAG,CAAC,GAAG,CAAC,CAAC;QACjB,MAAM,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,GAAG,aAAa,CAAC,QAAQ,CAAC,EAAE,CAAC,CAAC,GAAG,CAAC,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;QAChE,OAAO,EAAE,SAAS,EAAE,CAAC,EAAE,EAAE,EAAE,MAAM,EAAE,CAAC,EAAE,CAAC;IACzC,CAAC;IAED,qDAAqD;IACrD,MAAM,CAAC,cAAc,CAAC,CAAW,EAAE,CAAW;QAC5C,MAAM,EAAE,GAAG,CAAC,CAAC,MAAM,CAAC;QACpB,MAAM,EAAE,GAAG,CAAC,CAAC,MAAM,CAAC;QACpB,MAAM,EAAE,GAAG,CAAC,CAAC,MAAM,CAAC,CAAC,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC,CAAC,GAAG,CAAC,EAAE,CAAC,CAAC,GAAG,EAAE,CAAC;QAC7C,MAAM,EAAE,GAAG,CAAC,CAAC,MAAM,CAAC,CAAC,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC,CAAC,GAAG,CAAC,EAAE,CAAC,CAAC,GAAG,EAAE,CAAC;QAC7C,MAAM,EAAE,GAAG,CAAC,CAAC,MAAM,CAAC,CAAC,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,EAAE,CAAC,IAAI,CAAC,EAAE,CAAC,CAAC,GAAG,CAAC,EAAE,GAAG,CAAC,CAAC,CAAC;QAC/D,MAAM,EAAE,GAAG,CAAC,CAAC,MAAM,CAAC,CAAC,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,EAAE,CAAC,IAAI,CAAC,EAAE,CAAC,CAAC,GAAG,CAAC,EAAE,GAAG,CAAC,CAAC,CAAC;QAC/D,MAAM,EAAE,GAAG,IAAI,CAAC,IAAI,CAAC,EAAE,GAAG,EAAE,GAAG,EAAE,GAAG,EAAE,CAAC,CAAC;QACxC,MAAM,CAAC,GAAG,CAAC,EAAE,GAAG,EAAE,CAAC,GAAG,EAAE,CAAC;QACzB,MAAM,EAAE,GAAG,CAAC,EAAE,GAAG,EAAE,GAAG,EAAE,GAAG,EAAE,CAAC,IAAI,CAAC,GAAG,CAAC,CAAC,EAAE,GAAG,EAAE,CAAC,IAAI,CAAC,GAAG,CAAC,EAAE,GAAG,CAAC,CAAC,GAAG,CAAC,EAAE,GAAG,EAAE,CAAC,IAAI,CAAC,GAAG,CAAC,EAAE,GAAG,CAAC,CAAC,CAAC,CAAC;QAC9F,MAAM,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,GAAG,aAAa,CAAC,QAAQ,CAAC,EAAE,CAAC,CAAC,GAAG,CAAC,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;QAChE,OAAO,EAAE,SAAS,EAAE,CAAC,EAAE,EAAE,EAAE,MAAM,EAAE,CAAC,EAAE,CAAC;IACzC,CAAC;IAED,+CAA+C;IAC/C,MAAM,CAAC,sBAAsB,CAAC,QAAkB,EAAE,QAAkB;QAClE,MAAM,IAAI,GAAG,QAAQ,CAAC,MAAM,CAAC,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,QAAQ,CAAC,CAAC,CAAC,CAAC,IAAI,CAAC,GAAG,QAAQ,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC;QACvF,MAAM,EAAE,GAAG,QAAQ,CAAC,MAAM,GAAG,CAAC,CAAC;QAC/B,MAAM,CAAC,GAAG,CAAC,GAAG,aAAa,CAAC,SAAS,CAAC,EAAE,CAAC,CAAC,GAAG,CAAC,IAAI,CAAC,CAAC;QACpD,OAAO,EAAE,SAAS,EAAE,IAAI,EAAE,EAAE,EAAE,MAAM,EAAE,CAAC,EAAE,CAAC;IAC5C,CAAC;IAED,iFAAiF;IACjF,MAAM,CAAC,sBAAsB,CAAC,MAAgB,EAAE,KAAK,GAAG,IAAI;QAC1D,MAAM,CAAC,GAAG,MAAM,CAAC,MAAM,CAAC;QACxB,MAAM,CAAC,GAAG,MAAM,CAAC,MAAM,CAAC,CAAC,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC,CAAC,GAAG,CAAC,EAAE,CAAC,CAAC,GAAG,CAAC,CAAC;QAChD,MAAM,QAAQ,GAAG,MAAM,CAAC,MAAM,CAAC,CAAC,GAAG,EAAE,CAAC,EAAE,EAAE,CAAC,GAAG,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,IAAI,CAAC,EAAE,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC;QAC5E,MAAM,EAAE,GAAG,IAAI,CAAC,IAAI,CAAC,QAAQ,GAAG,CAAC,CAAC,CAAC;QACnC,MAAM,KAAK,GAAG,CAAC,GAAG,KAAK,CAAC;QACxB,MAAM,KAAK,GAAG,SAAS,CAAC,CAAC,GAAG,KAAK,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,CAAC,CAAC;QAC9C,OAAO,CAAC,CAAC,GAAG,KAAK,GAAG,EAAE,EAAE,CAAC,GAAG,KAAK,GAAG,EAAE,CAAC,CAAC;IAC1C,CAAC;CACF;AAED,8EAA8E;AAE9E,sDAAsD;AACtD,SAAS,SAAS,CAAC,CAAS,EAAE,EAAU;IACtC,MAAM,IAAI,GAAG,aAAa,CAAC,QAAQ,CAAC,EAAE,CAAC,CAAC;IACxC,IAAI,EAAE,GAAG,CAAC,GAAG,CAAC;IACd,IAAI,EAAE,GAAG,GAAG,CAAC;IACb,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,GAAG,EAAE,CAAC,EAAE,EAAE,CAAC;QAC7B,MAAM,GAAG,GAAG,CAAC,EAAE,GAAG,EAAE,CAAC,GAAG,CAAC,CAAC;QAC1B,IAAI,IAAI,CAAC,GAAG,CAAC,GAAG,CAAC,GAAG,CAAC;YAAE,EAAE,GAAG,GAAG,CAAC;;YAC3B,EAAE,GAAG,GAAG,CAAC;IAChB,CAAC;IACD,OAAO,CAAC,EAAE,GAAG,EAAE,CAAC,GAAG,CAAC,CAAC;AACvB,CAAC"}
|
|
1
|
+
{"version":3,"file":"Distributions.js","sourceRoot":"","sources":["../src/Distributions.ts"],"names":[],"mappings":"AAAA,OAAO,EAAE,gBAAgB,IAAI,EAAE,EAAE,MAAM,uBAAuB,CAAC;AAmD/D,sFAAsF;AACtF,SAAS,UAAU,CACjB,CAKC;IAED,MAAM,GAAG,GACP,CAAC,CAAC,GAAG;QACL,CAAC,CAAC,CAAS,EAAE,EAAE,CACb,gBAAgB,CAAC,CAAC,CAAC,GAAG,EAAE,CAAC,CAAC,EAAE,EAAE,CAAC,EAAE,CAAC,CAAC,KAAK,IAAI,MAAM,CAAC,iBAAiB,EAAE,CAAC,CAAC,KAAK,IAAI,MAAM,CAAC,iBAAiB,CAAC,CAAC,CAAC;IAChH,OAAO;QACL,GAAG,EAAE,CAAC,CAAC,GAAG;QACV,GAAG,EAAE,CAAC,CAAC,GAAG;QACV,EAAE,EAAE,CAAC,CAAC,EAAE;QACR,GAAG;QACH,QAAQ,EAAE,GAAG;QACb,IAAI,EAAE,CAAC,CAAC,IAAI;QACZ,QAAQ,EAAE,CAAC,CAAC,QAAQ;QACpB,GAAG,EAAE,GAAG,EAAE,CAAC,IAAI,CAAC,IAAI,CAAC,CAAC,CAAC,QAAQ,EAAE,CAAC;QAClC,MAAM,EAAE,CAAC,CAAC,MAAM;QAChB,GAAG,EAAE,CAAC,CAAC,EAAE,EAAE,CAAC,KAAK,CAAC,IAAI,CAAC,EAAE,MAAM,EAAE,CAAC,EAAE,EAAE,GAAG,EAAE,CAAC,CAAC,CAAC,MAAM,EAAE,CAAC;KACxD,CAAC;AACJ,CAAC;AAED,SAAS,QAAQ,CACf,CAAuG;IAEvG,MAAM,GAAG,GAAG,CAAC,CAAS,EAAU,EAAE;QAChC,IAAI,MAAM,CAAC,KAAK,CAAC,CAAC,CAAC,IAAI,CAAC,GAAG,CAAC,IAAI,CAAC,GAAG,CAAC;YAAE,OAAO,MAAM,CAAC,GAAG,CAAC;QACzD,IAAI,CAAC,KAAK,CAAC;YAAE,OAAO,CAAC,CAAC,CAAC;QACvB,MAAM,GAAG,GAAG,CAAC,CAAC,GAAG,IAAI,MAAM,CAAC,iBAAiB,CAAC;QAC9C,IAAI,CAAC,KAAK,CAAC;YAAE,OAAO,GAAG,CAAC;QACxB,IAAI,EAAE,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,EAAE,IAAI,CAAC,IAAI,CAAC,CAAC,CAAC,IAAI,EAAE,GAAG,CAAC,GAAG,IAAI,CAAC,IAAI,CAAC,CAAC,CAAC,QAAQ,EAAE,CAAC,CAAC,CAAC,CAAC;QACxE,OAAO,EAAE,GAAG,GAAG,IAAI,CAAC,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC;YAAE,EAAE,IAAI,CAAC,CAAC;QAC1C,EAAE,GAAG,IAAI,CAAC,GAAG,CAAC,EAAE,EAAE,GAAG,CAAC,CAAC;QACvB,IAAI,EAAE,GAAG,CAAC,CAAC,CAAC,CAAC,kBAAkB;QAC/B,OAAO,EAAE,GAAG,EAAE,GAAG,CAAC,EAAE,CAAC;YACnB,MAAM,GAAG,GAAG,IAAI,CAAC,KAAK,CAAC,CAAC,EAAE,GAAG,EAAE,CAAC,GAAG,CAAC,CAAC,CAAC;YACtC,IAAI,CAAC,CAAC,GAAG,CAAC,GAAG,CAAC,IAAI,CAAC;gBAAE,EAAE,GAAG,GAAG,CAAC;;gBACzB,EAAE,GAAG,GAAG,CAAC;QAChB,CAAC;QACD,OAAO,EAAE,CAAC;IACZ,CAAC,CAAC;IACF,OAAO;QACL,GAAG,EAAE,CAAC,CAAC,GAAG;QACV,GAAG,EAAE,CAAC,CAAC,GAAG;QACV,EAAE,EAAE,CAAC,CAAC,EAAE;QACR,GAAG;QACH,QAAQ,EAAE,GAAG;QACb,IAAI,EAAE,CAAC,CAAC,IAAI;QACZ,QAAQ,EAAE,CAAC,CAAC,QAAQ;QACpB,GAAG,EAAE,GAAG,EAAE,CAAC,IAAI,CAAC,IAAI,CAAC,CAAC,CAAC,QAAQ,EAAE,CAAC;QAClC,MAAM,EAAE,CAAC,CAAC,MAAM;QAChB,GAAG,EAAE,CAAC,CAAC,EAAE,EAAE,CAAC,KAAK,CAAC,IAAI,CAAC,EAAE,MAAM,EAAE,CAAC,EAAE,EAAE,GAAG,EAAE,CAAC,CAAC,CAAC,MAAM,EAAE,CAAC;KACxD,CAAC;AACJ,CAAC;AAED;;;GAGG;AACH,SAAS,gBAAgB,CACvB,GAA0B,EAC1B,EAAyB,EACzB,CAAS,EACT,KAAa,EACb,KAAa;IAEb,IAAI,MAAM,CAAC,KAAK,CAAC,CAAC,CAAC,IAAI,CAAC,GAAG,CAAC,IAAI,CAAC,GAAG,CAAC;QAAE,OAAO,MAAM,CAAC,GAAG,CAAC;IACzD,IAAI,CAAC,KAAK,CAAC;QAAE,OAAO,KAAK,CAAC;IAC1B,IAAI,CAAC,KAAK,CAAC;QAAE,OAAO,KAAK,CAAC;IAC1B,MAAM,KAAK,GAAG,CAAC,IAAI,GAAG,CAAC,CAAC,CAAC,CAAC,CAAS,EAAE,EAAE,CAAC,GAAG,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,CAAS,EAAE,EAAE,CAAC,EAAE,CAAC,CAAC,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC;IAClF,IAAI,EAAE,GAAG,KAAK,CAAC;IACf,IAAI,CAAC,MAAM,CAAC,QAAQ,CAAC,EAAE,CAAC,EAAE,CAAC;QACzB,EAAE,GAAG,CAAC,CAAC,CAAC;QACR,KAAK,IAAI,IAAI,GAAG,CAAC,EAAE,CAAC,KAAK,CAAC,EAAE,CAAC,EAAE,IAAI,IAAI,CAAC;YAAE,EAAE,IAAI,IAAI,CAAC;IACvD,CAAC;IACD,IAAI,EAAE,GAAG,MAAM,CAAC,QAAQ,CAAC,KAAK,CAAC,CAAC,CAAC,CAAC,KAAK,CAAC,CAAC,CAAC,EAAE,GAAG,CAAC,CAAC;IACjD,IAAI,CAAC,MAAM,CAAC,QAAQ,CAAC,KAAK,CAAC,EAAE,CAAC;QAC5B,KAAK,IAAI,IAAI,GAAG,CAAC,EAAE,KAAK,CAAC,EAAE,CAAC,EAAE,IAAI,IAAI,CAAC,EAAE,CAAC;YACxC,EAAE,GAAG,EAAE,CAAC;YACR,EAAE,IAAI,IAAI,CAAC;QACb,CAAC;IACH,CAAC;IACD,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,GAAG,EAAE,CAAC,EAAE,EAAE,CAAC;QAC7B,MAAM,GAAG,GAAG,EAAE,GAAG,CAAC,EAAE,GAAG,EAAE,CAAC,GAAG,CAAC,CAAC;QAC/B,IAAI,GAAG,KAAK,EAAE,IAAI,GAAG,KAAK,EAAE;YAAE,MAAM;QACpC,IAAI,KAAK,CAAC,GAAG,CAAC;YAAE,EAAE,GAAG,GAAG,CAAC;;YACpB,EAAE,GAAG,GAAG,CAAC;IAChB,CAAC;IACD,OAAO,EAAE,CAAC;AACZ,CAAC;AAED,yEAAyE;AACzE,SAAS,OAAO,CAAC,CAAS;IACxB,IAAI,CAAC,IAAI,CAAC;QAAE,OAAO,MAAM,CAAC,iBAAiB,CAAC;IAC5C,IAAI,CAAC,IAAI,CAAC;QAAE,OAAO,MAAM,CAAC,iBAAiB,CAAC;IAC5C,MAAM,CAAC,GAAG;QACR,CAAC,mBAAmB,EAAE,mBAAmB,EAAE,CAAC,mBAAmB,EAAE,kBAAkB,EAAE,CAAC,mBAAmB;QACzG,iBAAiB;KAClB,CAAC;IACF,MAAM,CAAC,GAAG;QACR,CAAC,mBAAmB,EAAE,mBAAmB,EAAE,CAAC,mBAAmB,EAAE,mBAAmB,EAAE,CAAC,mBAAmB;KAC3G,CAAC;IACF,MAAM,CAAC,GAAG;QACR,CAAC,oBAAoB,EAAE,CAAC,oBAAoB,EAAE,CAAC,iBAAiB,EAAE,CAAC,iBAAiB,EAAE,iBAAiB;QACvG,iBAAiB;KAClB,CAAC;IACF,MAAM,CAAC,GAAG,CAAC,oBAAoB,EAAE,oBAAoB,EAAE,iBAAiB,EAAE,iBAAiB,CAAC,CAAC;IAC7F,MAAM,IAAI,GAAG,OAAO,CAAC;IACrB,MAAM,KAAK,GAAG,CAAC,GAAG,IAAI,CAAC;IACvB,IAAI,CAAC,GAAG,IAAI,EAAE,CAAC;QACb,MAAM,CAAC,GAAG,IAAI,CAAC,IAAI,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC;QACtC,OAAO,CACL,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC;YACrE,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,CACtD,CAAC;IACJ,CAAC;IACD,IAAI,CAAC,IAAI,KAAK,EAAE,CAAC;QACf,MAAM,CAAC,GAAG,CAAC,GAAG,GAAG,CAAC;QAClB,MAAM,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC;QAChB,OAAO,CACL,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC;YAC3E,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,CACnE,CAAC;IACJ,CAAC;IACD,MAAM,CAAC,GAAG,IAAI,CAAC,IAAI,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC;IAC1C,OAAO,CACL,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC;QACtE,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,CACtD,CAAC;AACJ,CAAC;AAED,MAAM,QAAQ,GAAG,IAAI,CAAC,IAAI,CAAC,CAAC,GAAG,IAAI,CAAC,EAAE,CAAC,CAAC;AAExC,MAAM,OAAO,aAAa;IACxB,2GAA2G;IAC3G,MAAM,CAAC,cAAc,CAAC,CAAS,EAAE,IAAI,GAAG,CAAC,EAAE,EAAE,GAAG,CAAC;QAC/C,IAAI,MAAM,CAAC,KAAK,CAAC,CAAC,CAAC,IAAI,CAAC,GAAG,CAAC,IAAI,CAAC,GAAG,CAAC;YAAE,OAAO,MAAM,CAAC,GAAG,CAAC;QACzD,IAAI,CAAC,GAAG,OAAO,CAAC,CAAC,CAAC,CAAC;QACnB,IAAI,MAAM,CAAC,QAAQ,CAAC,CAAC,CAAC,EAAE,CAAC;YACvB,MAAM,CAAC,GAAG,GAAG,GAAG,EAAE,CAAC,IAAI,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,KAAK,CAAC,GAAG,CAAC,CAAC;YAC7C,MAAM,CAAC,GAAG,CAAC,GAAG,QAAQ,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC;YAC/C,CAAC,IAAI,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC;QAC7B,CAAC;QACD,OAAO,IAAI,GAAG,EAAE,GAAG,CAAC,CAAC;IACvB,CAAC;IAED,sEAAsE;IACtE,MAAM,CAAC,MAAM,CAAC,IAAI,GAAG,CAAC,EAAE,EAAE,GAAG,CAAC,EAAE,MAAW,IAAI,CAAC,MAAM;QACpD,OAAO,UAAU,CAAC;YAChB,GAAG,EAAE,CAAC,CAAC,EAAE,EAAE,CAAC,IAAI,CAAC,GAAG,CAAC,CAAC,GAAG,GAAG,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,GAAG,EAAE,CAAC,IAAI,CAAC,CAAC,GAAG,CAAC,EAAE,GAAG,QAAQ,CAAC;YACrE,GAAG,EAAE,CAAC,CAAC,EAAE,EAAE,CAAC,GAAG,GAAG,EAAE,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,EAAE,GAAG,IAAI,CAAC,KAAK,CAAC,CAAC;YAC1D,EAAE,EAAE,CAAC,CAAC,EAAE,EAAE,CAAC,GAAG,GAAG,EAAE,CAAC,IAAI,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,EAAE,GAAG,IAAI,CAAC,KAAK,CAAC,CAAC;YACxD,GAAG,EAAE,CAAC,CAAC,EAAE,EAAE,CAAC,aAAa,CAAC,cAAc,CAAC,CAAC,EAAE,IAAI,EAAE,EAAE,CAAC;YACrD,IAAI,EAAE,GAAG,EAAE,CAAC,IAAI;YAChB,QAAQ,EAAE,GAAG,EAAE,CAAC,EAAE,GAAG,EAAE;YACvB,MAAM,EAAE,GAAG,EAAE,CAAC,IAAI,GAAG,EAAE,GAAG,IAAI,CAAC,IAAI,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,GAAG,GAAG,EAAE,CAAC,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,GAAG,IAAI,CAAC,EAAE,GAAG,GAAG,EAAE,CAAC;SAC9F,CAAC,CAAC;IACL,CAAC;IAED,8EAA8E;IAC9E,MAAM,CAAC,WAAW,CAAC,IAAI,GAAG,CAAC,EAAE,MAAW,IAAI,CAAC,MAAM;QACjD,OAAO,UAAU,CAAC;YAChB,GAAG,EAAE,CAAC,CAAC,EAAE,EAAE,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,IAAI,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,IAAI,GAAG,CAAC,CAAC,CAAC;YACpD,GAAG,EAAE,CAAC,CAAC,EAAE,EAAE,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,IAAI,CAAC,KAAK,CAAC,CAAC,IAAI,GAAG,CAAC,CAAC,CAAC;YAChD,EAAE,EAAE,CAAC,CAAC,EAAE,EAAE,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,IAAI,CAAC,GAAG,CAAC,CAAC,IAAI,GAAG,CAAC,CAAC,CAAC;YAC5C,GAAG,EAAE,CAAC,CAAC,EAAE,EAAE,CAAC,CAAC,CAAC,GAAG,CAAC,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,MAAM,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,IAAI,CAAC,KAAK,CAAC,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC;YAClE,IAAI,EAAE,GAAG,EAAE,CAAC,CAAC,GAAG,IAAI;YACpB,QAAQ,EAAE,GAAG,EAAE,CAAC,CAAC,GAAG,CAAC,IAAI,GAAG,IAAI,CAAC;YACjC,MAAM,EAAE,GAAG,EAAE,CAAC,CAAC,IAAI,CAAC,GAAG,CAAC,CAAC,GAAG,GAAG,EAAE,CAAC,GAAG,IAAI;SAC1C,CAAC,CAAC;IACL,CAAC;IAED,sEAAsE;IACtE,MAAM,CAAC,OAAO,CAAC,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,EAAE,MAAW,IAAI,CAAC,MAAM;QACjD,OAAO,UAAU,CAAC;YAChB,GAAG,EAAE,CAAC,CAAC,EAAE,EAAE,CAAC,CAAC,CAAC,IAAI,CAAC,IAAI,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;YAChD,GAAG,EAAE,CAAC,CAAC,EAAE,EAAE,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC;YACvD,EAAE,EAAE,CAAC,CAAC,EAAE,EAAE,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC;YACtD,GAAG,EAAE,CAAC,CAAC,EAAE,EAAE,CAAC,CAAC,CAAC,GAAG,CAAC,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,MAAM,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC;YAC3D,IAAI,EAAE,GAAG,EAAE,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC;YACvB,QAAQ,EAAE,GAAG,EAAE,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,IAAI,CAAC,GAAG,EAAE;YACjC,MAAM,EAAE,GAAG,EAAE,CAAC,CAAC,GAAG,GAAG,EAAE,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC;SAClC,CAAC,CAAC;IACL,CAAC;IAED,mFAAmF;IACnF,MAAM,CAAC,KAAK,CAAC,KAAa,EAAE,KAAK,GAAG,CAAC,EAAE,MAAW,IAAI,CAAC,MAAM;QAC3D,OAAO,UAAU,CAAC;YAChB,GAAG,EAAE,CAAC,CAAC,EAAE,EAAE;gBACT,IAAI,CAAC,GAAG,CAAC;oBAAE,OAAO,CAAC,CAAC;gBACpB,IAAI,CAAC,KAAK,CAAC;oBAAE,OAAO,KAAK,GAAG,CAAC,CAAC,CAAC,CAAC,MAAM,CAAC,iBAAiB,CAAC,CAAC,CAAC,KAAK,KAAK,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,KAAK,CAAC,CAAC,CAAC,CAAC,CAAC;gBACvF,OAAO,IAAI,CAAC,GAAG,CAAC,CAAC,KAAK,GAAG,CAAC,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,GAAG,CAAC,GAAG,KAAK,GAAG,EAAE,CAAC,OAAO,CAAC,KAAK,CAAC,GAAG,KAAK,GAAG,IAAI,CAAC,GAAG,CAAC,KAAK,CAAC,CAAC,CAAC;YACvG,CAAC;YACD,GAAG,EAAE,CAAC,CAAC,EAAE,EAAE,CAAC,CAAC,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,EAAE,CAAC,iBAAiB,CAAC,KAAK,EAAE,CAAC,GAAG,KAAK,CAAC,CAAC;YACjE,EAAE,EAAE,CAAC,CAAC,EAAE,EAAE,CAAC,CAAC,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,EAAE,CAAC,iBAAiB,CAAC,KAAK,EAAE,CAAC,GAAG,KAAK,CAAC,CAAC;YAChE,KAAK,EAAE,CAAC;YACR,IAAI,EAAE,GAAG,EAAE,CAAC,KAAK,GAAG,KAAK;YACzB,QAAQ,EAAE,GAAG,EAAE,CAAC,KAAK,GAAG,KAAK,GAAG,KAAK;YACrC,MAAM,EAAE,GAAG,EAAE,CAAC,aAAa,CAAC,WAAW,CAAC,KAAK,EAAE,GAAG,CAAC,GAAG,KAAK;SAC5D,CAAC,CAAC;IACL,CAAC;IAED,yEAAyE;IACzE,MAAM,CAAC,SAAS,CAAC,EAAU,EAAE,MAAW,IAAI,CAAC,MAAM;QACjD,MAAM,CAAC,GAAG,aAAa,CAAC,KAAK,CAAC,EAAE,GAAG,CAAC,EAAE,CAAC,EAAE,GAAG,CAAC,CAAC;QAC9C,OAAO,EAAE,GAAG,CAAC,EAAE,IAAI,EAAE,GAAG,EAAE,CAAC,EAAE,EAAE,QAAQ,EAAE,GAAG,EAAE,CAAC,CAAC,GAAG,EAAE,EAAE,GAAG,EAAE,GAAG,EAAE,CAAC,IAAI,CAAC,IAAI,CAAC,CAAC,GAAG,EAAE,CAAC,EAAE,CAAC;IACxF,CAAC;IAED,uEAAuE;IACvE,MAAM,CAAC,QAAQ,CAAC,EAAU,EAAE,MAAW,IAAI,CAAC,MAAM;QAChD,MAAM,EAAE,GAAG,EAAE,CAAC,OAAO,CAAC,CAAC,EAAE,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,EAAE,CAAC,OAAO,CAAC,EAAE,GAAG,CAAC,CAAC,GAAG,GAAG,GAAG,IAAI,CAAC,GAAG,CAAC,EAAE,GAAG,IAAI,CAAC,EAAE,CAAC,CAAC;QACxF,MAAM,GAAG,GAAG,CAAC,CAAS,EAAU,EAAE;YAChC,IAAI,CAAC,KAAK,MAAM,CAAC,iBAAiB;gBAAE,OAAO,CAAC,CAAC;YAC7C,IAAI,CAAC,KAAK,MAAM,CAAC,iBAAiB;gBAAE,OAAO,CAAC,CAAC;YAC7C,MAAM,EAAE,GAAG,EAAE,CAAC,yBAAyB,CAAC,EAAE,GAAG,CAAC,EAAE,GAAG,CAAC,GAAG,CAAC,CAAC,EAAE,EAAE,GAAG,CAAC,EAAE,GAAG,CAAC,CAAC;YACxE,OAAO,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,GAAG,GAAG,EAAE,CAAC,CAAC,CAAC,GAAG,GAAG,EAAE,CAAC;QAC1C,CAAC,CAAC;QACF,OAAO,UAAU,CAAC;YAChB,GAAG,EAAE,CAAC,CAAC,EAAE,EAAE,CAAC,IAAI,CAAC,GAAG,CAAC,EAAE,GAAG,CAAC,CAAC,EAAE,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,IAAI,CAAC,KAAK,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,EAAE,CAAC,CAAC;YACpE,GAAG;YACH,EAAE,EAAE,CAAC,CAAC,EAAE,EAAE,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC;YAClB,IAAI,EAAE,GAAG,EAAE,CAAC,CAAC,EAAE,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,MAAM,CAAC,GAAG,CAAC;YACrC,QAAQ,EAAE,GAAG,EAAE,CAAC,CAAC,EAAE,GAAG,CAAC,CAAC,CAAC,CAAC,EAAE,GAAG,CAAC,EAAE,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,EAAE,GAAG,CAAC,CAAC,CAAC,CAAC,MAAM,CAAC,iBAAiB,CAAC,CAAC,CAAC,MAAM,CAAC,GAAG,CAAC;YACzF,MAAM,EAAE,GAAG,EAAE;gBACX,MAAM,CAAC,GAAG,aAAa,CAAC,MAAM,CAAC,CAAC,EAAE,CAAC,EAAE,GAAG,CAAC,CAAC,MAAM,EAAE,CAAC;gBACnD,MAAM,GAAG,GAAG,aAAa,CAAC,SAAS,CAAC,EAAE,EAAE,GAAG,CAAC,CAAC,MAAM,EAAE,CAAC;gBACtD,OAAO,CAAC,GAAG,IAAI,CAAC,IAAI,CAAC,GAAG,GAAG,EAAE,CAAC,CAAC;YACjC,CAAC;SACF,CAAC,CAAC;IACL,CAAC;IAED,gHAAgH;IAChH,MAAM,CAAC,OAAO,CAAC,EAAU,EAAE,EAAU,EAAE,MAAW,IAAI,CAAC,MAAM;QAC3D,MAAM,KAAK,GAAG,EAAE,CAAC,OAAO,CAAC,EAAE,GAAG,CAAC,CAAC,GAAG,EAAE,CAAC,OAAO,CAAC,EAAE,GAAG,CAAC,CAAC,GAAG,EAAE,CAAC,OAAO,CAAC,CAAC,EAAE,GAAG,EAAE,CAAC,GAAG,CAAC,CAAC,CAAC;QAClF,OAAO,UAAU,CAAC;YAChB,GAAG,EAAE,CAAC,CAAC,EAAE,EAAE;gBACT,IAAI,CAAC,GAAG,CAAC;oBAAE,OAAO,CAAC,CAAC;gBACpB,IAAI,CAAC,KAAK,CAAC;oBAAE,OAAO,EAAE,GAAG,CAAC,CAAC,CAAC,CAAC,MAAM,CAAC,iBAAiB,CAAC,CAAC,CAAC,EAAE,KAAK,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;gBACzE,OAAO,IAAI,CAAC,GAAG,CACb,GAAG,GAAG,EAAE,GAAG,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC;oBACrB,GAAG,GAAG,EAAE,GAAG,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC;oBACvB,CAAC,EAAE,GAAG,CAAC,GAAG,CAAC,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC;oBAC1B,CAAC,CAAC,EAAE,GAAG,EAAE,CAAC,GAAG,CAAC,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,EAAE,GAAG,EAAE,GAAG,CAAC,CAAC;oBACvC,KAAK,CACR,CAAC;YACJ,CAAC;YACD,GAAG,EAAE,CAAC,CAAC,EAAE,EAAE,CAAC,CAAC,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,EAAE,CAAC,yBAAyB,CAAC,CAAC,EAAE,GAAG,CAAC,CAAC,GAAG,CAAC,EAAE,GAAG,CAAC,GAAG,EAAE,CAAC,EAAE,EAAE,GAAG,CAAC,EAAE,EAAE,GAAG,CAAC,CAAC,CAAC;YACjG,EAAE,EAAE,CAAC,CAAC,EAAE,EAAE,CAAC,CAAC,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,EAAE,CAAC,yBAAyB,CAAC,EAAE,GAAG,CAAC,EAAE,GAAG,EAAE,GAAG,CAAC,CAAC,EAAE,EAAE,GAAG,CAAC,EAAE,EAAE,GAAG,CAAC,CAAC,CAAC;YAC1F,KAAK,EAAE,CAAC;YACR,IAAI,EAAE,GAAG,EAAE,CAAC,CAAC,EAAE,GAAG,CAAC,CAAC,CAAC,CAAC,EAAE,GAAG,CAAC,EAAE,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,MAAM,CAAC,GAAG,CAAC;YACjD,QAAQ,EAAE,GAAG,EAAE,CACb,EAAE,GAAG,CAAC;gBACJ,CAAC,CAAC,CAAC,CAAC,GAAG,EAAE,GAAG,EAAE,GAAG,CAAC,EAAE,GAAG,EAAE,GAAG,CAAC,CAAC,CAAC,GAAG,CAAC,EAAE,GAAG,CAAC,EAAE,GAAG,CAAC,CAAC,IAAI,CAAC,GAAG,CAAC,EAAE,GAAG,CAAC,CAAC,CAAC;gBACjE,CAAC,CAAC,EAAE,GAAG,CAAC;oBACN,CAAC,CAAC,MAAM,CAAC,iBAAiB;oBAC1B,CAAC,CAAC,MAAM,CAAC,GAAG;YAClB,MAAM,EAAE,GAAG,EAAE,CACX,CAAC,aAAa,CAAC,WAAW,CAAC,EAAE,GAAG,CAAC,EAAE,GAAG,CAAC,GAAG,CAAC,CAAC,GAAG,EAAE,GAAG,CAAC,CAAC,aAAa,CAAC,WAAW,CAAC,EAAE,GAAG,CAAC,EAAE,GAAG,CAAC,GAAG,CAAC,CAAC,GAAG,EAAE,CAAC;SAC1G,CAAC,CAAC;IACL,CAAC;IAED,uDAAuD;IACvD,MAAM,CAAC,QAAQ,CAAC,CAAS,EAAE,CAAS,EAAE,MAAW,IAAI,CAAC,MAAM;QAC1D,MAAM,SAAS,GAAG,CAAC,CAAS,EAAE,EAAE,CAAC,EAAE,CAAC,OAAO,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,EAAE,CAAC,OAAO,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,EAAE,CAAC,OAAO,CAAC,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,CAAC;QAC/F,MAAM,GAAG,GAAG,CAAC,CAAS,EAAU,EAAE;YAChC,IAAI,CAAC,GAAG,CAAC,IAAI,CAAC,GAAG,CAAC,IAAI,CAAC,MAAM,CAAC,SAAS,CAAC,CAAC,CAAC;gBAAE,OAAO,CAAC,CAAC;YACrD,IAAI,CAAC,KAAK,CAAC;gBAAE,OAAO,CAAC,KAAK,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;YACpC,IAAI,CAAC,KAAK,CAAC;gBAAE,OAAO,CAAC,KAAK,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;YACpC,OAAO,IAAI,CAAC,GAAG,CAAC,SAAS,CAAC,CAAC,CAAC,GAAG,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,IAAI,CAAC,KAAK,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;QAC7E,CAAC,CAAC;QACF,OAAO,QAAQ,CAAC;YACd,GAAG;YACH,GAAG,EAAE,CAAC,CAAC,EAAE,EAAE;gBACT,IAAI,CAAC,GAAG,CAAC;oBAAE,OAAO,CAAC,CAAC;gBACpB,MAAM,CAAC,GAAG,IAAI,CAAC,KAAK,CAAC,CAAC,CAAC,CAAC;gBACxB,OAAO,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,EAAE,CAAC,yBAAyB,CAAC,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,CAAC,CAAC;YACxE,CAAC;YACD,EAAE,EAAE,CAAC,CAAC,EAAE,EAAE;gBACR,IAAI,CAAC,GAAG,CAAC;oBAAE,OAAO,CAAC,CAAC;gBACpB,MAAM,CAAC,GAAG,IAAI,CAAC,KAAK,CAAC,CAAC,CAAC,CAAC;gBACxB,OAAO,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,EAAE,CAAC,yBAAyB,CAAC,CAAC,EAAE,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,CAAC,CAAC;YACpE,CAAC;YACD,GAAG,EAAE,CAAC;YACN,IAAI,EAAE,GAAG,EAAE,CAAC,CAAC,GAAG,CAAC;YACjB,QAAQ,EAAE,GAAG,EAAE,CAAC,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC;YAC/B,MAAM,EAAE,GAAG,EAAE;gBACX,IAAI,KAAK,GAAG,CAAC,CAAC;gBACd,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,EAAE,CAAC,EAAE;oBAAE,IAAI,GAAG,EAAE,GAAG,CAAC;wBAAE,KAAK,EAAE,CAAC;gBACnD,OAAO,KAAK,CAAC;YACf,CAAC;SACF,CAAC,CAAC;IACL,CAAC;IAED,6DAA6D;IAC7D,MAAM,CAAC,OAAO,CAAC,MAAc,EAAE,MAAW,IAAI,CAAC,MAAM;QACnD,OAAO,QAAQ,CAAC;YACd,GAAG,EAAE,CAAC,CAAC,EAAE,EAAE,CAAC,CAAC,CAAC,GAAG,CAAC,IAAI,CAAC,MAAM,CAAC,SAAS,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,IAAI,CAAC,GAAG,CAAC,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,MAAM,CAAC,GAAG,MAAM,GAAG,EAAE,CAAC,OAAO,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC;YAC7G,GAAG,EAAE,CAAC,CAAC,EAAE,EAAE,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,EAAE,CAAC,iBAAiB,CAAC,IAAI,CAAC,KAAK,CAAC,CAAC,CAAC,GAAG,CAAC,EAAE,MAAM,CAAC,CAAC;YACzE,EAAE,EAAE,CAAC,CAAC,EAAE,EAAE,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,EAAE,CAAC,iBAAiB,CAAC,IAAI,CAAC,KAAK,CAAC,CAAC,CAAC,GAAG,CAAC,EAAE,MAAM,CAAC,CAAC;YACxE,IAAI,EAAE,GAAG,EAAE,CAAC,MAAM;YAClB,QAAQ,EAAE,GAAG,EAAE,CAAC,MAAM;YACtB,MAAM,EAAE,GAAG,EAAE;gBACX,qEAAqE;gBACrE,mEAAmE;gBACnE,IAAI,SAAS,GAAG,MAAM,CAAC;gBACvB,IAAI,KAAK,GAAG,CAAC,CAAC;gBACd,OAAO,SAAS,GAAG,CAAC,EAAE,CAAC;oBACrB,MAAM,IAAI,GAAG,IAAI,CAAC,GAAG,CAAC,SAAS,EAAE,GAAG,CAAC,CAAC;oBACtC,SAAS,IAAI,IAAI,CAAC;oBAClB,MAAM,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,IAAI,CAAC,CAAC;oBAC1B,IAAI,CAAC,GAAG,CAAC,CAAC;oBACV,IAAI,IAAI,GAAG,CAAC,CAAC;oBACb,GAAG,CAAC;wBACF,CAAC,EAAE,CAAC;wBACJ,IAAI,IAAI,GAAG,EAAE,CAAC;oBAChB,CAAC,QAAQ,IAAI,GAAG,CAAC,EAAE;oBACnB,KAAK,IAAI,CAAC,GAAG,CAAC,CAAC;gBACjB,CAAC;gBACD,OAAO,KAAK,CAAC;YACf,CAAC;SACF,CAAC,CAAC;IACL,CAAC;IAED,iDAAiD;IACzC,MAAM,CAAC,WAAW,CAAC,KAAa,EAAE,MAAW,IAAI,CAAC,MAAM;QAC9D,IAAI,KAAK,GAAG,CAAC;YAAE,OAAO,aAAa,CAAC,WAAW,CAAC,KAAK,GAAG,CAAC,EAAE,GAAG,CAAC,GAAG,GAAG,EAAE,IAAI,CAAC,CAAC,GAAG,KAAK,CAAC,CAAC;QACvF,MAAM,CAAC,GAAG,KAAK,GAAG,CAAC,GAAG,CAAC,CAAC;QACxB,MAAM,CAAC,GAAG,CAAC,GAAG,IAAI,CAAC,IAAI,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC;QAC/B,SAAS,CAAC;YACR,IAAI,CAAS,CAAC;YACd,IAAI,CAAS,CAAC;YACd,GAAG,CAAC;gBACF,CAAC,GAAG,aAAa,CAAC,MAAM,CAAC,CAAC,EAAE,CAAC,EAAE,GAAG,CAAC,CAAC,MAAM,EAAE,CAAC;gBAC7C,CAAC,GAAG,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC;YAChB,CAAC,QAAQ,CAAC,IAAI,CAAC,EAAE;YACjB,CAAC,GAAG,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC;YACd,MAAM,CAAC,GAAG,GAAG,EAAE,CAAC;YAChB,IAAI,CAAC,GAAG,CAAC,GAAG,MAAM,GAAG,CAAC,IAAI,CAAC;gBAAE,OAAO,CAAC,GAAG,CAAC,CAAC;YAC1C,IAAI,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,GAAG,GAAG,GAAG,CAAC,GAAG,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC;gBAAE,OAAO,CAAC,GAAG,CAAC,CAAC;QAC1E,CAAC;IACH,CAAC;CACF;AA0BD,MAAM,KAAK,GAAG,CAAC,EAAqB,EAAE,EAAE,CAAC,EAAE,CAAC,MAAM,CAAC,CAAC,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC,CAAC,GAAG,CAAC,EAAE,CAAC,CAAC,CAAC;AACvE,MAAM,MAAM,GAAG,CAAC,EAAqB,EAAE,EAAE,CAAC,KAAK,CAAC,EAAE,CAAC,GAAG,EAAE,CAAC,MAAM,CAAC;AAChE,MAAM,SAAS,GAAG,CAAC,EAAqB,EAAE,EAAE;IAC1C,MAAM,CAAC,GAAG,MAAM,CAAC,EAAE,CAAC,CAAC;IACrB,OAAO,EAAE,CAAC,MAAM,CAAC,CAAC,GAAG,EAAE,CAAC,EAAE,EAAE,CAAC,GAAG,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,IAAI,CAAC,EAAE,CAAC,CAAC,GAAG,CAAC,EAAE,CAAC,MAAM,GAAG,CAAC,CAAC,CAAC;AACxE,CAAC,CAAC;AAEF,6DAA6D;AAC7D,SAAS,OAAO,CAAC,CAAS,EAAE,EAAU,EAAE,cAA2B,WAAW;IAC5E,MAAM,IAAI,GAAG,aAAa,CAAC,QAAQ,CAAC,EAAE,CAAC,CAAC;IACxC,IAAI,WAAW,KAAK,MAAM;QAAE,OAAO,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC;IAC/C,IAAI,WAAW,KAAK,SAAS;QAAE,OAAO,IAAI,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC;IACjD,OAAO,IAAI,CAAC,GAAG,CAAC,CAAC,EAAE,CAAC,GAAG,IAAI,CAAC,EAAE,CAAC,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;AAC/C,CAAC;AAED,MAAM,OAAO,eAAe;IAC1B,gGAAgG;IAChG,MAAM,CAAC,cAAc,CAAC,MAAgB,EAAE,GAAG,GAAG,CAAC,EAAE,UAAwB,EAAE;QACzE,MAAM,CAAC,GAAG,MAAM,CAAC,MAAM,CAAC;QACxB,MAAM,EAAE,GAAG,IAAI,CAAC,IAAI,CAAC,SAAS,CAAC,MAAM,CAAC,GAAG,CAAC,CAAC,CAAC;QAC5C,MAAM,CAAC,GAAG,CAAC,MAAM,CAAC,MAAM,CAAC,GAAG,GAAG,CAAC,GAAG,EAAE,CAAC;QACtC,MAAM,EAAE,GAAG,CAAC,GAAG,CAAC,CAAC;QACjB,OAAO,EAAE,SAAS,EAAE,CAAC,EAAE,EAAE,EAAE,MAAM,EAAE,OAAO,CAAC,CAAC,EAAE,EAAE,EAAE,OAAO,CAAC,WAAW,CAAC,EAAE,CAAC;IAC3E,CAAC;IAED;;;OAGG;IACH,MAAM,CAAC,cAAc,CAAC,CAAW,EAAE,CAAW,EAAE,UAAwB,EAAE;QACxE,OAAO,eAAe,CAAC,gBAAgB,CAAC,CAAC,EAAE,CAAC,EAAE,EAAE,GAAG,OAAO,EAAE,QAAQ,EAAE,KAAK,EAAE,CAAC,CAAC;IACjF,CAAC;IAED,gDAAgD;IAChD,MAAM,CAAC,KAAK,CAAC,CAAW,EAAE,CAAW,EAAE,UAAwB,EAAE;QAC/D,OAAO,eAAe,CAAC,cAAc,CAAC,CAAC,EAAE,CAAC,EAAE,OAAO,CAAC,CAAC;IACvD,CAAC;IAED,2HAA2H;IAC3H,MAAM,CAAC,gBAAgB,CAAC,CAAW,EAAE,CAAW,EAAE,UAAmC,EAAE;QACrF,MAAM,EAAE,QAAQ,GAAG,IAAI,EAAE,WAAW,EAAE,GAAG,OAAO,CAAC;QACjD,MAAM,EAAE,GAAG,CAAC,CAAC,MAAM,CAAC;QACpB,MAAM,EAAE,GAAG,CAAC,CAAC,MAAM,CAAC;QACpB,MAAM,EAAE,GAAG,SAAS,CAAC,CAAC,CAAC,CAAC;QACxB,MAAM,EAAE,GAAG,SAAS,CAAC,CAAC,CAAC,CAAC;QACxB,IAAI,CAAS,CAAC;QACd,IAAI,EAAU,CAAC;QACf,IAAI,QAAQ,EAAE,CAAC;YACb,EAAE,GAAG,EAAE,GAAG,EAAE,GAAG,CAAC,CAAC;YACjB,MAAM,MAAM,GAAG,CAAC,CAAC,EAAE,GAAG,CAAC,CAAC,GAAG,EAAE,GAAG,CAAC,EAAE,GAAG,CAAC,CAAC,GAAG,EAAE,CAAC,GAAG,EAAE,CAAC;YACpD,CAAC,GAAG,CAAC,MAAM,CAAC,CAAC,CAAC,GAAG,MAAM,CAAC,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,IAAI,CAAC,MAAM,GAAG,CAAC,CAAC,GAAG,EAAE,GAAG,CAAC,GAAG,EAAE,CAAC,CAAC,CAAC;QACtE,CAAC;aAAM,CAAC;YACN,MAAM,EAAE,GAAG,EAAE,GAAG,EAAE,CAAC;YACnB,MAAM,EAAE,GAAG,EAAE,GAAG,EAAE,CAAC;YACnB,CAAC,GAAG,CAAC,MAAM,CAAC,CAAC,CAAC,GAAG,MAAM,CAAC,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,IAAI,CAAC,EAAE,GAAG,EAAE,CAAC,CAAC;YACjD,EAAE,GAAG,CAAC,EAAE,GAAG,EAAE,CAAC,IAAI,CAAC,GAAG,CAAC,EAAE,IAAI,CAAC,GAAG,CAAC,EAAE,GAAG,CAAC,CAAC,GAAG,EAAE,IAAI,CAAC,GAAG,CAAC,EAAE,GAAG,CAAC,CAAC,CAAC,CAAC;QAClE,CAAC;QACD,OAAO,EAAE,SAAS,EAAE,CAAC,EAAE,EAAE,EAAE,MAAM,EAAE,OAAO,CAAC,CAAC,EAAE,EAAE,EAAE,WAAW,CAAC,EAAE,CAAC;IACnE,CAAC;IAED,mEAAmE;IACnE,MAAM,CAAC,WAAW,CAAC,CAAW,EAAE,CAAW,EAAE,UAAwB,EAAE;QACrE,IAAI,CAAC,CAAC,MAAM,KAAK,CAAC,CAAC,MAAM;YAAE,MAAM,IAAI,KAAK,CAAC,+CAA+C,CAAC,CAAC;QAC5F,OAAO,eAAe,CAAC,cAAc,CACnC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,EACzB,CAAC,EACD,OAAO,CACR,CAAC;IACJ,CAAC;IAED;;;;OAIG;IACH,MAAM,CAAC,sBAAsB,CAAC,QAAkB,EAAE,QAAmB,EAAE,IAAI,GAAG,CAAC;QAC7E,MAAM,CAAC,GAAG,QAAQ,CAAC,MAAM,CAAC;QAC1B,MAAM,GAAG,GAAG,QAAQ,IAAI,IAAI,KAAK,CAAS,CAAC,CAAC,CAAC,IAAI,CAAC,KAAK,CAAC,QAAQ,CAAC,GAAG,CAAC,CAAC,CAAC;QACvE,IAAI,GAAG,CAAC,MAAM,KAAK,CAAC;YAAE,MAAM,IAAI,KAAK,CAAC,kDAAkD,CAAC,CAAC;QAC1F,MAAM,EAAE,GAAG,KAAK,CAAC,QAAQ,CAAC,CAAC;QAC3B,MAAM,EAAE,GAAG,KAAK,CAAC,GAAG,CAAC,CAAC;QACtB,IAAI,IAAI,CAAC,GAAG,CAAC,EAAE,GAAG,EAAE,CAAC,GAAG,IAAI,GAAG,IAAI,CAAC,GAAG,CAAC,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,EAAE,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,CAAC,EAAE,CAAC;YACpE,MAAM,IAAI,KAAK,CAAC,aAAa,EAAE,mBAAmB,EAAE,2CAA2C,CAAC,CAAC;QACnG,CAAC;QACD,MAAM,IAAI,GAAG,QAAQ,CAAC,MAAM,CAAC,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,GAAG,CAAC,CAAC,CAAC,CAAC,IAAI,CAAC,GAAG,GAAG,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC;QAC7E,MAAM,EAAE,GAAG,CAAC,GAAG,CAAC,GAAG,IAAI,CAAC;QACxB,OAAO,EAAE,SAAS,EAAE,IAAI,EAAE,EAAE,EAAE,MAAM,EAAE,aAAa,CAAC,SAAS,CAAC,EAAE,CAAC,CAAC,EAAE,CAAC,IAAI,CAAC,EAAE,CAAC;IAC/E,CAAC;IAED;;;;OAIG;IACH,MAAM,CAAC,qBAAqB,CAAC,KAAiB,EAAE,UAAoC,EAAE;QACpF,MAAM,EAAE,UAAU,GAAG,IAAI,EAAE,GAAG,OAAO,CAAC;QACtC,MAAM,CAAC,GAAG,KAAK,CAAC,MAAM,CAAC;QACvB,MAAM,CAAC,GAAG,KAAK,CAAC,CAAC,CAAC,CAAC,MAAM,CAAC;QAC1B,MAAM,IAAI,GAAG,KAAK,CAAC,GAAG,CAAC,KAAK,CAAC,CAAC;QAC9B,MAAM,IAAI,GAAG,KAAK,CAAC,IAAI,CAAC,EAAE,MAAM,EAAE,CAAC,EAAE,EAAE,CAAC,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC,KAAK,CAAC,KAAK,CAAC,GAAG,CAAC,CAAC,GAAG,EAAE,EAAE,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;QACpF,MAAM,KAAK,GAAG,KAAK,CAAC,IAAI,CAAC,CAAC;QAC1B,MAAM,QAAQ,GAAG,KAAK,CAAC,GAAG,CAAC,CAAC,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC,IAAI,CAAC,GAAG,CAAC,CAAC,EAAE,EAAE,EAAE,CAAC,CAAC,IAAI,CAAC,CAAC,CAAC,GAAG,EAAE,CAAC,GAAG,KAAK,CAAC,CAAC,CAAC;QAC/E,MAAM,EAAE,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC;QAC7B,IAAI,IAAI,GAAG,CAAC,CAAC;QACb,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC;YAC3B,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC;gBAC3B,IAAI,IAAI,GAAG,IAAI,CAAC,GAAG,CAAC,KAAK,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,QAAQ,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;gBAClD,IAAI,UAAU,IAAI,EAAE,KAAK,CAAC;oBAAE,IAAI,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,EAAE,IAAI,GAAG,GAAG,CAAC,CAAC;gBAC3D,IAAI,IAAI,IAAI,IAAI,CAAC,GAAG,QAAQ,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;YACrC,CAAC;QACH,CAAC;QACD,OAAO,EAAE,SAAS,EAAE,IAAI,EAAE,EAAE,EAAE,MAAM,EAAE,EAAE,KAAK,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,aAAa,CAAC,SAAS,CAAC,EAAE,CAAC,CAAC,EAAE,CAAC,IAAI,CAAC,EAAE,QAAQ,EAAE,CAAC;IACxG,CAAC;IAED,4DAA4D;IAC5D,MAAM,CAAC,sBAAsB,CAAC,MAAgB,EAAE,KAAK,GAAG,IAAI;QAC1D,MAAM,CAAC,GAAG,MAAM,CAAC,MAAM,CAAC;QACxB,MAAM,CAAC,GAAG,MAAM,CAAC,MAAM,CAAC,CAAC;QACzB,MAAM,EAAE,GAAG,IAAI,CAAC,IAAI,CAAC,SAAS,CAAC,MAAM,CAAC,GAAG,CAAC,CAAC,CAAC;QAC5C,MAAM,KAAK,GAAG,aAAa,CAAC,QAAQ,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,KAAK,CAAC,GAAG,CAAC,CAAC,CAAC;QACrE,OAAO,CAAC,CAAC,GAAG,KAAK,GAAG,EAAE,EAAE,CAAC,GAAG,KAAK,GAAG,EAAE,CAAC,CAAC;IAC1C,CAAC;IAED,yDAAyD;IACzD,MAAM,CAAC,4BAA4B,CAAC,CAAW,EAAE,CAAW,EAAE,KAAK,GAAG,IAAI;QACxE,MAAM,EAAE,GAAG,SAAS,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,MAAM,CAAC;QACnC,MAAM,EAAE,GAAG,SAAS,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,MAAM,CAAC;QACnC,MAAM,EAAE,GAAG,CAAC,EAAE,GAAG,EAAE,CAAC,IAAI,CAAC,GAAG,CAAC,EAAE,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,MAAM,GAAG,CAAC,CAAC,GAAG,EAAE,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,MAAM,GAAG,CAAC,CAAC,CAAC,CAAC;QAClF,MAAM,CAAC,GAAG,MAAM,CAAC,CAAC,CAAC,GAAG,MAAM,CAAC,CAAC,CAAC,CAAC;QAChC,MAAM,IAAI,GAAG,aAAa,CAAC,QAAQ,CAAC,EAAE,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,KAAK,CAAC,GAAG,CAAC,CAAC,GAAG,IAAI,CAAC,IAAI,CAAC,EAAE,GAAG,EAAE,CAAC,CAAC;QACtF,OAAO,CAAC,CAAC,GAAG,IAAI,EAAE,CAAC,GAAG,IAAI,CAAC,CAAC;IAC9B,CAAC;IAED,wJAAwJ;IACxJ,MAAM,CAAC,4BAA4B,CAAC,SAAiB,EAAE,CAAS,EAAE,KAAK,GAAG,IAAI;QAC5E,MAAM,CAAC,GAAG,aAAa,CAAC,cAAc,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,KAAK,CAAC,GAAG,CAAC,CAAC,CAAC;QAC5D,MAAM,IAAI,GAAG,SAAS,GAAG,CAAC,CAAC;QAC3B,MAAM,KAAK,GAAG,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC;QAC9B,MAAM,MAAM,GAAG,CAAC,IAAI,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,GAAG,KAAK,CAAC;QAClD,MAAM,IAAI,GAAG,CAAC,CAAC,GAAG,IAAI,CAAC,IAAI,CAAC,CAAC,IAAI,GAAG,CAAC,CAAC,GAAG,IAAI,CAAC,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,GAAG,KAAK,CAAC;QACtF,OAAO,CAAC,IAAI,CAAC,GAAG,CAAC,CAAC,EAAE,MAAM,GAAG,IAAI,CAAC,EAAE,IAAI,CAAC,GAAG,CAAC,CAAC,EAAE,MAAM,GAAG,IAAI,CAAC,CAAC,CAAC;IAClE,CAAC;CACF"}
|
|
@@ -0,0 +1,22 @@
|
|
|
1
|
+
import { ComplexNumber } from "./ComplexNumber.ts";
|
|
2
|
+
/**
|
|
3
|
+
* General (non-symmetric) eigen solver: Householder reduction to upper
|
|
4
|
+
* Hessenberg form, then single-shift QR iteration (Wilkinson shifts, Givens
|
|
5
|
+
* rotations, deflation) in complex arithmetic until the matrix is the complex
|
|
6
|
+
* Schur form `A = Z·T·Zᴴ`; eigenvectors come from back-substitution on the
|
|
7
|
+
* triangular `T`. Real input with conjugate eigenvalue pairs works without a
|
|
8
|
+
* separate real-Schur code path, at the cost of doing the arithmetic in ℂ.
|
|
9
|
+
*
|
|
10
|
+
* Mirrors `scipy.linalg.eig` / `numpy.linalg.eig` (right eigenvectors,
|
|
11
|
+
* columns of `vectors`, unit 2-norm). Eigenvalues are sorted by descending
|
|
12
|
+
* real part, then descending imaginary part; each eigenvector's largest
|
|
13
|
+
* component is real and positive (LAPACK leaves that sign arbitrary).
|
|
14
|
+
*/
|
|
15
|
+
export interface GeneralEigenResult {
|
|
16
|
+
/** Eigenvalues, `values[k]` pairs with column `k` of `vectors`. */
|
|
17
|
+
values: ComplexNumber[];
|
|
18
|
+
/** `n×n` matrix whose column `k` is the unit-norm eigenvector for `values[k]`. */
|
|
19
|
+
vectors: ComplexNumber[][];
|
|
20
|
+
}
|
|
21
|
+
export declare function eigenGeneral(input: ReadonlyArray<ReadonlyArray<number | ComplexNumber>>, maxIterations?: number): GeneralEigenResult;
|
|
22
|
+
//# sourceMappingURL=EigenGeneral.d.ts.map
|
|
@@ -0,0 +1 @@
|
|
|
1
|
+
{"version":3,"file":"EigenGeneral.d.ts","sourceRoot":"","sources":["../src/EigenGeneral.ts"],"names":[],"mappings":"AAAA,OAAO,EAAE,aAAa,EAAE,MAAM,oBAAoB,CAAC;AAEnD;;;;;;;;;;;;GAYG;AAEH,MAAM,WAAW,kBAAkB;IACjC,mEAAmE;IACnE,MAAM,EAAE,aAAa,EAAE,CAAC;IACxB,kFAAkF;IAClF,OAAO,EAAE,aAAa,EAAE,EAAE,CAAC;CAC5B;AAMD,wBAAgB,YAAY,CAC1B,KAAK,EAAE,aAAa,CAAC,aAAa,CAAC,MAAM,GAAG,aAAa,CAAC,CAAC,EAC3D,aAAa,SAAK,GACjB,kBAAkB,CAsGpB"}
|