is_upper_triangular.hpp

Go to the documentation of this file.

// SPDX-License-Identifier: BSD-2-Clause
 
#ifndef pRC_CORE_TENSOR_OPERATOR_FUNCTIONS_IS_UPPER_TRIANGULAR_H
#define pRC_CORE_TENSOR_OPERATOR_FUNCTIONS_IS_UPPER_TRIANGULAR_H
 
#include <prc/core/functors/upper_triangular.hpp>
#include <prc/core/tensor/declarations.hpp>
#include <prc/core/value/limits.hpp>
 
namespace pRC
{
    template<class X, IsTensorish R = RemoveReference<X>, IsValue TT = Value<R>>
        requires IsInvocable<UpperTriangular, X>
    static inline constexpr auto isUpperTriangular(X &&a,
        TT const &tolerance = NumericLimits<TT>::tolerance())
    {
        return isApprox(upperTriangular(a), a, tolerance);
    }
}
#endif // pRC_CORE_TENSOR_OPERATOR_FUNCTIONS_IS_UPPER_TRIANGULAR_H