operator.cpp

Go to the documentation of this file.

#include "operator/dense.hpp"
#include "operator/sparse.hpp"
 
#include <stdexcept>
 
namespace num::operators {
 
void DenseOp::apply(const Vector& x, Vector& y) const {
    if (x.size() != A_.cols()) {
        throw std::invalid_argument("DenseOp::apply: input dimension mismatch");
    }
    if (y.size() != A_.rows()) {
        y = Vector(A_.rows());
    }
    matvec(A_, x, y, b_);
}
 
void SparseOp::apply(const Vector& x, Vector& y) const {
    if (x.size() != A_.n_cols()) {
        throw std::invalid_argument("SparseOp::apply: input dimension mismatch");
    }
    if (y.size() != A_.n_rows()) {
        y = Vector(A_.n_rows());
    }
    sparse_matvec(A_, x, y);
}
 
} // namespace num::operators