Thanks for sending an issue! Please fill in the following template to help quickly solve your problem.

Describe the current behavior / 问题描述 (Mandatory / 必填)

SIMT Bessel 函数族（j0f/j1f/jnf/y0f/y1f/ynf）在以下输入区间存在精度问题：

1. j0f/j1f/y0f/y1f 大 x（x ≳ 4e8）输出完全错误

原 >8 段用 __internal_trig_red_slowpath_f_fast_mode（uint64 两步 Cody-Waite 相位规约），x > ~4e8 后 a * (2/π) 的 float 舍入把相位低位吃掉，fmaf(q, -π/2_hi, a) 消差后相位全错，sin/cos 多项式喂入大残差 → 输出 1e22~1e28 量级垃圾值或 -inf。>1e13 直接 return 0 放弃精度。

典型数据（y0f，x ≈ 1e12）：

2. ynf 小 x / subnormal 段 ULP 超标

__internal_yn_case1 用 powf(2/x, n) 走 exp(n·log(2/x)) 链，subnormal x 下 log 精度有限 + n 次幂放大误差，且 (2/x)^n 易溢出 -inf。

典型数据（ynf）：

3. 特殊值语义不完整

nan 返回 INF/INF（payload 不固定，非标准传播）；+inf 走 larger8 的 isinf → 0，但与 golden（scipy 返回 nan）不一致；ynf 对 x<0 返回 NaN 正确但 ±0 语义未显式处理。

Environment / 环境信息 (Mandatory / 必填)

• 测试环境：Ascend NPU，CANN 9.2.0，--npu-arch=dav-3510

Steps to reproduce the issue / 重现步骤 (Mandatory / 必填)

同上典型数据

Describe the expected behavior / 预期结果 (Mandatory / 必填)

• j0f/j1f/y0f/y1f：|x|<8 走 9 ULP，|x|≥8 走 atol 2.2e-6
• ynf：|x|<n 走 ULP ceil(2+2.5n)，否则走 atol 2.2e-6

Related log / screenshot / 日志 / 截图 (Mandatory / 必填)

见上方典型数据表。

Special notes for this issue/备注 (Optional / 选填)

修复方案见 PR #5661：

• j0f/j1f/y0f/y1f：三段拆分（less8 / middle_range / huge_range）+ Payne-Hanek 大 x 相位规约
• ynf：case1 改前向递推 + Newton 精化 2/x + 激进归一化，case3 调 case1
• 共用 helper 改名（__internal_bessel_middle_* 等）