comparison: src/Util.jl
src/Util.jl
- changeset 28
- ffd693c381f2
- parent 19
- 63849571a046
- child 29
- 69242e8598eb
equal
deleted
inserted
replaced
| 27:62c62f451a41 | 28:ffd693c381f2 |
|---|---|
| 193 y′=reshape(y, (size(y, 1), prod(size(y)[2:end]))) | 193 y′=reshape(y, (size(y, 1), prod(size(y)[2:end]))) |
| 194 | 194 |
| 195 @inbounds @simd for i=1:size(y′, 2)# in CartesianIndices(size(y)[2:end]) | 195 @inbounds @simd for i=1:size(y′, 2)# in CartesianIndices(size(y)[2:end]) |
| 196 n² = norm₂²(@view(y′[:, i])) | 196 n² = norm₂²(@view(y′[:, i])) |
| 197 if n²>α² | 197 if n²>α² |
| 198 y′[:, i] .*= (α/√n²) | 198 @. @views y′[:, i] = y′[:, i]*(α/√n²) |
| 199 end | 199 end |
| 200 end | 200 end |
| 201 end | 201 end |
| 202 end | 202 end |
| 203 | 203 |
