DiagonalOperator.cpp

Go to the documentation of this file.

#include "MUQ/Modeling/LinearAlgebra/DiagonalOperator.h"
 
using namespace muq::Modeling;
 
 
DiagonalOperator::DiagonalOperator(Eigen::VectorXd const& diagIn) : LinearOperator(diagIn.rows(), diagIn.rows()),
                                                                    diag(diagIn)
{}
 
Eigen::MatrixXd DiagonalOperator::Apply(Eigen::Ref<const Eigen::MatrixXd> const& x)
{
    return diag.asDiagonal()*x;
}