Hello there, I’m trying to figure out how to do matrix operations for 2D slices of a 3D array whose shape is (3,3,N) where N is a large number (e.g. N=10^6). I’d like to write a vectorized code of C=A+B’ using loop fusion as explained in More Dots: Syntactic Loop Fusion in Julia for the first 1024 …

Since you’re working with such small slices, the cost of creating all those views (and doing the associated bookkeeping on transposes) is hard to amortize. Adding @inbounds speeds up the expanded loop even further, but if you can use StaticArrays, I think they give you the syntax you’re looking for: …

Efficient vectorized matrix operations for 2D matrix slices of a 3D array?

New to Julia

GregVernon May 1, 2019, 5:47pm 8

I asked a question a while ago, which has some similar themes to your question. You might find some of the answers I got to be useful.

1 Like