tests

Author: vpuri3

test/basic.jl

@@ -31,6 +31,7 @@ K = 12
     @test one(A) isa IdentityOperator{N}
     @test convert(AbstractMatrix, Id) == Matrix(I, N, N)
 
+    @test iscached(Id)
     @test size(Id) == (N, N)
     @test Id' isa IdentityOperator{N}
 
@@ -65,6 +66,7 @@ end
     @test zero(A) isa NullOperator{N}
     @test convert(AbstractMatrix, Z) == zeros(size(Z))
 
+    @test iscached(Z)
     @test size(Z) == (N, N)
     @test Z' isa NullOperator{N}
 
@@ -99,13 +101,15 @@ end
     op = ScaledOperator(α, MatrixOperator(A))
 
     @test op isa ScaledOperator
+    @test iscached(op)
 
     @test α * A * u ≈ op * u
     @test (β * op) * u ≈ β * α * A * u
 
     opF = factorize(op)
 
     @test opF isa ScaledOperator
+    @test iscached(opF)
 
     @test α * A  ≈ convert(AbstractMatrix, op) ≈ convert(AbstractMatrix, opF)
 
@@ -145,6 +149,8 @@ end
     end
 
     op = AddedOperator(A, B)
+    @test iscached(op)
+
     v=rand(N,K); @test mul!(v, op, u) ≈ (A+B) * u
     v=rand(N,K); w=copy(v); @test mul!(v, op, u, α, β) ≈ α*(A+B)*u + β*w
 end
@@ -172,7 +178,10 @@ end
     @test ABCmulu ≈ op * u
     @test ABCdivu ≈ op \ u ≈ opF \ u
 
+    @test !iscached(op)
     op = cache_operator(op, u)
+    @test iscached(op)
+
     v=rand(N,K); @test mul!(v, op, u) ≈ ABCmulu
     v=rand(N,K); w=copy(v); @test mul!(v, op, u, α, β) ≈ α*ABCmulu + β*w
 
@@ -181,7 +190,9 @@ end
     C = rand(N) |> Diagonal
 
     op = ∘(MatrixOperator.((A, B, C))...)
+    @test !iscached(op)
     op = cache_operator(op, u)
+    @test iscached(op)
     v=rand(N,K); @test ldiv!(v, op, u) ≈ (A * B * C) \ u
     v=copy(u);   @test ldiv!(op, u)    ≈ (A * B * C) \ v
 
@@ -195,7 +206,9 @@ end
     @test_throws MethodError inner_op * u
     # We can now test that caching does not rely on matmul
     op = inner_op * factorize(MatrixOperator(rand(N, N)))
+    @test !iscached(op)
     @test_nowarn op = cache_operator(op, rand(N)) 
+    @test iscached(op)
     u = rand(N)
     @test ldiv!(rand(N), op, u) ≈ op \ u
 end
@@ -245,7 +258,9 @@ end
     α  = rand()
     β  = rand()
 
+    @test !iscached(Di)
     Di = cache_operator(Di, u)
+    @test iscached(Di)
 
     @test Di * u ≈ u ./ s
     v=rand(N); @test mul!(v, Di, u) ≈ u ./ s

test/func.jl

@@ -62,7 +62,8 @@ K = 12
     @test has_ldiv(op2)
     @test has_ldiv!(op2)
 
-    op2 = cache_operator(op2, u)
+    @test iscached(op1)
+    @test iscached(op2)
 
     v = rand(N,K); @test A * u ≈ op1 * u ≈ mul!(v, op2, u)
     v = rand(N,K); @test A * u ≈ op1(u,p,t) ≈ op2(v,u,p,t)

test/matrix.jl

@@ -207,12 +207,24 @@ for square in [false, true] #for K in [1, K]
     @test AB  \ v2 ≈ opAB  \ v2 ≈ opAB_F  \ v2
     @test ABC \ v3 ≈ opABC \ v3 ≈ opABC_F \ v3
 
+    @test !iscached(opAB)
+    @test !iscached(opABC)
+
+    @test !iscached(opAB_F)
+    @test !iscached(opABC_F)
+
     opAB  = cache_operator(opAB,  u2)
     opABC = cache_operator(opABC, u3)
 
     opAB_F  = cache_operator(opAB_F,  u2)
     opABC_F = cache_operator(opABC_F, u3)
 
+    @test iscached(opAB)
+    @test iscached(opABC)
+
+    @test iscached(opAB_F)
+    @test iscached(opABC_F)
+
     N2 = n1*n2
     N3 = n1*n2*n3
     M2 = m1*m2

test/scalar.jl

@@ -14,6 +14,7 @@ K = 12
     u = rand(N,K)
 
     @test α isa ScalarOperator
+    @test iscached(α)
     @test convert(Number, α) isa Number
     @test convert(ScalarOperator, a) isa ScalarOperator