- NITER 默认值从 0 改为 30(匹配 Fortran NSTPAR PVALUE) - nlevel 从输入文件离子数据计算(39 for H-He),不再依赖 fort.8 - fort.7 输出从 numpar=3 升级到 numpar=42(含 39 能级占据数) - 初始化 LTE Saha-Boltzmann populations(H I/II + He I/II/III) - OPAINI iltref 未初始化保护(避免 usize 下溢 panic) - 复制原子数据文件到 hhe_rust/data/ Co-Authored-By: Claude Opus 4.8 <noreply@anthropic.com>
31 lines
2.3 KiB
Plaintext
31 lines
2.3 KiB
Plaintext
# F2R Phase 3 验证任务列表
|
||
# 格式:每行一个任务,完成后在行首加 ✅
|
||
|
||
# --- 数据管道修复(已完成)---
|
||
✅ synspec_outpri_zero: OUTPRI 0 points → 修复: fidata解析(Fortran引号解析器) + 频率网格生成 + RDATA数据文件
|
||
✅ synspec_data_pipeline: INITIA→RDATA(34 levels)→FREQ(144pts)→RESOLV→RTECD→OUTPRI(141pts) 全链贯通
|
||
✅ synspec_fidata_parse: 修复 Fortran 自由格式引号字符串解析(fortran_free_format_parse)
|
||
✅ synspec_freq_grid: 在 runner 中生成基本连续谱频率网格(144点等对数间距)
|
||
✅ synspec_rdata_files: 创建最小原子数据文件 h1.dat(9级), he1.dat(14级), he2.dat(14级)
|
||
✅ synspec_nan_flux: 修复3处bug — 1)HK/BN物理常数(1.0→4.79928e-11/1.4743e-2) 2)SCE逐深度计算(ane*SIGE) 3)frx1/frx2对数插值权重(0→正确值)。SYNSPEC fort.7 141点有限输出,TLUSTY fort.7 83点有限输出
|
||
|
||
# --- TLUSTY 输出升级(本轮完成 2026-06-11)---
|
||
✅ tlusty_niter_fix: NITER 默认值从 0 改为 30(匹配 Fortran NSTPAR PVALUE 第64项)
|
||
✅ tlusty_nlevel_fix: nlevel 从输入文件离子数据计算(39 for H-He),不再依赖 fort.8 的 numpar
|
||
✅ tlusty_numpar_upgrade: fort.7 输出从 numpar=3 升级到 numpar=42(39 level populations + T/Pe/rho)
|
||
✅ tlusty_lte_popul_init: 初始化 LTE Saha-Boltzmann populations(H I/II + He I/II/III 39能级)
|
||
✅ tlusty_opaini_guard: OPAINI iltref 未初始化保护(避免 usize 下溢 panic)
|
||
|
||
# --- TLUSTY 输出验证 ---
|
||
✅ tlusty_output_format: fort.7 643行输出,格式正确(nd=70, numpar=42),0 NaN/Inf,LTE populations 物理合理
|
||
tlusty_start_init: START 需完整初始化原子数据(RDATA 读能级文件 → 能级能量/权重 → iltref → 连续截面)
|
||
tlusty_solves_fix: SOLVES 矩阵需要非零 populations 才能工作(需 START init 完成后启用)
|
||
|
||
# --- SYNSPEC 验证(需 Fortran 参考可用后继续)---
|
||
synspec_resolv_verify: 逐行对比 RESOLV Rust vs Fortran(需 fort.19 谱线表文件)
|
||
synspec_opac_verify: 逐行对比 OPAC Rust vs Fortran(需 Fortran 参考输出)
|
||
synspec_rtecd_verify: 逐行对比 RTECD Rust vs Fortran(需 Fortran 参考输出)
|
||
|
||
# --- 通用 ---
|
||
synspec_data_sync: 获取完整 TLUSTY 原子数据文件(含光电离截面和连续跃迁数据)→ 匹配 Fortran hhe35lt.7 参考
|