matmul_ozaki — DRIVE ozaki op

• 資料種類:matrix × matrix → matrix

• 呼叫: import fullseye as fs; fs.ledger.matmul_ozaki(A, B, n_slices=None, scheme: 'str' = 'ozaki1', tol: 'float' = 1.1102230246251565e-16, on_insufficient: 'str' = 'raise', return_info: 'bool' = False)(要直接呼叫實作,import ozakimm; ozakimm.matmul_ozaki(A, B, n_slices=None, scheme: 'str' = 'ozaki1', tol: 'float' = 1.1102230246251565e-16, on_insufficient: 'str' = 'raise', return_info: 'bool' = False);從台帳取用則 opsdrive.get("matmul_ozaki"))

用法

由精確的低精度部分積組裝 FP64 精度的矩陣乘積 A @ B(Ozaki-I / Ozaki-II)。

scheme="ozaki1":切成 n_slices 片(乘積 n(n+1)/2 次)。n_slices=None 時,自動選擇使誤差保證上界不超過 tol · (|A||B|)_ij 的最少片數(1〜16 片)。達不到時,on_insufficient="raise" 引發 ValueError,"fallback" 回傳一般的 float64 乘積(只有這種情況下不再保證可重現性。可由 return_info=True 的 fallback 得知)。scheme="ozaki2":模數的個數必須以 n_slices(2〜20)指定(乘積 n 次)。沒有自動選擇。return_info=True 時回傳 (C, info)(n_slices、n_products、scheme、fallback)。

無論哪種方式,相同輸入都不受內部重排、BLAS 執行緒數影響而得到相同的位元。在 CPU 上比 DGEMM 慢 1〜2 個數量級(用法區分參見 module 的 docstring)。

參考(範例資料、文獻)

• 範例資料目錄(下載 URL / 授權) —— 2-D 用 skimage.data(BSD/公有領域)加合成圖,3-D 給出真實資料源(Stanford/PDS 等)的下載 URL。

• 運算子來歷與參考文獻 —— 該運算子族所依據的研究/方法出處。

• 演算法的正典(作者、年份)與用途見上面的族使用指南。

可執行的範例(實際呼叫該運算子並已驗證的樣例)

• poc_reproducible_icp — py -3.11 examples/poc_reproducible_icp.py

型別可銜接的下一個運算子(可接受 matrix 作為輸入)

world_camera · lidar_scan · relative_motion · foe_from_motion · flow_from_depth_motion · ttc_truth · world_materials · mesh_signed_volume

同類別(ozaki)

matmul_reproducible · ozaki_error_bound · cross_covariance_reproducible · kabsch_reproducible · fp64_emulation_probe

*Provenance: ozakimm.py — DRIVE 運算子登記表。本條目由 tools/opdocs.py md 自動產生(請勿手動編輯)。*

© 2026 Kazufumi Furuse — Fullseye operator documentation. Licensed under Apache-2.0.