is_unitary Function

public pure function is_unitary(A, rtol) result(res)

Uses

Determine whether a matrix is unitary w.r.t the Frobenius norm

Arguments

Type IntentOptional Attributes Name
class(*), intent(in) :: A(:,:)

The matrix
real(kind=dp), intent(in), optional :: rtol

The optional relative tolerance, default 1e-10

Return Value logical