DiagonalOperator.h

Go to the documentation of this file.

#ifndef DIAGONALOPERATOR_H
#define DIAGONALOPERATOR_H
 
#include "MUQ/Modeling/LinearAlgebra/LinearOperator.h"
 
 
namespace muq
{
namespace Modeling
{
 
    class DiagonalOperator : public LinearOperator
    {
    public:
        DiagonalOperator(Eigen::VectorXd const& diagIn);
 
        virtual ~DiagonalOperator(){};
        
        virtual Eigen::MatrixXd Apply(Eigen::Ref<const Eigen::MatrixXd> const& x) override;
  
        virtual Eigen::MatrixXd ApplyTranspose(Eigen::Ref<const Eigen::MatrixXd> const& x) override{return Apply(x);};
 
        virtual Eigen::MatrixXd GetMatrix() override{return diag.asDiagonal();};
        
    private:
        const Eigen::VectorXd diag;
 
    }; // class SumOperator
 
 
} // namespace Modeling
} // namespace muq
 
#endif // #ifndef DIAGONALOPERATOR_H