Specialize for adjoint/transpose

Author: jishnub

src/adjtrans.jl

@@ -572,3 +572,8 @@ _vecadjoint(A::Base.ReshapedArray{<:Any,1,<:AdjointAbsVec}) = adjoint(parent(A))
 
 diagview(A::Transpose, k::Integer = 0) = _vectranspose(diagview(parent(A), -k))
 diagview(A::Adjoint, k::Integer = 0) = _vecadjoint(diagview(parent(A), -k))
+
+function fillband!(A::AdjOrTrans, v, k1, k2)
+    fillband!(parent(A), wrapperop(A)(v), -k2, -k1)
+    return A
+end

test/adjtrans.jl

@@ -798,4 +798,16 @@ end
     end
 end
 
+@testset "fillband!" begin
+    for A in (rand(4, 4), rand(ComplexF64,4,4))
+        B = similar(A)
+        for op in (adjoint, transpose), k in -3:3
+            B .= op(A)
+            LinearAlgebra.fillband!(op(A), 1, k, k)
+            LinearAlgebra.fillband!(B, 1, k, k)
+            @test op(A) == B
+        end
+    end
+end
+
 end # module TestAdjointTranspose