精读笔记

Problem Setting

这篇论文不是在重新定义 ED degree，而是在解决“如何实际计算 ED degree”这个长期存在但常被低估的问题。ED degree 理论给出了 nearest-point critical equations，但研究中遇到的模型输入形态非常不统一：有的由隐式方程给出，有的是 complete intersection，有的是某个 complete intersection 的 component，有的只有参数化。每种情形下，直接构造临界理想都会遇到不同瓶颈。

真正困难点在于，ED degree 是一个代数计数问题，但计算路径会在 symbolic complexity 和 numerical reliability 之间摇摆。minors + saturation 在理论上干净，但 minors 数量和 saturation 成本很快爆炸；Lagrange multiplier formulation 更轻，但需要 complete intersection 假设；homotopy continuation 更 scalable，但会引入 path failure、infinity filtering 和低估风险；参数化模型若先隐式化，通常代价不可接受。关键矛盾是：几何定义统一，但可计算表示高度依赖模型结构。

Motivation

已有路线不够的地方不是没有算法，而是缺少一个研究者能直接使用的、覆盖多种 ED degree 计算场景的专门工具。AlgebraicOptimization 这类包更 general-purpose；ToricInvariants 只覆盖特定 toric 设置；Bertini、PHCpack、HomotopyContinuation.jl 能求 polynomial systems，但不会替用户决定如何形成 ED critical system。

作者的核心观察是，ED degree 计算的不同公式本质上都在表达同一个条件：数据点到 variety 上 critical point 的 displacement vector 位于法空间。于是 package 的价值在于将这个几何条件根据输入表示重新编码，而不是让用户每次手工搭建方程。关键缺口是 workflow-level 的：从 variety description 到 ED degree count 之间缺少可复用、可检查、可切换 symbolic / numerical 的中间层。

Core Idea

论文的核心思想是把 ED degree computation 组织成一组围绕法空间条件的计算后端，而不是绑定到单一公式。对于隐式 variety，augmented Jacobian 的 rank-drop 给出 critical locus；对于 complete intersection，引入 Lagrange multipliers 后可以避免显式 minors；对于重复改变 weights/data 的实验，用 homotopy continuation 将 generic start solve 和 target tracking 分离；对于参数化模型，直接在参数空间构造 conormal condition，避免隐式化。

这个思路理论上成立，因为 ED criticality 本质上是梯度正交于 tangent space，等价于 displacement vector 属于 conormal space。不同方法的区别只是如何表示 tangent / normal information。和 prior 的本质差异不是数学定理层面的新 ED characterization，而是把这些 characterization 变成一个 package-level computational abstraction，使用户可以在 exactness、scalability、输入表示之间切换。

Method

Minors 方法解决的是最直接的 exact computation：通过 augmented Jacobian 的 rank condition 构造 critical ideal，并用 saturation 去除奇异点贡献。它的必要性在于给出可验证、可导出的 symbolic object；核心变化是 ED degree 不只是一个数，也可以返回 critical ideal，便于其他系统复现。

Left-kernel 方法解决 large codimension 下 minors 爆炸的问题。它用 Lagrange multipliers 表达距离梯度落入 constraint gradients span 中，把 rank condition 改写成 square polynomial system。核心变化是从组合数量的 minors 转向一个数值可解的临界方程组，但代价是 complete intersection 假设变成正确性前提。

Homotopy 方法解决的是 repeated ED computation，尤其是不同 weights 或 data 下的实验。它先在 generic setting 下求 start system，再追踪到 target system。核心变化是把一次昂贵的 generic solve 复用为后续计算的起点，因此更适合探索 special weights/data，而不是每次从头求解。

Conormal 方法解决参数化 variety 的 ED degree 计算。它不先求隐式方程，而是从参数化 Jacobian 的左核生成法空间，再令 weighted displacement 落在这个法空间中。核心变化是把计算留在参数域，避免隐式化这个通常不可承受的中间步骤。

Key Insight / Why It Works

最关键的 insight 是：ED degree computation 的可扩展性主要取决于如何表示法空间，而不是取决于 ED degree 定义本身。minors、Lagrange multipliers、homotopy、conormal parameterization 都是在不同坐标系下表示同一个 conormal incidence relation。论文真正有价值的地方，是把这个关系做成多个 computational entry points。

最可能的核心贡献是 homotopy workflow 和 parameterized conormal workflow 的组合。前者把 repeated weighted/unit/generic experiments 的成本结构改变了：stage one 只做一次，stage two 做参数追踪；后者让 PNN、multiview 这类天然参数化模型不用先走隐式化。这里的增益不是新的 algebraic shortcut，而是避免选择错误表示导致的计算爆炸。

哪些部分只是辅助也很清楚：具体 Macaulay2 method 命名、TempDirectory、Bertini file inspection 等主要是 engineering usability。ReturnCriticalIdeal=>true 有实用价值，但它的理论内容仍是经典 critical ideal。数值方法的性能提升可能主要来自 homotopy continuation 和已有 Bertini / MonodromySolver 能力，package 自身更多是在组织 equations 和 workflow。

这不是 scaling law，也不是 data coverage，也不是 representation learning；它属于 computational algebraic geometry 中的 representation alignment：把同一个几何对象对齐到最适合求解器的方程表示。所谓 generality 来自多后端覆盖，而不是一个统一算法在所有输入上都更强。文中未充分说明在病态输入、非 complete intersection、奇异 component 或 ill-conditioned homotopy 下如何保证不低估。

Relation To Prior Work

最接近的路线是 ED degree 经典理论中的 critical ideal / conormal formulation，以及 algebraic optimization 软件中的 Fritz John 或 Lagrange multiplier 方程。minors 方法直接继承 Draisma et al. 的 ED degree critical equations；left-kernel 方法与 Fritz John / Lagrange multiplier 系统在精神上很接近；homotopy 部分依赖 numerical algebraic geometry 的 standard parameter homotopy；conormal 部分也是参数化模型中自然的法空间构造。

真正不同点在于工具整合和问题边界，而不是单个数学构造。相比 AlgebraicOptimization，它更专门面向 ED degree，并把 symbolic 与 numerical routes 都纳入同一接口；相比 ToricInvariants，它不局限于 toric/projective special structure；相比直接调用 Bertini 或 HomotopyContinuation.jl，它替用户完成 ED-specific equation formation 和 filtering。

看似新的地方，如 left-kernel 和 homotopy method，本质上是已有思想的工程化重组；实质创新更像 package architecture：按输入表示和计算目标选择合适的 conormal/critical formulation，并暴露足够中间对象让研究者复现或换 solver。

Dataset / Evaluation

evaluation 不是传统 benchmark，而是一组 algebraic examples 和 application demonstrations。覆盖了 symbolic exact examples、complete intersection numerical examples、homotopy reuse、参数化 variety、PNN neurovariety 和 multiview variety。任务范围比较贴近 applied algebraic geometry 研究实践，尤其是那些 ED degree 作为复杂度指标的场景。

但这些实验主要验证“能算”和“接口覆盖”，没有充分验证 robustness 或 scalability claim。论文没有系统比较不同方法在同一族问题上的 runtime、memory、failure rate、path count，也没有评估数值低估发生频率。PNN 和 multiview 的例子更像用 package 生成新现象或 conjectural evidence，而不是严格 benchmark。

因此 evaluation 支持的是：这个包能把多个已有计算范式统一进 Macaulay2，并可用于实际研究探索。它不充分支持更强的说法，例如相对已有软件显著更快、更稳定、或更自动化。

Limitation

核心限制首先是数学前提。left-kernel 方法要求 complete intersection；homotopy 方法要求 G 定义的 complete intersection 包含 V(F) 作为 irreducible component；参数化 conormal 方法要求参数化 Jacobian generically full rank，并且选择的 kernel generators 在一般点能张成正确维数。这些不是小注释，而是方法成立的边界。

第二个限制是数值可靠性。论文明确承认 numerical methods 可能因 path failure 或路径发散到无穷而低估 ED degree。这个问题没有被 package 消除，只是通过暴露 Bertini 配置和临时文件让有经验用户介入。换言之，困难从“构造方程”部分转移到了“诊断数值求解是否完整”部分。

第三个限制是 scalability 上限仍在。symbolic minors + saturation 对复杂 variety 仍然会爆；homotopy 的路径数本身可能等于或超过 ED degree 的复杂度；参数化方法避免隐式化，但会把复杂度转移到 kernel computation 和 critical equations。增益来源不清的地方在于，论文没有量化到底是 formulation 更好、solver 更强，还是 examples 本身规模可控。

第四，package 没有提供强自动策略：什么时候该用 minors、left-kernel、homotopy、conormal，文中更多靠用户判断。对于研究者这可以接受，但对于可靠批量计算仍不足。

Takeaway

• 第一，ED degree 计算的关键不是寻找一个万能公式，而是根据 variety 的表示选择法空间条件的编码方式。
• 这个 insight 可以迁移到其他 algebraic optimization degree 或 critical-point counting 问题。
• 第二，参数化模型应尽量在参数域构造 conormal / critical equations，隐式化往往是错误中间目标。
• PNN 和 multiview examples 的真正启发在这里。

一句话总结

这篇论文在 ED degree 方向中的位置是一篇工具型但有研究价值的 computational paper：它没有提出新的 ED 理论核心，而是把 conormal/critical-locus 思想系统地工程化为 Macaulay2 中可切换的 symbolic 与 numerical 计算框架。