# HG changeset patch # User Tuomo Valkonen # Date 1630071350 -10800 # Node ID ffd693c381f2a177d3ca6c8c4e308e7426fd36aa # Parent 62c62f451a410394aedde8056b0d68565f23e788 Julia stupid memory allocation bug workaround diff -r 62c62f451a41 -r ffd693c381f2 src/Util.jl --- 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