module somemodule using BenchmarkTools struct MyType A::AbstractMatrix B::AbstractMatrix end myinc(A::AbstractMatrix, B::AbstractMatrix) = A .+= B mytypeinc(t::MyType, B::AbstractMatrix) = t.A .+= B mytypeinc2(t::MyType, B::AbstractMatrix) = myinc(t.A, B) function perf() H = [1.0 0 ;…

Why does this matrix increment allocate?

lmiq February 20, 2022, 11:37am 2

Probably because these are abstract containers. Use Matrix{T}, or a full parametric type.

1 Like