https://github.com/JuliaLang/julia/pull/33130

README.md

@@ -114,6 +114,8 @@ Currently, the `@compat` macro supports the following syntaxes:
 
 ## New functions, macros, and methods
 
+* A new `dropdims(f, args...; dims, kwargs...)` method computes the reduction `f` over the region described by `dims` and then drops those dimensions from the result ([#33130]).
+
 * `eachrow`, `eachcol`, and `eachslice` to iterate over first, second, or given dimension
   of an array ([#29749]).
 

src/Compat.jl

@@ -1866,6 +1866,11 @@ if v"0.7.0" <= VERSION < v"1.1.0-DEV.594"
     Base.merge(a::NamedTuple) = a
 end
 
+# https://github.com/JuliaLang/julia/pull/33130
+if VERSION < v"1.4.0-DEV.79" # not in fact merged yet!
+    Base.dropdims(f, args...; dims, kwargs...) = _dropdims(f(args...; kwargs..., dims=dims), dims)
+end
+
 include("deprecated.jl")
 
 end # module Compat

test/runtests.jl

@@ -38,7 +38,7 @@ if VERSION >= v"0.7"
         @test collect(eachrow(M)) == collect(eachslice(M, dims = 1)) == [[1, 2, 3], [4, 5, 6], [7, 8, 9]]
         @test collect(eachcol(M)) == collect(eachslice(M, dims = 2)) == [[1, 4, 7], [2, 5, 8], [3, 6, 9]]
         @test_throws DimensionMismatch eachslice(M, dims = 4)
-        
+
         # Higher-dimensional case
         M = reshape([(1:16)...], 2, 2, 2, 2)
         @test_throws MethodError collect(eachrow(M))
@@ -1499,4 +1499,37 @@ end
     @test merge((a=1,), (b=2,), (c=3,)) == (a=1,b=2,c=3)
 end
 
+# https://github.com/JuliaLang/julia/pull/33130
+@testset "dropdims with reductions" begin # issue #16606
+     @test (@inferred(dropdims(sum, a, dims=1)) ==
+            @inferred(dropdims(sum, a, dims=(1,))) ==
+            reshape(sum(a, dims=1), (1, 8, 8, 1)))
+     @test (@inferred(dropdims(sum, a, dims=3)) ==
+            @inferred(dropdims(sum, a, dims=(3,))) ==
+            reshape(sum(a, dims=3), (1, 1, 8, 1)))
+     @test (@inferred(dropdims(sum, a, dims=4)) ==
+            @inferred(dropdims(sum, a, dims=(4,))) ==
+            reshape(sum(a, dims=4), (1, 1, 8, 1)))
+     @test (@inferred(dropdims(sum, a, dims=(1, 5))) ==
+            dropdims(sum, a, dims=(5, 1)) ==
+            reshape(sum(a, dims=(5, 1)), (1, 8, 8)))
+     @test (@inferred(dropdims(sum, a, dims=(1, 2, 5))) ==
+            dropdims(sum, a, dims=(5, 2, 1)) ==
+            reshape(sum(a, dims=(5, 2, 1)), (8, 8)))
+     @test (@inferred(dropdims(sum, abs2, a, dims=1)) ==
+            @inferred(dropdims(sum, abs2, a, dims=(1,))) ==
+            reshape(sum(abs2, a, dims=1), (1, 8, 8, 1)))
+     _sumplus(x; dims, plus) = sum(x; dims=dims) .+ plus  # reduction with keywords
+     @test (@inferred(dropdims(_sumplus, a, dims=4, plus=1)) ==
+            @inferred(dropdims(_sumplus, a, dims=(4,), plus=1)) ==
+            reshape(sum(a, dims=4) .+ 1, (1, 1, 8, 1)))
+     @test_throws UndefKeywordError dropdims(sum, a)
+     @test_throws UndefKeywordError dropdims(sum, a, 1)
+     @test_throws ArgumentError dropdims(sum, a, dims=0)
+     @test_throws ArgumentError dropdims(sum, a, dims=(1, 1))
+     @test_throws ArgumentError dropdims(sum, a, dims=(1, 2, 1))
+     @test_throws ArgumentError dropdims(sum, a, dims=(1, 1, 2))
+     @test_throws ArgumentError dropdims(sum, a, dims=6)
+end
+
 nothing