Mercurial > repos > AlgTools / changeset

--- a/src/Util.jl	Tue Apr 13 15:51:51 2021 -0500
+++ b/src/Util.jl	Fri Aug 27 16:35:50 2021 +0300
@@ -195,7 +195,7 @@
         @inbounds @simd for i=1:size(y′, 2)# in CartesianIndices(size(y)[2:end])
             n² = norm₂²(@view(y′[:, i]))
             if n²>α²
-                y′[:, i] .*= (α/√n²)
+                @. @views y′[:, i] = y′[:, i]*(α/√n²)
             end
         end
     end