lingshu-solver 1.0.4 → 1.0.6

package/index.html

@@ -5833,6 +5833,9 @@ function suan22(state) {
             _walk1(_eq);
         })();
         if (_lo < _hi) {
+            var _dense = false;  // 修3：根稠密/超时提前终止扫描（防宽域周期方程穷举挂起）
+            var _denseByCount = false;  // 修3：仅由"根数超阈值"触发的稠密标记（用于无限/截断判定，区别于超时触发）
+            var _denseNonPeriodic = false;  // 修3：非周期稠密截断标记（与 periodic→无限 分支互补）
             // ---- 周期感知步长（issue A 修复 2026-09-01）：高频/多周期方程用 < 周期/2 的网格，
             //      确保不漏根；非周期方程回退到 4096 点默认细网格。点数上限 _NMAX 护栏防宽域爆炸。 ----
             var _P = _detectPeriod1D(_eq, _vn);
@@ -5902,8 +5905,11 @@ function suan22(state) {
                             _recScanInterval(_eq, _vn, _gxP, _gx, state, allSolutions, resultVarNames, 0);
                         }
                     }
+                    if (allSolutions.length > 200) { _dense = true; _denseByCount = true; break; }
+                    if ((performance.now() - state.startTime) > 1500) { _dense = true; break; }
                     _gxP = _gx; _gfP = _gf;
                 }
+                if (_dense) break;
                 // 终止端点若有根也记录
                 if (_gfP !== null && Math.abs(_gfP) < 1e-7) {
                     var _fbV = reconstructSolution(state, [_B]).map(roundToGrid);
@@ -5951,6 +5957,31 @@ function suan22(state) {
             if (_dedupMap[_mi]) mergedSolutions.push(_dedupMap[_mi]);
         }
 
+        // 修3（2026-09-04）：稠密根防护——单变量周期方程在宽域内根稠密，
+        // 继续穷举既耗时（曾致 solve 挂起 >110s）又产出无用的海量解集。
+        // 含周期函数且根数超阈值 ⇒ 解集实际为无限（周期平移），如实标 rt=3 无限解集(代表解)；
+        // 非周期但稠密（如高次多项式）⇒ 不谎称无限，置 truncated 并后续告警，诚实返回部分代表根。
+        if (_denseByCount) {
+            if (_P !== null) {
+                state.done = true;
+                state.result = {
+                    solutions: mergedSolutions,
+                    error: null,
+                    message: "单变量周期方程在声明域内出现大量（>" + 200 + "）孤立根，呈周期性稠密分布，解集实际为无限（如 x = x₀ + k·" + _P.toFixed(4) + "）。已给出 " + mergedSolutions.length + " 个代表根；如需有限个解请收窄变量定义域。",
+                    executionPath: "单变量稠密周期根判定(无限解集)",
+                    timeMs: performance.now() - state.startTime,
+                    confidence: "high",
+                    varNames: resultVarNames,
+                    resultType: 3,
+                    resultTypeName: "无限解集(代表解)",
+                    resultTypeDesc: "周期方程在宽域内根稠密，解集为周期平移的无穷集合，已给出代表根"
+                };
+                return;
+            }
+            state.truncated = true;  // 非周期稠密：诚实截断（避免海量/不完整解集），下方 domainFiltered 产出部分代表根
+            _denseNonPeriodic = true;
+        }
+
         // 2) 原始方程回代验证：对每组候选解，代入全部原始方程重新计算残差
         //    二分查找得到的解（_bisect=true）跳过残差验证（区间收敛保证精度）
         var _origEqs = state.originalEquations || state.equations;
@@ -6082,6 +6113,12 @@ function suan22(state) {
                 }
             }
         }
+        // 修3：非周期稠密截断告警（与上方 periodic→无限 分支互补）
+        if (_denseNonPeriodic && state.result) {
+            state.result.warnings = (state.result.warnings || []).concat([
+                "方程在声明域内根数量过多（>" + 200 + "），已截断扫描；结果不完整，缩小变量范围或确认方程规模可获完整解集。"
+            ]);
+        }
     }
 }
 
@@ -6624,7 +6661,13 @@ function _affineEval(ast, env) {
     // 超越函数降级：仿射自变量 → 保守区间函数包围（复用 _rangeEval，无递归环）
     var tmp = { type: 'func', name: ast.name, arg: { type: 'var', name: '__a' } };
     var rv = _rangeEval(tmp, { __a: iv });
-    return rv ? _Aff((rv.min + rv.max) / 2, {}) : null;
+    if (!rv || !isFinite(rv.min) || !isFinite(rv.max)) return null; // 含 ±Inf：交还 _rangeEval 直接给出保守区间
+    // 关键修复（P0-2）：旧实现 return _Aff((rv.min+rv.max)/2, {}) 把合法区间宽度丢成 0，
+    // 致 intervalEval(sin x)@[-2,2] 塌成 [0,0]（unsound：真值须含 [-1,1]）。
+    // 正确做法：用半径 (rv.max-rv.min)/2 的噪声符号把 rv 完整编码为仿射，
+    // 使 _affToInterval 还原为 [rv.min, rv.max]（sound 且紧致），Krawczyk 雅可比随之可信。
+    var e = {}; e[_affSym++] = (rv.max - rv.min) / 2;
+    return _Aff((rv.min + rv.max) / 2, e);
   }
   if (ast.type === 'unary') {
     var u = _affineEval(ast.operand, env);
@@ -11378,7 +11421,11 @@ function _newtonRefine(eqs, vns, x0, box) {
         var F = []; for (var e = 0; e < n; e++) { var fe = evalAST(eqs[e], vmap); if (!isFinite(fe)) return null; F.push(fe); }
         var rms = 0; for (var i2 = 0; i2 < n; i2++) rms += F[i2] * F[i2]; rms = Math.sqrt(rms / n);
         if (rms < 1e-11) return x;
-        var J = []; for (var a = 0; a < n; a++) { J.push([]); for (var b = 0; b < n; b++) { var d = _diffAST(eqs[a], vns[b]); var dv = intervalEval(d, vmap); if (!dv) return null; J[a].push((dv.min + dv.max) / 2); } }
+        // 雅可比区间求导：必须喂退化的点区间 {min:x,max:x}（而非裸点值）。
+        // 旧实现 intervalEval(d, vmap) 把裸数当 iv，iv.min/iv.max 为 undefined → _buildAffEnv 算出 NaN 噪声
+        // → 雅可比整行 NaN → realMatInv 失败 → 本函数对一切良置问题恒返 null（P0-1）。
+        var pmap = {}; for (var i3 = 0; i3 < n; i3++) pmap[vns[i3]] = { min: x[i3], max: x[i3] };
+        var J = []; for (var a = 0; a < n; a++) { J.push([]); for (var b = 0; b < n; b++) { var d = _diffAST(eqs[a], vns[b]); var dv = intervalEval(d, pmap); if (!dv) return null; J[a].push((dv.min + dv.max) / 2); } }
         var Y = realMatInv(J); if (!Y) return null;
         var dX = rMatVec(Y, F);
         for (var j = 0; j < n; j++) x[j] = x[j] - dX[j];

package/package.json

@@ -1,6 +1,6 @@
 {
   "name": "lingshu-solver",
-  "version": "1.0.4",
+  "version": "1.0.6",
   "description": "灵数求解器 · 确定性实数方程组求解 MCP 工具（≤6变量）：区间收缩 + Krawczyk 认证，离线零依赖零数据，诚实标记 truncated。",
   "homepage": "https://genesis-plan.github.io/lingshu-solver/",
   "repository": {
@@ -50,6 +50,6 @@
   },
   "scripts": {
     "start": "node mcp-server.js",
-    "test": "node test/regression.js"
+    "test": "node test/p0_fix_regression.js"
   }
 }