concepts.hpp

Go to the documentation of this file.

/// @file operator/concepts.hpp
/// @brief Concepts for operator-oriented numerical algorithms.
#pragma once
 
#include "core/concepts.hpp"
#include "core/vector.hpp"
#include <concepts>
 
namespace num::operators {
 
template<class V>
concept VectorLike = num::VectorLike<V>;
 
template<class V>
concept MutableVectorLike = num::MutableVectorLike<V>;
 
/// @brief Compile-time contract for y = A*x.
template<class Op, class X = Vector, class Y = Vector>
concept LinearOperator =
    VectorLike<X> && MutableVectorLike<Y> && requires(const Op& A, const X& x, Y& y) {
        { A.rows() } -> std::convertible_to<idx>;
        { A.cols() } -> std::convertible_to<idx>;
        { A.apply(x, y) };
    };
 
} // namespace num::operators