Is it possible to do a
mapreduce with multiple arrays while broadcasting over some dimensions? Something like this
U = rand(5, 5, 5) Δ = rand(1, 1, 5) mapreduce((u, Δx) -> u / Δx, max, U, Δ)
where I really want to find the maximum of
U[i, j, k] / Δ[k].
Note that the original poster on Slack cannot see your response here on Discourse. Consider transcribing the appropriate answer back to Slack, or pinging the poster here on Discourse so they can follow this thread.
(Original message ) (More Info)