Fri, 27 Aug 2021 16:56:16 +0300
a bit less workaround
src/Util.jl | file | annotate | diff | comparison | revisions |
--- a/src/Util.jl Fri Aug 27 16:35:50 2021 +0300 +++ b/src/Util.jl Fri Aug 27 16:56:16 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²>α² - @. @views y′[:, i] = y′[:, i]*(α/√n²) + @views y′[:, i] .*= (α/√n²) end end end