diff-grok 1.0.10 → 1.2.0

package/src/solver-tools/ab5-method.js

@@ -0,0 +1,334 @@
+/* The Adams-Bashforth-Moulton predictor-corrector method of order 5 (AB5).
+
+   Uses the 5-step Adams-Bashforth formula as predictor and the 4-step
+   Adams-Moulton formula as corrector in PECE mode (Predict-Evaluate-
+   Correct-Evaluate), with Milne's device for local error estimation.
+   Bootstrapping is performed using the Runge-Kutta-Fehlberg 4(5) method.
+
+   References:
+   [1] E. Hairer, S.P. Norsett, G. Wanner, "Solving Ordinary Differential
+       Equations I: Nonstiff Problems", Springer, 2nd ed., 1993.
+   [2] L.F. Shampine, M.K. Gordon, "Computer Solution of Ordinary Differential
+       Equations: The Initial Value Problem", W.H. Freeman, 1975.
+*/
+import { max, abs, SAFETY, PSHRNK, PSGROW, REDUCE_COEF, GROW_COEF, ERR_CONTR, TINY, ERROR_MSG } from './solver-defs';
+// Adams-Bashforth 5-step coefficients (pre-divided by 720)
+// y_{n+1} = y_n + h * (AB0*f_n + AB1*f_{n-1} + AB2*f_{n-2} + AB3*f_{n-3} + AB4*f_{n-4})
+const AB0 = 1901.0 / 720.0;
+const AB1 = -2774.0 / 720.0;
+const AB2 = 2616.0 / 720.0;
+const AB3 = -1274.0 / 720.0;
+const AB4 = 251.0 / 720.0;
+// Adams-Moulton 4-step coefficients (pre-divided by 720)
+// y_{n+1} = y_n + h * (AM_NEW*f_{n+1} + AM0*f_n + AM1*f_{n-1} + AM2*f_{n-2} + AM3*f_{n-3})
+const AM_NEW = 251.0 / 720.0;
+const AM0 = 646.0 / 720.0;
+const AM1 = -264.0 / 720.0;
+const AM2 = 106.0 / 720.0;
+const AM3 = -19.0 / 720.0;
+// Milne error estimation coefficient for the AB5-AM4 pair: 2/27 = 38/513
+const MILNE = 2.0 / 27.0;
+// Number of uniformly-spaced history points required for AB5
+const HIST_LEN = 5;
+// Fehlberg 4(5) Butcher tableau coefficients (used during bootstrap phase)
+const C2 = 1.0 / 4.0;
+const C3 = 3.0 / 8.0;
+const C4 = 12.0 / 13.0;
+const C6 = 1.0 / 2.0;
+const A21 = 1.0 / 4.0;
+const A31 = 3.0 / 32.0;
+const A32 = 9.0 / 32.0;
+const A41 = 1932.0 / 2197.0;
+const A42 = -7200.0 / 2197.0;
+const A43 = 7296.0 / 2197.0;
+const A51 = 439.0 / 216.0;
+const A52 = -8.0;
+const A53 = 3680.0 / 513.0;
+const A54 = -845.0 / 4104.0;
+const A61 = -8.0 / 27.0;
+const A62 = 2.0;
+const A63 = -3544.0 / 2565.0;
+const A64 = 1859.0 / 4104.0;
+const A65 = -11.0 / 40.0;
+const B1 = 25.0 / 216.0;
+const B3 = 1408.0 / 2565.0;
+const B4 = 2197.0 / 4104.0;
+const B5 = -1.0 / 5.0;
+const RKE1 = 1.0 / 360.0;
+const RKE3 = -128.0 / 4275.0;
+const RKE4 = -2197.0 / 75240.0;
+const RKE5 = 1.0 / 50.0;
+const RKE6 = 2.0 / 55.0;
+/** Solve initial value problem using the Adams-Bashforth-Moulton predictor-corrector method of order 5
+ * @param odes initial value problem for ordinary differential equations
+ * @param callback computations control callback
+ * @returns solution of the problem
+ * @example
+ * // Example. Solve the following initial value problem using the AB5 method:
+ * //   dx/dt = x + y - t
+ * //   dy/dt = x * y + t
+ * //   x(0) = 1
+ * //   y(0) = -1
+ * // on [0, 2] with the step 0.01.
+ * import { ODEs, ab5 } from 'diff-grok';
+ *
+ * const task: ODEs = {
+ *   name: 'Example', // name of your model
+ *   arg: {
+ *       name: 't',  // name of the argument
+ *       start: 0,   // initial value of the argument
+ *       finish: 2,  // final value of the argument
+ *       step: 0.01, // solution grid step
+ *   },
+ *   initial: [1, -1], // initial values
+ *   func: (t: number, y: Float64Array, output: Float64Array) => { // right-hand side of the system
+ *     output[0] = y[0] + y[1] - t; // 1-st equation
+ *     output[1] = y[0] * y[1] + t; // 2-nd equation
+ *   },
+ *   tolerance: 1e-7, // tolerance
+ *   solutionColNames: ['x', 'y'], // names of solution functions
+ * };
+ *
+ * try {
+ *   // Solve the problem using the AB5 method
+ *   const solution = ab5(task);
+ *
+ *   // Print a table with the results
+ *   console.log(task.arg.name, task.solutionColNames[0], task.solutionColNames[1]);
+ *   const length = solution[0].length;
+ *   for (let i = 0; i < length; ++i)
+ *     console.log(solution[0][i], solution[1][i], solution[2][i]);
+ * } catch (err) {
+ *   console.log('Solver failed: ', err instanceof Error ? err.message : 'Unknown problem!');
+ * }
+*/
+export function ab5(odes, callback) {
+    const f = odes.func;
+    // operating variables
+    const t0 = odes.arg.start;
+    const t1 = odes.arg.finish;
+    let h = odes.arg.step;
+    const hDataframe = h;
+    const hMax = hDataframe * 10;
+    const tolerance = odes.tolerance;
+    const rowCount = Math.trunc((t1 - t0) / h) + 1;
+    const dim = odes.initial.length;
+    // result arrays
+    const tArr = new Float64Array(rowCount);
+    const yArrs = Array(dim);
+    for (let i = 0; i < dim; ++i)
+        yArrs[i] = new Float64Array(rowCount);
+    // method routine
+    let timeDataframe = t0 + hDataframe;
+    let t = t0;
+    let tPrev = t0;
+    let hNext = 0.0;
+    let flag = true;
+    let index = 1;
+    let errmax = 0;
+    let hTemp = 0;
+    let tNew = 0;
+    // working buffers
+    const y = new Float64Array(odes.initial);
+    const yPrev = new Float64Array(odes.initial);
+    const dydt = new Float64Array(dim);
+    const yScale = new Float64Array(dim);
+    const yTemp = new Float64Array(dim);
+    const yErr = new Float64Array(dim);
+    // AB5-AM4 history: fBuf[0] = f(t_n, y_n), fBuf[4] = f(t_{n-4}, y_{n-4})
+    const fBuf = [];
+    for (let i = 0; i < HIST_LEN; ++i)
+        fBuf[i] = new Float64Array(dim);
+    let histCount = 0;
+    let lastH = 0;
+    // AB5-AM4 corrector buffers
+    const yCorr = new Float64Array(dim);
+    const fPred = new Float64Array(dim);
+    // RKF45 stage buffers (bootstrap)
+    const k1 = new Float64Array(dim);
+    const k2 = new Float64Array(dim);
+    const k3 = new Float64Array(dim);
+    const k4 = new Float64Array(dim);
+    const k5 = new Float64Array(dim);
+    const k6 = new Float64Array(dim);
+    // store initial f-value
+    f(t, y, fBuf[0]);
+    histCount = 1;
+    // 1. SOLUTION AT THE POINT t0
+    tArr[0] = t0;
+    for (let i = 0; i < dim; ++i)
+        yArrs[i][0] = y[i];
+    // 2. COMPUTE NUMERICAL SOLUTION FOR THE POINTS FROM THE INTERVAL (t0, t1)
+    while (flag) {
+        // compute derivative
+        f(t, y, dydt);
+        // check whether to go on computations
+        if (callback)
+            callback.onIterationStart();
+        // compute scale vector
+        for (let i = 0; i < dim; ++i)
+            yScale[i] = abs(y[i]) + h * abs(dydt[i]) + TINY;
+        // check end point
+        if (t + h > t1) {
+            h = t1 - t;
+            flag = false;
+        }
+        // if step size changed, history is invalid (AB5 needs uniform spacing)
+        if (histCount > 1 && h !== lastH) {
+            histCount = 1;
+            f(t, y, fBuf[0]);
+        }
+        // adaptive step loop
+        while (true) {
+            if (histCount < HIST_LEN) {
+                // BOOTSTRAP PHASE: RKF45 single step to build up history
+                // Stage 1: k1 = f(t, y)
+                f(t, y, k1);
+                // Stage 2: k2 = f(t + C2*h, y + h*A21*k1)
+                for (let i = 0; i < dim; ++i)
+                    yTemp[i] = y[i] + h * A21 * k1[i];
+                f(t + C2 * h, yTemp, k2);
+                // Stage 3: k3 = f(t + C3*h, y + h*(A31*k1 + A32*k2))
+                for (let i = 0; i < dim; ++i)
+                    yTemp[i] = y[i] + h * (A31 * k1[i] + A32 * k2[i]);
+                f(t + C3 * h, yTemp, k3);
+                // Stage 4: k4 = f(t + C4*h, y + h*(A41*k1 + A42*k2 + A43*k3))
+                for (let i = 0; i < dim; ++i)
+                    yTemp[i] = y[i] + h * (A41 * k1[i] + A42 * k2[i] + A43 * k3[i]);
+                f(t + C4 * h, yTemp, k4);
+                // Stage 5: k5 = f(t + h, y + h*(A51*k1 + A52*k2 + A53*k3 + A54*k4))
+                for (let i = 0; i < dim; ++i)
+                    yTemp[i] = y[i] + h * (A51 * k1[i] + A52 * k2[i] + A53 * k3[i] + A54 * k4[i]);
+                f(t + h, yTemp, k5);
+                // Stage 6: k6 = f(t + C6*h, y + h*(A61*k1 + A62*k2 + A63*k3 + A64*k4 + A65*k5))
+                for (let i = 0; i < dim; ++i)
+                    yTemp[i] = y[i] + h * (A61 * k1[i] + A62 * k2[i] + A63 * k3[i] + A64 * k4[i] + A65 * k5[i]);
+                f(t + C6 * h, yTemp, k6);
+                // 4th-order solution
+                for (let i = 0; i < dim; ++i)
+                    yTemp[i] = y[i] + h * (B1 * k1[i] + B3 * k3[i] + B4 * k4[i] + B5 * k5[i]);
+                // Error estimate: difference between 4th and 5th order solutions
+                for (let i = 0; i < dim; ++i)
+                    yErr[i] = h * (RKE1 * k1[i] + RKE3 * k3[i] + RKE4 * k4[i] + RKE5 * k5[i] + RKE6 * k6[i]);
+                // estimating error
+                errmax = 0;
+                for (let i = 0; i < dim; ++i)
+                    errmax = max(errmax, abs(yErr[i] / yScale[i]));
+                errmax /= tolerance;
+                // processing the error obtained
+                if (errmax > 1) {
+                    hTemp = SAFETY * h * Math.pow(errmax, PSHRNK);
+                    h = max(hTemp, REDUCE_COEF * h);
+                    tNew = t + h;
+                    if (tNew == t)
+                        throw new Error(ERROR_MSG.AB5_FAILS);
+                    // step size changed, reset history if needed
+                    if (histCount > 1) {
+                        histCount = 1;
+                        f(t, y, fBuf[0]);
+                    }
+                }
+                else {
+                    // step accepted
+                    t = t + h;
+                    for (let i = 0; i < dim; ++i)
+                        y[i] = yTemp[i];
+                    // shift history and store new f-value
+                    for (let j = Math.min(histCount, HIST_LEN - 1); j > 0; --j)
+                        fBuf[j].set(fBuf[j - 1]);
+                    f(t, y, fBuf[0]);
+                    histCount++;
+                    lastH = h;
+                    // during bootstrap, keep h constant for uniform spacing
+                    hNext = h;
+                    if (hNext > hMax)
+                        hNext = hMax;
+                    break;
+                }
+            }
+            else {
+                // AB5-AM4 PHASE: Predictor-Corrector (PECE)
+                // Predict (Adams-Bashforth 5-step)
+                for (let i = 0; i < dim; ++i) {
+                    yTemp[i] = y[i] + h * (AB0 * fBuf[0][i] + AB1 * fBuf[1][i] +
+                        AB2 * fBuf[2][i] + AB3 * fBuf[3][i] + AB4 * fBuf[4][i]);
+                }
+                // Evaluate f at predicted point
+                f(t + h, yTemp, fPred);
+                // Correct (Adams-Moulton 4-step)
+                for (let i = 0; i < dim; ++i) {
+                    yCorr[i] = y[i] + h * (AM_NEW * fPred[i] + AM0 * fBuf[0][i] +
+                        AM1 * fBuf[1][i] + AM2 * fBuf[2][i] + AM3 * fBuf[3][i]);
+                }
+                // Error estimate (Milne's device)
+                for (let i = 0; i < dim; ++i)
+                    yErr[i] = MILNE * (yCorr[i] - yTemp[i]);
+                // estimating error
+                errmax = 0;
+                for (let i = 0; i < dim; ++i)
+                    errmax = max(errmax, abs(yErr[i] / yScale[i]));
+                errmax /= tolerance;
+                // processing the error obtained
+                if (errmax > 1) {
+                    // step rejected: shrink step and re-bootstrap
+                    hTemp = SAFETY * h * Math.pow(errmax, PSHRNK);
+                    h = max(hTemp, REDUCE_COEF * h);
+                    tNew = t + h;
+                    if (tNew == t)
+                        throw new Error(ERROR_MSG.AB5_FAILS);
+                    histCount = 1;
+                    f(t, y, fBuf[0]);
+                }
+                else {
+                    if (errmax > ERR_CONTR)
+                        hNext = SAFETY * h * Math.pow(errmax, PSGROW);
+                    else
+                        hNext = GROW_COEF * h;
+                    if (hNext > hMax)
+                        hNext = hMax;
+                    t = t + h;
+                    for (let i = 0; i < dim; ++i)
+                        y[i] = yCorr[i];
+                    // shift history and store f at corrected point (PECE final evaluate)
+                    const recycled = fBuf[HIST_LEN - 1];
+                    for (let j = HIST_LEN - 1; j > 0; --j)
+                        fBuf[j] = fBuf[j - 1];
+                    fBuf[0] = recycled;
+                    f(t, y, fBuf[0]);
+                    lastH = h;
+                    break;
+                }
+            }
+        } // while (true)
+        // compute linearly interpolated results and store them in dataframe
+        while (timeDataframe < t) {
+            const cLeft = (t - timeDataframe) / (t - tPrev);
+            const cRight = 1.0 - cLeft;
+            tArr[index] = timeDataframe;
+            for (let j = 0; j < dim; ++j)
+                yArrs[j][index] = cRight * y[j] + cLeft * yPrev[j];
+            timeDataframe += hDataframe;
+            ++index;
+        }
+        h = hNext;
+        tPrev = t;
+        for (let i = 0; i < dim; ++i)
+            yPrev[i] = y[i];
+    } // while (flag)
+    // perform final callback actions
+    if (callback)
+        callback.onComputationsCompleted();
+    // 3. solution at the point t1
+    tArr[rowCount - 1] = t1;
+    for (let i = 0; i < dim; ++i)
+        yArrs[i][rowCount - 1] = y[i];
+    // 4. prepare output
+    const solution = Array(dim + 1);
+    // independent variable
+    solution[0] = tArr;
+    // functions
+    for (let i = 0; i < dim; ++i)
+        solution[i + 1] = yArrs[i];
+    return solution;
+} // ab5
+//# sourceMappingURL=ab5-method.js.map

package/src/solver-tools/ab5-method.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"ab5-method.js","sourceRoot":"","sources":["ab5-method.ts"],"names":[],"mappings":"AAAA;;;;;;;;;;;;EAYE;AAEF,OAAO,EAAO,GAAG,EAAE,GAAG,EAAE,MAAM,EAAE,MAAM,EAAE,MAAM,EAAE,WAAW,EAAE,SAAS,EACpE,SAAS,EAAE,IAAI,EAAE,SAAS,EAAC,MAAM,eAAe,CAAC;AAGnD,2DAA2D;AAC3D,wFAAwF;AACxF,MAAM,GAAG,GAAG,MAAM,GAAG,KAAK,CAAC;AAC3B,MAAM,GAAG,GAAG,CAAC,MAAM,GAAG,KAAK,CAAC;AAC5B,MAAM,GAAG,GAAG,MAAM,GAAG,KAAK,CAAC;AAC3B,MAAM,GAAG,GAAG,CAAC,MAAM,GAAG,KAAK,CAAC;AAC5B,MAAM,GAAG,GAAG,KAAK,GAAG,KAAK,CAAC;AAE1B,yDAAyD;AACzD,2FAA2F;AAC3F,MAAM,MAAM,GAAG,KAAK,GAAG,KAAK,CAAC;AAC7B,MAAM,GAAG,GAAG,KAAK,GAAG,KAAK,CAAC;AAC1B,MAAM,GAAG,GAAG,CAAC,KAAK,GAAG,KAAK,CAAC;AAC3B,MAAM,GAAG,GAAG,KAAK,GAAG,KAAK,CAAC;AAC1B,MAAM,GAAG,GAAG,CAAC,IAAI,GAAG,KAAK,CAAC;AAE1B,yEAAyE;AACzE,MAAM,KAAK,GAAG,GAAG,GAAG,IAAI,CAAC;AAEzB,6DAA6D;AAC7D,MAAM,QAAQ,GAAG,CAAC,CAAC;AAEnB,2EAA2E;AAC3E,MAAM,EAAE,GAAG,GAAG,GAAG,GAAG,CAAC;AACrB,MAAM,EAAE,GAAG,GAAG,GAAG,GAAG,CAAC;AACrB,MAAM,EAAE,GAAG,IAAI,GAAG,IAAI,CAAC;AACvB,MAAM,EAAE,GAAG,GAAG,GAAG,GAAG,CAAC;AAErB,MAAM,GAAG,GAAG,GAAG,GAAG,GAAG,CAAC;AACtB,MAAM,GAAG,GAAG,GAAG,GAAG,IAAI,CAAC;AACvB,MAAM,GAAG,GAAG,GAAG,GAAG,IAAI,CAAC;AACvB,MAAM,GAAG,GAAG,MAAM,GAAG,MAAM,CAAC;AAC5B,MAAM,GAAG,GAAG,CAAC,MAAM,GAAG,MAAM,CAAC;AAC7B,MAAM,GAAG,GAAG,MAAM,GAAG,MAAM,CAAC;AAC5B,MAAM,GAAG,GAAG,KAAK,GAAG,KAAK,CAAC;AAC1B,MAAM,GAAG,GAAG,CAAC,GAAG,CAAC;AACjB,MAAM,GAAG,GAAG,MAAM,GAAG,KAAK,CAAC;AAC3B,MAAM,GAAG,GAAG,CAAC,KAAK,GAAG,MAAM,CAAC;AAC5B,MAAM,GAAG,GAAG,CAAC,GAAG,GAAG,IAAI,CAAC;AACxB,MAAM,GAAG,GAAG,GAAG,CAAC;AAChB,MAAM,GAAG,GAAG,CAAC,MAAM,GAAG,MAAM,CAAC;AAC7B,MAAM,GAAG,GAAG,MAAM,GAAG,MAAM,CAAC;AAC5B,MAAM,GAAG,GAAG,CAAC,IAAI,GAAG,IAAI,CAAC;AAEzB,MAAM,EAAE,GAAG,IAAI,GAAG,KAAK,CAAC;AACxB,MAAM,EAAE,GAAG,MAAM,GAAG,MAAM,CAAC;AAC3B,MAAM,EAAE,GAAG,MAAM,GAAG,MAAM,CAAC;AAC3B,MAAM,EAAE,GAAG,CAAC,GAAG,GAAG,GAAG,CAAC;AAEtB,MAAM,IAAI,GAAG,GAAG,GAAG,KAAK,CAAC;AACzB,MAAM,IAAI,GAAG,CAAC,KAAK,GAAG,MAAM,CAAC;AAC7B,MAAM,IAAI,GAAG,CAAC,MAAM,GAAG,OAAO,CAAC;AAC/B,MAAM,IAAI,GAAG,GAAG,GAAG,IAAI,CAAC;AACxB,MAAM,IAAI,GAAG,GAAG,GAAG,IAAI,CAAC;AAExB;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;EA0CE;AACF,MAAM,UAAU,GAAG,CAAC,IAAU,EAAE,QAAmB;IACjD,MAAM,CAAC,GAAG,IAAI,CAAC,IAAI,CAAC;IAEpB,sBAAsB;IACtB,MAAM,EAAE,GAAG,IAAI,CAAC,GAAG,CAAC,KAAK,CAAC;IAC1B,MAAM,EAAE,GAAG,IAAI,CAAC,GAAG,CAAC,MAAM,CAAC;IAC3B,IAAI,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,IAAI,CAAC;IACtB,MAAM,UAAU,GAAG,CAAC,CAAC;IACrB,MAAM,IAAI,GAAG,UAAU,GAAG,EAAE,CAAC;IAC7B,MAAM,SAAS,GAAG,IAAI,CAAC,SAAS,CAAC;IAEjC,MAAM,QAAQ,GAAG,IAAI,CAAC,KAAK,CAAC,CAAC,EAAE,GAAG,EAAE,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC;IAC/C,MAAM,GAAG,GAAG,IAAI,CAAC,OAAO,CAAC,MAAM,CAAC;IAEhC,gBAAgB;IAChB,MAAM,IAAI,GAAG,IAAI,YAAY,CAAC,QAAQ,CAAC,CAAC;IACxC,MAAM,KAAK,GAAG,KAAK,CAAe,GAAG,CAAC,CAAC;IAEvC,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,GAAG,EAAE,EAAE,CAAC;QAC1B,KAAK,CAAC,CAAC,CAAC,GAAG,IAAI,YAAY,CAAC,QAAQ,CAAC,CAAC;IAExC,iBAAiB;IACjB,IAAI,aAAa,GAAG,EAAE,GAAG,UAAU,CAAC;IACpC,IAAI,CAAC,GAAG,EAAE,CAAC;IACX,IAAI,KAAK,GAAG,EAAE,CAAC;IACf,IAAI,KAAK,GAAG,GAAG,CAAC;IAChB,IAAI,IAAI,GAAG,IAAI,CAAC;IAChB,IAAI,KAAK,GAAG,CAAC,CAAC;IACd,IAAI,MAAM,GAAG,CAAC,CAAC;IACf,IAAI,KAAK,GAAG,CAAC,CAAC;IACd,IAAI,IAAI,GAAG,CAAC,CAAC;IAEb,kBAAkB;IAClB,MAAM,CAAC,GAAG,IAAI,YAAY,CAAC,IAAI,CAAC,OAAO,CAAC,CAAC;IACzC,MAAM,KAAK,GAAG,IAAI,YAAY,CAAC,IAAI,CAAC,OAAO,CAAC,CAAC;IAC7C,MAAM,IAAI,GAAG,IAAI,YAAY,CAAC,GAAG,CAAC,CAAC;IACnC,MAAM,MAAM,GAAG,IAAI,YAAY,CAAC,GAAG,CAAC,CAAC;IACrC,MAAM,KAAK,GAAG,IAAI,YAAY,CAAC,GAAG,CAAC,CAAC;IACpC,MAAM,IAAI,GAAG,IAAI,YAAY,CAAC,GAAG,CAAC,CAAC;IAEnC,wEAAwE;IACxE,MAAM,IAAI,GAAmB,EAAE,CAAC;IAChC,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,QAAQ,EAAE,EAAE,CAAC;QAC/B,IAAI,CAAC,CAAC,CAAC,GAAG,IAAI,YAAY,CAAC,GAAG,CAAC,CAAC;IAClC,IAAI,SAAS,GAAG,CAAC,CAAC;IAClB,IAAI,KAAK,GAAG,CAAC,CAAC;IAEd,4BAA4B;IAC5B,MAAM,KAAK,GAAG,IAAI,YAAY,CAAC,GAAG,CAAC,CAAC;IACpC,MAAM,KAAK,GAAG,IAAI,YAAY,CAAC,GAAG,CAAC,CAAC;IAEpC,kCAAkC;IAClC,MAAM,EAAE,GAAG,IAAI,YAAY,CAAC,GAAG,CAAC,CAAC;IACjC,MAAM,EAAE,GAAG,IAAI,YAAY,CAAC,GAAG,CAAC,CAAC;IACjC,MAAM,EAAE,GAAG,IAAI,YAAY,CAAC,GAAG,CAAC,CAAC;IACjC,MAAM,EAAE,GAAG,IAAI,YAAY,CAAC,GAAG,CAAC,CAAC;IACjC,MAAM,EAAE,GAAG,IAAI,YAAY,CAAC,GAAG,CAAC,CAAC;IACjC,MAAM,EAAE,GAAG,IAAI,YAAY,CAAC,GAAG,CAAC,CAAC;IAEjC,wBAAwB;IACxB,CAAC,CAAC,CAAC,EAAE,CAAC,EAAE,IAAI,CAAC,CAAC,CAAC,CAAC,CAAC;IACjB,SAAS,GAAG,CAAC,CAAC;IAEd,8BAA8B;IAC9B,IAAI,CAAC,CAAC,CAAC,GAAG,EAAE,CAAC;IACb,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,GAAG,EAAE,EAAE,CAAC;QAC1B,KAAK,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC;IAErB,0EAA0E;IAC1E,OAAO,IAAI,EAAE,CAAC;QACZ,qBAAqB;QACrB,CAAC,CAAC,CAAC,EAAE,CAAC,EAAE,IAAI,CAAC,CAAC;QAEd,sCAAsC;QACtC,IAAI,QAAQ;YACV,QAAQ,CAAC,gBAAgB,EAAE,CAAC;QAE9B,uBAAuB;QACvB,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,GAAG,EAAE,EAAE,CAAC;YAC1B,MAAM,CAAC,CAAC,CAAC,GAAG,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,GAAG,GAAG,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC;QAElD,kBAAkB;QAClB,IAAI,CAAC,GAAG,CAAC,GAAG,EAAE,EAAE,CAAC;YACf,CAAC,GAAG,EAAE,GAAG,CAAC,CAAC;YACX,IAAI,GAAG,KAAK,CAAC;QACf,CAAC;QAED,uEAAuE;QACvE,IAAI,SAAS,GAAG,CAAC,IAAI,CAAC,KAAK,KAAK,EAAE,CAAC;YACjC,SAAS,GAAG,CAAC,CAAC;YACd,CAAC,CAAC,CAAC,EAAE,CAAC,EAAE,IAAI,CAAC,CAAC,CAAC,CAAC,CAAC;QACnB,CAAC;QAED,qBAAqB;QACrB,OAAO,IAAI,EAAE,CAAC;YACZ,IAAI,SAAS,GAAG,QAAQ,EAAE,CAAC;gBACzB,yDAAyD;gBAEzD,wBAAwB;gBACxB,CAAC,CAAC,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC,CAAC;gBAEZ,0CAA0C;gBAC1C,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,GAAG,EAAE,EAAE,CAAC;oBAC1B,KAAK,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,GAAG,GAAG,GAAG,EAAE,CAAC,CAAC,CAAC,CAAC;gBACpC,CAAC,CAAC,CAAC,GAAG,EAAE,GAAG,CAAC,EAAE,KAAK,EAAE,EAAE,CAAC,CAAC;gBAEzB,qDAAqD;gBACrD,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,GAAG,EAAE,EAAE,CAAC;oBAC1B,KAAK,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,GAAG,CAAC,GAAG,GAAG,EAAE,CAAC,CAAC,CAAC,GAAG,GAAG,GAAG,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC;gBACpD,CAAC,CAAC,CAAC,GAAG,EAAE,GAAG,CAAC,EAAE,KAAK,EAAE,EAAE,CAAC,CAAC;gBAEzB,8DAA8D;gBAC9D,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,GAAG,EAAE,EAAE,CAAC;oBAC1B,KAAK,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,GAAG,CAAC,GAAG,GAAG,EAAE,CAAC,CAAC,CAAC,GAAG,GAAG,GAAG,EAAE,CAAC,CAAC,CAAC,GAAG,GAAG,GAAG,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC;gBAClE,CAAC,CAAC,CAAC,GAAG,EAAE,GAAG,CAAC,EAAE,KAAK,EAAE,EAAE,CAAC,CAAC;gBAEzB,oEAAoE;gBACpE,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,GAAG,EAAE,EAAE,CAAC;oBAC1B,KAAK,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,GAAG,CAAC,GAAG,GAAG,EAAE,CAAC,CAAC,CAAC,GAAG,GAAG,GAAG,EAAE,CAAC,CAAC,CAAC,GAAG,GAAG,GAAG,EAAE,CAAC,CAAC,CAAC,GAAG,GAAG,GAAG,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC;gBAChF,CAAC,CAAC,CAAC,GAAG,CAAC,EAAE,KAAK,EAAE,EAAE,CAAC,CAAC;gBAEpB,gFAAgF;gBAChF,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,GAAG,EAAE,EAAE,CAAC;oBAC1B,KAAK,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,GAAG,CAAC,GAAG,GAAG,EAAE,CAAC,CAAC,CAAC,GAAG,GAAG,GAAG,EAAE,CAAC,CAAC,CAAC,GAAG,GAAG,GAAG,EAAE,CAAC,CAAC,CAAC,GAAG,GAAG,GAAG,EAAE,CAAC,CAAC,CAAC,GAAG,GAAG,GAAG,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC;gBAC9F,CAAC,CAAC,CAAC,GAAG,EAAE,GAAG,CAAC,EAAE,KAAK,EAAE,EAAE,CAAC,CAAC;gBAEzB,qBAAqB;gBACrB,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,GAAG,EAAE,EAAE,CAAC;oBAC1B,KAAK,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,GAAG,CAAC,EAAE,GAAG,EAAE,CAAC,CAAC,CAAC,GAAG,EAAE,GAAG,EAAE,CAAC,CAAC,CAAC,GAAG,EAAE,GAAG,EAAE,CAAC,CAAC,CAAC,GAAG,EAAE,GAAG,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC;gBAE5E,iEAAiE;gBACjE,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,GAAG,EAAE,EAAE,CAAC;oBAC1B,IAAI,CAAC,CAAC,CAAC,GAAG,CAAC,GAAG,CAAC,IAAI,GAAG,EAAE,CAAC,CAAC,CAAC,GAAG,IAAI,GAAG,EAAE,CAAC,CAAC,CAAC,GAAG,IAAI,GAAG,EAAE,CAAC,CAAC,CAAC,GAAG,IAAI,GAAG,EAAE,CAAC,CAAC,CAAC,GAAG,IAAI,GAAG,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC;gBAE3F,mBAAmB;gBACnB,MAAM,GAAG,CAAC,CAAC;gBACX,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,GAAG,EAAE,EAAE,CAAC;oBAC1B,MAAM,GAAG,GAAG,CAAC,MAAM,EAAE,GAAG,CAAC,IAAI,CAAC,CAAC,CAAC,GAAG,MAAM,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;gBACjD,MAAM,IAAI,SAAS,CAAC;gBAEpB,gCAAgC;gBAChC,IAAI,MAAM,GAAG,CAAC,EAAE,CAAC;oBACf,KAAK,GAAG,MAAM,GAAG,CAAC,GAAG,SAAA,MAAM,EAAE,MAAM,CAAA,CAAC;oBACpC,CAAC,GAAG,GAAG,CAAC,KAAK,EAAE,WAAW,GAAG,CAAC,CAAC,CAAC;oBAChC,IAAI,GAAG,CAAC,GAAG,CAAC,CAAC;oBACb,IAAI,IAAI,IAAI,CAAC;wBACX,MAAM,IAAI,KAAK,CAAC,SAAS,CAAC,SAAS,CAAC,CAAC;oBACvC,6CAA6C;oBAC7C,IAAI,SAAS,GAAG,CAAC,EAAE,CAAC;wBAClB,SAAS,GAAG,CAAC,CAAC;wBACd,CAAC,CAAC,CAAC,EAAE,CAAC,EAAE,IAAI,CAAC,CAAC,CAAC,CAAC,CAAC;oBACnB,CAAC;gBACH,CAAC;qBAAM,CAAC;oBACN,gBAAgB;oBAChB,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC;oBAEV,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,GAAG,EAAE,EAAE,CAAC;wBAC1B,CAAC,CAAC,CAAC,CAAC,GAAG,KAAK,CAAC,CAAC,CAAC,CAAC;oBAElB,sCAAsC;oBACtC,KAAK,IAAI,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,SAAS,EAAE,QAAQ,GAAG,CAAC,CAAC,EAAE,CAAC,GAAG,CAAC,EAAE,EAAE,CAAC;wBACxD,IAAI,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,IAAI,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC;oBAC3B,CAAC,CAAC,CAAC,EAAE,CAAC,EAAE,IAAI,CAAC,CAAC,CAAC,CAAC,CAAC;oBACjB,SAAS,EAAE,CAAC;oBACZ,KAAK,GAAG,CAAC,CAAC;oBAEV,wDAAwD;oBACxD,KAAK,GAAG,CAAC,CAAC;oBACV,IAAI,KAAK,GAAG,IAAI;wBACd,KAAK,GAAG,IAAI,CAAC;oBAEf,MAAM;gBACR,CAAC;YACH,CAAC;iBAAM,CAAC;gBACN,4CAA4C;gBAE5C,mCAAmC;gBACnC,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,GAAG,EAAE,EAAE,CAAC,EAAE,CAAC;oBAC7B,KAAK,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,GAAG,CAAC,GAAG,GAAG,IAAI,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,GAAG,GAAG,IAAI,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;wBACxD,GAAG,GAAG,IAAI,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,GAAG,GAAG,IAAI,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,GAAG,GAAG,IAAI,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;gBAC5D,CAAC;gBAED,gCAAgC;gBAChC,CAAC,CAAC,CAAC,GAAG,CAAC,EAAE,KAAK,EAAE,KAAK,CAAC,CAAC;gBAEvB,iCAAiC;gBACjC,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,GAAG,EAAE,EAAE,CAAC,EAAE,CAAC;oBAC7B,KAAK,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,GAAG,CAAC,MAAM,GAAG,KAAK,CAAC,CAAC,CAAC,GAAG,GAAG,GAAG,IAAI,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;wBACzD,GAAG,GAAG,IAAI,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,GAAG,GAAG,IAAI,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,GAAG,GAAG,GAAG,IAAI,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;gBAC5D,CAAC;gBAED,kCAAkC;gBAClC,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,GAAG,EAAE,EAAE,CAAC;oBAC1B,IAAI,CAAC,CAAC,CAAC,GAAG,KAAK,GAAG,CAAC,KAAK,CAAC,CAAC,CAAC,GAAG,KAAK,CAAC,CAAC,CAAC,CAAC,CAAC;gBAE1C,mBAAmB;gBACnB,MAAM,GAAG,CAAC,CAAC;gBACX,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,GAAG,EAAE,EAAE,CAAC;oBAC1B,MAAM,GAAG,GAAG,CAAC,MAAM,EAAE,GAAG,CAAC,IAAI,CAAC,CAAC,CAAC,GAAG,MAAM,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;gBACjD,MAAM,IAAI,SAAS,CAAC;gBAEpB,gCAAgC;gBAChC,IAAI,MAAM,GAAG,CAAC,EAAE,CAAC;oBACf,8CAA8C;oBAC9C,KAAK,GAAG,MAAM,GAAG,CAAC,GAAG,SAAA,MAAM,EAAE,MAAM,CAAA,CAAC;oBACpC,CAAC,GAAG,GAAG,CAAC,KAAK,EAAE,WAAW,GAAG,CAAC,CAAC,CAAC;oBAChC,IAAI,GAAG,CAAC,GAAG,CAAC,CAAC;oBACb,IAAI,IAAI,IAAI,CAAC;wBACX,MAAM,IAAI,KAAK,CAAC,SAAS,CAAC,SAAS,CAAC,CAAC;oBACvC,SAAS,GAAG,CAAC,CAAC;oBACd,CAAC,CAAC,CAAC,EAAE,CAAC,EAAE,IAAI,CAAC,CAAC,CAAC,CAAC,CAAC;gBACnB,CAAC;qBAAM,CAAC;oBACN,IAAI,MAAM,GAAG,SAAS;wBACpB,KAAK,GAAG,MAAM,GAAG,CAAC,GAAG,SAAA,MAAM,EAAE,MAAM,CAAA,CAAC;;wBAEpC,KAAK,GAAG,SAAS,GAAG,CAAC,CAAC;oBAExB,IAAI,KAAK,GAAG,IAAI;wBACd,KAAK,GAAG,IAAI,CAAC;oBAEf,CAAC,GAAG,CAAC,GAAG,CAAC,CAAC;oBAEV,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,GAAG,EAAE,EAAE,CAAC;wBAC1B,CAAC,CAAC,CAAC,CAAC,GAAG,KAAK,CAAC,CAAC,CAAC,CAAC;oBAElB,qEAAqE;oBACrE,MAAM,QAAQ,GAAG,IAAI,CAAC,QAAQ,GAAG,CAAC,CAAC,CAAC;oBACpC,KAAK,IAAI,CAAC,GAAG,QAAQ,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,EAAE,EAAE,CAAC;wBACnC,IAAI,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC;oBACxB,IAAI,CAAC,CAAC,CAAC,GAAG,QAAQ,CAAC;oBACnB,CAAC,CAAC,CAAC,EAAE,CAAC,EAAE,IAAI,CAAC,CAAC,CAAC,CAAC,CAAC;oBACjB,KAAK,GAAG,CAAC,CAAC;oBAEV,MAAM;gBACR,CAAC;YACH,CAAC;QACH,CAAC,CAAC,eAAe;QAEjB,oEAAoE;QACpE,OAAO,aAAa,GAAG,CAAC,EAAE,CAAC;YACzB,MAAM,KAAK,GAAG,CAAC,CAAC,GAAG,aAAa,CAAC,GAAG,CAAC,CAAC,GAAG,KAAK,CAAC,CAAC;YAChD,MAAM,MAAM,GAAG,GAAG,GAAG,KAAK,CAAC;YAE3B,IAAI,CAAC,KAAK,CAAC,GAAG,aAAa,CAAC;YAE5B,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,GAAG,EAAE,EAAE,CAAC;gBAC1B,KAAK,CAAC,CAAC,CAAC,CAAC,KAAK,CAAC,GAAG,MAAM,GAAG,CAAC,CAAC,CAAC,CAAC,GAAG,KAAK,GAAG,KAAK,CAAC,CAAC,CAAC,CAAC;YAErD,aAAa,IAAI,UAAU,CAAC;YAC5B,EAAE,KAAK,CAAC;QACV,CAAC;QAED,CAAC,GAAG,KAAK,CAAC;QACV,KAAK,GAAG,CAAC,CAAC;QAEV,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,GAAG,EAAE,EAAE,CAAC;YAC1B,KAAK,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC;IACpB,CAAC,CAAC,eAAe;IAEjB,iCAAiC;IACjC,IAAI,QAAQ;QACV,QAAQ,CAAC,uBAAuB,EAAE,CAAC;IAErC,8BAA8B;IAC9B,IAAI,CAAC,QAAQ,GAAG,CAAC,CAAC,GAAG,EAAE,CAAC;IAExB,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,GAAG,EAAE,EAAE,CAAC;QAC1B,KAAK,CAAC,CAAC,CAAC,CAAC,QAAQ,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC;IAEhC,oBAAoB;IACpB,MAAM,QAAQ,GAAG,KAAK,CAAe,GAAG,GAAG,CAAC,CAAC,CAAC;IAE9C,uBAAuB;IACvB,QAAQ,CAAC,CAAC,CAAC,GAAG,IAAI,CAAC;IAEnB,YAAY;IACZ,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,GAAG,EAAE,EAAE,CAAC;QAC1B,QAAQ,CAAC,CAAC,GAAG,CAAC,CAAC,GAAG,KAAK,CAAC,CAAC,CAAC,CAAC;IAE7B,OAAO,QAAQ,CAAC;AAClB,CAAC,CAAC,MAAM"}