Added missing impl of output_size for ComposedFunction (#296)

* add `output_size` for `ComposedFunction`

* add tests for `output_size` for `ComposedFunction`

* bump patch version

* Update test/bijectors/stacked.jl

Co-authored-by: github-actions[bot] <41898282+github-actions[bot]@users.noreply.github.com>

---------

Co-authored-by: github-actions[bot] <41898282+github-actions[bot]@users.noreply.github.com>

Project.toml

@@ -1,6 +1,6 @@
 name = "Bijectors"
 uuid = "76274a88-744f-5084-9051-94815aaf08c4"
-version = "0.13.7"
+version = "0.13.8"
 
 [deps]
 ArgCheck = "dce04be8-c92d-5529-be00-80e4d2c0e197"

src/interface.jl

@@ -42,6 +42,7 @@ with_logabsdet_jacobian(f::Columnwise, x::AbstractMatrix) = (f(x), logabsdetjac(
 Returns the output size of `f` given the input size `sz`.
 """
 output_size(f, sz) = sz
+output_size(f::ComposedFunction, sz) = output_size(f.outer, output_size(f.inner, sz))
 
 """
     output_length(f, len::Int)

test/bijectors/stacked.jl

@@ -34,4 +34,44 @@ end
         @test y == [exp(1.0), 2.0]
         @test binv(y) == [1.0, 2.0, 0.0]
     end
+
+    @testset "composition" begin
+        # Composition with one dimension reduction.
+        b = Stacked((elementwise(exp), ProjectionBijector() ∘ identity), [1:1, 2:3])
+        binv = inverse(b)
+        x = [1.0, 2.0, 3.0]
+        y = b(x)
+        x_ = binv(y)
+
+        # Are the values of correct size?
+        @test size(y) == (2,)
+        @test size(x_) == (3,)
+        # Can we determine the sizes correctly?
+        @test Bijectors.output_size(b, size(x)) == (2,)
+        @test Bijectors.output_size(binv, size(y)) == (3,)
+
+        # Are values correct?
+        @test y == [exp(1.0), 2.0]
+        @test binv(y) == [1.0, 2.0, 0.0]
+
+        # Composition with two dimension reductions.
+        b = Stacked(
+            (elementwise(exp), ProjectionBijector() ∘ ProjectionBijector()), [1:1, 2:4]
+        )
+        binv = inverse(b)
+        x = [1.0, 2.0, 3.0, 4.0]
+        y = b(x)
+        x_ = binv(y)
+
+        # Are the values of correct size?
+        @test size(y) == (2,)
+        @test size(x_) == (4,)
+        # Can we determine the sizes correctly?
+        @test Bijectors.output_size(b, size(x)) == (2,)
+        @test Bijectors.output_size(binv, size(y)) == (4,)
+
+        # Are values correct?
+        @test y == [exp(1.0), 2.0]
+        @test binv(y) == [1.0, 2.0, 0.0, 0.0]
+    end
 end