tests

Author: vpuri3

test/basic.jl

@@ -27,6 +27,8 @@ K = 12
     β = rand()
     Id = IdentityOperator{N}()
 
+    @test issquare(Id)
+    @test islinear(Id)
     @test IdentityOperator(u) isa IdentityOperator{N}
     @test one(A) isa IdentityOperator{N}
     @test convert(AbstractMatrix, Id) == Matrix(I, N, N)
@@ -61,6 +63,8 @@ end
     β = rand()
     Z = NullOperator{N}()
 
+    @test issquare(Z)
+    @test islinear(Z)
     @test NullOperator(u) isa NullOperator{N}
     @test zero(A) isa NullOperator{N}
     @test convert(AbstractMatrix, Z) == zeros(size(Z))
@@ -99,6 +103,8 @@ end
     op = ScaledOperator(α, MatrixOperator(A))
 
     @test op isa ScaledOperator
+    @test issquare(op)
+    @test islinear(op)
 
     @test α * A * u ≈ op * u
     @test (β * op) * u ≈ β * α * A * u
@@ -164,10 +170,14 @@ end
 
     @test op isa ComposedOperator
     @test *(op.ops...) isa ComposedOperator
+    @test issquare(op)
+    @test islinear(op)
 
     opF = factorize(op)
 
     @test opF isa ComposedOperator
+    @test issquare(opF)
+    @test islinear(opF)
 
     @test ABCmulu ≈ op * u
     @test ABCdivu ≈ op \ u ≈ opF \ u
@@ -202,10 +212,14 @@ end
 
 @testset "Adjoint, Transpose" begin
 
-    for (op, LType, VType) in (
-                               (adjoint,   AdjointOperator,    AbstractAdjointVecOrMat   ),
-                               (transpose, TransposedOperator, AbstractTransposedVecOrMat),
-                              )
+    for (op,
+         LType,
+         VType
+        ) in (
+              (adjoint,   AdjointOperator,    AbstractAdjointVecOrMat   ),
+              (transpose, TransposedOperator, AbstractTransposedVecOrMat),
+             )
+
         A = rand(N,N)
         D = Bidiagonal(rand(N,N), :L)
         u = rand(N,K)
@@ -215,7 +229,13 @@ end
         b = rand()
 
         At = op(A)
-        Dt = op(A)
+        Dt = op(D)
+
+        @test issquare(At)
+        @test issquare(Dt)
+
+        @test islinear(At)
+        @test islinear(Dt)
 
         AA = MatrixOperator(A)
         DD = MatrixOperator(D)
@@ -247,6 +267,9 @@ end
 
     Di = cache_operator(Di, u)
 
+    @test issquare(Di)
+    @test islinear(Di)
+
     @test Di * u ≈ u ./ s
     v=rand(N); @test mul!(v, Di, u) ≈ u ./ s
     v=rand(N); w=copy(v); @test mul!(v, Di, u, α, β) ≈ α *(u ./ s) + β*w

test/func.jl

@@ -29,6 +29,7 @@ K = 12
 
                            op_inverse=f1i,
 
+                           islinear=true,
                            opnorm=true,
                            issymmetric=true,
                            ishermitian=true,
@@ -40,11 +41,17 @@ K = 12
 
                            op_inverse=f2i,
 
+                           islinear=true,
                            opnorm=true,
                            issymmetric=true,
                            ishermitian=true,
                            isposdef=true,
                           )
+    @test issquare(op1)
+    @test issquare(op2)
+
+    @test islinear(op1)
+    @test islinear(op2)
 
     @test op1' === op1
 

test/matrix.jl

@@ -24,6 +24,9 @@ K = 19
     @test AA  isa MatrixOperator
     @test AAt isa MatrixOperator
 
+    @test issquare(AA)
+    @test islinear(AA)
+
     FF  = factorize(AA)
     FFt = FF'
 
@@ -73,6 +76,9 @@ end
                          update_func= (diag,u,p,t) -> (diag .= p*t; nothing)
                         )
 
+    @test issquare(D)
+    @test islinear(D)
+
     ans = Diagonal(p*t) * u
     @test D(u,p,t) ≈ ans
     v=copy(u); @test D(v,u,p,t) ≈ ans
@@ -86,6 +92,9 @@ end
 
     L = DiagonalOperator(d)
 
+    @test issquare(L)
+    @test islinear(L)
+
     @test L * u ≈ d .* u
     v=rand(N,K); @test mul!(v, L, u) ≈ d .* u
     v=rand(N,K); w=copy(v); @test mul!(v, L, u, α, β) ≈ α*(d .* u) + β*w
@@ -105,17 +114,25 @@ end
     β = rand()
 
     L = AffineOperator(MatrixOperator(A), MatrixOperator(B), b)
+    @test issquare(L)
+    @test !islinear(L)
 
     @test L * u ≈ A * u + B*b
     v=rand(N,K); @test mul!(v, L, u) ≈ A*u + B*b
     v=rand(N,K); w=copy(v); @test mul!(v, L, u, α, β) ≈ α*(A*u + B*b) + β*w
 
     L = AffineOperator(MatrixOperator(D), MatrixOperator(B), b)
+    @test issquare(L)
+    @test !islinear(L)
+
     @test L \ u ≈ D \ (u - B * b)
     v=rand(N,K); @test ldiv!(v, L, u) ≈ D \ (u-B*b)
     v=copy(u); @test ldiv!(L, u) ≈ D \ (v-B*b)
 
     L = AddVector(b)
+    @test issquare(L)
+    @test !islinear(L)
+
     @test L * u ≈ u + b
     @test L \ u ≈ u - b
     v=rand(N,K); @test mul!(v, L, u) ≈  u + b
@@ -124,6 +141,9 @@ end
     v=copy(u); @test ldiv!(L, u) ≈  v - b
 
     L = AddVector(MatrixOperator(B), b)
+    @test issquare(L)
+    @test !islinear(L)
+
     @test L * u ≈ u + B * b
     @test L \ u ≈ u - B * b
     v=rand(N,K); @test mul!(v, L, u) ≈  u + B * b
@@ -188,6 +208,17 @@ for square in [false, true] #for K in [1, K]
     @test opAB  isa TensorProductOperator
     @test opABC isa TensorProductOperator
 
+    if square
+        @test issquare(opAB)
+        @test issquare(opABC)
+    else
+        @test !issquare(opAB)
+        @test !issquare(opABC)
+    end
+
+    @test islinear(opAB)
+    @test islinear(opABC)
+
     @test AB ≈ convert(AbstractMatrix, opAB)
     @test ABC ≈ convert(AbstractMatrix, opABC)
 

test/scalar.jl

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

test/total.jl

@@ -27,6 +27,8 @@ K = 12
                            op_adjoint = (du,u,p,t) -> ldiv!(du, tr, u),
                            op_inverse = (du,u,p,t) -> ldiv!(du, tr, u),
                            op_adjoint_inverse = (du,u,p,t) -> ldiv!(du, tr, u),
+
+                           islinear=true,
                           )
 
     # derivative test