check_symmetric Function

    pure logical function check_symmetric(matrix)
        ! function to check if a given matrix is symmetric
        ! for a square matrix!
        real(dp), intent(in) :: matrix(:, :)
        debug_function_name("check_symmetric")
        ASSERT(size(matrix, 1) == size(matrix, 2))
        check_symmetric = near_zero(sum(abs(matrix - transpose(matrix))))
    end function check_symmetric