@johnhenry/math 0.2.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.
Files changed (52) hide show
  1. package/README.md +4 -3
  2. package/dist/Distributions.d.ts +86 -16
  3. package/dist/Distributions.d.ts.map +1 -1
  4. package/dist/Distributions.js +340 -89
  5. package/dist/Distributions.js.map +1 -1
  6. package/dist/EigenGeneral.d.ts +22 -0
  7. package/dist/EigenGeneral.d.ts.map +1 -0
  8. package/dist/EigenGeneral.js +327 -0
  9. package/dist/EigenGeneral.js.map +1 -0
  10. package/dist/Frobenius.d.ts +60 -0
  11. package/dist/Frobenius.d.ts.map +1 -0
  12. package/dist/Frobenius.js +162 -0
  13. package/dist/Frobenius.js.map +1 -0
  14. package/dist/Graph.d.ts +26 -0
  15. package/dist/Graph.d.ts.map +1 -1
  16. package/dist/Graph.js +75 -0
  17. package/dist/Graph.js.map +1 -1
  18. package/dist/MatrixMath.d.ts +10 -0
  19. package/dist/MatrixMath.d.ts.map +1 -1
  20. package/dist/MatrixMath.js +11 -0
  21. package/dist/MatrixMath.js.map +1 -1
  22. package/dist/MonteCarlo.d.ts +106 -0
  23. package/dist/MonteCarlo.d.ts.map +1 -0
  24. package/dist/MonteCarlo.js +183 -0
  25. package/dist/MonteCarlo.js.map +1 -0
  26. package/dist/Random.d.ts +31 -0
  27. package/dist/Random.d.ts.map +1 -0
  28. package/dist/Random.js +95 -0
  29. package/dist/Random.js.map +1 -0
  30. package/dist/Regression.d.ts +114 -0
  31. package/dist/Regression.d.ts.map +1 -0
  32. package/dist/Regression.js +338 -0
  33. package/dist/Regression.js.map +1 -0
  34. package/dist/SpecialFunctions.d.ts +3 -3
  35. package/dist/SpecialFunctions.d.ts.map +1 -1
  36. package/dist/SpecialFunctions.js +61 -51
  37. package/dist/SpecialFunctions.js.map +1 -1
  38. package/dist/index.d.ts +6 -1
  39. package/dist/index.d.ts.map +1 -1
  40. package/dist/index.js +4 -0
  41. package/dist/index.js.map +1 -1
  42. package/package.json +1 -1
  43. package/src/Distributions.ts +391 -94
  44. package/src/EigenGeneral.ts +342 -0
  45. package/src/Frobenius.ts +203 -0
  46. package/src/Graph.ts +68 -0
  47. package/src/MatrixMath.ts +13 -0
  48. package/src/MonteCarlo.ts +257 -0
  49. package/src/Random.ts +104 -0
  50. package/src/Regression.ts +423 -0
  51. package/src/SpecialFunctions.ts +52 -48
  52. package/src/index.ts +33 -0
@@ -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
- return mean + sd * normInv(p);
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 * Math.sqrt(2 * Math.PI)),
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 : 1 - Math.exp(-rate * x)),
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) => (x <= 0 ? 0 : (x ** (shape - 1) * Math.exp(-x / scale)) / (SF.gamma(shape) * scale ** shape)),
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 c = SF.gamma((df + 1) / 2) / (Math.sqrt(df * Math.PI) * SF.gamma(df / 2));
84
- return {
85
- pdf: (x) => c * (1 + (x * x) / df) ** (-(df + 1) / 2),
86
- cdf: (x) => {
87
- const ib = SF.regularizedIncompleteBeta(df / (df + x * x), df / 2, 0.5);
88
- return x >= 0 ? 1 - 0.5 * ib : 0.5 * ib;
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) => k < 0 || k > n || !Number.isInteger(k) ? 0 : Math.exp(logChoose(k) + k * Math.log(p) + (n - k) * Math.log(1 - p));
102
- return {
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
- // Reuse the single pmf closure above instead of rebuilding a whole new
106
- // binomial distribution (and its logChoose closure) on every step.
107
- let s = 0;
108
- for (let i = 0; i <= Math.floor(k); i++)
109
- s += pmf(i);
110
- return Math.min(1, s);
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
- const L = Math.exp(-lambda);
131
- let k = 0;
132
- let prod = 1;
133
- do {
134
- k++;
135
- prod *= rng();
136
- } while (prod > L);
137
- return k - 1;
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 two-tailed Student t-test against a hypothesised mean `mu0`. */
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 mean = sample.reduce((a, b) => a + b, 0) / n;
168
- const variance = sample.reduce((acc, x) => acc + (x - mean) ** 2, 0) / (n - 1);
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
- const p = 2 * (1 - Distributions.studentT(df).cdf(Math.abs(t)));
173
- return { statistic: t, df, pValue: p };
361
+ return { statistic: t, df, pValue: tPValue(t, df, options.alternative) };
174
362
  }
175
- /** Welch's two-sample t-test (unequal variances). */
176
- static tTestTwoSample(a, b) {
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 ma = a.reduce((s, x) => s + x, 0) / na;
180
- const mb = b.reduce((s, x) => s + x, 0) / nb;
181
- const va = a.reduce((s, x) => s + (x - ma) ** 2, 0) / (na - 1);
182
- const vb = b.reduce((s, x) => s + (x - mb) ** 2, 0) / (nb - 1);
183
- const se = Math.sqrt(va / na + vb / nb);
184
- const t = (ma - mb) / se;
185
- const df = (va / na + vb / nb) ** 2 / ((va / na) ** 2 / (na - 1) + (vb / nb) ** 2 / (nb - 1));
186
- const p = 2 * (1 - Distributions.studentT(df).cdf(Math.abs(t)));
187
- return { statistic: t, df, pValue: p };
188
- }
189
- /** Pearson chi-square goodness-of-fit test. */
190
- static chiSquareGoodnessOfFit(observed, expected) {
191
- const chi2 = observed.reduce((s, o, i) => s + (o - expected[i]) ** 2 / expected[i], 0);
192
- const df = observed.length - 1;
193
- const p = 1 - Distributions.chiSquare(df).cdf(chi2);
194
- return { statistic: chi2, df, pValue: p };
195
- }
196
- /** Confidence interval for the mean of a sample at the given level (t-based). */
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.reduce((a, b) => a + b, 0) / n;
200
- const variance = sample.reduce((acc, x) => acc + (x - m) ** 2, 0) / (n - 1);
201
- const se = Math.sqrt(variance / n);
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
- // -- small internal helpers -------------------------------------------------
208
- /** Quantile of Student-t via bisection on its CDF. */
209
- function tQuantile(p, df) {
210
- const dist = Distributions.studentT(df);
211
- let lo = -100;
212
- let hi = 100;
213
- for (let i = 0; i < 200; i++) {
214
- const mid = (lo + hi) / 2;
215
- if (dist.cdf(mid) < p)
216
- lo = mid;
217
- else
218
- hi = mid;
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"}