fmq/DCTS

Files

T

fmq b058e66722 feat(all): NaN 伪收敛否决与 nl_tight 能量回退、emflux 检验 4π 修正、nst 行宽与种子边界修复、gfATO 谱线表接通与网格加密 9216 点、阶段分项统计与跳板机部署

物理修复（400 失败点归因，见 docs/failed400_nan_pseudo_convergence_2026_08_17.md）：
- runner: 阶段 converged 后复查 fort.7，含 NaN/Inf 即否决（fort.9 全零伪收敛，
  曾致 261 点误跳过 nl_direct 回退）；否决阶段不产出种子，阻断污染传播
- runner: nl_tight 回退——仅能量边际失败时以 CHMAX 收紧 10× 从自身模型续迭代，
  残差降幅 ~10×；execute_tlusty_stage 抽取供主链与回退共用
- conv_check: fort.14 为 Eddington 通量 Hλ，积分需乘 4π 再比 σTeff⁴
  （旧版 ratio 稳定 0.0796=1/4π，全点系统性假阳性）+ 回归测试
- nst_writer: 单行超 80 字符被 TLUSTY 静默截断，IFALI/JALI/TRAD 等从未生效；
  按 75 字符自动换行
- seed_finder: Teff 容忍度改含边界 <=，相邻 5000K 档恢复互为种子 + 回归测试

谱线表与网格：
- 默认线表 gfVIS99 → gfATO（全波段 18-23000Å），TaskSpec.linelist 支持工作流
  级覆盖，节点按需下载（进程互斥锁防并发重复下载 238MB）
- sdB_cno Teff 加密至 5000K 步长，432 → 9216 点；tlusty/synspec 静态二进制更新

统计与部署：
- grid 汇总改按 tlusty_status/synspec_status 分项计数，新增 tlusty_failed/
  synspec_failed/synspec_pending，前端详情页双视图适配
- deploy/fetch_results 支持跳板机 ProxyJump 与 SSH 主连接复用，fetch 新增 --force；
- Docker 构建支持 CARGO_MIRROR/USE_MIRRORS 国内镜像参数；移除 tools/ 拷贝
- 新增 tlusty-synspec-test skill 与 6 篇根因分析/验证文档

2026-08-17 23:55:26 +08:00

2.5 KiB

Raw Permalink Blame History

{test_id} / {type} / {grid} / {content}

目的

一句话说清这次测试想验证什么。例：复现生产 t60000 冷启动 nl 阶段发散，对比 baseline 与 itek0/iacc0/dpsilg15 修复组合。

配置

二进制: dcts/assets/tlusty_static（或 synspec_static）
网格: t{teff}_g{logg}_he{loghe}（如 t60000_g5.0_he-2，Teff=60000K, logg=5.0, He/H=1e-2）
链: 冷启动 lte→nc→nl / 种子 seed_nc→nl / SYNSPEC
变体: 与 baseline 不同的关键参数（ITEK/IACC/DPSILG/ORELAX/NITER）

输入

inputs/ 生成方式: 代码生成（gen_input5.rs/nst_writer.rs）/ 模板修改（留 diff）/ 复制已验证输入
关键参数: 例 nc NITER=10, nl NITER=100, DPSILG=1.5

执行

命令: bash scripts/run.sh
各阶段: lte(rc=0) → nc(rc=0) → nl(rc=?)
耗时: 例 lte 8s + nc 312s + nl 421s

结果（TLUSTY）

阶段	迭代	末 max_relc	STOP	温度结构
lte	0	—	—	—
nc	10	1.4e0	无	物理
nl	19	1.2e16	STOP iter19	物理(停止时)

max_relc 来自 check_fort9.py；温度结构来自 check_temperature.py（DCTS 权威判据）。

结果（SYNSPEC，若适用）

指标	值
采样点数	N
波长范围	ALAM0 – ALAST Å
NaN/Inf/坏行	0
分段 NaN (EUV/UV/可见/NIR)	各段计数

审查结论

数值: ❌/✅（max_relc 是否 < chmax=0.001；有无 STOP/雪崩）
物理: ❌/✅（温度结构三项：表层 < 3×Teff、各深度 ∈ [10,1e8]K、无 NaN）
结论: 一句话。区分数值收敛与物理收敛——只有后者才是真正成功。

例：

数值: ❌ nc 雪崩(9.7→2.3e15), nl STOP(1.23e16)
物理: ✅ 停止时温度结构物理（表层 0.85×Teff，无 NaN）
结论: 精确复现生产失败；修复组合(itek0/iacc0/dpsilg15)使 nl max_relc 降至 48.6 但仍 STOP，未达数值收敛。

复现

cd dcts/test/{test_id}/{type}/{grid}/{content}
bash scripts/run.sh
# 事后审查（TLUSTY）:
python3 ../../../../../../.agents/skills/tlusty-synspec-test/scripts/check_fort9.py outputs/nl.9 --teff 60000
python3 ../../../../../../.agents/skills/tlusty-synspec-test/scripts/check_temperature.py outputs/nl.7 --teff 60000

产物清单（outputs/）

文件	内容
`lte.6/.7/.9/.14`	lte 阶段产物
`nc.6/.7/.9/.14`	nc 阶段产物
`nl.6/.7/.9/.14`	nl 阶段产物（最终模型 nl.7）
`{name}.spec/.cont/.iden`	SYNSPEC 光谱产物（若适用）

2.5 KiB Raw Permalink Blame History Unescape Escape

{test_id} / {type} / {grid} / {content}

目的

配置

输入

执行

结果（TLUSTY）

结果（SYNSPEC，若适用）

审查结论

复现

产物清单（outputs/）

2.5 KiB

Raw Permalink Blame History