Skip to content

Drizzle_PRECISE

fujiaze edited this page Aug 5, 2026 · 4 revisions

PRECISE Drizzle

2026-08-05 冻结取代:本文档的精度/支持域/验收语义已由 Phase1_Drizzle_Freeze_SpecPhase1_Drizzle_Acceptance_Spec 取代 (内部关键几何统一 FP64,数据面支持 FP32/FP64,支持 Sphere→Plane 双向底层)。 本文保留为历史实现参考。

PRECISE是当前唯一生产Drizzle路径。FAST暂停。

科学语义不变

[ F_p=\sum_j L_jrac{a_{jp}}{A_{j,drop}} ]

support继续由真实球面重叠面积计算。

性能重构原则

本轮不采用降低精度的FAST近似,而是消除参考实现中的重复计算与低效数据结构:

  1. 以NESTED父Tile为累加单位;
  2. 每线程处理连续源图像条带,保持空间局部性;
  3. thread-local结构只以parent tile为key,不以每个leaf ipix为hash key;
  4. Tile内部用连续typed数组按local_ipix直接寻址;
  5. merge按Tile执行,不逐个hash key合并;
  6. release路径仅保留signal和support所需累计;
  7. pixfrac=1提供共享顶点行精确快速路径;
  8. 其他pixfrac使用批量WCS,不在内环分配vector;
  9. 候选生成与球面重叠分别剖析;
  10. HISS按Tile直接流式写入,不先构造全局6159万项map。

精度

FP32和FP64必须是真正两套Scalar实例。不得继续共享固定double几何。

Clone this wiki locally