fix/restore matrixproperties

Author: Jutho

src/common/matrixproperties.jl

@@ -144,79 +144,3 @@ function _ishermitian_exact_offdiag(Al, Au, ::Val{anti}) where {anti}
     end
     return true
 end
-
-"""
-    ishermitian(A; isapprox_kwargs...)
-
-Test whether a linear map is Hermitian, i.e. `A = A'`.
-The `isapprox_kwargs` can be used to control the tolerances of the equality.
-"""
-function ishermitian(A; atol::Real = 0, rtol::Real = 0, norm = LinearAlgebra.norm, kwargs...)
-    if iszero(atol) && iszero(rtol)
-        return ishermitian_exact(A; kwargs...)
-    else
-        return 2 * norm(project_antihermitian(A; kwargs...)) ≤ max(atol, rtol * norm(A))
-    end
-end
-function ishermitian_exact(A)
-    return A == A'
-end
-function ishermitian_exact(A::AbstractMatrix; kwargs...)
-    return _ishermitian_exact(A, Val(false); kwargs...)
-end
-
-"""
-    isantihermitian(A; isapprox_kwargs...)
-
-Test whether a linear map is anti-Hermitian, i.e. `A = -A'`.
-The `isapprox_kwargs` can be used to control the tolerances of the equality.
-"""
-function isantihermitian(A; atol::Real = 0, rtol::Real = 0, norm = LinearAlgebra.norm, kwargs...)
-    if iszero(atol) && iszero(rtol)
-        return isantihermitian_exact(A; kwargs...)
-    else
-        return 2 * norm(project_hermitian(A; kwargs...)) ≤ max(atol, rtol * norm(A))
-    end
-end
-function isantihermitian_exact(A)
-    return A == -A'
-end
-function isantihermitian_exact(A::AbstractMatrix; kwargs...)
-    return _ishermitian_exact(A, Val(true); kwargs...)
-end
-
-# block implementation of exact checks
-function _ishermitian_exact(A::AbstractMatrix, anti::Val; blocksize = 32)
-    n = size(A, 1)
-    for j in 1:blocksize:n
-        jb = min(blocksize, n - j + 1)
-        _ishermitian_exact_diag(view(A, j:(j + jb - 1), j:(j + jb - 1)), anti) || return false
-        for i in 1:blocksize:(j - 1)
-            ib = blocksize
-            _ishermitian_exact_offdiag(
-                view(A, i:(i + ib - 1), j:(j + jb - 1)),
-                view(A, j:(j + jb - 1), i:(i + ib - 1)),
-                anti
-            ) || return false
-        end
-    end
-    return true
-end
-function _ishermitian_exact_diag(A, ::Val{anti}) where {anti}
-    n = size(A, 1)
-    @inbounds for j in 1:n
-        @simd for i in 1:j
-            A[i, j] == (anti ? -adjoint(A[j, i]) : adjoint(A[j, i])) || return false
-        end
-    end
-    return true
-end
-function _ishermitian_exact_offdiag(Al, Au, ::Val{anti}) where {anti}
-    m, n = size(Al) # == reverse(size(Al))
-    @inbounds for j in 1:n
-        @simd for i in 1:m
-            Al[i, j] == (anti ? -adjoint(Au[j, i]) : adjoint(Au[j, i])) || return false
-        end
-    end
-    return true
-end