@@ -801,38 +801,38 @@ function _triscale!(A::LowerOrUnitLowerTriangular, c::Number, B::UnitLowerTriang
801
801
return A
802
802
end
803
803
804
- function _trirdiv! (A:: UpperTriangular , B:: UpperTriangular , c:: Number )
804
+ function _trirdiv! (A:: UpperTriangular , B:: UpperOrUnitUpperTriangular , c:: Number )
805
805
checksize1 (A, B)
806
806
for j in axes (B,2 )
807
807
for i in firstindex (B,1 ): j
808
- @inbounds A. data [i, j] = B. data [i, j] / c
808
+ @inbounds A[i, j] = B[i, j] / c
809
809
end
810
810
end
811
811
return A
812
812
end
813
- function _trirdiv! (A:: LowerTriangular , B:: LowerTriangular , c:: Number )
813
+ function _trirdiv! (A:: LowerTriangular , B:: LowerOrUnitLowerTriangular , c:: Number )
814
814
checksize1 (A, B)
815
815
for j in axes (B,2 )
816
816
for i in j: lastindex (B,1 )
817
- @inbounds A. data [i, j] = B. data [i, j] / c
817
+ @inbounds A[i, j] = B[i, j] / c
818
818
end
819
819
end
820
820
return A
821
821
end
822
- function _trildiv! (A:: UpperTriangular , c:: Number , B:: UpperTriangular )
822
+ function _trildiv! (A:: UpperTriangular , c:: Number , B:: UpperOrUnitUpperTriangular )
823
823
checksize1 (A, B)
824
824
for j in axes (B,2 )
825
825
for i in firstindex (B,1 ): j
826
- @inbounds A. data [i, j] = c \ B. data [i, j]
826
+ @inbounds A[i, j] = c \ B[i, j]
827
827
end
828
828
end
829
829
return A
830
830
end
831
- function _trildiv! (A:: LowerTriangular , c:: Number , B:: LowerTriangular )
831
+ function _trildiv! (A:: LowerTriangular , c:: Number , B:: LowerOrUnitLowerTriangular )
832
832
checksize1 (A, B)
833
833
for j in axes (B,2 )
834
834
for i in j: lastindex (B,1 )
835
- @inbounds A. data [i, j] = c \ B. data [i, j]
835
+ @inbounds A[i, j] = c \ B[i, j]
836
836
end
837
837
end
838
838
return A
0 commit comments