universal-physics-tensor 0.44.1 → 0.44.3
This diff represents the content of publicly available package versions that have been released to one of the supported registries. The information contained in this diff is provided for informational purposes only and reflects changes between package versions as they appear in their respective public registries.
- package/README.md +51 -164
- package/dist/bridges/bridge-equations.d.ts +28 -6
- package/dist/bridges/bridge-equations.d.ts.map +1 -1
- package/dist/bridges/bridge-equations.js +30 -8
- package/dist/bridges/bridge-equations.js.map +1 -1
- package/dist/cli/commands/index.d.ts +1 -0
- package/dist/cli/commands/index.d.ts.map +1 -1
- package/dist/cli/commands/index.js +1 -0
- package/dist/cli/commands/index.js.map +1 -1
- package/dist/cli/commands/probe.d.ts +3 -0
- package/dist/cli/commands/probe.d.ts.map +1 -0
- package/dist/cli/commands/probe.js +169 -0
- package/dist/cli/commands/probe.js.map +1 -0
- package/dist/cli/main.d.ts.map +1 -1
- package/dist/cli/main.js +5 -0
- package/dist/cli/main.js.map +1 -1
- package/dist/cli-api.d.ts +1 -0
- package/dist/cli-api.d.ts.map +1 -1
- package/dist/cli-api.js +3 -0
- package/dist/cli-api.js.map +1 -1
- package/dist/composition/probe/backend-protocol.d.ts +31 -0
- package/dist/composition/probe/backend-protocol.d.ts.map +1 -0
- package/dist/composition/probe/backend-protocol.js +83 -0
- package/dist/composition/probe/backend-protocol.js.map +1 -0
- package/dist/composition/probe/candidate-store.d.ts +33 -0
- package/dist/composition/probe/candidate-store.d.ts.map +1 -0
- package/dist/composition/probe/candidate-store.js +110 -0
- package/dist/composition/probe/candidate-store.js.map +1 -0
- package/dist/composition/probe/corpus.d.ts +23 -0
- package/dist/composition/probe/corpus.d.ts.map +1 -0
- package/dist/composition/probe/corpus.js +51 -0
- package/dist/composition/probe/corpus.js.map +1 -0
- package/dist/composition/probe/dataset.d.ts +35 -0
- package/dist/composition/probe/dataset.d.ts.map +1 -0
- package/dist/composition/probe/dataset.js +154 -0
- package/dist/composition/probe/dataset.js.map +1 -0
- package/dist/composition/probe/experiment-design.d.ts +32 -0
- package/dist/composition/probe/experiment-design.d.ts.map +1 -0
- package/dist/composition/probe/experiment-design.js +86 -0
- package/dist/composition/probe/experiment-design.js.map +1 -0
- package/dist/composition/probe/falsify.d.ts +29 -0
- package/dist/composition/probe/falsify.d.ts.map +1 -0
- package/dist/composition/probe/falsify.js +142 -0
- package/dist/composition/probe/falsify.js.map +1 -0
- package/dist/composition/probe/fingerprint.d.ts +20 -0
- package/dist/composition/probe/fingerprint.d.ts.map +1 -0
- package/dist/composition/probe/fingerprint.js +99 -0
- package/dist/composition/probe/fingerprint.js.map +1 -0
- package/dist/composition/probe/fit.d.ts +23 -0
- package/dist/composition/probe/fit.d.ts.map +1 -0
- package/dist/composition/probe/fit.js +92 -0
- package/dist/composition/probe/fit.js.map +1 -0
- package/dist/composition/probe/frontier.d.ts +30 -0
- package/dist/composition/probe/frontier.d.ts.map +1 -0
- package/dist/composition/probe/frontier.js +169 -0
- package/dist/composition/probe/frontier.js.map +1 -0
- package/dist/composition/probe/generator.d.ts +26 -0
- package/dist/composition/probe/generator.d.ts.map +1 -0
- package/dist/composition/probe/generator.js +123 -0
- package/dist/composition/probe/generator.js.map +1 -0
- package/dist/composition/probe/index.d.ts +49 -0
- package/dist/composition/probe/index.d.ts.map +1 -0
- package/dist/composition/probe/index.js +35 -0
- package/dist/composition/probe/index.js.map +1 -0
- package/dist/composition/probe/limits.d.ts +27 -0
- package/dist/composition/probe/limits.d.ts.map +1 -0
- package/dist/composition/probe/limits.js +71 -0
- package/dist/composition/probe/limits.js.map +1 -0
- package/dist/composition/probe/metadata.d.ts +19 -0
- package/dist/composition/probe/metadata.d.ts.map +1 -0
- package/dist/composition/probe/metadata.js +30 -0
- package/dist/composition/probe/metadata.js.map +1 -0
- package/dist/composition/probe/pipeline.d.ts +42 -0
- package/dist/composition/probe/pipeline.d.ts.map +1 -0
- package/dist/composition/probe/pipeline.js +307 -0
- package/dist/composition/probe/pipeline.js.map +1 -0
- package/dist/composition/probe/problem.d.ts +40 -0
- package/dist/composition/probe/problem.d.ts.map +1 -0
- package/dist/composition/probe/problem.js +101 -0
- package/dist/composition/probe/problem.js.map +1 -0
- package/dist/composition/probe/report.d.ts +15 -0
- package/dist/composition/probe/report.d.ts.map +1 -0
- package/dist/composition/probe/report.js +91 -0
- package/dist/composition/probe/report.js.map +1 -0
- package/dist/composition/probe/residual.d.ts +20 -0
- package/dist/composition/probe/residual.d.ts.map +1 -0
- package/dist/composition/probe/residual.js +72 -0
- package/dist/composition/probe/residual.js.map +1 -0
- package/dist/composition/probe/run-manifest.d.ts +24 -0
- package/dist/composition/probe/run-manifest.d.ts.map +1 -0
- package/dist/composition/probe/run-manifest.js +48 -0
- package/dist/composition/probe/run-manifest.js.map +1 -0
- package/dist/composition/probe/scoring.d.ts +19 -0
- package/dist/composition/probe/scoring.d.ts.map +1 -0
- package/dist/composition/probe/scoring.js +55 -0
- package/dist/composition/probe/scoring.js.map +1 -0
- package/dist/composition/probe/search-budget.d.ts +26 -0
- package/dist/composition/probe/search-budget.d.ts.map +1 -0
- package/dist/composition/probe/search-budget.js +36 -0
- package/dist/composition/probe/search-budget.js.map +1 -0
- package/dist/composition/probe/serialize.d.ts +33 -0
- package/dist/composition/probe/serialize.d.ts.map +1 -0
- package/dist/composition/probe/serialize.js +61 -0
- package/dist/composition/probe/serialize.js.map +1 -0
- package/dist/composition/probe/structure.d.ts +51 -0
- package/dist/composition/probe/structure.d.ts.map +1 -0
- package/dist/composition/probe/structure.js +93 -0
- package/dist/composition/probe/structure.js.map +1 -0
- package/dist/composition/probe/types.d.ts +281 -0
- package/dist/composition/probe/types.d.ts.map +1 -0
- package/dist/composition/probe/types.js +23 -0
- package/dist/composition/probe/types.js.map +1 -0
- package/dist/diff/bridge-ast-gradient.d.ts.map +1 -1
- package/dist/diff/bridge-ast-gradient.js +22 -2
- package/dist/diff/bridge-ast-gradient.js.map +1 -1
- package/dist/numerical/float64-engine.d.ts.map +1 -1
- package/dist/numerical/float64-engine.js +22 -8
- package/dist/numerical/float64-engine.js.map +1 -1
- package/dist/numerical/geodesic-integrator.d.ts.map +1 -1
- package/dist/numerical/geodesic-integrator.js +85 -14
- package/dist/numerical/geodesic-integrator.js.map +1 -1
- package/dist/numerical/gl4-integrator.d.ts +2 -0
- package/dist/numerical/gl4-integrator.d.ts.map +1 -1
- package/dist/numerical/gl4-integrator.js +345 -156
- package/dist/numerical/gl4-integrator.js.map +1 -1
- package/dist/numerical/weyl-lowering.d.ts.map +1 -1
- package/dist/numerical/weyl-lowering.js +158 -48
- package/dist/numerical/weyl-lowering.js.map +1 -1
- package/package.json +15 -5
|
@@ -108,6 +108,10 @@ export function solveGL4Stage(state, h, gInverseFn, dgInverseFn, opts) {
|
|
|
108
108
|
// Pre-allocate arrays to hoist dxStage and dpStage out of the `i` loop
|
|
109
109
|
const dxStageArr = [new Float64Array(dim), new Float64Array(dim)];
|
|
110
110
|
const dpStageArr = [new Float64Array(dim), new Float64Array(dim)];
|
|
111
|
+
const dxStageArr0 = dxStageArr[0];
|
|
112
|
+
const dxStageArr1 = dxStageArr[1];
|
|
113
|
+
const dpStageArr0 = dpStageArr[0];
|
|
114
|
+
const dpStageArr1 = dpStageArr[1];
|
|
111
115
|
const hA00 = h * GL4_A[0][0];
|
|
112
116
|
const hA01 = h * GL4_A[0][1];
|
|
113
117
|
const hA10 = h * GL4_A[1][0];
|
|
@@ -125,57 +129,284 @@ export function solveGL4Stage(state, h, gInverseFn, dgInverseFn, opts) {
|
|
|
125
129
|
// 2. Precompute dxStage and dpStage for all j and mu (only 8 combinations)
|
|
126
130
|
const p0_st = P[0];
|
|
127
131
|
const p1_st = P[1];
|
|
128
|
-
|
|
129
|
-
|
|
130
|
-
|
|
131
|
-
|
|
132
|
-
const
|
|
133
|
-
|
|
134
|
-
|
|
135
|
-
|
|
136
|
-
const
|
|
137
|
-
|
|
138
|
-
|
|
139
|
-
|
|
140
|
-
//
|
|
141
|
-
|
|
142
|
-
if (
|
|
143
|
-
|
|
144
|
-
|
|
145
|
-
|
|
146
|
-
|
|
147
|
-
|
|
148
|
-
|
|
149
|
-
|
|
150
|
-
|
|
151
|
-
|
|
152
|
-
|
|
153
|
-
|
|
154
|
-
|
|
155
|
-
|
|
156
|
-
|
|
157
|
-
|
|
132
|
+
if (dim === 4) {
|
|
133
|
+
// Bolt: Pre-cache momentum elements into local scalars to eliminate array lookups.
|
|
134
|
+
// This combined with manual loop unrolling completely eliminates loop overhead
|
|
135
|
+
// for 4D spacetime integrations without losing `g !== 0` bailout sparseness benefits.
|
|
136
|
+
const p00 = p0_st[0], p01 = p0_st[1], p02 = p0_st[2], p03 = p0_st[3];
|
|
137
|
+
const p10 = p1_st[0], p11 = p1_st[1], p12 = p1_st[2], p13 = p1_st[3];
|
|
138
|
+
for (let mu = 0; mu < 4; mu++) {
|
|
139
|
+
let dx0 = 0, dp0 = 0, dx1 = 0, dp1 = 0;
|
|
140
|
+
const mu4 = mu * 4;
|
|
141
|
+
let g0, g1, dg0, dg1;
|
|
142
|
+
let pDotTerm0, pDotTerm1;
|
|
143
|
+
let offset = mu * 16;
|
|
144
|
+
// nu = 0
|
|
145
|
+
g0 = gInvAtX0[mu4];
|
|
146
|
+
if (g0 !== 0)
|
|
147
|
+
dx0 += g0 * p00;
|
|
148
|
+
g1 = gInvAtX1[mu4];
|
|
149
|
+
if (g1 !== 0)
|
|
150
|
+
dx1 += g1 * p10;
|
|
151
|
+
dg0 = dgInvAtX0[offset];
|
|
152
|
+
dg1 = dgInvAtX1[offset];
|
|
153
|
+
// Guarded like every other rho term, and like master. Assigning
|
|
154
|
+
// unguarded here would compute 0 * p when dg is exactly 0, which is
|
|
155
|
+
// NaN for a non-finite momentum -- manufacturing a divergence signal
|
|
156
|
+
// out of a term that contributes nothing.
|
|
157
|
+
pDotTerm0 = 0;
|
|
158
|
+
if (dg0 !== 0) {
|
|
159
|
+
pDotTerm0 = dg0 * p00;
|
|
160
|
+
}
|
|
161
|
+
pDotTerm1 = 0;
|
|
162
|
+
if (dg1 !== 0) {
|
|
163
|
+
pDotTerm1 = dg1 * p10;
|
|
164
|
+
}
|
|
165
|
+
dg0 = dgInvAtX0[offset + 1];
|
|
166
|
+
dg1 = dgInvAtX1[offset + 1];
|
|
167
|
+
if (dg0 !== 0) {
|
|
168
|
+
pDotTerm0 += dg0 * p01;
|
|
169
|
+
}
|
|
170
|
+
if (dg1 !== 0) {
|
|
171
|
+
pDotTerm1 += dg1 * p11;
|
|
172
|
+
}
|
|
173
|
+
dg0 = dgInvAtX0[offset + 2];
|
|
174
|
+
dg1 = dgInvAtX1[offset + 2];
|
|
175
|
+
if (dg0 !== 0) {
|
|
176
|
+
pDotTerm0 += dg0 * p02;
|
|
177
|
+
}
|
|
178
|
+
if (dg1 !== 0) {
|
|
179
|
+
pDotTerm1 += dg1 * p12;
|
|
180
|
+
}
|
|
181
|
+
dg0 = dgInvAtX0[offset + 3];
|
|
182
|
+
dg1 = dgInvAtX1[offset + 3];
|
|
183
|
+
if (dg0 !== 0) {
|
|
184
|
+
pDotTerm0 += dg0 * p03;
|
|
185
|
+
}
|
|
186
|
+
if (dg1 !== 0) {
|
|
187
|
+
pDotTerm1 += dg1 * p13;
|
|
188
|
+
}
|
|
189
|
+
if (pDotTerm0 !== 0) {
|
|
190
|
+
dp0 += pDotTerm0 * p00;
|
|
191
|
+
}
|
|
192
|
+
if (pDotTerm1 !== 0) {
|
|
193
|
+
dp1 += pDotTerm1 * p10;
|
|
194
|
+
}
|
|
195
|
+
// nu = 1
|
|
196
|
+
offset += 4;
|
|
197
|
+
g0 = gInvAtX0[mu4 + 1];
|
|
198
|
+
if (g0 !== 0)
|
|
199
|
+
dx0 += g0 * p01;
|
|
200
|
+
g1 = gInvAtX1[mu4 + 1];
|
|
201
|
+
if (g1 !== 0)
|
|
202
|
+
dx1 += g1 * p11;
|
|
203
|
+
dg0 = dgInvAtX0[offset];
|
|
204
|
+
dg1 = dgInvAtX1[offset];
|
|
205
|
+
// Guarded like every other rho term, and like master. Assigning
|
|
206
|
+
// unguarded here would compute 0 * p when dg is exactly 0, which is
|
|
207
|
+
// NaN for a non-finite momentum -- manufacturing a divergence signal
|
|
208
|
+
// out of a term that contributes nothing.
|
|
209
|
+
pDotTerm0 = 0;
|
|
210
|
+
if (dg0 !== 0) {
|
|
211
|
+
pDotTerm0 = dg0 * p00;
|
|
212
|
+
}
|
|
213
|
+
pDotTerm1 = 0;
|
|
214
|
+
if (dg1 !== 0) {
|
|
215
|
+
pDotTerm1 = dg1 * p10;
|
|
216
|
+
}
|
|
217
|
+
dg0 = dgInvAtX0[offset + 1];
|
|
218
|
+
dg1 = dgInvAtX1[offset + 1];
|
|
219
|
+
if (dg0 !== 0) {
|
|
220
|
+
pDotTerm0 += dg0 * p01;
|
|
221
|
+
}
|
|
222
|
+
if (dg1 !== 0) {
|
|
223
|
+
pDotTerm1 += dg1 * p11;
|
|
224
|
+
}
|
|
225
|
+
dg0 = dgInvAtX0[offset + 2];
|
|
226
|
+
dg1 = dgInvAtX1[offset + 2];
|
|
227
|
+
if (dg0 !== 0) {
|
|
228
|
+
pDotTerm0 += dg0 * p02;
|
|
229
|
+
}
|
|
230
|
+
if (dg1 !== 0) {
|
|
231
|
+
pDotTerm1 += dg1 * p12;
|
|
232
|
+
}
|
|
233
|
+
dg0 = dgInvAtX0[offset + 3];
|
|
234
|
+
dg1 = dgInvAtX1[offset + 3];
|
|
235
|
+
if (dg0 !== 0) {
|
|
236
|
+
pDotTerm0 += dg0 * p03;
|
|
237
|
+
}
|
|
238
|
+
if (dg1 !== 0) {
|
|
239
|
+
pDotTerm1 += dg1 * p13;
|
|
158
240
|
}
|
|
159
241
|
if (pDotTerm0 !== 0) {
|
|
160
|
-
dp0 += pDotTerm0 *
|
|
242
|
+
dp0 += pDotTerm0 * p01;
|
|
161
243
|
}
|
|
162
244
|
if (pDotTerm1 !== 0) {
|
|
163
|
-
dp1 += pDotTerm1 *
|
|
245
|
+
dp1 += pDotTerm1 * p11;
|
|
164
246
|
}
|
|
247
|
+
// nu = 2
|
|
248
|
+
offset += 4;
|
|
249
|
+
g0 = gInvAtX0[mu4 + 2];
|
|
250
|
+
if (g0 !== 0)
|
|
251
|
+
dx0 += g0 * p02;
|
|
252
|
+
g1 = gInvAtX1[mu4 + 2];
|
|
253
|
+
if (g1 !== 0)
|
|
254
|
+
dx1 += g1 * p12;
|
|
255
|
+
dg0 = dgInvAtX0[offset];
|
|
256
|
+
dg1 = dgInvAtX1[offset];
|
|
257
|
+
// Guarded like every other rho term, and like master. Assigning
|
|
258
|
+
// unguarded here would compute 0 * p when dg is exactly 0, which is
|
|
259
|
+
// NaN for a non-finite momentum -- manufacturing a divergence signal
|
|
260
|
+
// out of a term that contributes nothing.
|
|
261
|
+
pDotTerm0 = 0;
|
|
262
|
+
if (dg0 !== 0) {
|
|
263
|
+
pDotTerm0 = dg0 * p00;
|
|
264
|
+
}
|
|
265
|
+
pDotTerm1 = 0;
|
|
266
|
+
if (dg1 !== 0) {
|
|
267
|
+
pDotTerm1 = dg1 * p10;
|
|
268
|
+
}
|
|
269
|
+
dg0 = dgInvAtX0[offset + 1];
|
|
270
|
+
dg1 = dgInvAtX1[offset + 1];
|
|
271
|
+
if (dg0 !== 0) {
|
|
272
|
+
pDotTerm0 += dg0 * p01;
|
|
273
|
+
}
|
|
274
|
+
if (dg1 !== 0) {
|
|
275
|
+
pDotTerm1 += dg1 * p11;
|
|
276
|
+
}
|
|
277
|
+
dg0 = dgInvAtX0[offset + 2];
|
|
278
|
+
dg1 = dgInvAtX1[offset + 2];
|
|
279
|
+
if (dg0 !== 0) {
|
|
280
|
+
pDotTerm0 += dg0 * p02;
|
|
281
|
+
}
|
|
282
|
+
if (dg1 !== 0) {
|
|
283
|
+
pDotTerm1 += dg1 * p12;
|
|
284
|
+
}
|
|
285
|
+
dg0 = dgInvAtX0[offset + 3];
|
|
286
|
+
dg1 = dgInvAtX1[offset + 3];
|
|
287
|
+
if (dg0 !== 0) {
|
|
288
|
+
pDotTerm0 += dg0 * p03;
|
|
289
|
+
}
|
|
290
|
+
if (dg1 !== 0) {
|
|
291
|
+
pDotTerm1 += dg1 * p13;
|
|
292
|
+
}
|
|
293
|
+
if (pDotTerm0 !== 0) {
|
|
294
|
+
dp0 += pDotTerm0 * p02;
|
|
295
|
+
}
|
|
296
|
+
if (pDotTerm1 !== 0) {
|
|
297
|
+
dp1 += pDotTerm1 * p12;
|
|
298
|
+
}
|
|
299
|
+
// nu = 3
|
|
300
|
+
offset += 4;
|
|
301
|
+
g0 = gInvAtX0[mu4 + 3];
|
|
302
|
+
if (g0 !== 0)
|
|
303
|
+
dx0 += g0 * p03;
|
|
304
|
+
g1 = gInvAtX1[mu4 + 3];
|
|
305
|
+
if (g1 !== 0)
|
|
306
|
+
dx1 += g1 * p13;
|
|
307
|
+
dg0 = dgInvAtX0[offset];
|
|
308
|
+
dg1 = dgInvAtX1[offset];
|
|
309
|
+
// Guarded like every other rho term, and like master. Assigning
|
|
310
|
+
// unguarded here would compute 0 * p when dg is exactly 0, which is
|
|
311
|
+
// NaN for a non-finite momentum -- manufacturing a divergence signal
|
|
312
|
+
// out of a term that contributes nothing.
|
|
313
|
+
pDotTerm0 = 0;
|
|
314
|
+
if (dg0 !== 0) {
|
|
315
|
+
pDotTerm0 = dg0 * p00;
|
|
316
|
+
}
|
|
317
|
+
pDotTerm1 = 0;
|
|
318
|
+
if (dg1 !== 0) {
|
|
319
|
+
pDotTerm1 = dg1 * p10;
|
|
320
|
+
}
|
|
321
|
+
dg0 = dgInvAtX0[offset + 1];
|
|
322
|
+
dg1 = dgInvAtX1[offset + 1];
|
|
323
|
+
if (dg0 !== 0) {
|
|
324
|
+
pDotTerm0 += dg0 * p01;
|
|
325
|
+
}
|
|
326
|
+
if (dg1 !== 0) {
|
|
327
|
+
pDotTerm1 += dg1 * p11;
|
|
328
|
+
}
|
|
329
|
+
dg0 = dgInvAtX0[offset + 2];
|
|
330
|
+
dg1 = dgInvAtX1[offset + 2];
|
|
331
|
+
if (dg0 !== 0) {
|
|
332
|
+
pDotTerm0 += dg0 * p02;
|
|
333
|
+
}
|
|
334
|
+
if (dg1 !== 0) {
|
|
335
|
+
pDotTerm1 += dg1 * p12;
|
|
336
|
+
}
|
|
337
|
+
dg0 = dgInvAtX0[offset + 3];
|
|
338
|
+
dg1 = dgInvAtX1[offset + 3];
|
|
339
|
+
if (dg0 !== 0) {
|
|
340
|
+
pDotTerm0 += dg0 * p03;
|
|
341
|
+
}
|
|
342
|
+
if (dg1 !== 0) {
|
|
343
|
+
pDotTerm1 += dg1 * p13;
|
|
344
|
+
}
|
|
345
|
+
if (pDotTerm0 !== 0) {
|
|
346
|
+
dp0 += pDotTerm0 * p03;
|
|
347
|
+
}
|
|
348
|
+
if (pDotTerm1 !== 0) {
|
|
349
|
+
dp1 += pDotTerm1 * p13;
|
|
350
|
+
}
|
|
351
|
+
dxStageArr0[mu] = dx0;
|
|
352
|
+
dpStageArr0[mu] = dp0;
|
|
353
|
+
dxStageArr1[mu] = dx1;
|
|
354
|
+
dpStageArr1[mu] = dp1;
|
|
355
|
+
}
|
|
356
|
+
}
|
|
357
|
+
else {
|
|
358
|
+
for (let mu = 0; mu < dim; mu++) {
|
|
359
|
+
let dx0 = 0, dp0 = 0;
|
|
360
|
+
let dx1 = 0, dp1 = 0;
|
|
361
|
+
const mu_dim = mu * dim;
|
|
362
|
+
const mu_dim_dim = mu_dim * dim;
|
|
363
|
+
for (let nu = 0; nu < dim; nu++) {
|
|
364
|
+
const idx_g = mu_dim + nu;
|
|
365
|
+
// Stage 0
|
|
366
|
+
const g0 = gInvAtX0[idx_g];
|
|
367
|
+
if (g0 !== 0) {
|
|
368
|
+
dx0 += g0 * p0_st[nu];
|
|
369
|
+
}
|
|
370
|
+
// Stage 1
|
|
371
|
+
const g1 = gInvAtX1[idx_g];
|
|
372
|
+
if (g1 !== 0) {
|
|
373
|
+
dx1 += g1 * p1_st[nu];
|
|
374
|
+
}
|
|
375
|
+
let pDotTerm0 = 0;
|
|
376
|
+
let pDotTerm1 = 0;
|
|
377
|
+
const offset = mu_dim_dim + nu * dim;
|
|
378
|
+
for (let rho = 0; rho < dim; rho++) {
|
|
379
|
+
const dgIdx = offset + rho;
|
|
380
|
+
const dg0 = dgInvAtX0[dgIdx];
|
|
381
|
+
if (dg0 !== 0) {
|
|
382
|
+
pDotTerm0 += dg0 * p0_st[rho];
|
|
383
|
+
}
|
|
384
|
+
const dg1 = dgInvAtX1[dgIdx];
|
|
385
|
+
if (dg1 !== 0) {
|
|
386
|
+
pDotTerm1 += dg1 * p1_st[rho];
|
|
387
|
+
}
|
|
388
|
+
}
|
|
389
|
+
if (pDotTerm0 !== 0) {
|
|
390
|
+
dp0 += pDotTerm0 * p0_st[nu];
|
|
391
|
+
}
|
|
392
|
+
if (pDotTerm1 !== 0) {
|
|
393
|
+
dp1 += pDotTerm1 * p1_st[nu];
|
|
394
|
+
}
|
|
395
|
+
}
|
|
396
|
+
dxStageArr0[mu] = dx0;
|
|
397
|
+
dpStageArr0[mu] = dp0;
|
|
398
|
+
dxStageArr1[mu] = dx1;
|
|
399
|
+
dpStageArr1[mu] = dp1;
|
|
165
400
|
}
|
|
166
|
-
dxStageArr[0][mu] = dx0;
|
|
167
|
-
dpStageArr[0][mu] = dp0;
|
|
168
|
-
dxStageArr[1][mu] = dx1;
|
|
169
|
-
dpStageArr[1][mu] = dp1;
|
|
170
401
|
}
|
|
171
402
|
// 3. Accumulate for i and mu
|
|
172
403
|
for (let mu = 0; mu < dim; mu++) {
|
|
173
404
|
const x0 = state.x[mu];
|
|
174
405
|
const p0 = state.p[mu];
|
|
175
|
-
const dx0 =
|
|
176
|
-
const dx1 =
|
|
177
|
-
const dp0 =
|
|
178
|
-
const dp1 =
|
|
406
|
+
const dx0 = dxStageArr0[mu];
|
|
407
|
+
const dx1 = dxStageArr1[mu];
|
|
408
|
+
const dp0 = dpStageArr0[mu];
|
|
409
|
+
const dp1 = dpStageArr1[mu];
|
|
179
410
|
Xnew[0][mu] = x0 + hA00 * dx0 + hA01 * dx1;
|
|
180
411
|
Pnew[0][mu] = p0 - halfhA00 * dp0 - halfhA01 * dp1;
|
|
181
412
|
Xnew[1][mu] = x0 + hA10 * dx0 + hA11 * dx1;
|
|
@@ -195,22 +426,44 @@ export function solveGL4Stage(state, h, gInverseFn, dgInverseFn, opts) {
|
|
|
195
426
|
if (maxDelta < opts.picardTol) {
|
|
196
427
|
// Clone on return — caller may retain references and the next
|
|
197
428
|
// solveGL4Stage call will overwrite our internal buffers.
|
|
198
|
-
|
|
199
|
-
|
|
200
|
-
|
|
201
|
-
|
|
202
|
-
|
|
203
|
-
|
|
204
|
-
|
|
205
|
-
|
|
206
|
-
|
|
207
|
-
|
|
429
|
+
if (dim === 4) {
|
|
430
|
+
// Bolt: Array literals are significantly faster than allocating arrays and manually populating them in tight loops.
|
|
431
|
+
return {
|
|
432
|
+
stageX: [[X[0][0], X[0][1], X[0][2], X[0][3]], [X[1][0], X[1][1], X[1][2], X[1][3]]],
|
|
433
|
+
stageP: [[P[0][0], P[0][1], P[0][2], P[0][3]], [P[1][0], P[1][1], P[1][2], P[1][3]]],
|
|
434
|
+
stageDx: [[dxStageArr0[0], dxStageArr0[1], dxStageArr0[2], dxStageArr0[3]], [dxStageArr1[0], dxStageArr1[1], dxStageArr1[2], dxStageArr1[3]]],
|
|
435
|
+
stageDp: [[dpStageArr0[0], dpStageArr0[1], dpStageArr0[2], dpStageArr0[3]], [dpStageArr1[0], dpStageArr1[1], dpStageArr1[2], dpStageArr1[3]]],
|
|
436
|
+
iterations: k + 1,
|
|
437
|
+
};
|
|
438
|
+
}
|
|
439
|
+
else {
|
|
440
|
+
// Bolt: Manual loop is significantly faster than Array.from for TypedArrays in tight loops.
|
|
441
|
+
const sX0 = new Array(dim);
|
|
442
|
+
const sX1 = new Array(dim);
|
|
443
|
+
const sP0 = new Array(dim);
|
|
444
|
+
const sP1 = new Array(dim);
|
|
445
|
+
const sDx0 = new Array(dim);
|
|
446
|
+
const sDx1 = new Array(dim);
|
|
447
|
+
const sDp0 = new Array(dim);
|
|
448
|
+
const sDp1 = new Array(dim);
|
|
449
|
+
for (let m = 0; m < dim; m++) {
|
|
450
|
+
sX0[m] = X[0][m];
|
|
451
|
+
sX1[m] = X[1][m];
|
|
452
|
+
sP0[m] = P[0][m];
|
|
453
|
+
sP1[m] = P[1][m];
|
|
454
|
+
sDx0[m] = dxStageArr0[m];
|
|
455
|
+
sDx1[m] = dxStageArr1[m];
|
|
456
|
+
sDp0[m] = dpStageArr0[m];
|
|
457
|
+
sDp1[m] = dpStageArr1[m];
|
|
458
|
+
}
|
|
459
|
+
return {
|
|
460
|
+
stageX: [sX0, sX1],
|
|
461
|
+
stageP: [sP0, sP1],
|
|
462
|
+
stageDx: [sDx0, sDx1],
|
|
463
|
+
stageDp: [sDp0, sDp1],
|
|
464
|
+
iterations: k + 1,
|
|
465
|
+
};
|
|
208
466
|
}
|
|
209
|
-
return {
|
|
210
|
-
stageX: [sX0, sX1],
|
|
211
|
-
stageP: [sP0, sP1],
|
|
212
|
-
iterations: k + 1,
|
|
213
|
-
};
|
|
214
467
|
}
|
|
215
468
|
}
|
|
216
469
|
throw new GL4ConvergenceError(`Picard iteration did not converge in ${opts.picardMaxIter} iterations (maxDelta above picardTol=${opts.picardTol})`);
|
|
@@ -218,105 +471,6 @@ export function solveGL4Stage(state, h, gInverseFn, dgInverseFn, opts) {
|
|
|
218
471
|
// ---------------------------------------------------------------------------
|
|
219
472
|
// Integrator entry-point (Task 3, Phase 1a-iii)
|
|
220
473
|
// ---------------------------------------------------------------------------
|
|
221
|
-
/**
|
|
222
|
-
* Step update for x^μ given converged GL4 stage values.
|
|
223
|
-
*
|
|
224
|
-
* x^μ_{n+1} = x^μ_n + h · Σ_i b_i · (g^{μν}(X_i) P_{i,ν})
|
|
225
|
-
*
|
|
226
|
-
* Module-private — only `integrateGeodesicGL4` consumes it.
|
|
227
|
-
*
|
|
228
|
-
* @internal
|
|
229
|
-
*/
|
|
230
|
-
function updateFromStages(xPrev, h, stageX, stageP, gInverseFn) {
|
|
231
|
-
const dim = xPrev.length;
|
|
232
|
-
const x = xPrev.slice();
|
|
233
|
-
// Precompute metric inverses to avoid redundant evaluations inside the loops
|
|
234
|
-
const gInv0 = gInverseFn(stageX[0]);
|
|
235
|
-
const gInv1 = gInverseFn(stageX[1]);
|
|
236
|
-
const stageP0 = stageP[0];
|
|
237
|
-
const stageP1 = stageP[1];
|
|
238
|
-
const b0 = GL4_B[0];
|
|
239
|
-
const b1 = GL4_B[1];
|
|
240
|
-
for (let mu = 0; mu < dim; mu++) {
|
|
241
|
-
let xDot0 = 0;
|
|
242
|
-
let xDot1 = 0;
|
|
243
|
-
// Bolt: Hoist the mu * dim multiplication out of the inner loop to avoid redundant recalculation
|
|
244
|
-
const mu_dim = mu * dim;
|
|
245
|
-
for (let nu = 0; nu < dim; nu++) {
|
|
246
|
-
const idx = mu_dim + nu;
|
|
247
|
-
const g0 = gInv0[idx];
|
|
248
|
-
if (g0 !== 0) {
|
|
249
|
-
xDot0 += g0 * stageP0[nu];
|
|
250
|
-
}
|
|
251
|
-
const g1 = gInv1[idx];
|
|
252
|
-
if (g1 !== 0) {
|
|
253
|
-
xDot1 += g1 * stageP1[nu];
|
|
254
|
-
}
|
|
255
|
-
}
|
|
256
|
-
// Bolt: Manually unroll the outer loop across stages (i=0,1).
|
|
257
|
-
x[mu] += h * (b0 * xDot0 + b1 * xDot1);
|
|
258
|
-
}
|
|
259
|
-
return x;
|
|
260
|
-
}
|
|
261
|
-
/**
|
|
262
|
-
* Step update for p_μ given converged GL4 stage values.
|
|
263
|
-
*
|
|
264
|
-
* p_μ_{n+1} = p_μ_n − ½ h · Σ_i b_i · (∂_μ g^{νρ})(X_i) P_{i,ν} P_{i,ρ}
|
|
265
|
-
*
|
|
266
|
-
* I2 pin: `dgInverseFn(x)[mu][nu][rho] = ∂_mu g^{nu rho}` — `mu` is the
|
|
267
|
-
* differentiation axis (`λ` in the pinned order). Consistent with the
|
|
268
|
-
* Picard solver's stage update.
|
|
269
|
-
*
|
|
270
|
-
* Module-private.
|
|
271
|
-
*
|
|
272
|
-
* @internal
|
|
273
|
-
*/
|
|
274
|
-
function updateMomentumFromStages(pPrev, h, stageX, stageP, dgInverseFn) {
|
|
275
|
-
const dim = pPrev.length;
|
|
276
|
-
const p = pPrev.slice();
|
|
277
|
-
// Precompute metric derivatives to avoid redundant evaluations inside the loops
|
|
278
|
-
const dg0 = dgInverseFn(stageX[0]);
|
|
279
|
-
const dg1 = dgInverseFn(stageX[1]);
|
|
280
|
-
const stageP0 = stageP[0];
|
|
281
|
-
const stageP1 = stageP[1];
|
|
282
|
-
const b0 = GL4_B[0];
|
|
283
|
-
const b1 = GL4_B[1];
|
|
284
|
-
for (let mu = 0; mu < dim; mu++) {
|
|
285
|
-
// Bolt: Hoist the invariant mu * dim * dim multiplication outside the nu, rho loops
|
|
286
|
-
const mu_dim_dim = mu * dim * dim;
|
|
287
|
-
let pDot0 = 0;
|
|
288
|
-
let pDot1 = 0;
|
|
289
|
-
for (let nu = 0; nu < dim; nu++) {
|
|
290
|
-
// Bolt: Factor out stageP0[nu] and stageP1[nu] from the inner rho loop
|
|
291
|
-
// Original logic was pDot += dg[...] * stageP[i][nu] * stageP[i][rho].
|
|
292
|
-
// This algebraic refactoring reduces total multiplications from ~64 to ~20 per update loop.
|
|
293
|
-
let pDotTerm0 = 0;
|
|
294
|
-
let pDotTerm1 = 0;
|
|
295
|
-
const offset = mu_dim_dim + nu * dim;
|
|
296
|
-
for (let rho = 0; rho < dim; rho++) {
|
|
297
|
-
const idx = offset + rho;
|
|
298
|
-
const dgVal0 = dg0[idx];
|
|
299
|
-
if (dgVal0 !== 0) {
|
|
300
|
-
pDotTerm0 += dgVal0 * stageP0[rho];
|
|
301
|
-
}
|
|
302
|
-
const dgVal1 = dg1[idx];
|
|
303
|
-
if (dgVal1 !== 0) {
|
|
304
|
-
pDotTerm1 += dgVal1 * stageP1[rho];
|
|
305
|
-
}
|
|
306
|
-
}
|
|
307
|
-
if (pDotTerm0 !== 0) {
|
|
308
|
-
pDot0 += pDotTerm0 * stageP0[nu];
|
|
309
|
-
}
|
|
310
|
-
if (pDotTerm1 !== 0) {
|
|
311
|
-
pDot1 += pDotTerm1 * stageP1[nu];
|
|
312
|
-
}
|
|
313
|
-
}
|
|
314
|
-
// Bolt: Manually unroll the outer loop across stages (i=0,1).
|
|
315
|
-
// This avoids inner-loop iterator, function calls and array indexing overhead.
|
|
316
|
-
p[mu] += h * (b0 * (-0.5 * pDot0) + b1 * (-0.5 * pDot1));
|
|
317
|
-
}
|
|
318
|
-
return p;
|
|
319
|
-
}
|
|
320
474
|
/**
|
|
321
475
|
* GL4 symplectic integrator on the canonical (x, p) geodesic Hamiltonian.
|
|
322
476
|
*
|
|
@@ -422,8 +576,43 @@ export function integrateGeodesicGL4(initialState, options) {
|
|
|
422
576
|
if (!stepSucceeded || stages === undefined) {
|
|
423
577
|
throw new GL4ConvergenceError(`GL4 integrator: Picard iteration did not converge even at h_min=${hFloor} (step ${n}). Diagnose step-size or metric singularity.`);
|
|
424
578
|
}
|
|
425
|
-
|
|
426
|
-
|
|
579
|
+
// State is cloned before modification. This avoids modifying the user-provided
|
|
580
|
+
// initial state referenced from the first loop iterations, and ensures that
|
|
581
|
+
// step-halving retries don't cumulatively corrupt `x` and `p`.
|
|
582
|
+
const newX = x.slice();
|
|
583
|
+
const newP = p.slice();
|
|
584
|
+
// Update state using the converged derivatives from the Picard solver.
|
|
585
|
+
// This algebraically factors out redundant closures and tensor contractions.
|
|
586
|
+
// Note: `stages.stageDp` does NOT have the -0.5 factor (see solveGL4Stage logic).
|
|
587
|
+
const b0 = GL4_B[0];
|
|
588
|
+
const b1 = GL4_B[1];
|
|
589
|
+
const dim = newX.length;
|
|
590
|
+
const dx0 = stages.stageDx[0];
|
|
591
|
+
const dx1 = stages.stageDx[1];
|
|
592
|
+
const dp0 = stages.stageDp[0];
|
|
593
|
+
const dp1 = stages.stageDp[1];
|
|
594
|
+
const hb0 = subH * b0;
|
|
595
|
+
const hb1 = subH * b1;
|
|
596
|
+
const mhb0 = hb0 * -0.5;
|
|
597
|
+
const mhb1 = hb1 * -0.5;
|
|
598
|
+
if (dim === 4) {
|
|
599
|
+
newX[0] += hb0 * dx0[0] + hb1 * dx1[0];
|
|
600
|
+
newP[0] += mhb0 * dp0[0] + mhb1 * dp1[0];
|
|
601
|
+
newX[1] += hb0 * dx0[1] + hb1 * dx1[1];
|
|
602
|
+
newP[1] += mhb0 * dp0[1] + mhb1 * dp1[1];
|
|
603
|
+
newX[2] += hb0 * dx0[2] + hb1 * dx1[2];
|
|
604
|
+
newP[2] += mhb0 * dp0[2] + mhb1 * dp1[2];
|
|
605
|
+
newX[3] += hb0 * dx0[3] + hb1 * dx1[3];
|
|
606
|
+
newP[3] += mhb0 * dp0[3] + mhb1 * dp1[3];
|
|
607
|
+
}
|
|
608
|
+
else {
|
|
609
|
+
for (let mu = 0; mu < dim; mu++) {
|
|
610
|
+
newX[mu] += hb0 * dx0[mu] + hb1 * dx1[mu];
|
|
611
|
+
newP[mu] += mhb0 * dp0[mu] + mhb1 * dp1[mu];
|
|
612
|
+
}
|
|
613
|
+
}
|
|
614
|
+
x = newX;
|
|
615
|
+
p = newP;
|
|
427
616
|
remaining -= subH;
|
|
428
617
|
stepH = subH; // keep the converging size for the rest of this macro-step
|
|
429
618
|
lastIterations = stages.iterations;
|
|
@@ -1 +1 @@
|
|
|
1
|
-
{"version":3,"file":"gl4-integrator.js","sourceRoot":"","sources":["../../src/numerical/gl4-integrator.ts"],"names":[],"mappings":"AAAA;;;;;;;;;;;;;;;;;;;;;;;;GAwBG;AACH,OAAO,EAAE,mBAAmB,EAAE,qBAAqB,EAAE,MAAM,aAAa,CAAC;AAEzE,MAAM,YAAY,GAAG,IAAI,CAAC,IAAI,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC;AAEtC;;;;;;;;GAQG;AACH,MAAM,CAAC,MAAM,KAAK,GAA8B,CAAC,GAAG,GAAG,YAAY,EAAE,GAAG,GAAG,YAAY,CAAC,CAAC;AAEzF;;;;;;;;;;;GAWG;AACH,MAAM,CAAC,MAAM,KAAK,GAAoE;IACpF,CAAC,IAAI,EAAE,IAAI,GAAG,YAAY,CAAC;IAC3B,CAAC,IAAI,GAAG,YAAY,EAAE,IAAI,CAAC;CAC5B,CAAC;AAEF;;;;;;GAMG;AACH,MAAM,CAAC,MAAM,KAAK,GAA8B,CAAC,GAAG,EAAE,GAAG,CAAC,CAAC;AAkG3D;;;;;;;;;;;;;;;;;;;;;;;;;GAyBG;AACH,MAAM,UAAU,aAAa,CAC3B,KAAe,EACf,CAAS,EACT,UAAkD,EAClD,WAAmD,EACnD,IAAkD;IAElD,MAAM,GAAG,GAAG,KAAK,CAAC,CAAC,CAAC,MAAM,CAAC;IAC3B,kEAAkE;IAClE,kEAAkE;IAClE,kEAAkE;IAClE,gEAAgE;IAChE,yBAAyB;IACzB,MAAM,KAAK,GAAmB,CAAC,IAAI,YAAY,CAAC,GAAG,CAAC,EAAE,IAAI,YAAY,CAAC,GAAG,CAAC,CAAC,CAAC;IAC7E,MAAM,KAAK,GAAmB,CAAC,IAAI,YAAY,CAAC,GAAG,CAAC,EAAE,IAAI,YAAY,CAAC,GAAG,CAAC,CAAC,CAAC;IAC7E,MAAM,KAAK,GAAmB,CAAC,IAAI,YAAY,CAAC,GAAG,CAAC,EAAE,IAAI,YAAY,CAAC,GAAG,CAAC,CAAC,CAAC;IAC7E,MAAM,KAAK,GAAmB,CAAC,IAAI,YAAY,CAAC,GAAG,CAAC,EAAE,IAAI,YAAY,CAAC,GAAG,CAAC,CAAC,CAAC;IAE7E,6EAA6E;IAC7E,KAAK,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,KAAK,CAAC,CAAC,CAAC,CAAC;IAAC,KAAK,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,KAAK,CAAC,CAAC,CAAC,CAAC;IAC7C,KAAK,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,KAAK,CAAC,CAAC,CAAC,CAAC;IAAC,KAAK,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,KAAK,CAAC,CAAC,CAAC,CAAC;IAC7C,IAAI,CAAC,GAAmB,KAAK,CAAC;IAC9B,IAAI,CAAC,GAAmB,KAAK,CAAC;IAC9B,IAAI,IAAI,GAAmB,KAAK,CAAC;IACjC,IAAI,IAAI,GAAmB,KAAK,CAAC;IAEjC,uEAAuE;IACvE,MAAM,UAAU,GAAmB,CAAC,IAAI,YAAY,CAAC,GAAG,CAAC,EAAE,IAAI,YAAY,CAAC,GAAG,CAAC,CAAC,CAAC;IAClF,MAAM,UAAU,GAAmB,CAAC,IAAI,YAAY,CAAC,GAAG,CAAC,EAAE,IAAI,YAAY,CAAC,GAAG,CAAC,CAAC,CAAC;IAElF,MAAM,IAAI,GAAG,CAAC,GAAG,KAAK,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;IAC7B,MAAM,IAAI,GAAG,CAAC,GAAG,KAAK,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;IAC7B,MAAM,IAAI,GAAG,CAAC,GAAG,KAAK,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;IAC7B,MAAM,IAAI,GAAG,CAAC,GAAG,KAAK,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;IAC7B,MAAM,QAAQ,GAAG,GAAG,GAAG,IAAI,CAAC;IAC5B,MAAM,QAAQ,GAAG,GAAG,GAAG,IAAI,CAAC;IAC5B,MAAM,QAAQ,GAAG,GAAG,GAAG,IAAI,CAAC;IAC5B,MAAM,QAAQ,GAAG,GAAG,GAAG,IAAI,CAAC;IAE5B,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,IAAI,CAAC,aAAa,EAAE,CAAC,EAAE,EAAE,CAAC;QAC5C,8DAA8D;QAC9D,MAAM,QAAQ,GAAG,UAAU,CAAC,CAAC,CAAC,CAAC,CAAiC,CAAC,CAAC;QAClE,MAAM,SAAS,GAAG,WAAW,CAAC,CAAC,CAAC,CAAC,CAAiC,CAAC,CAAC;QACpE,MAAM,QAAQ,GAAG,UAAU,CAAC,CAAC,CAAC,CAAC,CAAiC,CAAC,CAAC;QAClE,MAAM,SAAS,GAAG,WAAW,CAAC,CAAC,CAAC,CAAC,CAAiC,CAAC,CAAC;QAEpE,2EAA2E;QAC3E,MAAM,KAAK,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC;QACnB,MAAM,KAAK,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC;QAEnB,KAAK,IAAI,EAAE,GAAG,CAAC,EAAE,EAAE,GAAG,GAAG,EAAE,EAAE,EAAE,EAAE,CAAC;YAChC,IAAI,GAAG,GAAG,CAAC,EAAE,GAAG,GAAG,CAAC,CAAC;YACrB,IAAI,GAAG,GAAG,CAAC,EAAE,GAAG,GAAG,CAAC,CAAC;YACrB,MAAM,MAAM,GAAG,EAAE,GAAG,GAAG,CAAC;YACxB,MAAM,UAAU,GAAG,MAAM,GAAG,GAAG,CAAC;YAEhC,KAAK,IAAI,EAAE,GAAG,CAAC,EAAE,EAAE,GAAG,GAAG,EAAE,EAAE,EAAE,EAAE,CAAC;gBAChC,MAAM,KAAK,GAAG,MAAM,GAAG,EAAE,CAAC;gBAE1B,UAAU;gBACV,MAAM,EAAE,GAAG,QAAQ,CAAC,KAAK,CAAC,CAAC;gBAC3B,IAAI,EAAE,KAAK,CAAC,EAAE,CAAC;oBACb,GAAG,IAAI,EAAE,GAAG,KAAK,CAAC,EAAE,CAAC,CAAC;gBACxB,CAAC;gBAED,UAAU;gBACV,MAAM,EAAE,GAAG,QAAQ,CAAC,KAAK,CAAC,CAAC;gBAC3B,IAAI,EAAE,KAAK,CAAC,EAAE,CAAC;oBACb,GAAG,IAAI,EAAE,GAAG,KAAK,CAAC,EAAE,CAAC,CAAC;gBACxB,CAAC;gBAED,IAAI,SAAS,GAAG,CAAC,CAAC;gBAClB,IAAI,SAAS,GAAG,CAAC,CAAC;gBAClB,MAAM,MAAM,GAAG,UAAU,GAAG,EAAE,GAAG,GAAG,CAAC;gBAErC,KAAK,IAAI,GAAG,GAAG,CAAC,EAAE,GAAG,GAAG,GAAG,EAAE,GAAG,EAAE,EAAE,CAAC;oBACnC,MAAM,KAAK,GAAG,MAAM,GAAG,GAAG,CAAC;oBAE3B,MAAM,GAAG,GAAG,SAAS,CAAC,KAAK,CAAC,CAAC;oBAC7B,IAAI,GAAG,KAAK,CAAC,EAAE,CAAC;wBACd,SAAS,IAAI,GAAG,GAAG,KAAK,CAAC,GAAG,CAAC,CAAC;oBAChC,CAAC;oBAED,MAAM,GAAG,GAAG,SAAS,CAAC,KAAK,CAAC,CAAC;oBAC7B,IAAI,GAAG,KAAK,CAAC,EAAE,CAAC;wBACd,SAAS,IAAI,GAAG,GAAG,KAAK,CAAC,GAAG,CAAC,CAAC;oBAChC,CAAC;gBACH,CAAC;gBAED,IAAI,SAAS,KAAK,CAAC,EAAE,CAAC;oBACpB,GAAG,IAAI,SAAS,GAAG,KAAK,CAAC,EAAE,CAAC,CAAC;gBAC/B,CAAC;gBACD,IAAI,SAAS,KAAK,CAAC,EAAE,CAAC;oBACpB,GAAG,IAAI,SAAS,GAAG,KAAK,CAAC,EAAE,CAAC,CAAC;gBAC/B,CAAC;YACH,CAAC;YAED,UAAU,CAAC,CAAC,CAAC,CAAC,EAAE,CAAC,GAAG,GAAG,CAAC;YACxB,UAAU,CAAC,CAAC,CAAC,CAAC,EAAE,CAAC,GAAG,GAAG,CAAC;YACxB,UAAU,CAAC,CAAC,CAAC,CAAC,EAAE,CAAC,GAAG,GAAG,CAAC;YACxB,UAAU,CAAC,CAAC,CAAC,CAAC,EAAE,CAAC,GAAG,GAAG,CAAC;QAC1B,CAAC;QAED,6BAA6B;QAC7B,KAAK,IAAI,EAAE,GAAG,CAAC,EAAE,EAAE,GAAG,GAAG,EAAE,EAAE,EAAE,EAAE,CAAC;YAChC,MAAM,EAAE,GAAG,KAAK,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC;YACvB,MAAM,EAAE,GAAG,KAAK,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC;YACvB,MAAM,GAAG,GAAG,UAAU,CAAC,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC;YAC9B,MAAM,GAAG,GAAG,UAAU,CAAC,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC;YAC9B,MAAM,GAAG,GAAG,UAAU,CAAC,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC;YAC9B,MAAM,GAAG,GAAG,UAAU,CAAC,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC;YAE9B,IAAI,CAAC,CAAC,CAAC,CAAC,EAAE,CAAC,GAAG,EAAE,GAAG,IAAI,GAAG,GAAG,GAAG,IAAI,GAAG,GAAG,CAAC;YAC3C,IAAI,CAAC,CAAC,CAAC,CAAC,EAAE,CAAC,GAAG,EAAE,GAAG,QAAQ,GAAG,GAAG,GAAG,QAAQ,GAAG,GAAG,CAAC;YAEnD,IAAI,CAAC,CAAC,CAAC,CAAC,EAAE,CAAC,GAAG,EAAE,GAAG,IAAI,GAAG,GAAG,GAAG,IAAI,GAAG,GAAG,CAAC;YAC3C,IAAI,CAAC,CAAC,CAAC,CAAC,EAAE,CAAC,GAAG,EAAE,GAAG,QAAQ,GAAG,GAAG,GAAG,QAAQ,GAAG,GAAG,CAAC;QACrD,CAAC;QAED,kCAAkC;QAClC,IAAI,QAAQ,GAAG,CAAC,CAAC;QACjB,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC;YAC3B,KAAK,IAAI,EAAE,GAAG,CAAC,EAAE,EAAE,GAAG,GAAG,EAAE,EAAE,EAAE,EAAE,CAAC;gBAChC,QAAQ,GAAG,IAAI,CAAC,GAAG,CAAC,QAAQ,EAAE,IAAI,CAAC,GAAG,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC,EAAE,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC;gBAChE,QAAQ,GAAG,IAAI,CAAC,GAAG,CAAC,QAAQ,EAAE,IAAI,CAAC,GAAG,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC,EAAE,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC;YAClE,CAAC;QACH,CAAC;QAED,6EAA6E;QAC7E,CAAC,CAAC,EAAE,IAAI,CAAC,GAAG,CAAC,IAAI,EAAE,CAAC,CAAC,CAAC;QACtB,CAAC,CAAC,EAAE,IAAI,CAAC,GAAG,CAAC,IAAI,EAAE,CAAC,CAAC,CAAC;QAEtB,IAAI,QAAQ,GAAG,IAAI,CAAC,SAAS,EAAE,CAAC;YAC9B,8DAA8D;YAC9D,0DAA0D;YAC1D,4FAA4F;YAC5F,MAAM,GAAG,GAAG,IAAI,KAAK,CAAS,GAAG,CAAC,CAAC;YACnC,MAAM,GAAG,GAAG,IAAI,KAAK,CAAS,GAAG,CAAC,CAAC;YACnC,MAAM,GAAG,GAAG,IAAI,KAAK,CAAS,GAAG,CAAC,CAAC;YACnC,MAAM,GAAG,GAAG,IAAI,KAAK,CAAS,GAAG,CAAC,CAAC;YACnC,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,GAAG,EAAE,CAAC,EAAE,EAAE,CAAC;gBAC7B,GAAG,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;gBACjB,GAAG,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;gBACjB,GAAG,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;gBACjB,GAAG,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;YACnB,CAAC;YACD,OAAO;gBACL,MAAM,EAAE,CAAC,GAAG,EAAE,GAAG,CAAC;gBAClB,MAAM,EAAE,CAAC,GAAG,EAAE,GAAG,CAAC;gBAClB,UAAU,EAAE,CAAC,GAAG,CAAC;aAClB,CAAC;QACJ,CAAC;IACH,CAAC;IAED,MAAM,IAAI,mBAAmB,CAC3B,wCAAwC,IAAI,CAAC,aAAa,yCAAyC,IAAI,CAAC,SAAS,GAAG,CACrH,CAAC;AACJ,CAAC;AAED,8EAA8E;AAC9E,gDAAgD;AAChD,8EAA8E;AAE9E;;;;;;;;GAQG;AACH,SAAS,gBAAgB,CACvB,KAAwB,EACxB,CAAS,EACT,MAAuD,EACvD,MAAuD,EACvD,UAAkD;IAElD,MAAM,GAAG,GAAG,KAAK,CAAC,MAAM,CAAC;IACzB,MAAM,CAAC,GAAG,KAAK,CAAC,KAAK,EAAc,CAAC;IAEpC,6EAA6E;IAC7E,MAAM,KAAK,GAAG,UAAU,CAAC,MAAM,CAAC,CAAC,CAAC,CAAC,CAAC;IACpC,MAAM,KAAK,GAAG,UAAU,CAAC,MAAM,CAAC,CAAC,CAAC,CAAC,CAAC;IAEpC,MAAM,OAAO,GAAG,MAAM,CAAC,CAAC,CAAC,CAAC;IAC1B,MAAM,OAAO,GAAG,MAAM,CAAC,CAAC,CAAC,CAAC;IAC1B,MAAM,EAAE,GAAG,KAAK,CAAC,CAAC,CAAC,CAAC;IACpB,MAAM,EAAE,GAAG,KAAK,CAAC,CAAC,CAAC,CAAC;IAEpB,KAAK,IAAI,EAAE,GAAG,CAAC,EAAE,EAAE,GAAG,GAAG,EAAE,EAAE,EAAE,EAAE,CAAC;QAChC,IAAI,KAAK,GAAG,CAAC,CAAC;QACd,IAAI,KAAK,GAAG,CAAC,CAAC;QACd,iGAAiG;QACjG,MAAM,MAAM,GAAG,EAAE,GAAG,GAAG,CAAC;QAExB,KAAK,IAAI,EAAE,GAAG,CAAC,EAAE,EAAE,GAAG,GAAG,EAAE,EAAE,EAAE,EAAE,CAAC;YAChC,MAAM,GAAG,GAAG,MAAM,GAAG,EAAE,CAAC;YACxB,MAAM,EAAE,GAAG,KAAK,CAAC,GAAG,CAAC,CAAC;YACtB,IAAI,EAAE,KAAK,CAAC,EAAE,CAAC;gBACb,KAAK,IAAI,EAAE,GAAG,OAAO,CAAC,EAAE,CAAC,CAAC;YAC5B,CAAC;YACD,MAAM,EAAE,GAAG,KAAK,CAAC,GAAG,CAAC,CAAC;YACtB,IAAI,EAAE,KAAK,CAAC,EAAE,CAAC;gBACb,KAAK,IAAI,EAAE,GAAG,OAAO,CAAC,EAAE,CAAC,CAAC;YAC5B,CAAC;QACH,CAAC;QACD,8DAA8D;QAC9D,CAAC,CAAC,EAAE,CAAC,IAAI,CAAC,GAAG,CAAC,EAAE,GAAG,KAAK,GAAG,EAAE,GAAG,KAAK,CAAC,CAAC;IACzC,CAAC;IACD,OAAO,CAAC,CAAC;AACX,CAAC;AAED;;;;;;;;;;;;GAYG;AACH,SAAS,wBAAwB,CAC/B,KAAwB,EACxB,CAAS,EACT,MAAuD,EACvD,MAAuD,EACvD,WAAmD;IAEnD,MAAM,GAAG,GAAG,KAAK,CAAC,MAAM,CAAC;IACzB,MAAM,CAAC,GAAG,KAAK,CAAC,KAAK,EAAc,CAAC;IAEpC,gFAAgF;IAChF,MAAM,GAAG,GAAG,WAAW,CAAC,MAAM,CAAC,CAAC,CAAC,CAAC,CAAC;IACnC,MAAM,GAAG,GAAG,WAAW,CAAC,MAAM,CAAC,CAAC,CAAC,CAAC,CAAC;IAEnC,MAAM,OAAO,GAAG,MAAM,CAAC,CAAC,CAAC,CAAC;IAC1B,MAAM,OAAO,GAAG,MAAM,CAAC,CAAC,CAAC,CAAC;IAC1B,MAAM,EAAE,GAAG,KAAK,CAAC,CAAC,CAAC,CAAC;IACpB,MAAM,EAAE,GAAG,KAAK,CAAC,CAAC,CAAC,CAAC;IAEpB,KAAK,IAAI,EAAE,GAAG,CAAC,EAAE,EAAE,GAAG,GAAG,EAAE,EAAE,EAAE,EAAE,CAAC;QAChC,oFAAoF;QACpF,MAAM,UAAU,GAAG,EAAE,GAAG,GAAG,GAAG,GAAG,CAAC;QAClC,IAAI,KAAK,GAAG,CAAC,CAAC;QACd,IAAI,KAAK,GAAG,CAAC,CAAC;QAEd,KAAK,IAAI,EAAE,GAAG,CAAC,EAAE,EAAE,GAAG,GAAG,EAAE,EAAE,EAAE,EAAE,CAAC;YAChC,uEAAuE;YACvE,uEAAuE;YACvE,4FAA4F;YAC5F,IAAI,SAAS,GAAG,CAAC,CAAC;YAClB,IAAI,SAAS,GAAG,CAAC,CAAC;YAClB,MAAM,MAAM,GAAG,UAAU,GAAG,EAAE,GAAG,GAAG,CAAC;YAErC,KAAK,IAAI,GAAG,GAAG,CAAC,EAAE,GAAG,GAAG,GAAG,EAAE,GAAG,EAAE,EAAE,CAAC;gBACnC,MAAM,GAAG,GAAG,MAAM,GAAG,GAAG,CAAC;gBACzB,MAAM,MAAM,GAAG,GAAG,CAAC,GAAG,CAAC,CAAC;gBACxB,IAAI,MAAM,KAAK,CAAC,EAAE,CAAC;oBACjB,SAAS,IAAI,MAAM,GAAG,OAAO,CAAC,GAAG,CAAC,CAAC;gBACrC,CAAC;gBACD,MAAM,MAAM,GAAG,GAAG,CAAC,GAAG,CAAC,CAAC;gBACxB,IAAI,MAAM,KAAK,CAAC,EAAE,CAAC;oBACjB,SAAS,IAAI,MAAM,GAAG,OAAO,CAAC,GAAG,CAAC,CAAC;gBACrC,CAAC;YACH,CAAC;YACD,IAAI,SAAS,KAAK,CAAC,EAAE,CAAC;gBACpB,KAAK,IAAI,SAAS,GAAG,OAAO,CAAC,EAAE,CAAC,CAAC;YACnC,CAAC;YACD,IAAI,SAAS,KAAK,CAAC,EAAE,CAAC;gBACpB,KAAK,IAAI,SAAS,GAAG,OAAO,CAAC,EAAE,CAAC,CAAC;YACnC,CAAC;QACH,CAAC;QAED,8DAA8D;QAC9D,+EAA+E;QAC/E,CAAC,CAAC,EAAE,CAAC,IAAI,CAAC,GAAG,CAAC,EAAE,GAAG,CAAC,CAAC,GAAG,GAAG,KAAK,CAAC,GAAG,EAAE,GAAG,CAAC,CAAC,GAAG,GAAG,KAAK,CAAC,CAAC,CAAC;IAC3D,CAAC;IACD,OAAO,CAAC,CAAC;AACX,CAAC;AAED;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;GAwDG;AACH,MAAM,UAAU,oBAAoB,CAClC,YAAsB,EACtB,OAAmB;IAEnB,MAAM,EACJ,KAAK,EACL,MAAM,EACN,UAAU,EACV,WAAW,EACX,SAAS,GAAG,KAAK,EACjB,aAAa,GAAG,EAAE,EAClB,IAAI,EACJ,eAAe,EACf,MAAM,GACP,GAAG,OAAO,CAAC;IAEZ,IAAI,eAAe,KAAK,SAAS,IAAI,YAAY,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,eAAe,EAAE,CAAC;QACzE,MAAM,IAAI,qBAAqB,CAC7B,6BAA6B,YAAY,CAAC,CAAC,CAAC,CAAC,CAAC,sBAAsB,eAAe,qBAAqB,CACzG,CAAC;IACJ,CAAC;IAED,MAAM,CAAC,GAAG,MAAM,GAAG,KAAK,CAAC;IACzB,MAAM,MAAM,GAAG,IAAI,IAAI,CAAC,GAAG,IAAI,CAAC;IAChC,MAAM,SAAS,GAAkB;QAC/B,EAAE,GAAG,EAAE,CAAC,EAAE,CAAC,EAAE,YAAY,CAAC,CAAC,CAAC,KAAK,EAAc,EAAE,CAAC,EAAE,YAAY,CAAC,CAAC,CAAC,KAAK,EAAc,EAAE;KACzF,CAAC;IACF,IAAI,CAAC,GAAG,YAAY,CAAC,CAAC,CAAC,KAAK,EAAc,CAAC;IAC3C,IAAI,CAAC,GAAG,YAAY,CAAC,CAAC,CAAC,KAAK,EAAc,CAAC;IAE3C,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,KAAK,EAAE,CAAC,EAAE,EAAE,CAAC;QAC/B,8DAA8D;QAC9D,uEAAuE;QACvE,yEAAyE;QACzE,uEAAuE;QACvE,wEAAwE;QACxE,sEAAsE;QACtE,4DAA4D;QAC5D,IAAI,SAAS,GAAG,CAAC,CAAC;QAClB,IAAI,KAAK,GAAG,CAAC,CAAC;QACd,IAAI,QAAQ,GAAG,CAAC,CAAC;QACjB,IAAI,cAAc,GAAG,CAAC,CAAC;QACvB,+EAA+E;QAC/E,MAAM,SAAS,GAAG,CAAC,GAAG,KAAK,CAAC;QAC5B,OAAO,SAAS,GAAG,SAAS,EAAE,CAAC;YAC7B,MAAM,MAAM,GAAG,IAAI,CAAC,GAAG,CAAC,KAAK,EAAE,SAAS,CAAC,CAAC;YAC1C,IAAI,MAAoC,CAAC;YACzC,IAAI,aAAa,GAAG,KAAK,CAAC;YAC1B,IAAI,IAAI,GAAG,MAAM,CAAC;YAClB,OAAO,IAAI,IAAI,MAAM,EAAE,CAAC;gBACtB,IAAI,CAAC;oBACH,MAAM,GAAG,aAAa,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,EAAE,IAAI,EAAE,UAAU,EAAE,WAAW,EAAE;wBAC9D,SAAS;wBACT,aAAa;qBACd,CAAC,CAAC;oBACH,aAAa,GAAG,IAAI,CAAC;oBACrB,MAAM;gBACR,CAAC;gBAAC,MAAM,CAAC;oBACP,IAAI,IAAI,CAAC,CAAC;oBACV,QAAQ,EAAE,CAAC;gBACb,CAAC;YACH,CAAC;YACD,IAAI,CAAC,aAAa,IAAI,MAAM,KAAK,SAAS,EAAE,CAAC;gBAC3C,MAAM,IAAI,mBAAmB,CAC3B,mEAAmE,MAAM,UAAU,CAAC,8CAA8C,CACnI,CAAC;YACJ,CAAC;YACD,CAAC,GAAG,gBAAgB,CAAC,CAAC,EAAE,IAAI,EAAE,MAAM,CAAC,MAAM,EAAE,MAAM,CAAC,MAAM,EAAE,UAAU,CAAC,CAAC;YACxE,CAAC,GAAG,wBAAwB,CAAC,CAAC,EAAE,IAAI,EAAE,MAAM,CAAC,MAAM,EAAE,MAAM,CAAC,MAAM,EAAE,WAAW,CAAC,CAAC;YACjF,SAAS,IAAI,IAAI,CAAC;YAClB,KAAK,GAAG,IAAI,CAAC,CAAC,2DAA2D;YACzE,cAAc,GAAG,MAAM,CAAC,UAAU,CAAC;QACrC,CAAC;QACD,SAAS,CAAC,IAAI,CAAC,EAAE,GAAG,EAAE,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC,KAAK,EAAE,EAAE,CAAC,EAAE,CAAC,CAAC,KAAK,EAAE,EAAE,CAAC,CAAC;QACjE,IAAI,MAAM,KAAK,SAAS,EAAE,CAAC;YACzB,MAAM,CAAC,EAAE,IAAI,EAAE,CAAC,EAAE,UAAU,EAAE,cAAc,EAAE,QAAQ,EAAE,CAAC,CAAC;QAC5D,CAAC;IACH,CAAC;IAED,OAAO,SAAS,CAAC;AACnB,CAAC"}
|
|
1
|
+
{"version":3,"file":"gl4-integrator.js","sourceRoot":"","sources":["../../src/numerical/gl4-integrator.ts"],"names":[],"mappings":"AAAA;;;;;;;;;;;;;;;;;;;;;;;;GAwBG;AACH,OAAO,EAAE,mBAAmB,EAAE,qBAAqB,EAAE,MAAM,aAAa,CAAC;AAEzE,MAAM,YAAY,GAAG,IAAI,CAAC,IAAI,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC;AAEtC;;;;;;;;GAQG;AACH,MAAM,CAAC,MAAM,KAAK,GAA8B,CAAC,GAAG,GAAG,YAAY,EAAE,GAAG,GAAG,YAAY,CAAC,CAAC;AAEzF;;;;;;;;;;;GAWG;AACH,MAAM,CAAC,MAAM,KAAK,GAAoE;IACpF,CAAC,IAAI,EAAE,IAAI,GAAG,YAAY,CAAC;IAC3B,CAAC,IAAI,GAAG,YAAY,EAAE,IAAI,CAAC;CAC5B,CAAC;AAEF;;;;;;GAMG;AACH,MAAM,CAAC,MAAM,KAAK,GAA8B,CAAC,GAAG,EAAE,GAAG,CAAC,CAAC;AAoG3D;;;;;;;;;;;;;;;;;;;;;;;;;GAyBG;AACH,MAAM,UAAU,aAAa,CAC3B,KAAe,EACf,CAAS,EACT,UAAkD,EAClD,WAAmD,EACnD,IAAkD;IAElD,MAAM,GAAG,GAAG,KAAK,CAAC,CAAC,CAAC,MAAM,CAAC;IAC3B,kEAAkE;IAClE,kEAAkE;IAClE,kEAAkE;IAClE,gEAAgE;IAChE,yBAAyB;IACzB,MAAM,KAAK,GAAmB,CAAC,IAAI,YAAY,CAAC,GAAG,CAAC,EAAE,IAAI,YAAY,CAAC,GAAG,CAAC,CAAC,CAAC;IAC7E,MAAM,KAAK,GAAmB,CAAC,IAAI,YAAY,CAAC,GAAG,CAAC,EAAE,IAAI,YAAY,CAAC,GAAG,CAAC,CAAC,CAAC;IAC7E,MAAM,KAAK,GAAmB,CAAC,IAAI,YAAY,CAAC,GAAG,CAAC,EAAE,IAAI,YAAY,CAAC,GAAG,CAAC,CAAC,CAAC;IAC7E,MAAM,KAAK,GAAmB,CAAC,IAAI,YAAY,CAAC,GAAG,CAAC,EAAE,IAAI,YAAY,CAAC,GAAG,CAAC,CAAC,CAAC;IAE7E,6EAA6E;IAC7E,KAAK,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,KAAK,CAAC,CAAC,CAAC,CAAC;IAAC,KAAK,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,KAAK,CAAC,CAAC,CAAC,CAAC;IAC7C,KAAK,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,KAAK,CAAC,CAAC,CAAC,CAAC;IAAC,KAAK,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,KAAK,CAAC,CAAC,CAAC,CAAC;IAC7C,IAAI,CAAC,GAAmB,KAAK,CAAC;IAC9B,IAAI,CAAC,GAAmB,KAAK,CAAC;IAC9B,IAAI,IAAI,GAAmB,KAAK,CAAC;IACjC,IAAI,IAAI,GAAmB,KAAK,CAAC;IAEjC,uEAAuE;IACvE,MAAM,UAAU,GAAmB,CAAC,IAAI,YAAY,CAAC,GAAG,CAAC,EAAE,IAAI,YAAY,CAAC,GAAG,CAAC,CAAC,CAAC;IAClF,MAAM,UAAU,GAAmB,CAAC,IAAI,YAAY,CAAC,GAAG,CAAC,EAAE,IAAI,YAAY,CAAC,GAAG,CAAC,CAAC,CAAC;IAClF,MAAM,WAAW,GAAG,UAAU,CAAC,CAAC,CAAC,CAAC;IAClC,MAAM,WAAW,GAAG,UAAU,CAAC,CAAC,CAAC,CAAC;IAClC,MAAM,WAAW,GAAG,UAAU,CAAC,CAAC,CAAC,CAAC;IAClC,MAAM,WAAW,GAAG,UAAU,CAAC,CAAC,CAAC,CAAC;IAElC,MAAM,IAAI,GAAG,CAAC,GAAG,KAAK,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;IAC7B,MAAM,IAAI,GAAG,CAAC,GAAG,KAAK,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;IAC7B,MAAM,IAAI,GAAG,CAAC,GAAG,KAAK,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;IAC7B,MAAM,IAAI,GAAG,CAAC,GAAG,KAAK,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;IAC7B,MAAM,QAAQ,GAAG,GAAG,GAAG,IAAI,CAAC;IAC5B,MAAM,QAAQ,GAAG,GAAG,GAAG,IAAI,CAAC;IAC5B,MAAM,QAAQ,GAAG,GAAG,GAAG,IAAI,CAAC;IAC5B,MAAM,QAAQ,GAAG,GAAG,GAAG,IAAI,CAAC;IAE5B,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,IAAI,CAAC,aAAa,EAAE,CAAC,EAAE,EAAE,CAAC;QAC5C,8DAA8D;QAC9D,MAAM,QAAQ,GAAG,UAAU,CAAC,CAAC,CAAC,CAAC,CAAiC,CAAC,CAAC;QAClE,MAAM,SAAS,GAAG,WAAW,CAAC,CAAC,CAAC,CAAC,CAAiC,CAAC,CAAC;QACpE,MAAM,QAAQ,GAAG,UAAU,CAAC,CAAC,CAAC,CAAC,CAAiC,CAAC,CAAC;QAClE,MAAM,SAAS,GAAG,WAAW,CAAC,CAAC,CAAC,CAAC,CAAiC,CAAC,CAAC;QAEpE,2EAA2E;QAC3E,MAAM,KAAK,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC;QACnB,MAAM,KAAK,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC;QAEnB,IAAI,GAAG,KAAK,CAAC,EAAE,CAAC;YACd,mFAAmF;YACnF,+EAA+E;YAC/E,sFAAsF;YACtF,MAAM,GAAG,GAAG,KAAK,CAAC,CAAC,CAAC,EAAE,GAAG,GAAG,KAAK,CAAC,CAAC,CAAC,EAAE,GAAG,GAAG,KAAK,CAAC,CAAC,CAAC,EAAE,GAAG,GAAG,KAAK,CAAC,CAAC,CAAC,CAAC;YACrE,MAAM,GAAG,GAAG,KAAK,CAAC,CAAC,CAAC,EAAE,GAAG,GAAG,KAAK,CAAC,CAAC,CAAC,EAAE,GAAG,GAAG,KAAK,CAAC,CAAC,CAAC,EAAE,GAAG,GAAG,KAAK,CAAC,CAAC,CAAC,CAAC;YAErE,KAAK,IAAI,EAAE,GAAG,CAAC,EAAE,EAAE,GAAG,CAAC,EAAE,EAAE,EAAE,EAAE,CAAC;gBAC9B,IAAI,GAAG,GAAG,CAAC,EAAE,GAAG,GAAG,CAAC,EAAE,GAAG,GAAG,CAAC,EAAE,GAAG,GAAG,CAAC,CAAC;gBACvC,MAAM,GAAG,GAAG,EAAE,GAAG,CAAC,CAAC;gBAEnB,IAAI,EAAE,EAAE,EAAE,EAAE,GAAG,EAAE,GAAG,CAAC;gBACrB,IAAI,SAAS,EAAE,SAAS,CAAC;gBAEzB,IAAI,MAAM,GAAG,EAAE,GAAG,EAAE,CAAC;gBAErB,SAAS;gBACT,EAAE,GAAG,QAAQ,CAAC,GAAG,CAAC,CAAC;gBAAC,IAAI,EAAE,KAAK,CAAC;oBAAE,GAAG,IAAI,EAAE,GAAG,GAAG,CAAC;gBAClD,EAAE,GAAG,QAAQ,CAAC,GAAG,CAAC,CAAC;gBAAC,IAAI,EAAE,KAAK,CAAC;oBAAE,GAAG,IAAI,EAAE,GAAG,GAAG,CAAC;gBAClD,GAAG,GAAG,SAAS,CAAC,MAAM,CAAC,CAAC;gBAAC,GAAG,GAAG,SAAS,CAAC,MAAM,CAAC,CAAC;gBACjD,gEAAgE;gBAChE,oEAAoE;gBACpE,qEAAqE;gBACrE,0CAA0C;gBAC1C,SAAS,GAAG,CAAC,CAAC;gBAAC,IAAI,GAAG,KAAK,CAAC,EAAE,CAAC;oBAAC,SAAS,GAAG,GAAG,GAAG,GAAG,CAAC;gBAAC,CAAC;gBACxD,SAAS,GAAG,CAAC,CAAC;gBAAC,IAAI,GAAG,KAAK,CAAC,EAAE,CAAC;oBAAC,SAAS,GAAG,GAAG,GAAG,GAAG,CAAC;gBAAC,CAAC;gBACxD,GAAG,GAAG,SAAS,CAAC,MAAM,GAAC,CAAC,CAAC,CAAC;gBAAC,GAAG,GAAG,SAAS,CAAC,MAAM,GAAC,CAAC,CAAC,CAAC;gBACrD,IAAI,GAAG,KAAK,CAAC,EAAE,CAAC;oBAAC,SAAS,IAAI,GAAG,GAAG,GAAG,CAAC;gBAAC,CAAC;gBAC1C,IAAI,GAAG,KAAK,CAAC,EAAE,CAAC;oBAAC,SAAS,IAAI,GAAG,GAAG,GAAG,CAAC;gBAAC,CAAC;gBAC1C,GAAG,GAAG,SAAS,CAAC,MAAM,GAAC,CAAC,CAAC,CAAC;gBAAC,GAAG,GAAG,SAAS,CAAC,MAAM,GAAC,CAAC,CAAC,CAAC;gBACrD,IAAI,GAAG,KAAK,CAAC,EAAE,CAAC;oBAAC,SAAS,IAAI,GAAG,GAAG,GAAG,CAAC;gBAAC,CAAC;gBAC1C,IAAI,GAAG,KAAK,CAAC,EAAE,CAAC;oBAAC,SAAS,IAAI,GAAG,GAAG,GAAG,CAAC;gBAAC,CAAC;gBAC1C,GAAG,GAAG,SAAS,CAAC,MAAM,GAAC,CAAC,CAAC,CAAC;gBAAC,GAAG,GAAG,SAAS,CAAC,MAAM,GAAC,CAAC,CAAC,CAAC;gBACrD,IAAI,GAAG,KAAK,CAAC,EAAE,CAAC;oBAAC,SAAS,IAAI,GAAG,GAAG,GAAG,CAAC;gBAAC,CAAC;gBAC1C,IAAI,GAAG,KAAK,CAAC,EAAE,CAAC;oBAAC,SAAS,IAAI,GAAG,GAAG,GAAG,CAAC;gBAAC,CAAC;gBAC1C,IAAI,SAAS,KAAK,CAAC,EAAE,CAAC;oBAAC,GAAG,IAAI,SAAS,GAAG,GAAG,CAAC;gBAAC,CAAC;gBAChD,IAAI,SAAS,KAAK,CAAC,EAAE,CAAC;oBAAC,GAAG,IAAI,SAAS,GAAG,GAAG,CAAC;gBAAC,CAAC;gBAEhD,SAAS;gBACT,MAAM,IAAI,CAAC,CAAC;gBACZ,EAAE,GAAG,QAAQ,CAAC,GAAG,GAAC,CAAC,CAAC,CAAC;gBAAC,IAAI,EAAE,KAAK,CAAC;oBAAE,GAAG,IAAI,EAAE,GAAG,GAAG,CAAC;gBACpD,EAAE,GAAG,QAAQ,CAAC,GAAG,GAAC,CAAC,CAAC,CAAC;gBAAC,IAAI,EAAE,KAAK,CAAC;oBAAE,GAAG,IAAI,EAAE,GAAG,GAAG,CAAC;gBACpD,GAAG,GAAG,SAAS,CAAC,MAAM,CAAC,CAAC;gBAAC,GAAG,GAAG,SAAS,CAAC,MAAM,CAAC,CAAC;gBACjD,gEAAgE;gBAChE,oEAAoE;gBACpE,qEAAqE;gBACrE,0CAA0C;gBAC1C,SAAS,GAAG,CAAC,CAAC;gBAAC,IAAI,GAAG,KAAK,CAAC,EAAE,CAAC;oBAAC,SAAS,GAAG,GAAG,GAAG,GAAG,CAAC;gBAAC,CAAC;gBACxD,SAAS,GAAG,CAAC,CAAC;gBAAC,IAAI,GAAG,KAAK,CAAC,EAAE,CAAC;oBAAC,SAAS,GAAG,GAAG,GAAG,GAAG,CAAC;gBAAC,CAAC;gBACxD,GAAG,GAAG,SAAS,CAAC,MAAM,GAAC,CAAC,CAAC,CAAC;gBAAC,GAAG,GAAG,SAAS,CAAC,MAAM,GAAC,CAAC,CAAC,CAAC;gBACrD,IAAI,GAAG,KAAK,CAAC,EAAE,CAAC;oBAAC,SAAS,IAAI,GAAG,GAAG,GAAG,CAAC;gBAAC,CAAC;gBAC1C,IAAI,GAAG,KAAK,CAAC,EAAE,CAAC;oBAAC,SAAS,IAAI,GAAG,GAAG,GAAG,CAAC;gBAAC,CAAC;gBAC1C,GAAG,GAAG,SAAS,CAAC,MAAM,GAAC,CAAC,CAAC,CAAC;gBAAC,GAAG,GAAG,SAAS,CAAC,MAAM,GAAC,CAAC,CAAC,CAAC;gBACrD,IAAI,GAAG,KAAK,CAAC,EAAE,CAAC;oBAAC,SAAS,IAAI,GAAG,GAAG,GAAG,CAAC;gBAAC,CAAC;gBAC1C,IAAI,GAAG,KAAK,CAAC,EAAE,CAAC;oBAAC,SAAS,IAAI,GAAG,GAAG,GAAG,CAAC;gBAAC,CAAC;gBAC1C,GAAG,GAAG,SAAS,CAAC,MAAM,GAAC,CAAC,CAAC,CAAC;gBAAC,GAAG,GAAG,SAAS,CAAC,MAAM,GAAC,CAAC,CAAC,CAAC;gBACrD,IAAI,GAAG,KAAK,CAAC,EAAE,CAAC;oBAAC,SAAS,IAAI,GAAG,GAAG,GAAG,CAAC;gBAAC,CAAC;gBAC1C,IAAI,GAAG,KAAK,CAAC,EAAE,CAAC;oBAAC,SAAS,IAAI,GAAG,GAAG,GAAG,CAAC;gBAAC,CAAC;gBAC1C,IAAI,SAAS,KAAK,CAAC,EAAE,CAAC;oBAAC,GAAG,IAAI,SAAS,GAAG,GAAG,CAAC;gBAAC,CAAC;gBAChD,IAAI,SAAS,KAAK,CAAC,EAAE,CAAC;oBAAC,GAAG,IAAI,SAAS,GAAG,GAAG,CAAC;gBAAC,CAAC;gBAEhD,SAAS;gBACT,MAAM,IAAI,CAAC,CAAC;gBACZ,EAAE,GAAG,QAAQ,CAAC,GAAG,GAAC,CAAC,CAAC,CAAC;gBAAC,IAAI,EAAE,KAAK,CAAC;oBAAE,GAAG,IAAI,EAAE,GAAG,GAAG,CAAC;gBACpD,EAAE,GAAG,QAAQ,CAAC,GAAG,GAAC,CAAC,CAAC,CAAC;gBAAC,IAAI,EAAE,KAAK,CAAC;oBAAE,GAAG,IAAI,EAAE,GAAG,GAAG,CAAC;gBACpD,GAAG,GAAG,SAAS,CAAC,MAAM,CAAC,CAAC;gBAAC,GAAG,GAAG,SAAS,CAAC,MAAM,CAAC,CAAC;gBACjD,gEAAgE;gBAChE,oEAAoE;gBACpE,qEAAqE;gBACrE,0CAA0C;gBAC1C,SAAS,GAAG,CAAC,CAAC;gBAAC,IAAI,GAAG,KAAK,CAAC,EAAE,CAAC;oBAAC,SAAS,GAAG,GAAG,GAAG,GAAG,CAAC;gBAAC,CAAC;gBACxD,SAAS,GAAG,CAAC,CAAC;gBAAC,IAAI,GAAG,KAAK,CAAC,EAAE,CAAC;oBAAC,SAAS,GAAG,GAAG,GAAG,GAAG,CAAC;gBAAC,CAAC;gBACxD,GAAG,GAAG,SAAS,CAAC,MAAM,GAAC,CAAC,CAAC,CAAC;gBAAC,GAAG,GAAG,SAAS,CAAC,MAAM,GAAC,CAAC,CAAC,CAAC;gBACrD,IAAI,GAAG,KAAK,CAAC,EAAE,CAAC;oBAAC,SAAS,IAAI,GAAG,GAAG,GAAG,CAAC;gBAAC,CAAC;gBAC1C,IAAI,GAAG,KAAK,CAAC,EAAE,CAAC;oBAAC,SAAS,IAAI,GAAG,GAAG,GAAG,CAAC;gBAAC,CAAC;gBAC1C,GAAG,GAAG,SAAS,CAAC,MAAM,GAAC,CAAC,CAAC,CAAC;gBAAC,GAAG,GAAG,SAAS,CAAC,MAAM,GAAC,CAAC,CAAC,CAAC;gBACrD,IAAI,GAAG,KAAK,CAAC,EAAE,CAAC;oBAAC,SAAS,IAAI,GAAG,GAAG,GAAG,CAAC;gBAAC,CAAC;gBAC1C,IAAI,GAAG,KAAK,CAAC,EAAE,CAAC;oBAAC,SAAS,IAAI,GAAG,GAAG,GAAG,CAAC;gBAAC,CAAC;gBAC1C,GAAG,GAAG,SAAS,CAAC,MAAM,GAAC,CAAC,CAAC,CAAC;gBAAC,GAAG,GAAG,SAAS,CAAC,MAAM,GAAC,CAAC,CAAC,CAAC;gBACrD,IAAI,GAAG,KAAK,CAAC,EAAE,CAAC;oBAAC,SAAS,IAAI,GAAG,GAAG,GAAG,CAAC;gBAAC,CAAC;gBAC1C,IAAI,GAAG,KAAK,CAAC,EAAE,CAAC;oBAAC,SAAS,IAAI,GAAG,GAAG,GAAG,CAAC;gBAAC,CAAC;gBAC1C,IAAI,SAAS,KAAK,CAAC,EAAE,CAAC;oBAAC,GAAG,IAAI,SAAS,GAAG,GAAG,CAAC;gBAAC,CAAC;gBAChD,IAAI,SAAS,KAAK,CAAC,EAAE,CAAC;oBAAC,GAAG,IAAI,SAAS,GAAG,GAAG,CAAC;gBAAC,CAAC;gBAEhD,SAAS;gBACT,MAAM,IAAI,CAAC,CAAC;gBACZ,EAAE,GAAG,QAAQ,CAAC,GAAG,GAAC,CAAC,CAAC,CAAC;gBAAC,IAAI,EAAE,KAAK,CAAC;oBAAE,GAAG,IAAI,EAAE,GAAG,GAAG,CAAC;gBACpD,EAAE,GAAG,QAAQ,CAAC,GAAG,GAAC,CAAC,CAAC,CAAC;gBAAC,IAAI,EAAE,KAAK,CAAC;oBAAE,GAAG,IAAI,EAAE,GAAG,GAAG,CAAC;gBACpD,GAAG,GAAG,SAAS,CAAC,MAAM,CAAC,CAAC;gBAAC,GAAG,GAAG,SAAS,CAAC,MAAM,CAAC,CAAC;gBACjD,gEAAgE;gBAChE,oEAAoE;gBACpE,qEAAqE;gBACrE,0CAA0C;gBAC1C,SAAS,GAAG,CAAC,CAAC;gBAAC,IAAI,GAAG,KAAK,CAAC,EAAE,CAAC;oBAAC,SAAS,GAAG,GAAG,GAAG,GAAG,CAAC;gBAAC,CAAC;gBACxD,SAAS,GAAG,CAAC,CAAC;gBAAC,IAAI,GAAG,KAAK,CAAC,EAAE,CAAC;oBAAC,SAAS,GAAG,GAAG,GAAG,GAAG,CAAC;gBAAC,CAAC;gBACxD,GAAG,GAAG,SAAS,CAAC,MAAM,GAAC,CAAC,CAAC,CAAC;gBAAC,GAAG,GAAG,SAAS,CAAC,MAAM,GAAC,CAAC,CAAC,CAAC;gBACrD,IAAI,GAAG,KAAK,CAAC,EAAE,CAAC;oBAAC,SAAS,IAAI,GAAG,GAAG,GAAG,CAAC;gBAAC,CAAC;gBAC1C,IAAI,GAAG,KAAK,CAAC,EAAE,CAAC;oBAAC,SAAS,IAAI,GAAG,GAAG,GAAG,CAAC;gBAAC,CAAC;gBAC1C,GAAG,GAAG,SAAS,CAAC,MAAM,GAAC,CAAC,CAAC,CAAC;gBAAC,GAAG,GAAG,SAAS,CAAC,MAAM,GAAC,CAAC,CAAC,CAAC;gBACrD,IAAI,GAAG,KAAK,CAAC,EAAE,CAAC;oBAAC,SAAS,IAAI,GAAG,GAAG,GAAG,CAAC;gBAAC,CAAC;gBAC1C,IAAI,GAAG,KAAK,CAAC,EAAE,CAAC;oBAAC,SAAS,IAAI,GAAG,GAAG,GAAG,CAAC;gBAAC,CAAC;gBAC1C,GAAG,GAAG,SAAS,CAAC,MAAM,GAAC,CAAC,CAAC,CAAC;gBAAC,GAAG,GAAG,SAAS,CAAC,MAAM,GAAC,CAAC,CAAC,CAAC;gBACrD,IAAI,GAAG,KAAK,CAAC,EAAE,CAAC;oBAAC,SAAS,IAAI,GAAG,GAAG,GAAG,CAAC;gBAAC,CAAC;gBAC1C,IAAI,GAAG,KAAK,CAAC,EAAE,CAAC;oBAAC,SAAS,IAAI,GAAG,GAAG,GAAG,CAAC;gBAAC,CAAC;gBAC1C,IAAI,SAAS,KAAK,CAAC,EAAE,CAAC;oBAAC,GAAG,IAAI,SAAS,GAAG,GAAG,CAAC;gBAAC,CAAC;gBAChD,IAAI,SAAS,KAAK,CAAC,EAAE,CAAC;oBAAC,GAAG,IAAI,SAAS,GAAG,GAAG,CAAC;gBAAC,CAAC;gBAEhD,WAAW,CAAC,EAAE,CAAC,GAAG,GAAG,CAAC;gBACtB,WAAW,CAAC,EAAE,CAAC,GAAG,GAAG,CAAC;gBACtB,WAAW,CAAC,EAAE,CAAC,GAAG,GAAG,CAAC;gBACtB,WAAW,CAAC,EAAE,CAAC,GAAG,GAAG,CAAC;YACxB,CAAC;QACH,CAAC;aAAM,CAAC;YACN,KAAK,IAAI,EAAE,GAAG,CAAC,EAAE,EAAE,GAAG,GAAG,EAAE,EAAE,EAAE,EAAE,CAAC;gBAChC,IAAI,GAAG,GAAG,CAAC,EAAE,GAAG,GAAG,CAAC,CAAC;gBACrB,IAAI,GAAG,GAAG,CAAC,EAAE,GAAG,GAAG,CAAC,CAAC;gBACrB,MAAM,MAAM,GAAG,EAAE,GAAG,GAAG,CAAC;gBACxB,MAAM,UAAU,GAAG,MAAM,GAAG,GAAG,CAAC;gBAEhC,KAAK,IAAI,EAAE,GAAG,CAAC,EAAE,EAAE,GAAG,GAAG,EAAE,EAAE,EAAE,EAAE,CAAC;oBAChC,MAAM,KAAK,GAAG,MAAM,GAAG,EAAE,CAAC;oBAE1B,UAAU;oBACV,MAAM,EAAE,GAAG,QAAQ,CAAC,KAAK,CAAC,CAAC;oBAC3B,IAAI,EAAE,KAAK,CAAC,EAAE,CAAC;wBACb,GAAG,IAAI,EAAE,GAAG,KAAK,CAAC,EAAE,CAAC,CAAC;oBACxB,CAAC;oBAED,UAAU;oBACV,MAAM,EAAE,GAAG,QAAQ,CAAC,KAAK,CAAC,CAAC;oBAC3B,IAAI,EAAE,KAAK,CAAC,EAAE,CAAC;wBACb,GAAG,IAAI,EAAE,GAAG,KAAK,CAAC,EAAE,CAAC,CAAC;oBACxB,CAAC;oBAED,IAAI,SAAS,GAAG,CAAC,CAAC;oBAClB,IAAI,SAAS,GAAG,CAAC,CAAC;oBAClB,MAAM,MAAM,GAAG,UAAU,GAAG,EAAE,GAAG,GAAG,CAAC;oBAErC,KAAK,IAAI,GAAG,GAAG,CAAC,EAAE,GAAG,GAAG,GAAG,EAAE,GAAG,EAAE,EAAE,CAAC;wBACnC,MAAM,KAAK,GAAG,MAAM,GAAG,GAAG,CAAC;wBAE3B,MAAM,GAAG,GAAG,SAAS,CAAC,KAAK,CAAC,CAAC;wBAC7B,IAAI,GAAG,KAAK,CAAC,EAAE,CAAC;4BACd,SAAS,IAAI,GAAG,GAAG,KAAK,CAAC,GAAG,CAAC,CAAC;wBAChC,CAAC;wBAED,MAAM,GAAG,GAAG,SAAS,CAAC,KAAK,CAAC,CAAC;wBAC7B,IAAI,GAAG,KAAK,CAAC,EAAE,CAAC;4BACd,SAAS,IAAI,GAAG,GAAG,KAAK,CAAC,GAAG,CAAC,CAAC;wBAChC,CAAC;oBACH,CAAC;oBAED,IAAI,SAAS,KAAK,CAAC,EAAE,CAAC;wBACpB,GAAG,IAAI,SAAS,GAAG,KAAK,CAAC,EAAE,CAAC,CAAC;oBAC/B,CAAC;oBACD,IAAI,SAAS,KAAK,CAAC,EAAE,CAAC;wBACpB,GAAG,IAAI,SAAS,GAAG,KAAK,CAAC,EAAE,CAAC,CAAC;oBAC/B,CAAC;gBACH,CAAC;gBAED,WAAW,CAAC,EAAE,CAAC,GAAG,GAAG,CAAC;gBACtB,WAAW,CAAC,EAAE,CAAC,GAAG,GAAG,CAAC;gBACtB,WAAW,CAAC,EAAE,CAAC,GAAG,GAAG,CAAC;gBACtB,WAAW,CAAC,EAAE,CAAC,GAAG,GAAG,CAAC;YACxB,CAAC;QACH,CAAC;QAED,6BAA6B;QAC7B,KAAK,IAAI,EAAE,GAAG,CAAC,EAAE,EAAE,GAAG,GAAG,EAAE,EAAE,EAAE,EAAE,CAAC;YAChC,MAAM,EAAE,GAAG,KAAK,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC;YACvB,MAAM,EAAE,GAAG,KAAK,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC;YACvB,MAAM,GAAG,GAAG,WAAW,CAAC,EAAE,CAAC,CAAC;YAC5B,MAAM,GAAG,GAAG,WAAW,CAAC,EAAE,CAAC,CAAC;YAC5B,MAAM,GAAG,GAAG,WAAW,CAAC,EAAE,CAAC,CAAC;YAC5B,MAAM,GAAG,GAAG,WAAW,CAAC,EAAE,CAAC,CAAC;YAE5B,IAAI,CAAC,CAAC,CAAC,CAAC,EAAE,CAAC,GAAG,EAAE,GAAG,IAAI,GAAG,GAAG,GAAG,IAAI,GAAG,GAAG,CAAC;YAC3C,IAAI,CAAC,CAAC,CAAC,CAAC,EAAE,CAAC,GAAG,EAAE,GAAG,QAAQ,GAAG,GAAG,GAAG,QAAQ,GAAG,GAAG,CAAC;YAEnD,IAAI,CAAC,CAAC,CAAC,CAAC,EAAE,CAAC,GAAG,EAAE,GAAG,IAAI,GAAG,GAAG,GAAG,IAAI,GAAG,GAAG,CAAC;YAC3C,IAAI,CAAC,CAAC,CAAC,CAAC,EAAE,CAAC,GAAG,EAAE,GAAG,QAAQ,GAAG,GAAG,GAAG,QAAQ,GAAG,GAAG,CAAC;QACrD,CAAC;QAED,kCAAkC;QAClC,IAAI,QAAQ,GAAG,CAAC,CAAC;QACjB,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC;YAC3B,KAAK,IAAI,EAAE,GAAG,CAAC,EAAE,EAAE,GAAG,GAAG,EAAE,EAAE,EAAE,EAAE,CAAC;gBAChC,QAAQ,GAAG,IAAI,CAAC,GAAG,CAAC,QAAQ,EAAE,IAAI,CAAC,GAAG,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC,EAAE,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC;gBAChE,QAAQ,GAAG,IAAI,CAAC,GAAG,CAAC,QAAQ,EAAE,IAAI,CAAC,GAAG,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC,EAAE,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC;YAClE,CAAC;QACH,CAAC;QAED,6EAA6E;QAC7E,CAAC,CAAC,EAAE,IAAI,CAAC,GAAG,CAAC,IAAI,EAAE,CAAC,CAAC,CAAC;QACtB,CAAC,CAAC,EAAE,IAAI,CAAC,GAAG,CAAC,IAAI,EAAE,CAAC,CAAC,CAAC;QAEtB,IAAI,QAAQ,GAAG,IAAI,CAAC,SAAS,EAAE,CAAC;YAC9B,8DAA8D;YAC9D,0DAA0D;YAC1D,IAAI,GAAG,KAAK,CAAC,EAAE,CAAC;gBACd,oHAAoH;gBACpH,OAAO;oBACL,MAAM,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;oBACpF,MAAM,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;oBACpF,OAAO,EAAE,CAAC,CAAC,WAAW,CAAC,CAAC,CAAC,EAAE,WAAW,CAAC,CAAC,CAAC,EAAE,WAAW,CAAC,CAAC,CAAC,EAAE,WAAW,CAAC,CAAC,CAAC,CAAC,EAAE,CAAC,WAAW,CAAC,CAAC,CAAC,EAAE,WAAW,CAAC,CAAC,CAAC,EAAE,WAAW,CAAC,CAAC,CAAC,EAAE,WAAW,CAAC,CAAC,CAAC,CAAC,CAAC;oBAC7I,OAAO,EAAE,CAAC,CAAC,WAAW,CAAC,CAAC,CAAC,EAAE,WAAW,CAAC,CAAC,CAAC,EAAE,WAAW,CAAC,CAAC,CAAC,EAAE,WAAW,CAAC,CAAC,CAAC,CAAC,EAAE,CAAC,WAAW,CAAC,CAAC,CAAC,EAAE,WAAW,CAAC,CAAC,CAAC,EAAE,WAAW,CAAC,CAAC,CAAC,EAAE,WAAW,CAAC,CAAC,CAAC,CAAC,CAAC;oBAC7I,UAAU,EAAE,CAAC,GAAG,CAAC;iBAClB,CAAC;YACJ,CAAC;iBAAM,CAAC;gBACN,4FAA4F;gBAC5F,MAAM,GAAG,GAAG,IAAI,KAAK,CAAS,GAAG,CAAC,CAAC;gBACnC,MAAM,GAAG,GAAG,IAAI,KAAK,CAAS,GAAG,CAAC,CAAC;gBACnC,MAAM,GAAG,GAAG,IAAI,KAAK,CAAS,GAAG,CAAC,CAAC;gBACnC,MAAM,GAAG,GAAG,IAAI,KAAK,CAAS,GAAG,CAAC,CAAC;gBACnC,MAAM,IAAI,GAAG,IAAI,KAAK,CAAS,GAAG,CAAC,CAAC;gBACpC,MAAM,IAAI,GAAG,IAAI,KAAK,CAAS,GAAG,CAAC,CAAC;gBACpC,MAAM,IAAI,GAAG,IAAI,KAAK,CAAS,GAAG,CAAC,CAAC;gBACpC,MAAM,IAAI,GAAG,IAAI,KAAK,CAAS,GAAG,CAAC,CAAC;gBACpC,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,GAAG,EAAE,CAAC,EAAE,EAAE,CAAC;oBAC7B,GAAG,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;oBACjB,GAAG,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;oBACjB,GAAG,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;oBACjB,GAAG,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;oBACjB,IAAI,CAAC,CAAC,CAAC,GAAG,WAAW,CAAC,CAAC,CAAC,CAAC;oBACzB,IAAI,CAAC,CAAC,CAAC,GAAG,WAAW,CAAC,CAAC,CAAC,CAAC;oBACzB,IAAI,CAAC,CAAC,CAAC,GAAG,WAAW,CAAC,CAAC,CAAC,CAAC;oBACzB,IAAI,CAAC,CAAC,CAAC,GAAG,WAAW,CAAC,CAAC,CAAC,CAAC;gBAC3B,CAAC;gBACD,OAAO;oBACL,MAAM,EAAE,CAAC,GAAG,EAAE,GAAG,CAAC;oBAClB,MAAM,EAAE,CAAC,GAAG,EAAE,GAAG,CAAC;oBAClB,OAAO,EAAE,CAAC,IAAI,EAAE,IAAI,CAAC;oBACrB,OAAO,EAAE,CAAC,IAAI,EAAE,IAAI,CAAC;oBACrB,UAAU,EAAE,CAAC,GAAG,CAAC;iBAClB,CAAC;YACJ,CAAC;QACH,CAAC;IACH,CAAC;IAED,MAAM,IAAI,mBAAmB,CAC3B,wCAAwC,IAAI,CAAC,aAAa,yCAAyC,IAAI,CAAC,SAAS,GAAG,CACrH,CAAC;AACJ,CAAC;AAED,8EAA8E;AAC9E,gDAAgD;AAChD,8EAA8E;AAE9E;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;GAwDG;AACH,MAAM,UAAU,oBAAoB,CAClC,YAAsB,EACtB,OAAmB;IAEnB,MAAM,EACJ,KAAK,EACL,MAAM,EACN,UAAU,EACV,WAAW,EACX,SAAS,GAAG,KAAK,EACjB,aAAa,GAAG,EAAE,EAClB,IAAI,EACJ,eAAe,EACf,MAAM,GACP,GAAG,OAAO,CAAC;IAEZ,IAAI,eAAe,KAAK,SAAS,IAAI,YAAY,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,eAAe,EAAE,CAAC;QACzE,MAAM,IAAI,qBAAqB,CAC7B,6BAA6B,YAAY,CAAC,CAAC,CAAC,CAAC,CAAC,sBAAsB,eAAe,qBAAqB,CACzG,CAAC;IACJ,CAAC;IAED,MAAM,CAAC,GAAG,MAAM,GAAG,KAAK,CAAC;IACzB,MAAM,MAAM,GAAG,IAAI,IAAI,CAAC,GAAG,IAAI,CAAC;IAChC,MAAM,SAAS,GAAkB;QAC/B,EAAE,GAAG,EAAE,CAAC,EAAE,CAAC,EAAE,YAAY,CAAC,CAAC,CAAC,KAAK,EAAc,EAAE,CAAC,EAAE,YAAY,CAAC,CAAC,CAAC,KAAK,EAAc,EAAE;KACzF,CAAC;IACF,IAAI,CAAC,GAAG,YAAY,CAAC,CAAC,CAAC,KAAK,EAAc,CAAC;IAC3C,IAAI,CAAC,GAAG,YAAY,CAAC,CAAC,CAAC,KAAK,EAAc,CAAC;IAE3C,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,KAAK,EAAE,CAAC,EAAE,EAAE,CAAC;QAC/B,8DAA8D;QAC9D,uEAAuE;QACvE,yEAAyE;QACzE,uEAAuE;QACvE,wEAAwE;QACxE,sEAAsE;QACtE,4DAA4D;QAC5D,IAAI,SAAS,GAAG,CAAC,CAAC;QAClB,IAAI,KAAK,GAAG,CAAC,CAAC;QACd,IAAI,QAAQ,GAAG,CAAC,CAAC;QACjB,IAAI,cAAc,GAAG,CAAC,CAAC;QACvB,+EAA+E;QAC/E,MAAM,SAAS,GAAG,CAAC,GAAG,KAAK,CAAC;QAC5B,OAAO,SAAS,GAAG,SAAS,EAAE,CAAC;YAC7B,MAAM,MAAM,GAAG,IAAI,CAAC,GAAG,CAAC,KAAK,EAAE,SAAS,CAAC,CAAC;YAC1C,IAAI,MAAoC,CAAC;YACzC,IAAI,aAAa,GAAG,KAAK,CAAC;YAC1B,IAAI,IAAI,GAAG,MAAM,CAAC;YAClB,OAAO,IAAI,IAAI,MAAM,EAAE,CAAC;gBACtB,IAAI,CAAC;oBACH,MAAM,GAAG,aAAa,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,EAAE,IAAI,EAAE,UAAU,EAAE,WAAW,EAAE;wBAC9D,SAAS;wBACT,aAAa;qBACd,CAAC,CAAC;oBACH,aAAa,GAAG,IAAI,CAAC;oBACrB,MAAM;gBACR,CAAC;gBAAC,MAAM,CAAC;oBACP,IAAI,IAAI,CAAC,CAAC;oBACV,QAAQ,EAAE,CAAC;gBACb,CAAC;YACH,CAAC;YACD,IAAI,CAAC,aAAa,IAAI,MAAM,KAAK,SAAS,EAAE,CAAC;gBAC3C,MAAM,IAAI,mBAAmB,CAC3B,mEAAmE,MAAM,UAAU,CAAC,8CAA8C,CACnI,CAAC;YACJ,CAAC;YACD,+EAA+E;YAC/E,4EAA4E;YAC5E,+DAA+D;YAC/D,MAAM,IAAI,GAAG,CAAC,CAAC,KAAK,EAAE,CAAC;YACvB,MAAM,IAAI,GAAG,CAAC,CAAC,KAAK,EAAE,CAAC;YAEvB,uEAAuE;YACvE,6EAA6E;YAC7E,kFAAkF;YAClF,MAAM,EAAE,GAAG,KAAK,CAAC,CAAC,CAAC,CAAC;YACpB,MAAM,EAAE,GAAG,KAAK,CAAC,CAAC,CAAC,CAAC;YACpB,MAAM,GAAG,GAAG,IAAI,CAAC,MAAM,CAAC;YACxB,MAAM,GAAG,GAAG,MAAM,CAAC,OAAO,CAAC,CAAC,CAAC,CAAC;YAC9B,MAAM,GAAG,GAAG,MAAM,CAAC,OAAO,CAAC,CAAC,CAAC,CAAC;YAC9B,MAAM,GAAG,GAAG,MAAM,CAAC,OAAO,CAAC,CAAC,CAAC,CAAC;YAC9B,MAAM,GAAG,GAAG,MAAM,CAAC,OAAO,CAAC,CAAC,CAAC,CAAC;YAC9B,MAAM,GAAG,GAAG,IAAI,GAAG,EAAE,CAAC;YACtB,MAAM,GAAG,GAAG,IAAI,GAAG,EAAE,CAAC;YACtB,MAAM,IAAI,GAAG,GAAG,GAAG,CAAC,GAAG,CAAC;YACxB,MAAM,IAAI,GAAG,GAAG,GAAG,CAAC,GAAG,CAAC;YACxB,IAAI,GAAG,KAAK,CAAC,EAAE,CAAC;gBACd,IAAI,CAAC,CAAC,CAAC,IAAI,GAAG,GAAG,GAAG,CAAC,CAAC,CAAC,GAAG,GAAG,GAAG,GAAG,CAAC,CAAC,CAAC,CAAC;gBAAC,IAAI,CAAC,CAAC,CAAC,IAAI,IAAI,GAAG,GAAG,CAAC,CAAC,CAAC,GAAG,IAAI,GAAG,GAAG,CAAC,CAAC,CAAC,CAAC;gBACjF,IAAI,CAAC,CAAC,CAAC,IAAI,GAAG,GAAG,GAAG,CAAC,CAAC,CAAC,GAAG,GAAG,GAAG,GAAG,CAAC,CAAC,CAAC,CAAC;gBAAC,IAAI,CAAC,CAAC,CAAC,IAAI,IAAI,GAAG,GAAG,CAAC,CAAC,CAAC,GAAG,IAAI,GAAG,GAAG,CAAC,CAAC,CAAC,CAAC;gBACjF,IAAI,CAAC,CAAC,CAAC,IAAI,GAAG,GAAG,GAAG,CAAC,CAAC,CAAC,GAAG,GAAG,GAAG,GAAG,CAAC,CAAC,CAAC,CAAC;gBAAC,IAAI,CAAC,CAAC,CAAC,IAAI,IAAI,GAAG,GAAG,CAAC,CAAC,CAAC,GAAG,IAAI,GAAG,GAAG,CAAC,CAAC,CAAC,CAAC;gBACjF,IAAI,CAAC,CAAC,CAAC,IAAI,GAAG,GAAG,GAAG,CAAC,CAAC,CAAC,GAAG,GAAG,GAAG,GAAG,CAAC,CAAC,CAAC,CAAC;gBAAC,IAAI,CAAC,CAAC,CAAC,IAAI,IAAI,GAAG,GAAG,CAAC,CAAC,CAAC,GAAG,IAAI,GAAG,GAAG,CAAC,CAAC,CAAC,CAAC;YACnF,CAAC;iBAAM,CAAC;gBACN,KAAK,IAAI,EAAE,GAAG,CAAC,EAAE,EAAE,GAAG,GAAG,EAAE,EAAE,EAAE,EAAE,CAAC;oBAChC,IAAI,CAAC,EAAE,CAAC,IAAI,GAAG,GAAG,GAAG,CAAC,EAAE,CAAC,GAAG,GAAG,GAAG,GAAG,CAAC,EAAE,CAAC,CAAC;oBAC1C,IAAI,CAAC,EAAE,CAAC,IAAI,IAAI,GAAG,GAAG,CAAC,EAAE,CAAC,GAAG,IAAI,GAAG,GAAG,CAAC,EAAE,CAAC,CAAC;gBAC9C,CAAC;YACH,CAAC;YACD,CAAC,GAAG,IAAI,CAAC;YACT,CAAC,GAAG,IAAI,CAAC;YAET,SAAS,IAAI,IAAI,CAAC;YAClB,KAAK,GAAG,IAAI,CAAC,CAAC,2DAA2D;YACzE,cAAc,GAAG,MAAM,CAAC,UAAU,CAAC;QACrC,CAAC;QACD,SAAS,CAAC,IAAI,CAAC,EAAE,GAAG,EAAE,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC,KAAK,EAAE,EAAE,CAAC,EAAE,CAAC,CAAC,KAAK,EAAE,EAAE,CAAC,CAAC;QACjE,IAAI,MAAM,KAAK,SAAS,EAAE,CAAC;YACzB,MAAM,CAAC,EAAE,IAAI,EAAE,CAAC,EAAE,UAAU,EAAE,cAAc,EAAE,QAAQ,EAAE,CAAC,CAAC;QAC5D,CAAC;IACH,CAAC;IAED,OAAO,SAAS,CAAC;AACnB,CAAC"}
|
|
@@ -1 +1 @@
|
|
|
1
|
-
{"version":3,"file":"weyl-lowering.d.ts","sourceRoot":"","sources":["../../src/numerical/weyl-lowering.ts"],"names":[],"mappings":"AAAA;;;;;;;;;;;;;;;;;;;;;;;;;GAyBG;AAMH;;;;;;GAMG;AAMH,UAAU,UAAU;IAClB,6CAA6C;IAC7C,OAAO,EAAE,MAAM,EAAE,EAAE,EAAE,EAAE,CAAC;IACxB,+CAA+C;IAC/C,KAAK,EAAE,MAAM,EAAE,EAAE,CAAC;IAClB,sCAAsC;IACtC,WAAW,EAAE,MAAM,CAAC;IACpB,4EAA4E;IAC5E,MAAM,EAAE,MAAM,EAAE,EAAE,GAAG,YAAY,CAAC;IAClC,0EAA0E;IAC1E,aAAa,EAAE,MAAM,EAAE,EAAE,GAAG,YAAY,CAAC;CAC1C;
|
|
1
|
+
{"version":3,"file":"weyl-lowering.d.ts","sourceRoot":"","sources":["../../src/numerical/weyl-lowering.ts"],"names":[],"mappings":"AAAA;;;;;;;;;;;;;;;;;;;;;;;;;GAyBG;AAMH;;;;;;GAMG;AAMH,UAAU,UAAU;IAClB,6CAA6C;IAC7C,OAAO,EAAE,MAAM,EAAE,EAAE,EAAE,EAAE,CAAC;IACxB,+CAA+C;IAC/C,KAAK,EAAE,MAAM,EAAE,EAAE,CAAC;IAClB,sCAAsC;IACtC,WAAW,EAAE,MAAM,CAAC;IACpB,4EAA4E;IAC5E,MAAM,EAAE,MAAM,EAAE,EAAE,GAAG,YAAY,CAAC;IAClC,0EAA0E;IAC1E,aAAa,EAAE,MAAM,EAAE,EAAE,GAAG,YAAY,CAAC;CAC1C;AAsED;;;;;;;;;;;GAWG;AACH,wBAAgB,iBAAiB,CAAC,KAAK,EAAE,UAAU,GAAG,MAAM,EAAE,EAAE,EAAE,EAAE,CA4FnE"}
|