diagzero in tridiagonal indexing

Author: jishnub

src/tridiag.jl

@@ -470,7 +470,7 @@ end
     elseif i + 1 == j
         return @inbounds A.ev[i]
     else
-        return zero(T)
+        return diagzero(A, i, j)
     end
 end
 
@@ -714,7 +714,7 @@ end
     elseif i + 1 == j
         return @inbounds A.du[i]
     else
-        return zero(T)
+        return diagzero(A, i, j)
     end
 end
 
@@ -727,7 +727,7 @@ end
     elseif b.band == 1
         return @inbounds A.du[b.index]
     else
-        return zero(T)
+        return diagzero(A, Tuple(_cartinds(b))...)
     end
 end