¿Cuál es la forma recomendada de iterar una matriz sobre filas?

La solución que usted mismo enumeró, así como mapslices, ambos funcionan bien. Pero si por " recomendado "lo que realmente quieres decir es" alto rendimiento", entonces la mejor respuesta es: no iteres sobre filas.

El problema es que, dado que los arrays se almacenan en orden mayor de columna, para cualquier cosa que no sea una matriz pequeña, terminará con una pobre relación de aciertos de caché si recorre el array en orden mayor de fila.

Como se señala en un excelente entrada de blog , si desea resumir filas, su mejor apuesta es hacer algo como esto:

msum = zeros(eltype(m), size(m, 1))
for j = 1:size(m,2)
    for i = 1:size(m,1)
        msum[i] += m[i,j]
    end
end

Atravesamos tanto m como msum en su orden de almacenamiento nativo, por lo que cada vez que cargamos una línea de caché usamos todos los valores, lo que produce una relación de aciertos de caché de 1. Podría pensar ingenuamente que es mejor recorrerlo en orden de fila mayor y acumular el resultado a una variable tmp, pero en cualquier máquina moderna la falta de caché es mucho más cara que la búsqueda msum[i].

Muchos de los algoritmos internos de Julia que toman un region parámetro, como sum(m, 2), manejar esto para usted.