Optimizing matrix multiplication code and being compatible with autodiff

General Usage Performance
,
6

I see. My experience with Turing models has been that I can allocate an Array{T} and update its elements each iteration, with reversediff. But the same model fails with zygote.