Dyadic mapslices

[mschauer]

It would be a nice feature if mapslices could have a dyadic counterpart (like map has map(f, a, b)) which takes two matrices A and B and two slices and applies the two argument function f to pairs of slices of A and B.

Many operations can be expressed in that form, as a silly example left and right matrix multiplication by a diagonal matrix diagm(d) would correspond to mapslices(.*, A, 1, d, 1) respective mapslices(.*, A, 2, d, 1)

[mschauer]

I took a second look at mapslices, I can translate that into the dyadic case.

[mschauer]

For monadic mapslices, the images of the slices get stored along the dimensions sliced over. For dyadic mapslices one would rather break that convention (how should that work when slicing over different dimensions in both arguments...) and reorder dimensions so that the sliced dimensions come after the dimensions iterated over. Having a second version of mapslices with different behaviour in base is not so nice. So I close this for now.