sbplib_julia: test/testLazyTensors.jl comparison

Start implementing TensorComposition

author	Jonatan Werpers <jonatan@werpers.com>
date	Thu, 15 Oct 2020 22:05:22 +0200
parents	3b4b1758a8ad
children	814865d40f48

comparison

equal deleted inserted replaced

-:26e186b565b3
+:d94891b8dfca
 @test_throws BoundsError LazyFunctionArray((i,j)->i*j, (3,2))[4,2]
 @test_throws BoundsError LazyFunctionArray((i,j)->i*j, (3,2))[2,3]
 end
+@testset "TensorMappingComposition" begin
+A = rand(2,3)
+B = rand(3,4)
+Ã = LazyLinearMap(A, (1,), (2,))
+B̃ = LazyLinearMap(B, (1,), (2,))
+@test Ã∘B̃ isa TensorMappingComposition
+@test range_size(Ã∘B̃) == (2,)
+@test domain_size(Ã∘B̃) == (4,)
+# @test_throws DimensionMismatch B̃∘Ã
+v = rand(4)
+@test Ã∘B̃*v ≈ A*B*v
+v = rand(2)
+@test_broken (Ã∘B̃)'*v ≈ B'*A'*v
+end
 @testset "LazyLinearMap" begin
 # Test a standard matrix-vector product
 # mapping vectors of size 4 to vectors of size 3.
 A = rand(3,4)
 Ã = LazyLinearMap(A, (1,), (2,))

Mercurial > repos > public > sbplib_julia