Files
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

66 lines
2.5 KiB
Markdown
Raw Permalink Blame History

This file contains ambiguous Unicode characters
This file contains Unicode characters that might be confused with other characters. If you think that this is intentional, you can safely ignore this warning. Use the Escape button to reveal them.
# {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 | NaN | 温度结构 |
|---|---|---|---|---|---|
| lte | 0 | — | — | 0 | — |
| nc | 10 | 1.4e0 | 无 | 0 | 物理 |
| nl | 19 | 1.2e16 | STOP iter19 | 0 | 物理(停止时) |
> 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,未达数值收敛。
## 复现
```bash
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 光谱产物(若适用) |