sparse.hpp

Go to the documentation of this file.

/// @file operator/sparse.hpp
/// @brief SparseMatrix adapter for the operator protocol.
#pragma once
 
#include "linalg/sparse/sparse.hpp"
#include "operator/concepts.hpp"
 
namespace num::operators {
 
/// @brief Adapt a SparseMatrix to the operator protocol.
struct SparseOp final {
    explicit SparseOp(const SparseMatrix& A)
        : A_(A) {}
 
    void apply(const Vector& x, Vector& y) const;
    [[nodiscard]] idx rows() const noexcept { return A_.n_rows(); }
    [[nodiscard]] idx cols() const noexcept { return A_.n_cols(); }
 
  private:
    const SparseMatrix& A_;
};
 
static_assert(LinearOperator<SparseOp>);
 
} // namespace num::operators