Consider the following code doing a simple computation: function test_vectorized(a, a_old) @. a[:] = a_old[:] * 2 end function test_vectorized2(a, a_old) @. a = a_old * 2 end function test_devectorized(a, a_old) nx = size(a, 1) @inbounds @simd for i = 1:nx a[i] = a_old[i…

@views function test_vectorized3(a, a_old) @. a[:] = a_old[:] * 2 end julia> @time test_vectorized3(a2, a); 0.126813 seconds (5 allocations: 208 bytes) I think views will also work for your stencils.

Performance drop for explicit indexing in fused loops

General Usage

dpsanders June 26, 2017, 1:38pm 5

Don’t benchmark in global scope.

CSV Reading (rewrite in C?)

Topic		Replies	Views
Comprehensions and vectorization New to Julia	13	4789	September 20, 2017
Large performance regression of vectorized functions on master General Usage	4	609	July 25, 2018
Performance difference between two code New to Julia	3	778	January 18, 2017
Broadcast vs. scalar loop, can Julia vectorize better? Internals & Design	8	1889	February 15, 2020
What is idiomatic way of vectorizing longer operations in Julia? Performance performance , loopvectorization	5	409	September 21, 2023

Performance drop for explicit indexing in fused loops

Related topics