Error of the fixed-point complex product
ProvedIntMul.HvdH.mulZ_errinteger-multiplication
Let with , and let with . The fixed-point product (both parts rounded towards zero) approximates with
This is the fixed-point multiplication error estimate (Lemma 2.4 and Corollary 2.10 of Harvey–van der Hoeven) in the form used for twiddle factors, chirps and pointwise products.
Preamble
import Definitions.Def_IntMul_HvdH_ConvC
Formal statement
namespace IntMul.HvdH
theorem mulZ_err (p : ℕ) (a b : ℤ × ℤ) (x y : ℂ) (ha : ‖gC a‖ ≤ 2 ^ p) (hy : ‖y‖ ≤ 1) :
‖gC (mulZ p a b) - 2 ^ p * (x * y)‖ ≤
‖gC a - 2 ^ p * x‖ + ‖gC b - 2 ^ p * y‖ + Real.sqrt 2 := by sorry
end IntMul.HvdHSource
Harvey, van der Hoeven, Integer multiplication in time O(n log n), Annals of Math. 193 (2021), https://hal.science/hal-02070778, Lemma 2.4, Corollary 2.10