We read every piece of feedback, and take your input very seriously.
To see all available qualifiers, see our documentation.
1 parent 1e36ef1 commit dc17ad7Copy full SHA for dc17ad7
src/common/matrixproperties.jl
@@ -61,13 +61,7 @@ The `isapprox_kwargs` can be used to control the tolerances of the equality.
61
62
See also [`isisometric`](@ref) and [`MatrixAlgebraKit.is_left_isometric`](@ref).
63
""" is_right_isometric
64
-
65
-function is_right_isometric(A::AbstractMatrix; atol::Real = 0, rtol::Real = defaulttol(A), norm = LinearAlgebra.norm)
66
- P = A * A'
67
- nP = norm(P) # isapprox would use `rtol * max(norm(P), norm(I))`
68
- diagview(P) .-= 1
69
- return norm(P) <= max(atol, rtol * nP) # assume that the norm of I is `sqrt(n)`
70
-end
+is_right_isometric(A; kwargs...) = is_left_isometric(A'; kwargs...)
71
72
"""
73
ishermitian(A; isapprox_kwargs...)
0 commit comments