fix: Normalize the column once more

The column may not be normalized if the `factor` is on a scale of 1e-40.
Possibly, f32 just runs out of precision.

There is likely a better solution to the problem.

src/linalg/householder.rs

@@ -34,6 +34,7 @@ pub fn reflection_axis_mut<T: ComplexField, D: Dim, S: StorageMut<T, D>>(
 
     if !factor.is_zero() {
         column.unscale_mut(factor.sqrt());
+        let _ = column.normalize_mut();
         (-signed_norm, true)
     } else {
         // TODO: not sure why we don't have a - sign here.