@@ -694,28 +694,11 @@ mul!(C::AbstractMatrix, A::AbstractTriangular, adjB::Adjoint{<:Any,<:AbstractVec
694
694
mul! (C:: AbstractVecOrMat , A:: AbstractTriangular , adjB:: Adjoint{<:Any,<:AbstractVecOrMat} ) =
695
695
(B = adjB. parent; lmul! (A, adjoint! (C, B)))
696
696
697
- # The three methods for each op are neceesary to avoid ambiguities with definitions in matmul.jl
697
+ # The three methods are neceesary to avoid ambiguities with definitions in matmul.jl
698
698
mul! (C:: AbstractVector , A:: AbstractTriangular , B:: AbstractVector ) = lmul! (A, copyto! (C, B))
699
699
mul! (C:: AbstractMatrix , A:: AbstractTriangular , B:: AbstractVecOrMat ) = lmul! (A, copyto! (C, B))
700
700
mul! (C:: AbstractVecOrMat , A:: AbstractTriangular , B:: AbstractVecOrMat ) = lmul! (A, copyto! (C, B))
701
- # function mul!(C::AbstractVector, adjA::Adjoint{<:Any,<:AbstractTriangular}, B::AbstractVector)
702
- # return lmul!(adjA, copyto!(C, B))
703
- # end
704
- # function mul!(C::AbstractMatrix, adjA::Adjoint{<:Any,<:AbstractTriangular}, B::AbstractVecOrMat)
705
- # return lmul!(adjA, copyto!(C, B))
706
- # end
707
- # function mul!(C::AbstractVecOrMat, adjA::Adjoint{<:Any,<:AbstractTriangular}, B::AbstractVecOrMat)
708
- # return lmul!(adjA, copyto!(C, B))
709
- # end
710
- # function mul!(C::AbstractVector, transA::Transpose{<:Any,<:AbstractTriangular}, B::AbstractVector)
711
- # return lmul!(transA, copyto!(C, B))
712
- # end
713
- # function mul!(C::AbstractMatrix, transA::Transpose{<:Any,<:AbstractTriangular}, B::AbstractVecOrMat)
714
- # return lmul!(transA, copyto!(C, B))
715
- # end
716
- # function mul!(C::AbstractVecOrMat, transA::Transpose{<:Any,<:AbstractTriangular}, B::AbstractVecOrMat)
717
- # return lmul!(transA, copyto!(C, B))
718
- # end
701
+
719
702
@inline mul! (C:: AbstractMatrix , A:: AbstractTriangular , B:: Adjoint{<:Any,<:AbstractVecOrMat} , alpha:: Number , beta:: Number ) =
720
703
mul! (C, A, copy (B), alpha, beta)
721
704
@inline mul! (C:: AbstractMatrix , A:: AbstractTriangular , B:: Transpose{<:Any,<:AbstractVecOrMat} , alpha:: Number , beta:: Number ) =
@@ -1742,7 +1725,7 @@ end
1742
1725
powm (A:: LowerTriangular , p:: Real ) = copy (transpose (powm! (copy (transpose (A)), p:: Real )))
1743
1726
1744
1727
# Complex matrix logarithm for the upper triangular factor, see:
1745
- # Al-Mohy and Higham, "Improved inverse scaling and squaring algorithms for
1728
+ # Al-Mohy and Higham, "Improved inverse scaling and squaring algorithms for
1746
1729
# the matrix logarithm", SIAM J. Sci. Comput., 34(4), (2012), pp. C153–C169.
1747
1730
# Al-Mohy, Higham and Relton, "Computing the Frechet derivative of the matrix
1748
1731
# logarithm and estimating the condition number", SIAM J. Sci. Comput.,
0 commit comments