Is there a shorter way to write `for m₁ in axes(M, 1), m₂ in axes(M, 2)` in the co

[tim.holy]

for m in CartesianIndices(M) and then use m[1] and m[2] in your code.