HvdH Prop. 4.7(i): with
ProvedIntMul.HvdH.prop_4_7_ifourier-analysisinteger-multiplication
Gaussian resampling reduces a DFT of length to one of length (Harvey–van der Hoeven, Proposition 4.7(i), one-dimensional case).
Let be coprime positive integers and let be an integer with , where . Write for the normalised DFT (Mathlib's ZMod.dft divided by ). Then there are linear maps and , both of operator norm at most for the sup norm, with
In the paper, is the Gaussian resampling map and is built from the permutations , the row selection , the diagonal and the inverse of with . The proof formalises the resampling identity (Theorem 4.2) via Poisson summation for Gaussians and Lemma 4.6.
This is a node in the integer-multiplication-kappa campaign: it is the analytic core of the HvdH algorithm.
Preamble
import Mathlib
Formal statement
namespace IntMul.HvdH
theorem prop_4_7_i (s t α : ℕ) [NeZero s] [NeZero t] (hst : s < t) (hcop : Nat.Coprime s t)
(hα : 2 ≤ α) (hθ : 1 ≤ (α : ℝ) ^ 2 * ((t : ℝ) / s - 1)) :
∃ (A : (ZMod s → ℂ) →L[ℂ] (ZMod t → ℂ)) (B : (ZMod t → ℂ) →L[ℂ] (ZMod s → ℂ)),
‖A‖ ≤ 1 ∧ ‖B‖ ≤ 1 ∧
∀ u, (s : ℂ)⁻¹ • ZMod.dft u = (2 : ℂ) ^ (2 * α ^ 2) • B ((t : ℂ)⁻¹ • ZMod.dft (A u)) := 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 4.3, Proposition 4.7(i) (one-dimensional case; alpha taken integral), using Theorem 4.2 and Lemmas 4.5, 4.6