diff-grok 1.2.0 → 1.3.0

package/src/solver-tools/cvode/cvode_nls.ts

@@ -1,259 +0,0 @@
-/**
- * TypeScript port of SUNDIALS CVODE v7.5.0 (https://github.com/LLNL/sundials)
- *
- * Copyright (c) 2025, Lawrence Livermore National Security, University of
- *   Maryland Baltimore County, and the SUNDIALS contributors.
- * Copyright (c) 2013-2025, Lawrence Livermore National Security and
- *   Southern Methodist University.
- * Copyright (c) 2002-2013, Lawrence Livermore National Security.
- * Copyright (c) 2026 Datagrok.
- *
- * Licensed under the BSD-3-Clause License. See THIRD_PARTY_LICENSES
- * for the full license text and provenance chain of the original code.
- */
-
-import {CvodeMem, wrmsNorm, vLinearSum, vScale, vConst,
-  CV_SUCCESS, CV_BDF, CV_ADAMS, CV_RHSFUNC_FAIL, CV_LSETUP_FAIL, CV_LSOLVE_FAIL,
-  RHSFUNC_RECVR, SUN_NLS_CONV_RECVR, SUN_NLS_CONTINUE,
-  CV_NO_FAILURES, CV_FAIL_BAD_J, CV_FAIL_OTHER,
-  FIRST_CALL, PREV_CONV_FAIL, PREV_ERR_FAIL, DO_ERROR_TEST,
-  NLS_MAXCOR, CRDOWN, RDIV, ONEPSM} from './common';
-
-/**
- * cvNls - Solve the nonlinear system for one step.
- *
- * For BDF: Newton iteration solving G(ycor) = 0 where
- *   G(ycor) = rl1*zn[1] + ycor - gamma*f(tn, zn[0]+ycor)
- * The Newton update is: (I - gamma*J) * delta = -G(ycor)
- *
- * For Adams: Fixed-point iteration
- *   ycor_new = h * rl1 * f(tn, zn[0]+ycor) - rl1*zn[1]
- */
-export function cvNls(mem: CvodeMem, nflag: number): number {
-  const N = mem.cv_N;
-  let callSetup: boolean;
-
-  // Decide whether to call linear solver setup
-  if (mem.cv_lsetup) {
-    mem.convfail = (nflag === FIRST_CALL || nflag === PREV_ERR_FAIL) ?
-      CV_NO_FAILURES :
-      CV_FAIL_OTHER;
-
-    callSetup = (nflag === PREV_CONV_FAIL) || (nflag === PREV_ERR_FAIL) ||
-                (mem.cv_nst === 0) || (mem.first_step_after_resize) ||
-                (mem.cv_nst >= mem.cv_nstlp + mem.cv_msbp) ||
-                (Math.abs(mem.cv_gamrat - 1.0) > mem.cv_dgmax_lsetup);
-  } else {
-    mem.cv_crate = 1.0;
-    callSetup = false;
-  }
-
-  // Initial guess for correction: acor = 0
-  vConst(0.0, mem.cv_acor, N);
-
-  // Set acnrmcur to false
-  mem.cv_acnrmcur = false;
-
-  // Dispatch to Newton or fixed-point
-  let flag: number;
-  if (mem.cv_lmm === CV_BDF)
-    flag = cvNewtonIteration(mem, callSetup);
-  else
-    flag = cvFixedPointIteration(mem);
-
-
-  if (flag !== CV_SUCCESS) return flag;
-
-  // Solve successful - update y = zn[0] + acor
-  vLinearSum(1.0, mem.cv_zn[0], 1.0, mem.cv_acor, mem.cv_y, N);
-
-  // Compute acnrm if not already done
-  if (!mem.cv_acnrmcur)
-    mem.cv_acnrm = wrmsNorm(N, mem.cv_acor, mem.cv_ewt);
-
-
-  // Update Jacobian status
-  mem.cv_jcur = false;
-
-  return CV_SUCCESS;
-}
-
-/**
- * Newton iteration for BDF method.
- * Solves: (I - gamma*J) * delta = -G(ycor)
- * where G(ycor) = rl1*zn[1] + ycor - gamma*f(tn, zn[0]+ycor)
- */
-function cvNewtonIteration(mem: CvodeMem, callSetup: boolean): number {
-  const N = mem.cv_N;
-  let retval: number;
-
-  // Call linear solver setup if needed
-  if (callSetup && mem.cv_lsetup) {
-    // y = zn[0] + acor (current approximation)
-    vLinearSum(1.0, mem.cv_zn[0], 1.0, mem.cv_acor, mem.cv_y, N);
-
-    // Evaluate f at current approximation
-    retval = mem.cv_f!(mem.cv_tn, mem.cv_y, mem.cv_ftemp, mem.cv_user_data);
-    mem.cv_nfe++;
-    if (retval < 0) return CV_RHSFUNC_FAIL;
-    if (retval > 0) return RHSFUNC_RECVR;
-
-    // Call lsetup
-    const jcurPtr = {value: false};
-    retval = mem.cv_lsetup!(mem, mem.convfail, mem.cv_y, mem.cv_ftemp, jcurPtr,
-      mem.cv_vtemp1, mem.cv_vtemp2, mem.cv_vtemp3);
-    mem.cv_nsetups++;
-
-    mem.cv_jcur = jcurPtr.value;
-    mem.cv_gamrat = 1.0;
-    mem.cv_gammap = mem.cv_gamma;
-    mem.cv_crate = 1.0;
-    mem.cv_nstlp = mem.cv_nst;
-
-    if (retval < 0) return CV_LSETUP_FAIL;
-    if (retval > 0) return SUN_NLS_CONV_RECVR;
-  }
-
-  // Newton iteration loop
-  for (let m = 0; m < NLS_MAXCOR; m++) {
-    mem.cv_nni++;
-
-    // y = zn[0] + acor
-    vLinearSum(1.0, mem.cv_zn[0], 1.0, mem.cv_acor, mem.cv_y, N);
-
-    // Evaluate f(tn, y)
-    retval = mem.cv_f!(mem.cv_tn, mem.cv_y, mem.cv_ftemp, mem.cv_user_data);
-    mem.cv_nfe++;
-    if (retval < 0) return CV_RHSFUNC_FAIL;
-    if (retval > 0) return RHSFUNC_RECVR;
-
-    // Compute residual: b = rl1*zn[1] + acor - gamma*ftemp
-    // Store in tempv as the RHS for the linear solve
-    const b = mem.cv_tempv;
-    for (let i = 0; i < N; i++)
-      b[i] = mem.cv_rl1 * mem.cv_zn[1][i] + mem.cv_acor[i] - mem.cv_gamma * mem.cv_ftemp[i];
-
-
-    // Solve: (I - gamma*J) * delta = -b, so pass b (which is G, so we negate)
-    // Actually the convention is lsolve solves (I-gamma*J)*x = b
-    // and b here is -G so we negate
-    for (let i = 0; i < N; i++) b[i] = -b[i];
-
-    retval = mem.cv_lsolve!(mem, b, mem.cv_ewt, mem.cv_y, mem.cv_ftemp);
-    if (retval < 0) return CV_LSOLVE_FAIL;
-    if (retval > 0) return SUN_NLS_CONV_RECVR;
-
-    // Update acor += delta (delta is in b after solve)
-    for (let i = 0; i < N; i++) mem.cv_acor[i] += b[i];
-
-    // Convergence test
-    const del = wrmsNorm(N, b, mem.cv_ewt);
-    const convTestResult = cvNlsConvTest(mem, del, m, mem.cv_tq[4]);
-
-    if (convTestResult === CV_SUCCESS) return CV_SUCCESS;
-    if (convTestResult === SUN_NLS_CONV_RECVR) {
-      mem.cv_nnf++;
-      return SUN_NLS_CONV_RECVR;
-    }
-    // SUN_NLS_CONTINUE: continue iterating
-  }
-
-  // Max iterations reached without convergence
-  mem.cv_nnf++;
-  return SUN_NLS_CONV_RECVR;
-}
-
-/**
- * Fixed-point iteration for Adams method.
- * ycor = rl1 * (h * f(tn, zn[0]+ycor) - zn[1])
- */
-function cvFixedPointIteration(mem: CvodeMem): number {
-  const N = mem.cv_N;
-  let retval: number;
-
-  for (let m = 0; m < NLS_MAXCOR; m++) {
-    mem.cv_nni++;
-
-    // y = zn[0] + acor
-    vLinearSum(1.0, mem.cv_zn[0], 1.0, mem.cv_acor, mem.cv_y, N);
-
-    // Evaluate f(tn, y) -> store in tempv
-    retval = mem.cv_f!(mem.cv_tn, mem.cv_y, mem.cv_tempv, mem.cv_user_data);
-    mem.cv_nfe++;
-    if (retval < 0) return CV_RHSFUNC_FAIL;
-    if (retval > 0) return RHSFUNC_RECVR;
-
-    // Compute new ycor: rl1 * (h * f - zn[1])
-    // First: tempv = h * tempv - zn[1]
-    for (let i = 0; i < N; i++)
-      mem.cv_tempv[i] = mem.cv_h * mem.cv_tempv[i] - mem.cv_zn[1][i];
-
-    // Then scale by rl1
-    const newAcor = mem.cv_vtemp1;
-    vScale(mem.cv_rl1, mem.cv_tempv, newAcor, N);
-
-    // Compute delta = newAcor - acor
-    const delta = mem.cv_tempv;
-    for (let i = 0; i < N; i++)
-      delta[i] = newAcor[i] - mem.cv_acor[i];
-
-
-    // Update acor
-    for (let i = 0; i < N; i++) mem.cv_acor[i] = newAcor[i];
-
-    // Convergence test
-    const del = wrmsNorm(N, delta, mem.cv_ewt);
-    const convTestResult = cvNlsConvTest(mem, del, m, mem.cv_tq[4]);
-
-    if (convTestResult === CV_SUCCESS) return CV_SUCCESS;
-    if (convTestResult === SUN_NLS_CONV_RECVR) {
-      mem.cv_nnf++;
-      return SUN_NLS_CONV_RECVR;
-    }
-    // SUN_NLS_CONTINUE: continue iterating
-  }
-
-  // Max iterations reached
-  mem.cv_nnf++;
-  return SUN_NLS_CONV_RECVR;
-}
-
-/**
- * cvNlsConvTest - Test for convergence of the nonlinear iteration.
- *
- * del: WRMS norm of the current update
- * m: iteration count (0-based)
- * tol: tolerance (tq[4])
- *
- * Returns: CV_SUCCESS if converged, SUN_NLS_CONV_RECVR if diverging,
- *          SUN_NLS_CONTINUE if not yet converged
- */
-function cvNlsConvTest(mem: CvodeMem, del: number, m: number, tol: number): number {
-  // Update convergence rate estimate
-  if (m > 0)
-    mem.cv_crate = Math.max(CRDOWN * mem.cv_crate, del / mem.cv_delp);
-
-
-  // Test for convergence
-  const dcon = del * Math.min(1.0, mem.cv_crate) / tol;
-
-  if (dcon <= 1.0) {
-    // Converged
-    if (m === 0)
-      mem.cv_acnrm = del;
-    else
-      mem.cv_acnrm = wrmsNorm(mem.cv_N, mem.cv_acor, mem.cv_ewt);
-
-    mem.cv_acnrmcur = true;
-    return CV_SUCCESS;
-  }
-
-  // Check for divergence
-  if (m >= 1 && del > RDIV * mem.cv_delp)
-    return SUN_NLS_CONV_RECVR;
-
-
-  // Save del and continue
-  mem.cv_delp = del;
-  return SUN_NLS_CONTINUE;
-}

package/src/solver-tools/cvode/cvode_root.js.map

@@ -1 +0,0 @@
-{"version":3,"file":"cvode_root.js","sourceRoot":"","sources":["cvode_root.ts"],"names":[],"mappings":"AAAA;;;;;;;;;;;;GAYG;AAEH,6BAA6B;AAC7B,6EAA6E;AAE7E,OAAO,EAGL,WAAW,EACX,UAAU,EACV,MAAM,EACN,UAAU,EACV,cAAc,EACd,OAAO,EACP,OAAO,EACP,MAAM,EAEN,SAAS,EACT,WAAW,GACZ,MAAM,UAAU,CAAC;AAElB,+EAA+E;AAC/E,8CAA8C;AAC9C,+EAA+E;AAE/E,MAAM,IAAI,GAAG,GAAG,CAAC;AACjB,MAAM,IAAI,GAAG,GAAG,CAAC;AACjB,MAAM,GAAG,GAAG,GAAG,CAAC;AAChB,MAAM,GAAG,GAAG,GAAG,CAAC;AAChB,MAAM,IAAI,GAAG,GAAG,CAAC;AACjB,MAAM,OAAO,GAAG,KAAK,CAAC;AACtB,MAAM,GAAG,GAAG,GAAG,CAAC;AAEhB,+EAA+E;AAC/E,qBAAqB;AACrB,+EAA+E;AAE/E;;;;;;;;;;;;GAYG;AACH,MAAM,UAAU,aAAa,CAAC,GAAa,EAAE,KAAa,EAAE,CAAc;IACxE,qCAAqC;IACrC,IAAI,KAAK,IAAI,CAAC,EAAE,CAAC;QACf,GAAG,CAAC,QAAQ,GAAG,CAAC,CAAC;QACjB,GAAG,CAAC,OAAO,GAAG,IAAI,CAAC;QACnB,OAAO,UAAU,CAAC;IACpB,CAAC;IAED,GAAG,CAAC,QAAQ,GAAG,KAAK,CAAC;IACrB,GAAG,CAAC,OAAO,GAAG,CAAC,CAAC;IAEhB,sCAAsC;IACtC,GAAG,CAAC,MAAM,GAAG,IAAI,YAAY,CAAC,KAAK,CAAC,CAAC;IACrC,GAAG,CAAC,MAAM,GAAG,IAAI,YAAY,CAAC,KAAK,CAAC,CAAC;IACrC,GAAG,CAAC,QAAQ,GAAG,IAAI,YAAY,CAAC,KAAK,CAAC,CAAC;IAEvC,mDAAmD;IACnD,GAAG,CAAC,SAAS,GAAG,IAAI,UAAU,CAAC,KAAK,CAAC,CAAC;IACtC,GAAG,CAAC,UAAU,GAAG,IAAI,UAAU,CAAC,KAAK,CAAC,CAAC;IAEvC,uCAAuC;IACvC,GAAG,CAAC,UAAU,GAAG,IAAI,UAAU,CAAC,KAAK,CAAC,CAAC;IAEvC,8EAA8E;IAC9E,GAAG,CAAC,UAAU,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC;IACvB,GAAG,CAAC,UAAU,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC;IAEvB,OAAO,UAAU,CAAC;AACpB,CAAC;AAED;;;;;;;;;;GAUG;AACH,MAAM,UAAU,SAAS,CAAC,GAAa;IACrC,MAAM,KAAK,GAAG,GAAG,CAAC,QAAQ,CAAC;IAC3B,MAAM,CAAC,GAAG,GAAG,CAAC,IAAI,CAAC;IAEnB,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,KAAK,EAAE,CAAC,EAAE;QAC5B,GAAG,CAAC,SAAS,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC;IAGvB,GAAG,CAAC,MAAM,GAAG,GAAG,CAAC,KAAK,CAAC;IACvB,GAAG,CAAC,OAAO,GAAG,CAAC,IAAI,CAAC,GAAG,CAAC,GAAG,CAAC,KAAK,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,GAAG,CAAC,IAAI,CAAC,CAAC,GAAG,MAAM,GAAG,OAAO,CAAC;IAE5E,oDAAoD;IACpD,MAAM,MAAM,GAAG,GAAG,CAAC,OAAQ,CAAC,GAAG,CAAC,MAAM,EAAE,GAAG,CAAC,KAAK,CAAC,CAAC,CAAC,EAAE,GAAG,CAAC,MAAM,EAAE,GAAG,CAAC,YAAY,CAAC,CAAC;IACpF,GAAG,CAAC,MAAM,GAAG,CAAC,CAAC;IACf,IAAI,MAAM,KAAK,CAAC;QACd,OAAO,cAAc,CAAC;IAGxB,IAAI,KAAK,GAAG,KAAK,CAAC;IAClB,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,KAAK,EAAE,CAAC,EAAE,EAAE,CAAC;QAC/B,IAAI,IAAI,CAAC,GAAG,CAAC,GAAG,CAAC,MAAM,CAAC,CAAC,CAAC,CAAC,KAAK,IAAI,EAAE,CAAC;YACrC,KAAK,GAAG,IAAI,CAAC;YACb,GAAG,CAAC,UAAU,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,WAAW;QACpC,CAAC;IACH,CAAC;IACD,IAAI,CAAC,KAAK;QACR,OAAO,UAAU,CAAC;IAGpB,4DAA4D;IAC5D,MAAM,MAAM,GAAG,IAAI,CAAC,GAAG,CAAC,GAAG,CAAC,OAAO,GAAG,IAAI,CAAC,GAAG,CAAC,GAAG,CAAC,IAAI,CAAC,EAAE,GAAG,CAAC,CAAC;IAC/D,MAAM,MAAM,GAAG,MAAM,GAAG,GAAG,CAAC,IAAI,CAAC;IACjC,MAAM,KAAK,GAAG,GAAG,CAAC,MAAM,GAAG,MAAM,CAAC;IAElC,6BAA6B;IAC7B,UAAU,CAAC,GAAG,EAAE,GAAG,CAAC,KAAK,CAAC,CAAC,CAAC,EAAE,MAAM,EAAE,GAAG,CAAC,KAAK,CAAC,CAAC,CAAC,EAAE,GAAG,CAAC,IAAI,EAAE,CAAC,CAAC,CAAC;IAEjE,MAAM,OAAO,GAAG,GAAG,CAAC,OAAQ,CAAC,KAAK,EAAE,GAAG,CAAC,IAAI,EAAE,GAAG,CAAC,MAAM,EAAE,GAAG,CAAC,YAAY,CAAC,CAAC;IAC5E,GAAG,CAAC,MAAM,EAAE,CAAC;IACb,IAAI,OAAO,KAAK,CAAC;QACf,OAAO,cAAc,CAAC;IAGxB,sDAAsD;IACtD,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,KAAK,EAAE,CAAC,EAAE,EAAE,CAAC;QAC/B,IAAI,GAAG,CAAC,UAAU,CAAC,CAAC,CAAC,KAAK,CAAC,IAAI,IAAI,CAAC,GAAG,CAAC,GAAG,CAAC,MAAM,CAAC,CAAC,CAAC,CAAC,KAAK,IAAI,EAAE,CAAC;YAChE,GAAG,CAAC,UAAU,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,UAAU;YACjC,GAAG,CAAC,MAAM,CAAC,CAAC,CAAC,GAAG,GAAG,CAAC,MAAM,CAAC,CAAC,CAAC,CAAC;QAChC,CAAC;IACH,CAAC;IAED,OAAO,UAAU,CAAC;AACpB,CAAC;AAED;;;;;;;;;;GAUG;AACH,MAAM,UAAU,SAAS,CAAC,GAAa;IACrC,MAAM,KAAK,GAAG,GAAG,CAAC,QAAQ,CAAC;IAC3B,MAAM,CAAC,GAAG,GAAG,CAAC,IAAI,CAAC;IAEnB,IAAI,GAAG,CAAC,QAAQ,KAAK,CAAC;QACpB,OAAO,UAAU,CAAC;IAGpB,iCAAiC;IACjC,WAAW,CAAC,GAAG,EAAE,GAAG,CAAC,MAAM,EAAE,CAAC,EAAE,GAAG,CAAC,IAAI,CAAC,CAAC;IAE1C,qBAAqB;IACrB,IAAI,MAAM,GAAG,GAAG,CAAC,OAAQ,CAAC,GAAG,CAAC,MAAM,EAAE,GAAG,CAAC,IAAI,EAAE,GAAG,CAAC,MAAM,EAAE,GAAG,CAAC,YAAY,CAAC,CAAC;IAC9E,GAAG,CAAC,MAAM,EAAE,CAAC;IACb,IAAI,MAAM,KAAK,CAAC;QACd,OAAO,cAAc,CAAC;IAGxB,IAAI,KAAK,GAAG,KAAK,CAAC;IAClB,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,KAAK,EAAE,CAAC,EAAE;QAC5B,GAAG,CAAC,SAAS,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC;IAEvB,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,KAAK,EAAE,CAAC,EAAE,EAAE,CAAC;QAC/B,IAAI,GAAG,CAAC,UAAU,CAAC,CAAC,CAAC,KAAK,CAAC;YAAE,SAAS;QACtC,IAAI,IAAI,CAAC,GAAG,CAAC,GAAG,CAAC,MAAM,CAAC,CAAC,CAAC,CAAC,KAAK,IAAI,EAAE,CAAC;YACrC,KAAK,GAAG,IAAI,CAAC;YACb,GAAG,CAAC,SAAS,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC;QACvB,CAAC;IACH,CAAC;IACD,IAAI,CAAC,KAAK;QACR,OAAO,UAAU,CAAC;IAGpB,4DAA4D;IAC5D,GAAG,CAAC,OAAO,GAAG,CAAC,IAAI,CAAC,GAAG,CAAC,GAAG,CAAC,KAAK,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,GAAG,CAAC,IAAI,CAAC,CAAC,GAAG,MAAM,GAAG,OAAO,CAAC;IAC5E,MAAM,MAAM,GAAG,CAAC,GAAG,CAAC,IAAI,GAAG,IAAI,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,OAAO,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,OAAO,CAAC;IAC9D,MAAM,KAAK,GAAG,GAAG,CAAC,MAAM,GAAG,MAAM,CAAC;IAElC,IAAI,CAAC,KAAK,GAAG,GAAG,CAAC,KAAK,CAAC,GAAG,GAAG,CAAC,IAAI,IAAI,IAAI,EAAE,CAAC;QAC3C,MAAM,MAAM,GAAG,MAAM,GAAG,GAAG,CAAC,IAAI,CAAC;QACjC,yBAAyB;QACzB,UAAU,CAAC,GAAG,EAAE,GAAG,CAAC,IAAI,EAAE,MAAM,EAAE,GAAG,CAAC,KAAK,CAAC,CAAC,CAAC,EAAE,GAAG,CAAC,IAAI,EAAE,CAAC,CAAC,CAAC;IAC/D,CAAC;;QACC,WAAW,CAAC,GAAG,EAAE,KAAK,EAAE,CAAC,EAAE,GAAG,CAAC,IAAI,CAAC,CAAC;IAGvC,MAAM,GAAG,GAAG,CAAC,OAAQ,CAAC,KAAK,EAAE,GAAG,CAAC,IAAI,EAAE,GAAG,CAAC,MAAM,EAAE,GAAG,CAAC,YAAY,CAAC,CAAC;IACrE,GAAG,CAAC,MAAM,EAAE,CAAC;IACb,IAAI,MAAM,KAAK,CAAC;QACd,OAAO,cAAc,CAAC;IAGxB,gEAAgE;IAChE,6CAA6C;IAC7C,KAAK,GAAG,KAAK,CAAC;IACd,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,KAAK,EAAE,CAAC,EAAE,EAAE,CAAC;QAC/B,IAAI,GAAG,CAAC,UAAU,CAAC,CAAC,CAAC,KAAK,CAAC;YAAE,SAAS;QACtC,IAAI,IAAI,CAAC,GAAG,CAAC,GAAG,CAAC,MAAM,CAAC,CAAC,CAAC,CAAC,KAAK,IAAI,EAAE,CAAC;YACrC,IAAI,GAAG,CAAC,SAAS,CAAC,CAAC,CAAC,KAAK,CAAC;gBACxB,OAAO,OAAO,CAAC;YAEjB,KAAK,GAAG,IAAI,CAAC;YACb,GAAG,CAAC,SAAS,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC;QACvB,CAAC;aAAM,CAAC;YACN,IAAI,GAAG,CAAC,SAAS,CAAC,CAAC,CAAC,KAAK,CAAC;gBACxB,GAAG,CAAC,MAAM,CAAC,CAAC,CAAC,GAAG,GAAG,CAAC,MAAM,CAAC,CAAC,CAAC,CAAC;QAClC,CAAC;IACH,CAAC;IACD,IAAI,KAAK;QACP,OAAO,OAAO,CAAC;IAGjB,OAAO,UAAU,CAAC;AACpB,CAAC;AAED;;;;;;;;;;GAUG;AACH,MAAM,UAAU,SAAS,CAAC,GAAa;IACrC,MAAM,KAAK,GAAG,GAAG,CAAC,QAAQ,CAAC;IAC3B,MAAM,CAAC,GAAG,GAAG,CAAC,IAAI,CAAC;IAEnB,8DAA8D;IAC9D,IAAI,GAAG,CAAC,QAAQ,KAAK,WAAW,EAAE,CAAC;QACjC,GAAG,CAAC,MAAM,GAAG,GAAG,CAAC,KAAK,CAAC;QACvB,MAAM,CAAC,GAAG,EAAE,GAAG,CAAC,KAAK,CAAC,CAAC,CAAC,EAAE,GAAG,CAAC,IAAI,EAAE,CAAC,CAAC,CAAC;IACzC,CAAC;IACD,IAAI,GAAG,CAAC,QAAQ,KAAK,SAAS,EAAE,CAAC;QAC/B,IAAI,CAAC,GAAG,CAAC,QAAQ,GAAG,GAAG,CAAC,KAAK,CAAC,GAAG,GAAG,CAAC,IAAI,IAAI,IAAI,EAAE,CAAC;YAClD,GAAG,CAAC,MAAM,GAAG,GAAG,CAAC,KAAK,CAAC;YACvB,MAAM,CAAC,GAAG,EAAE,GAAG,CAAC,KAAK,CAAC,CAAC,CAAC,EAAE,GAAG,CAAC,IAAI,EAAE,CAAC,CAAC,CAAC;QACzC,CAAC;aAAM,CAAC;YACN,GAAG,CAAC,MAAM,GAAG,GAAG,CAAC,QAAQ,CAAC;YAC1B,WAAW,CAAC,GAAG,EAAE,GAAG,CAAC,MAAM,EAAE,CAAC,EAAE,GAAG,CAAC,IAAI,CAAC,CAAC;QAC5C,CAAC;IACH,CAAC;IAED,yCAAyC;IACzC,MAAM,MAAM,GAAG,GAAG,CAAC,OAAQ,CAAC,GAAG,CAAC,MAAM,EAAE,GAAG,CAAC,IAAI,EAAE,GAAG,CAAC,MAAM,EAAE,GAAG,CAAC,YAAY,CAAC,CAAC;IAChF,GAAG,CAAC,MAAM,EAAE,CAAC;IACb,IAAI,MAAM,KAAK,CAAC;QACd,OAAO,cAAc,CAAC;IAGxB,GAAG,CAAC,OAAO,GAAG,CAAC,IAAI,CAAC,GAAG,CAAC,GAAG,CAAC,KAAK,CAAC,GAAG,IAAI,CAAC,GAAG,CAAC,GAAG,CAAC,IAAI,CAAC,CAAC,GAAG,MAAM,GAAG,OAAO,CAAC;IAC5E,MAAM,GAAG,GAAG,UAAU,CAAC,GAAG,CAAC,CAAC;IAC5B,IAAI,GAAG,KAAK,cAAc;QACxB,OAAO,cAAc,CAAC;IAGxB,qDAAqD;IACrD,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,KAAK,EAAE,CAAC,EAAE,EAAE,CAAC;QAC/B,IAAI,GAAG,CAAC,UAAU,CAAC,CAAC,CAAC,KAAK,CAAC,IAAI,GAAG,CAAC,QAAQ,CAAC,CAAC,CAAC,KAAK,IAAI;YACrD,GAAG,CAAC,UAAU,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,UAAU;IACrC,CAAC;IAED,GAAG,CAAC,MAAM,GAAG,GAAG,CAAC,QAAQ,CAAC;IAC1B,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,KAAK,EAAE,CAAC,EAAE;QAC5B,GAAG,CAAC,MAAM,CAAC,CAAC,CAAC,GAAG,GAAG,CAAC,QAAQ,CAAC,CAAC,CAAC,CAAC;IAGlC,sCAAsC;IACtC,IAAI,GAAG,KAAK,UAAU;QACpB,OAAO,UAAU,CAAC;IAGpB,8DAA8D;IAC9D,WAAW,CAAC,GAAG,EAAE,GAAG,CAAC,QAAQ,EAAE,CAAC,EAAE,GAAG,CAAC,IAAI,CAAC,CAAC;IAC5C,OAAO,OAAO,CAAC;AACjB,CAAC;AAED,+EAA+E;AAC/E,oCAAoC;AACpC,+EAA+E;AAE/E;;;;;;;;;;;;GAYG;AACH,SAAS,UAAU,CAAC,GAAa;IAC/B,MAAM,KAAK,GAAG,GAAG,CAAC,QAAQ,CAAC;IAE3B,IAAI,IAAI,GAAG,CAAC,CAAC;IACb,IAAI,IAAY,CAAC;IACjB,IAAI,IAAY,CAAC;IACjB,IAAI,KAAa,CAAC;IAClB,IAAI,OAAe,CAAC;IACpB,IAAI,OAAe,CAAC;IACpB,IAAI,OAAe,CAAC;IAEpB,8DAA8D;IAC9D,OAAO,GAAG,IAAI,CAAC;IACf,IAAI,KAAK,GAAG,KAAK,CAAC;IAClB,IAAI,MAAM,GAAG,KAAK,CAAC;IAEnB,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,KAAK,EAAE,CAAC,EAAE,EAAE,CAAC;QAC/B,IAAI,GAAG,CAAC,UAAU,CAAC,CAAC,CAAC,KAAK,CAAC;YAAE,SAAS;QACtC,IAAI,IAAI,CAAC,GAAG,CAAC,GAAG,CAAC,MAAM,CAAC,CAAC,CAAC,CAAC,KAAK,IAAI,EAAE,CAAC;YACrC,IAAI,GAAG,CAAC,UAAU,CAAC,CAAC,CAAC,GAAG,GAAG,CAAC,MAAM,CAAC,CAAC,CAAC,IAAI,IAAI;gBAC3C,KAAK,GAAG,IAAI,CAAC;QACjB,CAAC;aAAM,CAAC;YACN,IAAI,CAAC,GAAG,CAAC,MAAM,CAAC,CAAC,CAAC,GAAG,GAAG,CAAC,MAAM,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC;gBACnC,CAAC,GAAG,CAAC,UAAU,CAAC,CAAC,CAAC,GAAG,GAAG,CAAC,MAAM,CAAC,CAAC,CAAC,IAAI,IAAI,CAAC,EAAE,CAAC;gBAChD,KAAK,GAAG,IAAI,CAAC,GAAG,CAAC,GAAG,CAAC,MAAM,CAAC,CAAC,CAAC,GAAG,CAAC,GAAG,CAAC,MAAM,CAAC,CAAC,CAAC,GAAG,GAAG,CAAC,MAAM,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;gBAClE,IAAI,KAAK,GAAG,OAAO,EAAE,CAAC;oBACpB,MAAM,GAAG,IAAI,CAAC;oBACd,OAAO,GAAG,KAAK,CAAC;oBAChB,IAAI,GAAG,CAAC,CAAC;gBACX,CAAC;YACH,CAAC;QACH,CAAC;IACH,CAAC;IAED,sDAAsD;IACtD,4EAA4E;IAC5E,IAAI,CAAC,MAAM,EAAE,CAAC;QACZ,GAAG,CAAC,QAAQ,GAAG,GAAG,CAAC,MAAM,CAAC;QAC1B,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,KAAK,EAAE,CAAC,EAAE;YAC5B,GAAG,CAAC,QAAQ,CAAC,CAAC,CAAC,GAAG,GAAG,CAAC,MAAM,CAAC,CAAC,CAAC,CAAC;QAElC,IAAI,CAAC,KAAK;YACR,OAAO,UAAU,CAAC;QAEpB,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,KAAK,EAAE,CAAC,EAAE,EAAE,CAAC;YAC/B,GAAG,CAAC,SAAS,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC;YACrB,IAAI,GAAG,CAAC,UAAU,CAAC,CAAC,CAAC,KAAK,CAAC;gBAAE,SAAS;YACtC,IAAI,CAAC,IAAI,CAAC,GAAG,CAAC,GAAG,CAAC,MAAM,CAAC,CAAC,CAAC,CAAC,KAAK,IAAI,CAAC;gBAClC,CAAC,GAAG,CAAC,UAAU,CAAC,CAAC,CAAC,GAAG,GAAG,CAAC,MAAM,CAAC,CAAC,CAAC,IAAI,IAAI,CAAC;gBAC7C,GAAG,CAAC,SAAS,CAAC,CAAC,CAAC,GAAG,GAAG,CAAC,MAAM,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;QAClD,CAAC;QACD,OAAO,OAAO,CAAC;IACjB,CAAC;IAED,kBAAkB;IAClB,IAAI,GAAG,GAAG,CAAC;IAEX,wDAAwD;IACxD,IAAI,IAAI,GAAG,CAAC,CAAC;IACb,IAAI,QAAQ,GAAG,CAAC,CAAC,CAAC;IAElB,SAAS,CAAC;QACR,+DAA+D;QAC/D,IAAI,IAAI,CAAC,GAAG,CAAC,GAAG,CAAC,MAAM,GAAG,GAAG,CAAC,MAAM,CAAC,IAAI,GAAG,CAAC,OAAO;YAClD,MAAM;QAGR,mBAAmB;QACnB,gEAAgE;QAChE,4EAA4E;QAC5E,IAAI,QAAQ,KAAK,IAAI;YACnB,IAAI,GAAG,CAAC,IAAI,KAAK,CAAC,CAAC,CAAC,CAAC,CAAC,IAAI,GAAG,GAAG,CAAC,CAAC,CAAC,IAAI,GAAG,IAAI,CAAC;;YAE/C,IAAI,GAAG,GAAG,CAAC;QAGb,oDAAoD;QACpD,wDAAwD;QACxD,IAAI,GAAG,GAAG,CAAC,MAAM;YACf,CAAC,GAAG,CAAC,MAAM,GAAG,GAAG,CAAC,MAAM,CAAC,GAAG,GAAG,CAAC,MAAM,CAAC,IAAI,CAAC;gBAC5C,CAAC,GAAG,CAAC,MAAM,CAAC,IAAI,CAAC,GAAG,IAAI,GAAG,GAAG,CAAC,MAAM,CAAC,IAAI,CAAC,CAAC,CAAC;QAE/C,IAAI,IAAI,CAAC,GAAG,CAAC,IAAI,GAAG,GAAG,CAAC,MAAM,CAAC,GAAG,IAAI,GAAG,GAAG,CAAC,OAAO,EAAE,CAAC;YACrD,OAAO,GAAG,IAAI,CAAC,GAAG,CAAC,GAAG,CAAC,MAAM,GAAG,GAAG,CAAC,MAAM,CAAC,GAAG,GAAG,CAAC,OAAO,CAAC;YAC1D,OAAO,GAAG,CAAC,OAAO,GAAG,IAAI,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,IAAI,GAAG,OAAO,CAAC;YAClD,IAAI,GAAG,GAAG,CAAC,MAAM,GAAG,OAAO,GAAG,CAAC,GAAG,CAAC,MAAM,GAAG,GAAG,CAAC,MAAM,CAAC,CAAC;QAC1D,CAAC;QACD,IAAI,IAAI,CAAC,GAAG,CAAC,GAAG,CAAC,MAAM,GAAG,IAAI,CAAC,GAAG,IAAI,GAAG,GAAG,CAAC,OAAO,EAAE,CAAC;YACrD,OAAO,GAAG,IAAI,CAAC,GAAG,CAAC,GAAG,CAAC,MAAM,GAAG,GAAG,CAAC,MAAM,CAAC,GAAG,GAAG,CAAC,OAAO,CAAC;YAC1D,OAAO,GAAG,CAAC,OAAO,GAAG,IAAI,CAAC,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,IAAI,GAAG,OAAO,CAAC;YAClD,IAAI,GAAG,GAAG,CAAC,MAAM,GAAG,OAAO,GAAG,CAAC,GAAG,CAAC,MAAM,GAAG,GAAG,CAAC,MAAM,CAAC,CAAC;QAC1D,CAAC;QAED,kDAAkD;QAClD,WAAW,CAAC,GAAG,EAAE,IAAI,EAAE,CAAC,EAAE,GAAG,CAAC,IAAI,CAAC,CAAC;QACpC,MAAM,MAAM,GAAG,GAAG,CAAC,OAAQ,CAAC,IAAI,EAAE,GAAG,CAAC,IAAI,EAAE,GAAG,CAAC,QAAQ,EAAE,GAAG,CAAC,YAAY,CAAC,CAAC;QAC5E,GAAG,CAAC,MAAM,EAAE,CAAC;QACb,IAAI,MAAM,KAAK,CAAC;YACd,OAAO,cAAc,CAAC;QAGxB,6DAA6D;QAC7D,oEAAoE;QACpE,OAAO,GAAG,IAAI,CAAC;QACf,KAAK,GAAG,KAAK,CAAC;QACd,MAAM,GAAG,KAAK,CAAC;QACf,QAAQ,GAAG,IAAI,CAAC;QAEhB,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,KAAK,EAAE,CAAC,EAAE,EAAE,CAAC;YAC/B,IAAI,GAAG,CAAC,UAAU,CAAC,CAAC,CAAC,KAAK,CAAC;gBAAE,SAAS;YACtC,IAAI,IAAI,CAAC,GAAG,CAAC,GAAG,CAAC,QAAQ,CAAC,CAAC,CAAC,CAAC,KAAK,IAAI,EAAE,CAAC;gBACvC,IAAI,GAAG,CAAC,UAAU,CAAC,CAAC,CAAC,GAAG,GAAG,CAAC,MAAM,CAAC,CAAC,CAAC,IAAI,IAAI;oBAC3C,KAAK,GAAG,IAAI,CAAC;YACjB,CAAC;iBAAM,CAAC;gBACN,IAAI,CAAC,GAAG,CAAC,MAAM,CAAC,CAAC,CAAC,GAAG,GAAG,CAAC,QAAQ,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC;oBACrC,CAAC,GAAG,CAAC,UAAU,CAAC,CAAC,CAAC,GAAG,GAAG,CAAC,MAAM,CAAC,CAAC,CAAC,IAAI,IAAI,CAAC,EAAE,CAAC;oBAChD,KAAK,GAAG,IAAI,CAAC,GAAG,CAAC,GAAG,CAAC,QAAQ,CAAC,CAAC,CAAC;wBAC9B,CAAC,GAAG,CAAC,QAAQ,CAAC,CAAC,CAAC,GAAG,GAAG,CAAC,MAAM,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;oBACrC,IAAI,KAAK,GAAG,OAAO,EAAE,CAAC;wBACpB,MAAM,GAAG,IAAI,CAAC;wBACd,OAAO,GAAG,KAAK,CAAC;wBAChB,IAAI,GAAG,CAAC,CAAC;oBACX,CAAC;gBACH,CAAC;YACH,CAAC;QACH,CAAC;QAED,IAAI,MAAM,EAAE,CAAC;YACX,2DAA2D;YAC3D,GAAG,CAAC,MAAM,GAAG,IAAI,CAAC;YAClB,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,KAAK,EAAE,CAAC,EAAE;gBAC5B,GAAG,CAAC,MAAM,CAAC,CAAC,CAAC,GAAG,GAAG,CAAC,QAAQ,CAAC,CAAC,CAAC,CAAC;YAElC,IAAI,GAAG,CAAC,CAAC;YACT,iDAAiD;YACjD,IAAI,IAAI,CAAC,GAAG,CAAC,GAAG,CAAC,MAAM,GAAG,GAAG,CAAC,MAAM,CAAC,IAAI,GAAG,CAAC,OAAO;gBAClD,MAAM;YAER,SAAS;QACX,CAAC;QAED,IAAI,KAAK,EAAE,CAAC;YACV,sEAAsE;YACtE,GAAG,CAAC,MAAM,GAAG,IAAI,CAAC;YAClB,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,KAAK,EAAE,CAAC,EAAE;gBAC5B,GAAG,CAAC,MAAM,CAAC,CAAC,CAAC,GAAG,GAAG,CAAC,QAAQ,CAAC,CAAC,CAAC,CAAC;YAElC,MAAM;QACR,CAAC;QAED,sDAAsD;QACtD,6DAA6D;QAC7D,GAAG,CAAC,MAAM,GAAG,IAAI,CAAC;QAClB,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,KAAK,EAAE,CAAC,EAAE;YAC5B,GAAG,CAAC,MAAM,CAAC,CAAC,CAAC,GAAG,GAAG,CAAC,QAAQ,CAAC,CAAC,CAAC,CAAC;QAElC,IAAI,GAAG,CAAC,CAAC;QACT,sDAAsD;QACtD,IAAI,IAAI,CAAC,GAAG,CAAC,GAAG,CAAC,MAAM,GAAG,GAAG,CAAC,MAAM,CAAC,IAAI,GAAG,CAAC,OAAO;YAClD,MAAM;IACV,CAAC,CAAC,0BAA0B;IAE5B,yDAAyD;IACzD,GAAG,CAAC,QAAQ,GAAG,GAAG,CAAC,MAAM,CAAC;IAC1B,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,KAAK,EAAE,CAAC,EAAE,EAAE,CAAC;QAC/B,GAAG,CAAC,QAAQ,CAAC,CAAC,CAAC,GAAG,GAAG,CAAC,MAAM,CAAC,CAAC,CAAC,CAAC;QAChC,GAAG,CAAC,SAAS,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC;QACrB,IAAI,GAAG,CAAC,UAAU,CAAC,CAAC,CAAC,KAAK,CAAC;YAAE,SAAS;QACtC,IAAI,CAAC,IAAI,CAAC,GAAG,CAAC,GAAG,CAAC,MAAM,CAAC,CAAC,CAAC,CAAC,KAAK,IAAI,CAAC;YAClC,CAAC,GAAG,CAAC,UAAU,CAAC,CAAC,CAAC,GAAG,GAAG,CAAC,MAAM,CAAC,CAAC,CAAC,IAAI,IAAI,CAAC;YAC7C,GAAG,CAAC,SAAS,CAAC,CAAC,CAAC,GAAG,GAAG,CAAC,MAAM,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;QAEhD,IAAI,CAAC,GAAG,CAAC,MAAM,CAAC,CAAC,CAAC,GAAG,GAAG,CAAC,MAAM,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC;YACnC,CAAC,GAAG,CAAC,UAAU,CAAC,CAAC,CAAC,GAAG,GAAG,CAAC,MAAM,CAAC,CAAC,CAAC,IAAI,IAAI,CAAC;YAC7C,GAAG,CAAC,SAAS,CAAC,CAAC,CAAC,GAAG,GAAG,CAAC,MAAM,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;IAClD,CAAC;IAED,OAAO,OAAO,CAAC;AACjB,CAAC"}