We read every piece of feedback, and take your input very seriously.
To see all available qualifiers, see our documentation.
There was an error while loading. Please reload this page.
1 parent 098966a commit c4ef85eCopy full SHA for c4ef85e
src/symmetric.jl
@@ -841,10 +841,11 @@ end
841
^(A::Hermitian, p::Integer) = sympow(A, p)
842
function sympow(A, p::Integer)
843
if p < 0
844
- return wrappertype(A)(Base.power_by_squaring(inv(A), -p))
+ retmat = Base.power_by_squaring(inv(A), -p)
845
else
846
- return wrappertype(A)(Base.power_by_squaring(A, p))
+ retmat = Base.power_by_squaring(A, p)
847
end
848
+ return wrappertype(A)(retmat)
849
850
function ^(A::SelfAdjoint, p::Real)
851
isinteger(p) && return integerpow(A, p)
0 commit comments