Fixed-point scaling error (Lemma 2.2):
ProvedIntMul.HvdH.scaleFP_boundinteger-multiplication
Let , let be an integer, and let with , so that lies in the unit disc. Let be any approximation of , and let be the output of the scaling procedure scaleFP (Lemma 2.2: take if it lies in the unit disc, otherwise renormalise by and round). Then , and the errors and satisfy
This is the error statement of Lemma 2.2 of Harvey–van der Hoeven, used for example in the final scaling step of an inverse FFT.
Preamble
import Definitions.Def_IntMul_HvdH_FixedPoint
Formal statement
namespace IntMul.HvdH
theorem scaleFP_bound (p c : ℕ) (hc : 1 ≤ c) (u ũ : ℂ) (hu : ‖u‖ ≤ (c : ℝ)⁻¹) :
‖scaleFP p c ũ‖ ≤ 1 ∧
(2 : ℝ) ^ p * ‖scaleFP p c ũ - c * u‖ < 2 * c * ((2 : ℝ) ^ p * ‖ũ - u‖) + 3 := 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, Section 2.2, Lemma 2.2 (error bound)