Back to empty array as sentinel

Author: timholy

src/abstract.jl

@@ -54,8 +54,8 @@ mutable struct LinearOperator{T, S, I <: Integer, F, Ft, Fct} <: AbstractLinearO
   nprod::I
   ntprod::I
   nctprod::I
-  Mv::Union{S, Nothing}
-  Mtu::Union{S, Nothing}
+  Mv::S
+  Mtu::S
 end
 
 function LinearOperator{T, S}(
@@ -81,8 +81,8 @@ function LinearOperator{T, S}(
     nprod,
     ntprod,
     nctprod,
-    nothing,
-    nothing,
+    S(undef, 0),
+    S(undef, 0),
   )
 end
 
@@ -161,7 +161,7 @@ No additional vectors are generated when using the 3-args `mul!`.
 """
 has_args5(op::AbstractLinearOperator) = get_nargs(op.prod!) == 4
 
-isallocated5(op::LinearOperator) = !(isnothing(op.Mv) || isnothing(op.Mtu))
+isallocated5(op::LinearOperator) = !(isempty(op.Mv) || isempty(op.Mtu))
 
 has_args5(op::AbstractMatrix) = true  # Needed for BlockDiagonalOperator
 

src/adjtrans.jl

@@ -112,7 +112,7 @@ function mul!(
     if hasmethod(ctprod!, Tuple{typeof(res), typeof(v), typeof(α), typeof(β)})
       return ctprod!(res, v, α, β)
     else
-      iszero(β) || (p.Mtu !== nothing) || allocate_vectors_args3!(p)
+      iszero(β) || !isempty(p.Mtu) || allocate_vectors_args3!(p)
       return prod3!(res, ctprod!, v, α, β, p.Mtu)
     end
   end
@@ -135,7 +135,7 @@ function mul!(
   if hasmethod(tprod!, Tuple{typeof(res), typeof(v), typeof(α), typeof(β)})
     tprod!(res, conj.(v), conj(α), conj(β))
   else
-    iszero(β) || (p.Mtu !== nothing) || allocate_vectors_args3!(p)
+    iszero(β) || !isempty(p.Mtu) || allocate_vectors_args3!(p)
     prod3!(res, tprod!, conj.(v), conj(α), conj(β), p.Mtu)
   end
   conj!(res)
@@ -179,7 +179,7 @@ function mul!(
     if hasmethod(tprod!, Tuple{typeof(res), typeof(v), typeof(α), typeof(β)})
       return tprod!(res, v, α, β)
     else
-      iszero(β) || (p.Mtu !== nothing) || allocate_vectors_args3!(p)
+      iszero(β) || !isempty(p.Mtu) || allocate_vectors_args3!(p)
       return prod3!(res, tprod!, v, α, β, p.Mtu)
     end
   end
@@ -202,7 +202,7 @@ function mul!(
   if hasmethod(ctprod!, Tuple{typeof(res), typeof(v), typeof(α), typeof(β)})
     ctprod!(res, conj.(v), conj(α), conj(β))
   else
-    iszero(β) || (p.Mtu !== nothing) || allocate_vectors_args3!(p)
+    iszero(β) || !isempty(p.Mtu) || allocate_vectors_args3!(p)
     prod3!(res, ctprod!, conj.(v), conj(α), conj(β), p.Mtu)
   end
   conj!(res)

src/operations.jl

@@ -26,7 +26,7 @@ function mul!(res::AbstractVector, op::AbstractLinearOperator{T}, v::AbstractVec
   if hasmethod(op.prod!, Tuple{typeof(res), typeof(v), typeof(α), typeof(β)})
     op.prod!(res, v, α, β)
   else
-    iszero(β) || (op.Mv !== nothing) || allocate_vectors_args3!(op)
+    iszero(β) || !isempty(op.Mv) || allocate_vectors_args3!(op)
     prod3!(res, op.prod!, v, α, β, op.Mv)
   end
 end