@@ -14,24 +14,24 @@ sizes_from_blocks(A::Diagonal, _) = (size.(A.diag, 1), size.(A.diag,2))
1414
1515function sizes_from_blocks (A:: Tridiagonal , _)
1616 for k = 1 : length (A. du)
17- size (A. du[k],1 ) == size (A. d[k],1 ) || throw (ArgumentError (" block sizes of upper diagonal inconsisent with diagonal" ))
18- size (A. du[k],2 ) == size (A. d[k+ 1 ],2 ) || throw (ArgumentError (" block sizes of upper diagonal inconsisent with diagonal" ))
19- size (A. dl[k],1 ) == size (A. d[k+ 1 ],1 ) || throw (ArgumentError (" block sizes of lower diagonal inconsisent with diagonal" ))
20- size (A. dl[k],2 ) == size (A. d[k],2 ) || throw (ArgumentError (" block sizes of lower diagonal inconsisent with diagonal" ))
17+ size (A. du[k],1 ) == size (A. d[k],1 ) || throw (ArgumentError (" block sizes of upper diagonal inconsistent with diagonal" ))
18+ size (A. du[k],2 ) == size (A. d[k+ 1 ],2 ) || throw (ArgumentError (" block sizes of upper diagonal inconsistent with diagonal" ))
19+ size (A. dl[k],1 ) == size (A. d[k+ 1 ],1 ) || throw (ArgumentError (" block sizes of lower diagonal inconsistent with diagonal" ))
20+ size (A. dl[k],2 ) == size (A. d[k],2 ) || throw (ArgumentError (" block sizes of lower diagonal inconsistent with diagonal" ))
2121 end
2222 (size .(A. d, 1 ), size .(A. d,2 ))
2323end
2424
2525function sizes_from_blocks (A:: Bidiagonal , _)
2626 if A. uplo == ' U'
2727 for k = 1 : length (A. ev)
28- size (A. ev[k],1 ) == size (A. dv[k],1 ) || throw (ArgumentError (" block sizes of upper diagonal inconsisent with diagonal" ))
29- size (A. ev[k],2 ) == size (A. dv[k+ 1 ],2 ) || throw (ArgumentError (" block sizes of upper diagonal inconsisent with diagonal" ))
28+ size (A. ev[k],1 ) == size (A. dv[k],1 ) || throw (ArgumentError (" block sizes of upper diagonal inconsistent with diagonal" ))
29+ size (A. ev[k],2 ) == size (A. dv[k+ 1 ],2 ) || throw (ArgumentError (" block sizes of upper diagonal inconsistent with diagonal" ))
3030 end
3131 else
3232 for k = 1 : length (A. ev)
33- size (A. ev[k],1 ) == size (A. dv[k+ 1 ],1 ) || throw (ArgumentError (" block sizes of lower diagonal inconsisent with diagonal" ))
34- size (A. ev[k],2 ) == size (A. dv[k],2 ) || throw (ArgumentError (" block sizes of lower diagonal inconsisent with diagonal" ))
33+ size (A. ev[k],1 ) == size (A. dv[k+ 1 ],1 ) || throw (ArgumentError (" block sizes of lower diagonal inconsistent with diagonal" ))
34+ size (A. ev[k],2 ) == size (A. dv[k],2 ) || throw (ArgumentError (" block sizes of lower diagonal inconsistent with diagonal" ))
3535 end
3636 end
3737 (size .(A. dv, 1 ), size .(A. dv,2 ))
@@ -89,4 +89,4 @@ for op in (:-, :+)
8989 mortar (Bidiagonal (broadcast ($ op, Ref (λ), A. blocks. dv), broadcast ($ op,A. blocks. ev), A. blocks. uplo))
9090 end
9191 end
92- end
92+ end
0 commit comments