CompanionMatrix.h

Go to the documentation of this file.

#ifndef COMPANIONMATRIX_H
#define COMPANIONMATRIX_H
 
#include "MUQ/Modeling/LinearAlgebra/LinearOperator.h"
#include <Eigen/Core>
 
namespace muq
{
namespace Modeling
{
 
    class CompanionMatrix : public LinearOperator
    {
    public:
        CompanionMatrix(Eigen::VectorXd const& lastRowIn) : LinearOperator(lastRowIn.size(), lastRowIn.size()), lastRow(lastRowIn){};
        
        virtual Eigen::MatrixXd Apply(Eigen::Ref<const Eigen::MatrixXd> const& x) override;
 
        virtual Eigen::MatrixXd ApplyTranspose(Eigen::Ref<const Eigen::MatrixXd> const& x) override;
 
        virtual Eigen::MatrixXd GetMatrix() override;
        
    private:
        Eigen::VectorXd lastRow;
    };
 
} // namespace Modeling
} // namespace muq
 
#endif // #ifndef COMPANIONMATRIX_H