Function l2_norm 

pub fn l2_norm(u: &[f64], dx: f64) -> f64

Computes the discrete L2 norm ||u||₂ = sqrt(dx * Σ u[i]²).