pub fn l2_norm(f: &[f64], dx: f64) -> f64
Compute the discrete L² norm ‖f‖ = √(⟨f, f⟩).
‖f‖ = √(⟨f, f⟩)