Multi-arguments mapslices

[fabiangans]

I think you can solve this using LazyArrays.jl:

using LazyArrays
reduce(+, @~(A.*B), dims=1,init=0.0)

For reference, I am linking a few related topics:

Is there something like broadcast_mapreduce?

Just in case someone stumbles upon this thread and to answer my own question. Today I thought about the problem again and I found a nice and fast way to do this by using the mutating broadcast! on a custom AbstractArray which does the reduction upon calling setindex! Here is my final implementation of broadcast_reduce which has no allocations and is faster than the allocating version: #Define abstract Array which consumes values on setindex! type ArrayReducer{T,N,F} <: AbstractArray{T,N} v:…

Non-allocating reduced broadcast?

There is https://github.com/JuliaLang/julia/pull/31020 which is blocked by https://github.com/JuliaLang/julia/issues/19198

Is it possible to do a `mapreduce` with multiple arrays while broadcasting over so

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 ) (…