diff-grok 1.0.9 → 1.1.0

This diff represents the content of publicly available package versions that have been released to one of the supported registries. The information contained in this diff is provided for informational purposes only and reflects changes between package versions as they appear in their respective public registries.
Files changed (200) hide show
  1. package/.claude/settings.local.json +11 -0
  2. package/.github/ISSUE_TEMPLATE/bug_report.md +29 -0
  3. package/.github/ISSUE_TEMPLATE/feature_request.md +20 -0
  4. package/.github/PULL_REQUEST_TEMPLATE.md +22 -0
  5. package/.github/workflows/deploy-docs.yml +15 -0
  6. package/CLAUDE.md +98 -0
  7. package/CONTRIBUTING.md +74 -0
  8. package/README.MD +272 -18
  9. package/images/E5.png +0 -0
  10. package/images/HIRES.png +0 -0
  11. package/images/OREGO.png +0 -0
  12. package/images/POLL-full.png +0 -0
  13. package/images/POLL.png +0 -0
  14. package/images/ROBER.png +0 -0
  15. package/images/VDP0L.png +0 -0
  16. package/images/VDPOL.png +0 -0
  17. package/index.js +2 -2
  18. package/index.js.map +1 -1
  19. package/index.ts +7 -5
  20. package/jest.config.js +1 -0
  21. package/package.json +2 -1
  22. package/src/examples/bioreactor.js +125 -0
  23. package/src/examples/bioreactor.js.map +1 -0
  24. package/src/examples/bioreactor.ts +132 -0
  25. package/src/examples/check-methods.js +23 -9
  26. package/src/examples/check-methods.js.map +1 -1
  27. package/src/examples/check-methods.ts +28 -9
  28. package/src/examples/corr-probs.js +3 -1
  29. package/src/examples/corr-probs.js.map +1 -1
  30. package/src/examples/corr-probs.ts +7 -3
  31. package/src/examples/cyclic-model.js +3 -1
  32. package/src/examples/cyclic-model.js.map +1 -1
  33. package/src/examples/cyclic-model.ts +3 -1
  34. package/src/examples/e5.js +6 -0
  35. package/src/examples/e5.js.map +1 -1
  36. package/src/examples/e5.ts +7 -0
  37. package/src/examples/hires.js +10 -0
  38. package/src/examples/hires.js.map +1 -1
  39. package/src/examples/hires.ts +11 -0
  40. package/src/examples/index.js +21 -7
  41. package/src/examples/index.js.map +1 -1
  42. package/src/examples/index.ts +23 -7
  43. package/src/examples/orego.js +5 -0
  44. package/src/examples/orego.js.map +1 -1
  45. package/src/examples/orego.ts +6 -0
  46. package/src/examples/pollution.js +22 -0
  47. package/src/examples/pollution.js.map +1 -1
  48. package/src/examples/pollution.ts +23 -0
  49. package/src/examples/print-benchmark.js +581 -0
  50. package/src/examples/print-benchmark.js.map +1 -0
  51. package/src/examples/print-benchmark.ts +594 -0
  52. package/src/examples/robertson.js +6 -1
  53. package/src/examples/robertson.js.map +1 -1
  54. package/src/examples/robertson.ts +7 -1
  55. package/src/examples/scripting.js +1 -1
  56. package/src/examples/scripting.js.map +1 -1
  57. package/src/examples/scripting.ts +1 -1
  58. package/src/examples/vdpol.js +4 -0
  59. package/src/examples/vdpol.js.map +1 -1
  60. package/src/examples/vdpol.ts +5 -0
  61. package/src/pipeline/basic-pipeline-creator.js +96 -1
  62. package/src/pipeline/basic-pipeline-creator.js.map +1 -1
  63. package/src/pipeline/basic-pipeline-creator.ts +97 -1
  64. package/src/pipeline/loops-pipeline-creator.js +115 -2
  65. package/src/pipeline/loops-pipeline-creator.js.map +1 -1
  66. package/src/pipeline/loops-pipeline-creator.ts +116 -2
  67. package/src/pipeline/output.js +7 -2
  68. package/src/pipeline/output.js.map +1 -1
  69. package/src/pipeline/output.ts +7 -2
  70. package/src/pipeline/pipeline-creator.js.map +1 -1
  71. package/src/pipeline/pipeline-creator.ts +4 -0
  72. package/src/pipeline/pipeline.js +193 -2
  73. package/src/pipeline/pipeline.js.map +1 -1
  74. package/src/pipeline/pipeline.ts +206 -4
  75. package/src/pipeline/updates-pipeline-creator.js +130 -2
  76. package/src/pipeline/updates-pipeline-creator.js.map +1 -1
  77. package/src/pipeline/updates-pipeline-creator.ts +131 -2
  78. package/src/scripting-tools/constants.js +18 -4
  79. package/src/scripting-tools/constants.js.map +1 -1
  80. package/src/scripting-tools/constants.ts +22 -4
  81. package/src/scripting-tools/index.js.map +1 -1
  82. package/src/scripting-tools/index.ts +3 -2
  83. package/src/scripting-tools/model-error.js +3 -1
  84. package/src/scripting-tools/model-error.js.map +1 -1
  85. package/src/scripting-tools/model-error.ts +3 -1
  86. package/src/scripting-tools/scripting-tools.js +139 -9
  87. package/src/scripting-tools/scripting-tools.js.map +1 -1
  88. package/src/scripting-tools/scripting-tools.ts +228 -21
  89. package/src/solver-tools/ab4-method.js +332 -0
  90. package/src/solver-tools/ab4-method.js.map +1 -0
  91. package/src/solver-tools/ab4-method.ts +395 -0
  92. package/src/solver-tools/ab5-method.js +334 -0
  93. package/src/solver-tools/ab5-method.js.map +1 -0
  94. package/src/solver-tools/ab5-method.ts +397 -0
  95. package/src/solver-tools/callbacks/callback-base.js +2 -0
  96. package/src/solver-tools/callbacks/callback-base.js.map +1 -1
  97. package/src/solver-tools/callbacks/callback-base.ts +4 -0
  98. package/src/solver-tools/callbacks/callback-tools.js +3 -1
  99. package/src/solver-tools/callbacks/callback-tools.js.map +1 -1
  100. package/src/solver-tools/callbacks/callback-tools.ts +3 -1
  101. package/src/solver-tools/callbacks/iter-checker-callback.js +4 -0
  102. package/src/solver-tools/callbacks/iter-checker-callback.js.map +1 -1
  103. package/src/solver-tools/callbacks/iter-checker-callback.ts +4 -0
  104. package/src/solver-tools/callbacks/time-checker-callback.js +4 -0
  105. package/src/solver-tools/callbacks/time-checker-callback.js.map +1 -1
  106. package/src/solver-tools/callbacks/time-checker-callback.ts +4 -0
  107. package/src/solver-tools/index.js +9 -0
  108. package/src/solver-tools/index.js.map +1 -1
  109. package/src/solver-tools/index.ts +9 -0
  110. package/src/solver-tools/lsoda/blas.js +246 -0
  111. package/src/solver-tools/lsoda/blas.js.map +1 -0
  112. package/src/solver-tools/lsoda/blas.ts +267 -0
  113. package/src/solver-tools/lsoda/cfode.js +73 -0
  114. package/src/solver-tools/lsoda/cfode.js.map +1 -0
  115. package/src/solver-tools/lsoda/cfode.ts +83 -0
  116. package/src/solver-tools/lsoda/common.js +94 -0
  117. package/src/solver-tools/lsoda/common.js.map +1 -0
  118. package/src/solver-tools/lsoda/common.ts +135 -0
  119. package/src/solver-tools/lsoda/corfailure.js +25 -0
  120. package/src/solver-tools/lsoda/corfailure.js.map +1 -0
  121. package/src/solver-tools/lsoda/corfailure.ts +30 -0
  122. package/src/solver-tools/lsoda/correction.js +101 -0
  123. package/src/solver-tools/lsoda/correction.js.map +1 -0
  124. package/src/solver-tools/lsoda/correction.ts +119 -0
  125. package/src/solver-tools/lsoda/dense.js +110 -0
  126. package/src/solver-tools/lsoda/dense.js.map +1 -0
  127. package/src/solver-tools/lsoda/dense.ts +129 -0
  128. package/src/solver-tools/lsoda/index.js +4 -0
  129. package/src/solver-tools/lsoda/index.js.map +1 -0
  130. package/src/solver-tools/lsoda/index.ts +5 -0
  131. package/src/solver-tools/lsoda/intdy.js +43 -0
  132. package/src/solver-tools/lsoda/intdy.js.map +1 -0
  133. package/src/solver-tools/lsoda/intdy.ts +50 -0
  134. package/src/solver-tools/lsoda/lsoda.js +567 -0
  135. package/src/solver-tools/lsoda/lsoda.js.map +1 -0
  136. package/src/solver-tools/lsoda/lsoda.ts +632 -0
  137. package/src/solver-tools/lsoda/methodswitch.js +97 -0
  138. package/src/solver-tools/lsoda/methodswitch.js.map +1 -0
  139. package/src/solver-tools/lsoda/methodswitch.ts +107 -0
  140. package/src/solver-tools/lsoda/orderswitch.js +86 -0
  141. package/src/solver-tools/lsoda/orderswitch.js.map +1 -0
  142. package/src/solver-tools/lsoda/orderswitch.ts +95 -0
  143. package/src/solver-tools/lsoda/prja.js +43 -0
  144. package/src/solver-tools/lsoda/prja.js.map +1 -0
  145. package/src/solver-tools/lsoda/prja.ts +49 -0
  146. package/src/solver-tools/lsoda/scaleh.js +30 -0
  147. package/src/solver-tools/lsoda/scaleh.js.map +1 -0
  148. package/src/solver-tools/lsoda/scaleh.ts +33 -0
  149. package/src/solver-tools/lsoda/solsy.js +14 -0
  150. package/src/solver-tools/lsoda/solsy.js.map +1 -0
  151. package/src/solver-tools/lsoda/solsy.ts +18 -0
  152. package/src/solver-tools/lsoda/stoda.js +259 -0
  153. package/src/solver-tools/lsoda/stoda.js.map +1 -0
  154. package/src/solver-tools/lsoda/stoda.ts +286 -0
  155. package/src/solver-tools/lsoda-method.js +120 -0
  156. package/src/solver-tools/lsoda-method.js.map +1 -0
  157. package/src/solver-tools/lsoda-method.ts +136 -0
  158. package/src/solver-tools/mrt-method.js +121 -19
  159. package/src/solver-tools/mrt-method.js.map +1 -1
  160. package/src/solver-tools/mrt-method.ts +151 -22
  161. package/src/solver-tools/rk3-method.js +210 -0
  162. package/src/solver-tools/rk3-method.js.map +1 -0
  163. package/src/solver-tools/rk3-method.ts +257 -0
  164. package/src/solver-tools/rk4-method.js +236 -0
  165. package/src/solver-tools/rk4-method.js.map +1 -0
  166. package/src/solver-tools/rk4-method.ts +286 -0
  167. package/src/solver-tools/rkdp-method.js +247 -0
  168. package/src/solver-tools/rkdp-method.js.map +1 -0
  169. package/src/solver-tools/rkdp-method.ts +300 -0
  170. package/src/solver-tools/ros34prw-method.js +40 -2
  171. package/src/solver-tools/ros34prw-method.js.map +1 -1
  172. package/src/solver-tools/ros34prw-method.ts +40 -2
  173. package/src/solver-tools/ros3prw-method.js +40 -2
  174. package/src/solver-tools/ros3prw-method.js.map +1 -1
  175. package/src/solver-tools/ros3prw-method.ts +40 -2
  176. package/src/solver-tools/solver-defs.js +15 -3
  177. package/src/solver-tools/solver-defs.js.map +1 -1
  178. package/src/solver-tools/solver-defs.ts +91 -7
  179. package/src/tests/correctness.test.js +6 -4
  180. package/src/tests/correctness.test.js.map +1 -1
  181. package/src/tests/correctness.test.ts +6 -4
  182. package/src/tests/lsoda-lib-correctness.test.js +223 -0
  183. package/src/tests/lsoda-lib-correctness.test.js.map +1 -0
  184. package/src/tests/lsoda-lib-correctness.test.ts +283 -0
  185. package/src/tests/lsoda-lib-performance.test.js +95 -0
  186. package/src/tests/lsoda-lib-performance.test.js.map +1 -0
  187. package/src/tests/lsoda-lib-performance.test.ts +126 -0
  188. package/src/tests/performance.test.js +10 -8
  189. package/src/tests/performance.test.js.map +1 -1
  190. package/src/tests/performance.test.ts +10 -8
  191. package/src/tests/test-defs.js +14 -1
  192. package/src/tests/test-defs.js.map +1 -1
  193. package/src/tests/test-defs.ts +15 -1
  194. package/src/worker-tools/scripting.js +95 -1
  195. package/src/worker-tools/scripting.js.map +1 -1
  196. package/src/worker-tools/scripting.ts +104 -2
  197. package/src/worker-tools/solving.js +18 -2
  198. package/src/worker-tools/solving.js.map +1 -1
  199. package/src/worker-tools/solving.ts +25 -2
  200. package/typedoc.json +8 -0
@@ -0,0 +1,567 @@
1
+ import { LsodaContext, LsodaCommon, ETA, } from './common';
2
+ import { vmnorm } from './blas';
3
+ import { intdy } from './intdy';
4
+ import { stoda } from './stoda';
5
+ import { DenseOutput } from './dense';
6
+ // ---- Internal helpers ----
7
+ function captureSnapshot(c, neq) {
8
+ const nq = c.nq;
9
+ const yh = new Array(nq + 1);
10
+ for (let j = 0; j <= nq; j++) {
11
+ const row = new Float64Array(neq);
12
+ for (let i = 0; i < neq; i++)
13
+ row[i] = c.yh[j + 1][i + 1]; // convert 1-based to 0-based
14
+ yh[j] = row;
15
+ }
16
+ return { tn: c.tn, h: c.h, hu: c.hu, nq, yh };
17
+ }
18
+ function ewset(ctx, ycur) {
19
+ const c = ctx.common;
20
+ const neq = ctx.neq;
21
+ const rtol = ctx.opt.rtol;
22
+ const atol = ctx.opt.atol;
23
+ for (let i = 1; i <= neq; i++)
24
+ c.ewt[i] = rtol[i] * Math.abs(ycur[i]) + atol[i];
25
+ for (let i = 1; i <= neq; i++)
26
+ c.ewt[i] = 1. / c.ewt[i];
27
+ }
28
+ function checkOpt(ctx, opt) {
29
+ const mxstp0 = 500;
30
+ const mord = [0, 12, 5];
31
+ if (ctx.state === 0)
32
+ ctx.state = 1;
33
+ if (ctx.state === 1) {
34
+ opt.h0 = 0.;
35
+ opt.mxordn = mord[1];
36
+ opt.mxords = mord[2];
37
+ }
38
+ if (ctx.neq <= 0) {
39
+ ctx.error = `[lsoda] neq = ${ctx.neq} is less than 1`;
40
+ return false;
41
+ }
42
+ // Check rtol and atol (1-indexed)
43
+ if (ctx.state === 1 || ctx.state === 3) {
44
+ for (let i = 1; i <= ctx.neq; i++) {
45
+ if (opt.rtol[i] < 0.) {
46
+ ctx.error = `[lsoda] rtol = ${opt.rtol[i]} is less than 0.`;
47
+ return false;
48
+ }
49
+ if (opt.atol[i] < 0.) {
50
+ ctx.error = `[lsoda] atol = ${opt.atol[i]} is less than 0.`;
51
+ return false;
52
+ }
53
+ }
54
+ }
55
+ if (opt.itask === 0)
56
+ opt.itask = 1;
57
+ if (opt.itask < 1 || opt.itask > 5) {
58
+ ctx.error = `[lsoda] illegal itask = ${opt.itask}`;
59
+ return false;
60
+ }
61
+ if (opt.ixpr < 0 || opt.ixpr > 1) {
62
+ ctx.error = `[lsoda] ixpr = ${opt.ixpr} is illegal`;
63
+ return false;
64
+ }
65
+ if (opt.mxstep < 0) {
66
+ ctx.error = '[lsoda] mxstep < 0';
67
+ return false;
68
+ }
69
+ if (opt.mxstep === 0)
70
+ opt.mxstep = mxstp0;
71
+ if (opt.mxhnil < 0) {
72
+ ctx.error = '[lsoda] mxhnil < 0';
73
+ return false;
74
+ }
75
+ if (ctx.state === 1) {
76
+ if (opt.mxordn < 0) {
77
+ ctx.error = `[lsoda] mxordn = ${opt.mxordn} is less than 0`;
78
+ return false;
79
+ }
80
+ if (opt.mxordn === 0)
81
+ opt.mxordn = 100;
82
+ opt.mxordn = Math.min(opt.mxordn, mord[1]);
83
+ if (opt.mxords < 0) {
84
+ ctx.error = `[lsoda] mxords = ${opt.mxords} is less than 0`;
85
+ return false;
86
+ }
87
+ if (opt.mxords === 0)
88
+ opt.mxords = 100;
89
+ opt.mxords = Math.min(opt.mxords, mord[2]);
90
+ }
91
+ if (opt.hmax < 0.) {
92
+ ctx.error = '[lsoda] hmax < 0.';
93
+ return false;
94
+ }
95
+ opt.hmxi = 0.;
96
+ if (opt.hmax > 0)
97
+ opt.hmxi = 1. / opt.hmax;
98
+ if (opt.hmin < 0.) {
99
+ ctx.error = '[lsoda] hmin < 0.';
100
+ return false;
101
+ }
102
+ return true;
103
+ }
104
+ function allocMem(ctx) {
105
+ const c = ctx.common;
106
+ const nyh = ctx.neq;
107
+ const lenyh = 1 + Math.max(ctx.opt.mxordn, ctx.opt.mxords);
108
+ // yh: array of (lenyh+1) Float64Arrays, each of size (nyh+1)
109
+ c.yh = new Array(lenyh + 1);
110
+ for (let i = 0; i <= lenyh; i++)
111
+ c.yh[i] = new Float64Array(nyh + 1);
112
+ // wm: array of (nyh+1) Float64Arrays, each of size (nyh+1)
113
+ c.wm = new Array(nyh + 1);
114
+ for (let i = 0; i <= nyh; i++)
115
+ c.wm[i] = new Float64Array(nyh + 1);
116
+ c.ewt = new Float64Array(nyh + 1);
117
+ c.savf = new Float64Array(nyh + 1);
118
+ c.acor = new Float64Array(nyh + 1);
119
+ c.ipvt = new Int32Array(nyh + 1);
120
+ // Allocate elco and tesco
121
+ c.elco = new Array(13);
122
+ for (let i = 0; i < 13; i++)
123
+ c.elco[i] = new Float64Array(14);
124
+ c.tesco = new Array(13);
125
+ for (let i = 0; i < 13; i++)
126
+ c.tesco[i] = new Float64Array(4);
127
+ return true;
128
+ }
129
+ // ---- Public API ----
130
+ /**
131
+ * Prepare a context with options. Must be called before lsoda().
132
+ * Returns true on success.
133
+ */
134
+ export function lsodaPrepare(ctx, opt) {
135
+ ctx.common = new LsodaCommon();
136
+ ctx.opt = opt;
137
+ if (!checkOpt(ctx, opt))
138
+ return false;
139
+ return allocMem(ctx);
140
+ }
141
+ /**
142
+ * Reset context for reuse with a different problem.
143
+ * All options are preserved, internal state is zeroed.
144
+ */
145
+ export function lsodaReset(ctx) {
146
+ const c = ctx.common;
147
+ // Zero all scalar state
148
+ c.h = 0;
149
+ c.hu = 0;
150
+ c.rc = 0;
151
+ c.tn = 0;
152
+ c.tsw = 0;
153
+ c.pdnorm = 0;
154
+ c.crate = 0;
155
+ c.hold = 0;
156
+ c.rmax = 0;
157
+ c.pdest = 0;
158
+ c.pdlast = 0;
159
+ c.ialth = 0;
160
+ c.ipup = 0;
161
+ c.nslp = 0;
162
+ c.icount = 0;
163
+ c.irflag = 0;
164
+ c.imxer = 0;
165
+ c.illin = 0;
166
+ c.nhnil = 0;
167
+ c.nslast = 0;
168
+ c.jcur = 0;
169
+ c.meth = 0;
170
+ c.mused = 0;
171
+ c.nq = 0;
172
+ c.nst = 0;
173
+ c.ncf = 0;
174
+ c.nfe = 0;
175
+ c.nje = 0;
176
+ c.nqu = 0;
177
+ c.miter = 0;
178
+ c.el.fill(0);
179
+ // Zero arrays
180
+ for (let i = 0; i < c.yh.length; i++)
181
+ c.yh[i].fill(0);
182
+ for (let i = 0; i < c.wm.length; i++)
183
+ c.wm[i].fill(0);
184
+ c.ewt.fill(0);
185
+ c.savf.fill(0);
186
+ c.acor.fill(0);
187
+ c.ipvt.fill(0);
188
+ for (let i = 0; i < c.elco.length; i++)
189
+ c.elco[i].fill(0);
190
+ for (let i = 0; i < c.tesco.length; i++)
191
+ c.tesco[i].fill(0);
192
+ ctx.state = 1;
193
+ ctx.error = null;
194
+ }
195
+ /**
196
+ * Free context resources. In JS/TS, this is mostly a no-op (GC handles it).
197
+ */
198
+ export function lsodaFree(ctx) {
199
+ if (ctx.error)
200
+ console.error(`unhandled error message: ${ctx.error}`);
201
+ ctx.common = null;
202
+ }
203
+ /**
204
+ * Main solver. Integrates from t towards tout.
205
+ * y: 1-indexed Float64Array (y[0] unused, y[1..neq] are state values)
206
+ * Returns { t, state } — updated time and solver state.
207
+ */
208
+ export function lsoda(ctx, y, tIn, tout) {
209
+ const c = ctx.common;
210
+ const opt = ctx.opt;
211
+ const neq = ctx.neq;
212
+ let t = tIn;
213
+ let kflag;
214
+ let jstart = 0;
215
+ let ihit = 0;
216
+ // Helper functions matching C macros
217
+ function hardfailure(msg) {
218
+ ctx.error = msg;
219
+ ctx.state = -3;
220
+ return { t, state: ctx.state };
221
+ }
222
+ function softfailure(code, msg) {
223
+ ctx.error = msg;
224
+ for (let i = 1; i <= neq; i++)
225
+ y[i] = c.yh[1][i];
226
+ t = c.tn;
227
+ ctx.state = code;
228
+ return { t, state: ctx.state };
229
+ }
230
+ function successreturn() {
231
+ for (let i = 1; i <= neq; i++)
232
+ y[i] = c.yh[1][i];
233
+ t = c.tn;
234
+ if (itask === 4 || itask === 5) {
235
+ if (ihit)
236
+ t = tcrit;
237
+ }
238
+ ctx.state = 2;
239
+ return { t, state: ctx.state };
240
+ }
241
+ function intdyreturn() {
242
+ const iflag = intdy(ctx, tout, 0, y);
243
+ if (iflag !== 0) {
244
+ ctx.error = `[lsoda] trouble from intdy, itask = ${itask}, tout = ${tout}`;
245
+ for (let i = 1; i <= neq; i++)
246
+ y[i] = c.yh[1][i];
247
+ t = c.tn;
248
+ }
249
+ t = tout;
250
+ ctx.state = 2;
251
+ return { t, state: ctx.state };
252
+ }
253
+ if (c === null)
254
+ return hardfailure('[lsoda] illegal common block did you call lsoda_prepare?');
255
+ // Block b: initial/continuation call setup
256
+ let h0 = 0;
257
+ let tcrit = 0;
258
+ const rtol = opt.rtol;
259
+ const atol = opt.atol;
260
+ if (ctx.state === 1 || ctx.state === 3) {
261
+ h0 = opt.h0;
262
+ if (ctx.state === 1) {
263
+ if ((tout - t) * h0 < 0.)
264
+ return hardfailure(`[lsoda] tout = ${tout} behind t = ${t}. integration direction is given by ${h0}`);
265
+ }
266
+ }
267
+ const itask = opt.itask;
268
+ if (ctx.state === 3)
269
+ jstart = -1;
270
+ // Block c: initial call only (state == 1)
271
+ if (ctx.state === 1) {
272
+ c.meth = 1;
273
+ c.tn = t;
274
+ c.tsw = t;
275
+ if (itask === 4 || itask === 5) {
276
+ tcrit = opt.tcrit;
277
+ if ((tcrit - tout) * (tout - t) < 0.)
278
+ return hardfailure('[lsoda] itask = 4 or 5 and tcrit behind tout');
279
+ if (h0 !== 0. && (t + h0 - tcrit) * h0 > 0.)
280
+ h0 = tcrit - t;
281
+ }
282
+ jstart = 0;
283
+ c.nq = 1;
284
+ // Initial call to f
285
+ ctx.func(t, y.subarray(1), c.yh[2].subarray(1), ctx.data);
286
+ c.nfe = 1;
287
+ // Load initial value vector
288
+ for (let i = 1; i <= neq; i++)
289
+ c.yh[1][i] = y[i];
290
+ // Load and invert ewt
291
+ ewset(ctx, y);
292
+ for (let i = 1; i <= neq; i++) {
293
+ if (c.ewt[i] <= 0.)
294
+ return hardfailure(`[lsoda] ewt[${i}] = ${c.ewt[i]} <= 0.`);
295
+ }
296
+ // Compute initial step size h0 if not given
297
+ if (h0 === 0.) {
298
+ const tdist = Math.abs(tout - t);
299
+ const w0 = Math.max(Math.abs(t), Math.abs(tout));
300
+ if (tdist < 2. * ETA * w0)
301
+ return hardfailure('[lsoda] tout too close to t to start integration');
302
+ let tol = 0.;
303
+ for (let i = 1; i <= neq; i++)
304
+ tol = Math.max(tol, rtol[i]);
305
+ if (tol <= 0.) {
306
+ for (let i = 1; i <= neq; i++) {
307
+ const atoli = atol[i];
308
+ const ayi = Math.abs(y[i]);
309
+ if (ayi !== 0.)
310
+ tol = Math.max(tol, atoli / ayi);
311
+ }
312
+ }
313
+ tol = Math.max(tol, 100. * ETA);
314
+ tol = Math.min(tol, 0.001);
315
+ const sum = vmnorm(neq, c.yh[2], c.ewt);
316
+ const sumSq = 1. / (tol * w0 * w0) + tol * sum * sum;
317
+ h0 = 1. / Math.sqrt(sumSq);
318
+ h0 = Math.min(h0, tdist);
319
+ h0 = h0 * ((tout - t >= 0.) ? 1. : -1.);
320
+ }
321
+ // Adjust h0 for hmax
322
+ const rh = Math.abs(h0) * opt.hmxi;
323
+ if (rh > 1.)
324
+ h0 /= rh;
325
+ // Load h and scale yh[2]
326
+ c.h = h0;
327
+ for (let i = 1; i <= neq; i++)
328
+ c.yh[2][i] *= h0;
329
+ }
330
+ // Block d: continuation calls (state == 2 or 3)
331
+ if (ctx.state === 2 || ctx.state === 3) {
332
+ jstart = 1;
333
+ c.nslast = c.nst;
334
+ switch (itask) {
335
+ case 1:
336
+ if ((c.tn - tout) * c.h >= 0.)
337
+ return intdyreturn();
338
+ break;
339
+ case 2:
340
+ break;
341
+ case 3: {
342
+ const tp = c.tn - c.hu * (1. + 100. * ETA);
343
+ if ((tp - tout) * c.h > 0.)
344
+ return hardfailure(`[lsoda] itask = ${itask} and tout behind tcur - hu`);
345
+ if ((c.tn - tout) * c.h < 0.)
346
+ break;
347
+ return successreturn();
348
+ }
349
+ case 4:
350
+ tcrit = opt.tcrit;
351
+ if ((c.tn - tcrit) * c.h > 0.)
352
+ return hardfailure('[lsoda] itask = 4 or 5 and tcrit behind tcur');
353
+ if ((tcrit - tout) * c.h < 0.)
354
+ return hardfailure('[lsoda] itask = 4 or 5 and tcrit behind tout');
355
+ if ((c.tn - tout) * c.h >= 0.)
356
+ return intdyreturn();
357
+ // fall through to case 5 logic
358
+ // eslint-disable-next-line no-fallthrough
359
+ case 5:
360
+ if (itask === 5) {
361
+ tcrit = opt.tcrit;
362
+ if ((c.tn - tcrit) * c.h > 0.)
363
+ return hardfailure('[lsoda] itask = 4 or 5 and tcrit behind tcur');
364
+ }
365
+ {
366
+ const hmx = Math.abs(c.tn) + Math.abs(c.h);
367
+ ihit = Math.abs(c.tn - tcrit) <= (100. * ETA * hmx) ? 1 : 0;
368
+ if (ihit) {
369
+ t = tcrit;
370
+ return successreturn();
371
+ }
372
+ const tnext = c.tn + c.h * (1. + 4. * ETA);
373
+ if ((tnext - tcrit) * c.h <= 0.)
374
+ break;
375
+ c.h = (tcrit - c.tn) * (1. - 4. * ETA);
376
+ if (ctx.state === 2)
377
+ jstart = -2;
378
+ }
379
+ break;
380
+ }
381
+ }
382
+ // Block e: main integration loop
383
+ while (true) {
384
+ if (ctx.state !== 1 || c.nst !== 0) {
385
+ if ((c.nst - c.nslast) >= opt.mxstep)
386
+ return softfailure(-1, `[lsoda] ${opt.mxstep} steps taken before reaching tout`);
387
+ ewset(ctx, c.yh[1]);
388
+ for (let i = 1; i <= neq; i++) {
389
+ if (c.ewt[i] <= 0.)
390
+ return softfailure(-6, `[lsoda] ewt[${i}] = ${c.ewt[i]} <= 0.`);
391
+ }
392
+ }
393
+ const tolsf = ETA * vmnorm(neq, c.yh[1], c.ewt);
394
+ if (tolsf > 0.01) {
395
+ const scaled = tolsf * 200.;
396
+ if (c.nst === 0) {
397
+ return hardfailure('lsoda -- at start of problem, too much accuracy requested' +
398
+ ` for precision of machine, suggested scaling factor = ${scaled}`);
399
+ }
400
+ return softfailure(-2, `lsoda -- at t = ${t}, too much accuracy requested` +
401
+ ` for precision of machine, suggested scaling factor = ${scaled}`);
402
+ }
403
+ if ((c.tn + c.h) === c.tn) {
404
+ c.nhnil++;
405
+ if (c.nhnil <= opt.mxhnil) {
406
+ console.error(`lsoda -- warning..internal t = ${c.tn} and h = ${c.h} are` +
407
+ ' such that t + h = t on the next step');
408
+ if (c.nhnil === opt.mxhnil) {
409
+ console.error(`lsoda -- above warning has been issued ${c.nhnil}` +
410
+ ' times, it will not be issued again for this problem');
411
+ }
412
+ }
413
+ }
414
+ // Call stoda
415
+ kflag = stoda(ctx, y, jstart);
416
+ if (kflag === 0) {
417
+ // Block f: successful return from stoda
418
+ // Dense output: capture snapshot
419
+ if (ctx.snapshots)
420
+ ctx.snapshots.push(captureSnapshot(c, neq));
421
+ jstart = 1;
422
+ if (c.meth !== c.mused) {
423
+ c.tsw = c.tn;
424
+ jstart = -1;
425
+ if (opt.ixpr) {
426
+ if (c.meth === 2) {
427
+ console.error('[lsoda] a switch to the stiff method has occurred' +
428
+ ` at t = ${c.tn}, tentative step size h = ${c.h}, step nst = ${c.nst}`);
429
+ }
430
+ if (c.meth === 1) {
431
+ console.error('[lsoda] a switch to the nonstiff method has occurred' +
432
+ ` at t = ${c.tn}, tentative step size h = ${c.h}, step nst = ${c.nst}`);
433
+ }
434
+ }
435
+ }
436
+ if (itask === 1) {
437
+ if ((c.tn - tout) * c.h < 0.)
438
+ continue;
439
+ return intdyreturn();
440
+ }
441
+ if (itask === 2)
442
+ return successreturn();
443
+ if (itask === 3) {
444
+ if ((c.tn - tout) * c.h >= 0.)
445
+ return successreturn();
446
+ continue;
447
+ }
448
+ if (itask === 4) {
449
+ tcrit = opt.tcrit;
450
+ if ((c.tn - tout) * c.h >= 0.)
451
+ return intdyreturn();
452
+ else {
453
+ const hmx = Math.abs(c.tn) + Math.abs(c.h);
454
+ ihit = Math.abs(c.tn - tcrit) <= (100. * ETA * hmx) ? 1 : 0;
455
+ if (ihit)
456
+ return successreturn();
457
+ const tnext = c.tn + c.h * (1. + 4. * ETA);
458
+ if ((tnext - tcrit) * c.h <= 0.)
459
+ continue;
460
+ c.h = (tcrit - c.tn) * (1. - 4. * ETA);
461
+ jstart = -2;
462
+ continue;
463
+ }
464
+ }
465
+ if (itask === 5) {
466
+ tcrit = opt.tcrit;
467
+ const hmx = Math.abs(c.tn) + Math.abs(c.h);
468
+ ihit = Math.abs(c.tn - tcrit) <= (100. * ETA * hmx) ? 1 : 0;
469
+ return successreturn();
470
+ }
471
+ }
472
+ // kflag == -1 or -2: error test / convergence failures
473
+ if (kflag === -1 || kflag === -2) {
474
+ let big = 0.;
475
+ c.imxer = 1;
476
+ for (let i = 1; i <= neq; i++) {
477
+ const size = Math.abs(c.acor[i]) * c.ewt[i];
478
+ if (big < size) {
479
+ big = size;
480
+ c.imxer = i;
481
+ }
482
+ }
483
+ if (kflag === -1) {
484
+ return softfailure(-4, `lsoda -- at t = ${c.tn} and step size h = ${c.h},` +
485
+ ' the error test failed repeatedly or with abs(h) = hmin');
486
+ }
487
+ if (kflag === -2) {
488
+ return softfailure(-5, `lsoda -- at t = ${c.tn} and step size h = ${c.h},` +
489
+ ' the corrector convergence failed repeatedly or with abs(h) = hmin');
490
+ }
491
+ }
492
+ }
493
+ }
494
+ /**
495
+ * High-level LSODA solver class with 0-based public interface.
496
+ */
497
+ export class Lsoda {
498
+ constructor(f, neq, opt, data) {
499
+ this.denseEnabled = false;
500
+ // Wrap user's 0-based function to internal 1-based convention
501
+ const wrappedFunc = (t, y, ydot, data) => {
502
+ return f(t, y, ydot, data);
503
+ };
504
+ this.ctx = new LsodaContext(wrappedFunc, neq, data);
505
+ this.internalY = new Float64Array(neq + 1); // 1-indexed
506
+ // Build full options with defaults
507
+ const fullOpt = Object.assign({ ixpr: 0, mxstep: 0, mxhnil: 0, mxordn: 0, mxords: 0, tcrit: 0, h0: 0, hmax: 0, hmin: 0, hmxi: 0, itask: 1, rtol: new Float64Array(neq + 1), atol: new Float64Array(neq + 1) }, opt);
508
+ // If rtol/atol provided as partial, make sure they're 1-indexed Float64Arrays
509
+ if (opt === null || opt === void 0 ? void 0 : opt.rtol) {
510
+ if (opt.rtol.length === neq) {
511
+ // User provided 0-based, convert to 1-based
512
+ const r = new Float64Array(neq + 1);
513
+ for (let i = 0; i < neq; i++)
514
+ r[i + 1] = opt.rtol[i];
515
+ fullOpt.rtol = r;
516
+ }
517
+ else
518
+ fullOpt.rtol = opt.rtol;
519
+ }
520
+ if (opt === null || opt === void 0 ? void 0 : opt.atol) {
521
+ if (opt.atol.length === neq) {
522
+ const a = new Float64Array(neq + 1);
523
+ for (let i = 0; i < neq; i++)
524
+ a[i + 1] = opt.atol[i];
525
+ fullOpt.atol = a;
526
+ }
527
+ else
528
+ fullOpt.atol = opt.atol;
529
+ }
530
+ lsodaPrepare(this.ctx, fullOpt);
531
+ if (opt === null || opt === void 0 ? void 0 : opt.dense) {
532
+ this.denseEnabled = true;
533
+ this.ctx.snapshots = [];
534
+ }
535
+ }
536
+ /**
537
+ * Integrate from current t to tout.
538
+ * y: 0-based array of state values.
539
+ * Returns updated { y, t }.
540
+ */
541
+ solve(y, t, tout) {
542
+ const neq = this.ctx.neq;
543
+ // Copy 0-based user y into 1-based internal array
544
+ for (let i = 0; i < neq; i++)
545
+ this.internalY[i + 1] = y[i];
546
+ const result = lsoda(this.ctx, this.internalY, t, tout);
547
+ // Copy back to 0-based
548
+ const out = new Float64Array(neq);
549
+ for (let i = 0; i < neq; i++)
550
+ out[i] = this.internalY[i + 1];
551
+ return { y: out, t: result.t };
552
+ }
553
+ get state() { return this.ctx.state; }
554
+ get error() { return this.ctx.error; }
555
+ /** Returns a DenseOutput interpolator from collected snapshots. */
556
+ getDenseOutput() {
557
+ if (!this.ctx.snapshots || this.ctx.snapshots.length === 0)
558
+ throw new Error('[Lsoda] dense output not enabled or no steps taken — pass { dense: true } in options');
559
+ return new DenseOutput(this.ctx.snapshots, this.ctx.neq);
560
+ }
561
+ reset() {
562
+ lsodaReset(this.ctx);
563
+ if (this.denseEnabled)
564
+ this.ctx.snapshots = [];
565
+ }
566
+ }
567
+ //# sourceMappingURL=lsoda.js.map
@@ -0,0 +1 @@
1
+ {"version":3,"file":"lsoda.js","sourceRoot":"","sources":["lsoda.ts"],"names":[],"mappings":"AAAA,OAAO,EACL,YAAY,EAAE,WAAW,EAA0C,GAAG,GACvE,MAAM,UAAU,CAAC;AAClB,OAAO,EAAC,MAAM,EAAC,MAAM,QAAQ,CAAC;AAC9B,OAAO,EAAC,KAAK,EAAC,MAAM,SAAS,CAAC;AAC9B,OAAO,EAAC,KAAK,EAAC,MAAM,SAAS,CAAC;AAC9B,OAAO,EAAC,WAAW,EAAC,MAAM,SAAS,CAAC;AAEpC,6BAA6B;AAE7B,SAAS,eAAe,CAAC,CAAc,EAAE,GAAW;IAClD,MAAM,EAAE,GAAG,CAAC,CAAC,EAAE,CAAC;IAChB,MAAM,EAAE,GAAmB,IAAI,KAAK,CAAC,EAAE,GAAG,CAAC,CAAC,CAAC;IAC7C,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,IAAI,EAAE,EAAE,CAAC,EAAE,EAAE,CAAC;QAC7B,MAAM,GAAG,GAAG,IAAI,YAAY,CAAC,GAAG,CAAC,CAAC;QAClC,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,GAAG,EAAE,CAAC,EAAE;YAC1B,GAAG,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,EAAE,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,6BAA6B;QAC5D,EAAE,CAAC,CAAC,CAAC,GAAG,GAAG,CAAC;IACd,CAAC;IACD,OAAO,EAAC,EAAE,EAAE,CAAC,CAAC,EAAE,EAAE,CAAC,EAAE,CAAC,CAAC,CAAC,EAAE,EAAE,EAAE,CAAC,CAAC,EAAE,EAAE,EAAE,EAAE,EAAE,EAAC,CAAC;AAC9C,CAAC;AAED,SAAS,KAAK,CAAC,GAAiB,EAAE,IAAkB;IAClD,MAAM,CAAC,GAAG,GAAG,CAAC,MAAO,CAAC;IACtB,MAAM,GAAG,GAAG,GAAG,CAAC,GAAG,CAAC;IACpB,MAAM,IAAI,GAAG,GAAG,CAAC,GAAI,CAAC,IAAI,CAAC;IAC3B,MAAM,IAAI,GAAG,GAAG,CAAC,GAAI,CAAC,IAAI,CAAC;IAC3B,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,IAAI,GAAG,EAAE,CAAC,EAAE;QAC3B,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,CAAC,CAAC,CAAC;IACnD,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,IAAI,GAAG,EAAE,CAAC,EAAE;QAC3B,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,GAAG,EAAE,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC;AAC7B,CAAC;AAED,SAAS,QAAQ,CAAC,GAAiB,EAAE,GAAa;IAChD,MAAM,MAAM,GAAG,GAAG,CAAC;IACnB,MAAM,IAAI,GAAG,CAAC,CAAC,EAAE,EAAE,EAAE,CAAC,CAAC,CAAC;IAExB,IAAI,GAAG,CAAC,KAAK,KAAK,CAAC;QAAE,GAAG,CAAC,KAAK,GAAG,CAAC,CAAC;IACnC,IAAI,GAAG,CAAC,KAAK,KAAK,CAAC,EAAE,CAAC;QACpB,GAAG,CAAC,EAAE,GAAG,EAAE,CAAC;QACZ,GAAG,CAAC,MAAM,GAAG,IAAI,CAAC,CAAC,CAAC,CAAC;QACrB,GAAG,CAAC,MAAM,GAAG,IAAI,CAAC,CAAC,CAAC,CAAC;IACvB,CAAC;IAED,IAAI,GAAG,CAAC,GAAG,IAAI,CAAC,EAAE,CAAC;QACjB,GAAG,CAAC,KAAK,GAAG,iBAAiB,GAAG,CAAC,GAAG,iBAAiB,CAAC;QACtD,OAAO,KAAK,CAAC;IACf,CAAC;IAED,kCAAkC;IAClC,IAAI,GAAG,CAAC,KAAK,KAAK,CAAC,IAAI,GAAG,CAAC,KAAK,KAAK,CAAC,EAAE,CAAC;QACvC,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,IAAI,GAAG,CAAC,GAAG,EAAE,CAAC,EAAE,EAAE,CAAC;YAClC,IAAI,GAAG,CAAC,IAAI,CAAC,CAAC,CAAC,GAAG,EAAE,EAAE,CAAC;gBACrB,GAAG,CAAC,KAAK,GAAG,kBAAkB,GAAG,CAAC,IAAI,CAAC,CAAC,CAAC,kBAAkB,CAAC;gBAC5D,OAAO,KAAK,CAAC;YACf,CAAC;YACD,IAAI,GAAG,CAAC,IAAI,CAAC,CAAC,CAAC,GAAG,EAAE,EAAE,CAAC;gBACrB,GAAG,CAAC,KAAK,GAAG,kBAAkB,GAAG,CAAC,IAAI,CAAC,CAAC,CAAC,kBAAkB,CAAC;gBAC5D,OAAO,KAAK,CAAC;YACf,CAAC;QACH,CAAC;IACH,CAAC;IAED,IAAI,GAAG,CAAC,KAAK,KAAK,CAAC;QAAE,GAAG,CAAC,KAAK,GAAG,CAAC,CAAC;IACnC,IAAI,GAAG,CAAC,KAAK,GAAG,CAAC,IAAI,GAAG,CAAC,KAAK,GAAG,CAAC,EAAE,CAAC;QACnC,GAAG,CAAC,KAAK,GAAG,2BAA2B,GAAG,CAAC,KAAK,EAAE,CAAC;QACnD,OAAO,KAAK,CAAC;IACf,CAAC;IACD,IAAI,GAAG,CAAC,IAAI,GAAG,CAAC,IAAI,GAAG,CAAC,IAAI,GAAG,CAAC,EAAE,CAAC;QACjC,GAAG,CAAC,KAAK,GAAG,kBAAkB,GAAG,CAAC,IAAI,aAAa,CAAC;QACpD,OAAO,KAAK,CAAC;IACf,CAAC;IACD,IAAI,GAAG,CAAC,MAAM,GAAG,CAAC,EAAE,CAAC;QACnB,GAAG,CAAC,KAAK,GAAG,oBAAoB,CAAC;QACjC,OAAO,KAAK,CAAC;IACf,CAAC;IACD,IAAI,GAAG,CAAC,MAAM,KAAK,CAAC;QAAE,GAAG,CAAC,MAAM,GAAG,MAAM,CAAC;IAC1C,IAAI,GAAG,CAAC,MAAM,GAAG,CAAC,EAAE,CAAC;QACnB,GAAG,CAAC,KAAK,GAAG,oBAAoB,CAAC;QACjC,OAAO,KAAK,CAAC;IACf,CAAC;IAED,IAAI,GAAG,CAAC,KAAK,KAAK,CAAC,EAAE,CAAC;QACpB,IAAI,GAAG,CAAC,MAAM,GAAG,CAAC,EAAE,CAAC;YACnB,GAAG,CAAC,KAAK,GAAG,oBAAoB,GAAG,CAAC,MAAM,iBAAiB,CAAC;YAC5D,OAAO,KAAK,CAAC;QACf,CAAC;QACD,IAAI,GAAG,CAAC,MAAM,KAAK,CAAC;YAAE,GAAG,CAAC,MAAM,GAAG,GAAG,CAAC;QACvC,GAAG,CAAC,MAAM,GAAG,IAAI,CAAC,GAAG,CAAC,GAAG,CAAC,MAAM,EAAE,IAAI,CAAC,CAAC,CAAC,CAAC,CAAC;QAC3C,IAAI,GAAG,CAAC,MAAM,GAAG,CAAC,EAAE,CAAC;YACnB,GAAG,CAAC,KAAK,GAAG,oBAAoB,GAAG,CAAC,MAAM,iBAAiB,CAAC;YAC5D,OAAO,KAAK,CAAC;QACf,CAAC;QACD,IAAI,GAAG,CAAC,MAAM,KAAK,CAAC;YAAE,GAAG,CAAC,MAAM,GAAG,GAAG,CAAC;QACvC,GAAG,CAAC,MAAM,GAAG,IAAI,CAAC,GAAG,CAAC,GAAG,CAAC,MAAM,EAAE,IAAI,CAAC,CAAC,CAAC,CAAC,CAAC;IAC7C,CAAC;IAED,IAAI,GAAG,CAAC,IAAI,GAAG,EAAE,EAAE,CAAC;QAClB,GAAG,CAAC,KAAK,GAAG,mBAAmB,CAAC;QAChC,OAAO,KAAK,CAAC;IACf,CAAC;IACD,GAAG,CAAC,IAAI,GAAG,EAAE,CAAC;IACd,IAAI,GAAG,CAAC,IAAI,GAAG,CAAC;QACd,GAAG,CAAC,IAAI,GAAG,EAAE,GAAG,GAAG,CAAC,IAAI,CAAC;IAC3B,IAAI,GAAG,CAAC,IAAI,GAAG,EAAE,EAAE,CAAC;QAClB,GAAG,CAAC,KAAK,GAAG,mBAAmB,CAAC;QAChC,OAAO,KAAK,CAAC;IACf,CAAC;IACD,OAAO,IAAI,CAAC;AACd,CAAC;AAED,SAAS,QAAQ,CAAC,GAAiB;IACjC,MAAM,CAAC,GAAG,GAAG,CAAC,MAAO,CAAC;IACtB,MAAM,GAAG,GAAG,GAAG,CAAC,GAAG,CAAC;IACpB,MAAM,KAAK,GAAG,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,GAAG,CAAC,GAAI,CAAC,MAAM,EAAE,GAAG,CAAC,GAAI,CAAC,MAAM,CAAC,CAAC;IAE7D,6DAA6D;IAC7D,CAAC,CAAC,EAAE,GAAG,IAAI,KAAK,CAAC,KAAK,GAAG,CAAC,CAAC,CAAC;IAC5B,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,IAAI,KAAK,EAAE,CAAC,EAAE;QAC7B,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,GAAG,IAAI,YAAY,CAAC,GAAG,GAAG,CAAC,CAAC,CAAC;IAEtC,2DAA2D;IAC3D,CAAC,CAAC,EAAE,GAAG,IAAI,KAAK,CAAC,GAAG,GAAG,CAAC,CAAC,CAAC;IAC1B,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,IAAI,GAAG,EAAE,CAAC,EAAE;QAC3B,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,GAAG,IAAI,YAAY,CAAC,GAAG,GAAG,CAAC,CAAC,CAAC;IAEtC,CAAC,CAAC,GAAG,GAAG,IAAI,YAAY,CAAC,GAAG,GAAG,CAAC,CAAC,CAAC;IAClC,CAAC,CAAC,IAAI,GAAG,IAAI,YAAY,CAAC,GAAG,GAAG,CAAC,CAAC,CAAC;IACnC,CAAC,CAAC,IAAI,GAAG,IAAI,YAAY,CAAC,GAAG,GAAG,CAAC,CAAC,CAAC;IACnC,CAAC,CAAC,IAAI,GAAG,IAAI,UAAU,CAAC,GAAG,GAAG,CAAC,CAAC,CAAC;IAEjC,0BAA0B;IAC1B,CAAC,CAAC,IAAI,GAAG,IAAI,KAAK,CAAC,EAAE,CAAC,CAAC;IACvB,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,EAAE,EAAE,CAAC,EAAE;QACzB,CAAC,CAAC,IAAI,CAAC,CAAC,CAAC,GAAG,IAAI,YAAY,CAAC,EAAE,CAAC,CAAC;IACnC,CAAC,CAAC,KAAK,GAAG,IAAI,KAAK,CAAC,EAAE,CAAC,CAAC;IACxB,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,EAAE,EAAE,CAAC,EAAE;QACzB,CAAC,CAAC,KAAK,CAAC,CAAC,CAAC,GAAG,IAAI,YAAY,CAAC,CAAC,CAAC,CAAC;IAEnC,OAAO,IAAI,CAAC;AACd,CAAC;AAED,uBAAuB;AAEvB;;;GAGG;AACH,MAAM,UAAU,YAAY,CAAC,GAAiB,EAAE,GAAa;IAC3D,GAAG,CAAC,MAAM,GAAG,IAAI,WAAW,EAAE,CAAC;IAC/B,GAAG,CAAC,GAAG,GAAG,GAAG,CAAC;IACd,IAAI,CAAC,QAAQ,CAAC,GAAG,EAAE,GAAG,CAAC;QAAE,OAAO,KAAK,CAAC;IACtC,OAAO,QAAQ,CAAC,GAAG,CAAC,CAAC;AACvB,CAAC;AAED;;;GAGG;AACH,MAAM,UAAU,UAAU,CAAC,GAAiB;IAC1C,MAAM,CAAC,GAAG,GAAG,CAAC,MAAO,CAAC;IACtB,wBAAwB;IACxB,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC;IAAC,CAAC,CAAC,EAAE,GAAG,CAAC,CAAC;IAAC,CAAC,CAAC,EAAE,GAAG,CAAC,CAAC;IAAC,CAAC,CAAC,EAAE,GAAG,CAAC,CAAC;IAAC,CAAC,CAAC,GAAG,GAAG,CAAC,CAAC;IAAC,CAAC,CAAC,MAAM,GAAG,CAAC,CAAC;IAC/D,CAAC,CAAC,KAAK,GAAG,CAAC,CAAC;IAAC,CAAC,CAAC,IAAI,GAAG,CAAC,CAAC;IAAC,CAAC,CAAC,IAAI,GAAG,CAAC,CAAC;IAAC,CAAC,CAAC,KAAK,GAAG,CAAC,CAAC;IAAC,CAAC,CAAC,MAAM,GAAG,CAAC,CAAC;IAC/D,CAAC,CAAC,KAAK,GAAG,CAAC,CAAC;IAAC,CAAC,CAAC,IAAI,GAAG,CAAC,CAAC;IAAC,CAAC,CAAC,IAAI,GAAG,CAAC,CAAC;IAAC,CAAC,CAAC,MAAM,GAAG,CAAC,CAAC;IAAC,CAAC,CAAC,MAAM,GAAG,CAAC,CAAC;IAChE,CAAC,CAAC,KAAK,GAAG,CAAC,CAAC;IAAC,CAAC,CAAC,KAAK,GAAG,CAAC,CAAC;IAAC,CAAC,CAAC,KAAK,GAAG,CAAC,CAAC;IAAC,CAAC,CAAC,MAAM,GAAG,CAAC,CAAC;IAAC,CAAC,CAAC,IAAI,GAAG,CAAC,CAAC;IAChE,CAAC,CAAC,IAAI,GAAG,CAAC,CAAC;IAAC,CAAC,CAAC,KAAK,GAAG,CAAC,CAAC;IAAC,CAAC,CAAC,EAAE,GAAG,CAAC,CAAC;IAAC,CAAC,CAAC,GAAG,GAAG,CAAC,CAAC;IAAC,CAAC,CAAC,GAAG,GAAG,CAAC,CAAC;IACxD,CAAC,CAAC,GAAG,GAAG,CAAC,CAAC;IAAC,CAAC,CAAC,GAAG,GAAG,CAAC,CAAC;IAAC,CAAC,CAAC,GAAG,GAAG,CAAC,CAAC;IAAC,CAAC,CAAC,KAAK,GAAG,CAAC,CAAC;IAC7C,CAAC,CAAC,EAAE,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC;IACb,cAAc;IACd,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,CAAC,EAAE,CAAC,MAAM,EAAE,CAAC,EAAE;QAAE,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC;IACtD,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,CAAC,EAAE,CAAC,MAAM,EAAE,CAAC,EAAE;QAAE,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC;IACtD,CAAC,CAAC,GAAG,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC;IAAC,CAAC,CAAC,IAAI,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC;IAAC,CAAC,CAAC,IAAI,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC;IAAC,CAAC,CAAC,IAAI,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC;IAC9D,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,CAAC,IAAI,CAAC,MAAM,EAAE,CAAC,EAAE;QAAE,CAAC,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC;IAC1D,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,CAAC,KAAK,CAAC,MAAM,EAAE,CAAC,EAAE;QAAE,CAAC,CAAC,KAAK,CAAC,CAAC,CAAC,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC;IAC5D,GAAG,CAAC,KAAK,GAAG,CAAC,CAAC;IACd,GAAG,CAAC,KAAK,GAAG,IAAI,CAAC;AACnB,CAAC;AAED;;GAEG;AACH,MAAM,UAAU,SAAS,CAAC,GAAiB;IACzC,IAAI,GAAG,CAAC,KAAK;QACX,OAAO,CAAC,KAAK,CAAC,4BAA4B,GAAG,CAAC,KAAK,EAAE,CAAC,CAAC;IAEzD,GAAG,CAAC,MAAM,GAAG,IAAI,CAAC;AACpB,CAAC;AAED;;;;GAIG;AACH,MAAM,UAAU,KAAK,CACnB,GAAiB,EAAE,CAAe,EAAE,GAAW,EAAE,IAAY;IAE7D,MAAM,CAAC,GAAG,GAAG,CAAC,MAAO,CAAC;IACtB,MAAM,GAAG,GAAG,GAAG,CAAC,GAAI,CAAC;IACrB,MAAM,GAAG,GAAG,GAAG,CAAC,GAAG,CAAC;IACpB,IAAI,CAAC,GAAG,GAAG,CAAC;IAEZ,IAAI,KAAa,CAAC;IAClB,IAAI,MAAM,GAAW,CAAC,CAAC;IACvB,IAAI,IAAI,GAAG,CAAC,CAAC;IAEb,qCAAqC;IACrC,SAAS,WAAW,CAAC,GAAW;QAC9B,GAAG,CAAC,KAAK,GAAG,GAAG,CAAC;QAChB,GAAG,CAAC,KAAK,GAAG,CAAC,CAAC,CAAC;QACf,OAAO,EAAC,CAAC,EAAE,KAAK,EAAE,GAAG,CAAC,KAAK,EAAC,CAAC;IAC/B,CAAC;IAED,SAAS,WAAW,CAAC,IAAY,EAAE,GAAW;QAC5C,GAAG,CAAC,KAAK,GAAG,GAAG,CAAC;QAChB,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,IAAI,GAAG,EAAE,CAAC,EAAE;YAC3B,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;QACpB,CAAC,GAAG,CAAC,CAAC,EAAE,CAAC;QACT,GAAG,CAAC,KAAK,GAAG,IAAI,CAAC;QACjB,OAAO,EAAC,CAAC,EAAE,KAAK,EAAE,GAAG,CAAC,KAAK,EAAC,CAAC;IAC/B,CAAC;IAED,SAAS,aAAa;QACpB,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,IAAI,GAAG,EAAE,CAAC,EAAE;YAC3B,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;QACpB,CAAC,GAAG,CAAC,CAAC,EAAE,CAAC;QACT,IAAI,KAAK,KAAK,CAAC,IAAI,KAAK,KAAK,CAAC,EAAE,CAAC;YAC/B,IAAI,IAAI;gBACN,CAAC,GAAG,KAAK,CAAC;QACd,CAAC;QACD,GAAG,CAAC,KAAK,GAAG,CAAC,CAAC;QACd,OAAO,EAAC,CAAC,EAAE,KAAK,EAAE,GAAG,CAAC,KAAK,EAAC,CAAC;IAC/B,CAAC;IAED,SAAS,WAAW;QAClB,MAAM,KAAK,GAAG,KAAK,CAAC,GAAG,EAAE,IAAI,EAAE,CAAC,EAAE,CAAC,CAAC,CAAC;QACrC,IAAI,KAAK,KAAK,CAAC,EAAE,CAAC;YAChB,GAAG,CAAC,KAAK,GAAG,uCAAuC,KAAK,YAAY,IAAI,EAAE,CAAC;YAC3E,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,IAAI,GAAG,EAAE,CAAC,EAAE;gBAC3B,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;YACpB,CAAC,GAAG,CAAC,CAAC,EAAE,CAAC;QACX,CAAC;QACD,CAAC,GAAG,IAAI,CAAC;QACT,GAAG,CAAC,KAAK,GAAG,CAAC,CAAC;QACd,OAAO,EAAC,CAAC,EAAE,KAAK,EAAE,GAAG,CAAC,KAAK,EAAC,CAAC;IAC/B,CAAC;IAED,IAAI,CAAC,KAAK,IAAI;QACZ,OAAO,WAAW,CAAC,0DAA0D,CAAC,CAAC;IAGjF,2CAA2C;IAC3C,IAAI,EAAE,GAAW,CAAC,CAAC;IACnB,IAAI,KAAK,GAAW,CAAC,CAAC;IACtB,MAAM,IAAI,GAAG,GAAG,CAAC,IAAI,CAAC;IACtB,MAAM,IAAI,GAAG,GAAG,CAAC,IAAI,CAAC;IAEtB,IAAI,GAAG,CAAC,KAAK,KAAK,CAAC,IAAI,GAAG,CAAC,KAAK,KAAK,CAAC,EAAE,CAAC;QACvC,EAAE,GAAG,GAAG,CAAC,EAAE,CAAC;QACZ,IAAI,GAAG,CAAC,KAAK,KAAK,CAAC,EAAE,CAAC;YACpB,IAAI,CAAC,IAAI,GAAG,CAAC,CAAC,GAAG,EAAE,GAAG,EAAE;gBACtB,OAAO,WAAW,CAAC,kBAAkB,IAAI,eAAe,CAAC,uCAAuC,EAAE,EAAE,CAAC,CAAC;QAC1G,CAAC;IACH,CAAC;IAED,MAAM,KAAK,GAAG,GAAG,CAAC,KAAK,CAAC;IAExB,IAAI,GAAG,CAAC,KAAK,KAAK,CAAC;QACjB,MAAM,GAAG,CAAC,CAAC,CAAC;IAGd,0CAA0C;IAC1C,IAAI,GAAG,CAAC,KAAK,KAAK,CAAC,EAAE,CAAC;QACpB,CAAC,CAAC,IAAI,GAAG,CAAC,CAAC;QACX,CAAC,CAAC,EAAE,GAAG,CAAC,CAAC;QACT,CAAC,CAAC,GAAG,GAAG,CAAC,CAAC;QAEV,IAAI,KAAK,KAAK,CAAC,IAAI,KAAK,KAAK,CAAC,EAAE,CAAC;YAC/B,KAAK,GAAG,GAAG,CAAC,KAAK,CAAC;YAClB,IAAI,CAAC,KAAK,GAAG,IAAI,CAAC,GAAG,CAAC,IAAI,GAAG,CAAC,CAAC,GAAG,EAAE;gBAClC,OAAO,WAAW,CAAC,8CAA8C,CAAC,CAAC;YAErE,IAAI,EAAE,KAAK,EAAE,IAAI,CAAC,CAAC,GAAG,EAAE,GAAG,KAAK,CAAC,GAAG,EAAE,GAAG,EAAE;gBACzC,EAAE,GAAG,KAAK,GAAG,CAAC,CAAC;QACnB,CAAC;QAED,MAAM,GAAG,CAAC,CAAC;QACX,CAAC,CAAC,EAAE,GAAG,CAAC,CAAC;QAET,oBAAoB;QACpB,GAAG,CAAC,IAAI,CAAC,CAAC,EAAE,CAAC,CAAC,QAAQ,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,QAAQ,CAAC,CAAC,CAAC,EAAE,GAAG,CAAC,IAAI,CAAC,CAAC;QAC1D,CAAC,CAAC,GAAG,GAAG,CAAC,CAAC;QAEV,4BAA4B;QAC5B,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,IAAI,GAAG,EAAE,CAAC,EAAE;YAC3B,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC;QAEpB,sBAAsB;QACtB,KAAK,CAAC,GAAG,EAAE,CAAC,CAAC,CAAC;QACd,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,IAAI,GAAG,EAAE,CAAC,EAAE,EAAE,CAAC;YAC9B,IAAI,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,IAAI,EAAE;gBAChB,OAAO,WAAW,CAAC,eAAe,CAAC,OAAO,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,QAAQ,CAAC,CAAC;QAChE,CAAC;QAED,4CAA4C;QAC5C,IAAI,EAAE,KAAK,EAAE,EAAE,CAAC;YACd,MAAM,KAAK,GAAG,IAAI,CAAC,GAAG,CAAC,IAAI,GAAG,CAAC,CAAC,CAAC;YACjC,MAAM,EAAE,GAAG,IAAI,CAAC,GAAG,CAAC,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,EAAE,IAAI,CAAC,GAAG,CAAC,IAAI,CAAC,CAAC,CAAC;YACjD,IAAI,KAAK,GAAG,EAAE,GAAG,GAAG,GAAG,EAAE;gBACvB,OAAO,WAAW,CAAC,kDAAkD,CAAC,CAAC;YAEzE,IAAI,GAAG,GAAG,EAAE,CAAC;YACb,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,IAAI,GAAG,EAAE,CAAC,EAAE;gBAC3B,GAAG,GAAG,IAAI,CAAC,GAAG,CAAC,GAAG,EAAE,IAAI,CAAC,CAAC,CAAC,CAAC,CAAC;YAC/B,IAAI,GAAG,IAAI,EAAE,EAAE,CAAC;gBACd,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,IAAI,GAAG,EAAE,CAAC,EAAE,EAAE,CAAC;oBAC9B,MAAM,KAAK,GAAG,IAAI,CAAC,CAAC,CAAC,CAAC;oBACtB,MAAM,GAAG,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;oBAC3B,IAAI,GAAG,KAAK,EAAE;wBACZ,GAAG,GAAG,IAAI,CAAC,GAAG,CAAC,GAAG,EAAE,KAAK,GAAG,GAAG,CAAC,CAAC;gBACrC,CAAC;YACH,CAAC;YACD,GAAG,GAAG,IAAI,CAAC,GAAG,CAAC,GAAG,EAAE,IAAI,GAAG,GAAG,CAAC,CAAC;YAChC,GAAG,GAAG,IAAI,CAAC,GAAG,CAAC,GAAG,EAAE,KAAK,CAAC,CAAC;YAC3B,MAAM,GAAG,GAAG,MAAM,CAAC,GAAG,EAAE,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,GAAG,CAAC,CAAC;YACxC,MAAM,KAAK,GAAG,EAAE,GAAG,CAAC,GAAG,GAAG,EAAE,GAAG,EAAE,CAAC,GAAG,GAAG,GAAG,GAAG,GAAG,GAAG,CAAC;YACrD,EAAE,GAAG,EAAE,GAAG,IAAI,CAAC,IAAI,CAAC,KAAK,CAAC,CAAC;YAC3B,EAAE,GAAG,IAAI,CAAC,GAAG,CAAC,EAAE,EAAE,KAAK,CAAC,CAAC;YACzB,EAAE,GAAG,EAAE,GAAG,CAAC,CAAC,IAAI,GAAG,CAAC,IAAI,EAAE,CAAC,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC;QAC1C,CAAC;QAED,qBAAqB;QACrB,MAAM,EAAE,GAAG,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,GAAG,CAAC,IAAI,CAAC;QACnC,IAAI,EAAE,GAAG,EAAE;YAAE,EAAE,IAAI,EAAE,CAAC;QAEtB,yBAAyB;QACzB,CAAC,CAAC,CAAC,GAAG,EAAE,CAAC;QACT,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,IAAI,GAAG,EAAE,CAAC,EAAE;YAC3B,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,IAAI,EAAE,CAAC;IACrB,CAAC;IAED,gDAAgD;IAChD,IAAI,GAAG,CAAC,KAAK,KAAK,CAAC,IAAI,GAAG,CAAC,KAAK,KAAK,CAAC,EAAE,CAAC;QACvC,MAAM,GAAG,CAAC,CAAC;QACX,CAAC,CAAC,MAAM,GAAG,CAAC,CAAC,GAAG,CAAC;QAEjB,QAAQ,KAAK,EAAE,CAAC;YAChB,KAAK,CAAC;gBACJ,IAAI,CAAC,CAAC,CAAC,EAAE,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,IAAI,EAAE;oBAC3B,OAAO,WAAW,EAAE,CAAC;gBAEvB,MAAM;YACR,KAAK,CAAC;gBACJ,MAAM;YACR,KAAK,CAAC,CAAC,CAAC,CAAC;gBACP,MAAM,EAAE,GAAG,CAAC,CAAC,EAAE,GAAG,CAAC,CAAC,EAAE,GAAG,CAAC,EAAE,GAAG,IAAI,GAAG,GAAG,CAAC,CAAC;gBAC3C,IAAI,CAAC,EAAE,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,GAAG,EAAE;oBACxB,OAAO,WAAW,CAAC,mBAAmB,KAAK,4BAA4B,CAAC,CAAC;gBAE3E,IAAI,CAAC,CAAC,CAAC,EAAE,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,GAAG,EAAE;oBAAE,MAAM;gBACpC,OAAO,aAAa,EAAE,CAAC;YACzB,CAAC;YACD,KAAK,CAAC;gBACJ,KAAK,GAAG,GAAG,CAAC,KAAK,CAAC;gBAClB,IAAI,CAAC,CAAC,CAAC,EAAE,GAAG,KAAK,CAAC,GAAG,CAAC,CAAC,CAAC,GAAG,EAAE;oBAC3B,OAAO,WAAW,CAAC,8CAA8C,CAAC,CAAC;gBAErE,IAAI,CAAC,KAAK,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,GAAG,EAAE;oBAC3B,OAAO,WAAW,CAAC,8CAA8C,CAAC,CAAC;gBAErE,IAAI,CAAC,CAAC,CAAC,EAAE,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,IAAI,EAAE;oBAC3B,OAAO,WAAW,EAAE,CAAC;YAEvB,+BAA+B;YAC/B,0CAA0C;YAC5C,KAAK,CAAC;gBACJ,IAAI,KAAK,KAAK,CAAC,EAAE,CAAC;oBAChB,KAAK,GAAG,GAAG,CAAC,KAAK,CAAC;oBAClB,IAAI,CAAC,CAAC,CAAC,EAAE,GAAG,KAAK,CAAC,GAAG,CAAC,CAAC,CAAC,GAAG,EAAE;wBAC3B,OAAO,WAAW,CAAC,8CAA8C,CAAC,CAAC;gBACvE,CAAC;gBACD,CAAC;oBACC,MAAM,GAAG,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,EAAE,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;oBAC3C,IAAI,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,EAAE,GAAG,KAAK,CAAC,IAAI,CAAC,IAAI,GAAG,GAAG,GAAG,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;oBAC5D,IAAI,IAAI,EAAE,CAAC;wBACT,CAAC,GAAG,KAAK,CAAC;wBACV,OAAO,aAAa,EAAE,CAAC;oBACzB,CAAC;oBACD,MAAM,KAAK,GAAG,CAAC,CAAC,EAAE,GAAG,CAAC,CAAC,CAAC,GAAG,CAAC,EAAE,GAAG,EAAE,GAAG,GAAG,CAAC,CAAC;oBAC3C,IAAI,CAAC,KAAK,GAAG,KAAK,CAAC,GAAG,CAAC,CAAC,CAAC,IAAI,EAAE;wBAAE,MAAM;oBACvC,CAAC,CAAC,CAAC,GAAG,CAAC,KAAK,GAAG,CAAC,CAAC,EAAE,CAAC,GAAG,CAAC,EAAE,GAAG,EAAE,GAAG,GAAG,CAAC,CAAC;oBACvC,IAAI,GAAG,CAAC,KAAK,KAAK,CAAC;wBAAE,MAAM,GAAG,CAAC,CAAC,CAAC;gBACnC,CAAC;gBACD,MAAM;QACR,CAAC;IACH,CAAC;IAED,iCAAiC;IACjC,OAAO,IAAI,EAAE,CAAC;QACZ,IAAI,GAAG,CAAC,KAAK,KAAK,CAAC,IAAI,CAAC,CAAC,GAAG,KAAK,CAAC,EAAE,CAAC;YACnC,IAAI,CAAC,CAAC,CAAC,GAAG,GAAG,CAAC,CAAC,MAAM,CAAC,IAAI,GAAG,CAAC,MAAM;gBAClC,OAAO,WAAW,CAAC,CAAC,CAAC,EAAE,WAAW,GAAG,CAAC,MAAM,mCAAmC,CAAC,CAAC;YAEnF,KAAK,CAAC,GAAG,EAAE,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC;YACpB,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,IAAI,GAAG,EAAE,CAAC,EAAE,EAAE,CAAC;gBAC9B,IAAI,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,IAAI,EAAE;oBAChB,OAAO,WAAW,CAAC,CAAC,CAAC,EAAE,eAAe,CAAC,OAAO,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,QAAQ,CAAC,CAAC;YACpE,CAAC;QACH,CAAC;QAED,MAAM,KAAK,GAAG,GAAG,GAAG,MAAM,CAAC,GAAG,EAAE,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,GAAG,CAAC,CAAC;QAChD,IAAI,KAAK,GAAG,IAAI,EAAE,CAAC;YACjB,MAAM,MAAM,GAAG,KAAK,GAAG,IAAI,CAAC;YAC5B,IAAI,CAAC,CAAC,GAAG,KAAK,CAAC,EAAE,CAAC;gBAChB,OAAO,WAAW,CAChB,2DAA2D;oBAC3D,yDAAyD,MAAM,EAAE,CAClE,CAAC;YACJ,CAAC;YACD,OAAO,WAAW,CAAC,CAAC,CAAC,EACnB,mBAAmB,CAAC,+BAA+B;gBACnD,yDAAyD,MAAM,EAAE,CAClE,CAAC;QACJ,CAAC;QAED,IAAI,CAAC,CAAC,CAAC,EAAE,GAAG,CAAC,CAAC,CAAC,CAAC,KAAK,CAAC,CAAC,EAAE,EAAE,CAAC;YAC1B,CAAC,CAAC,KAAK,EAAE,CAAC;YACV,IAAI,CAAC,CAAC,KAAK,IAAI,GAAG,CAAC,MAAM,EAAE,CAAC;gBAC1B,OAAO,CAAC,KAAK,CAAC,kCAAkC,CAAC,CAAC,EAAE,YAAY,CAAC,CAAC,CAAC,MAAM;oBACvE,uCAAuC,CAAC,CAAC;gBAC3C,IAAI,CAAC,CAAC,KAAK,KAAK,GAAG,CAAC,MAAM,EAAE,CAAC;oBAC3B,OAAO,CAAC,KAAK,CAAC,0CAA0C,CAAC,CAAC,KAAK,EAAE;wBAC/D,sDAAsD,CAAC,CAAC;gBAC5D,CAAC;YACH,CAAC;QACH,CAAC;QAED,aAAa;QACb,KAAK,GAAG,KAAK,CAAC,GAAG,EAAE,CAAC,EAAE,MAAM,CAAC,CAAC;QAE9B,IAAI,KAAK,KAAK,CAAC,EAAE,CAAC;YAChB,wCAAwC;YACxC,iCAAiC;YACjC,IAAI,GAAG,CAAC,SAAS;gBACf,GAAG,CAAC,SAAS,CAAC,IAAI,CAAC,eAAe,CAAC,CAAC,EAAE,GAAG,CAAC,CAAC,CAAC;YAE9C,MAAM,GAAG,CAAC,CAAC;YACX,IAAI,CAAC,CAAC,IAAI,KAAK,CAAC,CAAC,KAAK,EAAE,CAAC;gBACvB,CAAC,CAAC,GAAG,GAAG,CAAC,CAAC,EAAE,CAAC;gBACb,MAAM,GAAG,CAAC,CAAC,CAAC;gBACZ,IAAI,GAAG,CAAC,IAAI,EAAE,CAAC;oBACb,IAAI,CAAC,CAAC,IAAI,KAAK,CAAC,EAAE,CAAC;wBACjB,OAAO,CAAC,KAAK,CAAC,mDAAmD;4BAC/D,WAAW,CAAC,CAAC,EAAE,6BAA6B,CAAC,CAAC,CAAC,gBAAgB,CAAC,CAAC,GAAG,EAAE,CAAC,CAAC;oBAC5E,CAAC;oBACD,IAAI,CAAC,CAAC,IAAI,KAAK,CAAC,EAAE,CAAC;wBACjB,OAAO,CAAC,KAAK,CAAC,sDAAsD;4BAClE,WAAW,CAAC,CAAC,EAAE,6BAA6B,CAAC,CAAC,CAAC,gBAAgB,CAAC,CAAC,GAAG,EAAE,CAAC,CAAC;oBAC5E,CAAC;gBACH,CAAC;YACH,CAAC;YAED,IAAI,KAAK,KAAK,CAAC,EAAE,CAAC;gBAChB,IAAI,CAAC,CAAC,CAAC,EAAE,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,GAAG,EAAE;oBAAE,SAAS;gBACvC,OAAO,WAAW,EAAE,CAAC;YACvB,CAAC;YACD,IAAI,KAAK,KAAK,CAAC;gBACb,OAAO,aAAa,EAAE,CAAC;YAEzB,IAAI,KAAK,KAAK,CAAC,EAAE,CAAC;gBAChB,IAAI,CAAC,CAAC,CAAC,EAAE,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,IAAI,EAAE;oBAC3B,OAAO,aAAa,EAAE,CAAC;gBAEzB,SAAS;YACX,CAAC;YACD,IAAI,KAAK,KAAK,CAAC,EAAE,CAAC;gBAChB,KAAK,GAAG,GAAG,CAAC,KAAK,CAAC;gBAClB,IAAI,CAAC,CAAC,CAAC,EAAE,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,IAAI,EAAE;oBAC3B,OAAO,WAAW,EAAE,CAAC;qBAClB,CAAC;oBACJ,MAAM,GAAG,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,EAAE,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;oBAC3C,IAAI,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,EAAE,GAAG,KAAK,CAAC,IAAI,CAAC,IAAI,GAAG,GAAG,GAAG,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;oBAC5D,IAAI,IAAI;wBACN,OAAO,aAAa,EAAE,CAAC;oBAEzB,MAAM,KAAK,GAAG,CAAC,CAAC,EAAE,GAAG,CAAC,CAAC,CAAC,GAAG,CAAC,EAAE,GAAG,EAAE,GAAG,GAAG,CAAC,CAAC;oBAC3C,IAAI,CAAC,KAAK,GAAG,KAAK,CAAC,GAAG,CAAC,CAAC,CAAC,IAAI,EAAE;wBAAE,SAAS;oBAC1C,CAAC,CAAC,CAAC,GAAG,CAAC,KAAK,GAAG,CAAC,CAAC,EAAE,CAAC,GAAG,CAAC,EAAE,GAAG,EAAE,GAAG,GAAG,CAAC,CAAC;oBACvC,MAAM,GAAG,CAAC,CAAC,CAAC;oBACZ,SAAS;gBACX,CAAC;YACH,CAAC;YACD,IAAI,KAAK,KAAK,CAAC,EAAE,CAAC;gBAChB,KAAK,GAAG,GAAG,CAAC,KAAK,CAAC;gBAClB,MAAM,GAAG,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,EAAE,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;gBAC3C,IAAI,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,EAAE,GAAG,KAAK,CAAC,IAAI,CAAC,IAAI,GAAG,GAAG,GAAG,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;gBAC5D,OAAO,aAAa,EAAE,CAAC;YACzB,CAAC;QACH,CAAC;QAED,uDAAuD;QACvD,IAAI,KAAK,KAAK,CAAC,CAAC,IAAI,KAAK,KAAK,CAAC,CAAC,EAAE,CAAC;YACjC,IAAI,GAAG,GAAG,EAAE,CAAC;YACb,CAAC,CAAC,KAAK,GAAG,CAAC,CAAC;YACZ,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,IAAI,GAAG,EAAE,CAAC,EAAE,EAAE,CAAC;gBAC9B,MAAM,IAAI,GAAG,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC;gBAC5C,IAAI,GAAG,GAAG,IAAI,EAAE,CAAC;oBACf,GAAG,GAAG,IAAI,CAAC;oBACX,CAAC,CAAC,KAAK,GAAG,CAAC,CAAC;gBACd,CAAC;YACH,CAAC;YACD,IAAI,KAAK,KAAK,CAAC,CAAC,EAAE,CAAC;gBACjB,OAAO,WAAW,CAAC,CAAC,CAAC,EACnB,mBAAmB,CAAC,CAAC,EAAE,sBAAsB,CAAC,CAAC,CAAC,GAAG;oBACnD,yDAAyD,CAC1D,CAAC;YACJ,CAAC;YACD,IAAI,KAAK,KAAK,CAAC,CAAC,EAAE,CAAC;gBACjB,OAAO,WAAW,CAAC,CAAC,CAAC,EACnB,mBAAmB,CAAC,CAAC,EAAE,sBAAsB,CAAC,CAAC,CAAC,GAAG;oBACnD,oEAAoE,CACrE,CAAC;YACJ,CAAC;QACH,CAAC;IACH,CAAC;AACH,CAAC;AAcD;;GAEG;AACH,MAAM,OAAO,KAAK;IAKhB,YACE,CAAc,EAAE,GAAW,EAC3B,GAA4C,EAAE,IAAU;QAJlD,iBAAY,GAAY,KAAK,CAAC;QAMpC,8DAA8D;QAC9D,MAAM,WAAW,GAAc,CAAC,CAAC,EAAE,CAAC,EAAE,IAAI,EAAE,IAAI,EAAE,EAAE;YAClD,OAAO,CAAC,CAAC,CAAC,EAAE,CAAC,EAAE,IAAI,EAAE,IAAI,CAAC,CAAC;QAC7B,CAAC,CAAC;QAEF,IAAI,CAAC,GAAG,GAAG,IAAI,YAAY,CAAC,WAAW,EAAE,GAAG,EAAE,IAAI,CAAC,CAAC;QACpD,IAAI,CAAC,SAAS,GAAG,IAAI,YAAY,CAAC,GAAG,GAAG,CAAC,CAAC,CAAC,CAAC,YAAY;QAExD,mCAAmC;QACnC,MAAM,OAAO,mBACX,IAAI,EAAE,CAAC,EAAE,MAAM,EAAE,CAAC,EAAE,MAAM,EAAE,CAAC,EAAE,MAAM,EAAE,CAAC,EAAE,MAAM,EAAE,CAAC,EACnD,KAAK,EAAE,CAAC,EAAE,EAAE,EAAE,CAAC,EAAE,IAAI,EAAE,CAAC,EAAE,IAAI,EAAE,CAAC,EAAE,IAAI,EAAE,CAAC,EAAE,KAAK,EAAE,CAAC,EACpD,IAAI,EAAE,IAAI,YAAY,CAAC,GAAG,GAAG,CAAC,CAAC,EAC/B,IAAI,EAAE,IAAI,YAAY,CAAC,GAAG,GAAG,CAAC,CAAC,IAC5B,GAAG,CACP,CAAC;QAEF,8EAA8E;QAC9E,IAAI,GAAG,aAAH,GAAG,uBAAH,GAAG,CAAE,IAAI,EAAE,CAAC;YACd,IAAI,GAAG,CAAC,IAAI,CAAC,MAAM,KAAK,GAAG,EAAE,CAAC;gBAC5B,4CAA4C;gBAC5C,MAAM,CAAC,GAAG,IAAI,YAAY,CAAC,GAAG,GAAG,CAAC,CAAC,CAAC;gBACpC,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,GAAG,EAAE,CAAC,EAAE;oBAAE,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,GAAG,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC;gBACrD,OAAO,CAAC,IAAI,GAAG,CAAC,CAAC;YACnB,CAAC;;gBACC,OAAO,CAAC,IAAI,GAAG,GAAG,CAAC,IAAI,CAAC;QAC5B,CAAC;QACD,IAAI,GAAG,aAAH,GAAG,uBAAH,GAAG,CAAE,IAAI,EAAE,CAAC;YACd,IAAI,GAAG,CAAC,IAAI,CAAC,MAAM,KAAK,GAAG,EAAE,CAAC;gBAC5B,MAAM,CAAC,GAAG,IAAI,YAAY,CAAC,GAAG,GAAG,CAAC,CAAC,CAAC;gBACpC,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,GAAG,EAAE,CAAC,EAAE;oBAAE,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,GAAG,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC;gBACrD,OAAO,CAAC,IAAI,GAAG,CAAC,CAAC;YACnB,CAAC;;gBACC,OAAO,CAAC,IAAI,GAAG,GAAG,CAAC,IAAI,CAAC;QAC5B,CAAC;QAED,YAAY,CAAC,IAAI,CAAC,GAAG,EAAE,OAAO,CAAC,CAAC;QAEhC,IAAI,GAAG,aAAH,GAAG,uBAAH,GAAG,CAAE,KAAK,EAAE,CAAC;YACf,IAAI,CAAC,YAAY,GAAG,IAAI,CAAC;YACzB,IAAI,CAAC,GAAG,CAAC,SAAS,GAAG,EAAE,CAAC;QAC1B,CAAC;IACH,CAAC;IAED;;;;OAIG;IACH,KAAK,CAAC,CAAoB,EAAE,CAAS,EAAE,IAAY;QACjD,MAAM,GAAG,GAAG,IAAI,CAAC,GAAG,CAAC,GAAG,CAAC;QACzB,kDAAkD;QAClD,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,GAAG,EAAE,CAAC,EAAE;YAC1B,IAAI,CAAC,SAAS,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC;QAE/B,MAAM,MAAM,GAAG,KAAK,CAAC,IAAI,CAAC,GAAG,EAAE,IAAI,CAAC,SAAS,EAAE,CAAC,EAAE,IAAI,CAAC,CAAC;QAExD,uBAAuB;QACvB,MAAM,GAAG,GAAG,IAAI,YAAY,CAAC,GAAG,CAAC,CAAC;QAClC,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,GAAG,EAAE,CAAC,EAAE;YAC1B,GAAG,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,SAAS,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC;QAEjC,OAAO,EAAC,CAAC,EAAE,GAAG,EAAE,CAAC,EAAE,MAAM,CAAC,CAAC,EAAC,CAAC;IAC/B,CAAC;IAED,IAAI,KAAK,KAAY,OAAO,IAAI,CAAC,GAAG,CAAC,KAAK,CAAC,CAAA,CAAC;IAC5C,IAAI,KAAK,KAAmB,OAAO,IAAI,CAAC,GAAG,CAAC,KAAK,CAAC,CAAA,CAAC;IAEnD,mEAAmE;IACnE,cAAc;QACZ,IAAI,CAAC,IAAI,CAAC,GAAG,CAAC,SAAS,IAAI,IAAI,CAAC,GAAG,CAAC,SAAS,CAAC,MAAM,KAAK,CAAC;YACxD,MAAM,IAAI,KAAK,CAAC,sFAAsF,CAAC,CAAC;QAC1G,OAAO,IAAI,WAAW,CAAC,IAAI,CAAC,GAAG,CAAC,SAAS,EAAE,IAAI,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC;IAC3D,CAAC;IAED,KAAK;QACH,UAAU,CAAC,IAAI,CAAC,GAAG,CAAC,CAAC;QACrB,IAAI,IAAI,CAAC,YAAY;YACnB,IAAI,CAAC,GAAG,CAAC,SAAS,GAAG,EAAE,CAAC;IAC5B,CAAC;CACF"}