KroneckerProductOperator.h

Go to the documentation of this file.

#ifndef KRONECKERPRODUCTOPERATOR_H
#define KRONECKERPRODUCTOPERATOR_H
 
#include "MUQ/Modeling/LinearAlgebra/LinearOperator.h"
#include "MUQ/Modeling/LinearAlgebra/SumOperator.h"
 
namespace muq
{
namespace Modeling
{
 
    class KroneckerProductOperator : public LinearOperator
    {
 
    public:
        KroneckerProductOperator(std::shared_ptr<LinearOperator> Ain,
                                 std::shared_ptr<LinearOperator> Bin);
 
 
        virtual Eigen::MatrixXd Apply(Eigen::Ref<const Eigen::MatrixXd> const& x) override;
  
        virtual Eigen::MatrixXd ApplyTranspose(Eigen::Ref<const Eigen::MatrixXd> const& x) override;
 
    private:
        std::shared_ptr<LinearOperator> A;
        std::shared_ptr<LinearOperator> B;
 
    }; // class KroneckerProductOperator
 
 
 
    std::shared_ptr<LinearOperator> KroneckerSum(std::shared_ptr<LinearOperator> A,
                                                 std::shared_ptr<LinearOperator> B);
 
    Eigen::MatrixXd KroneckerProduct(Eigen::Ref<const Eigen::MatrixXd> const& A,
                                     Eigen::Ref<const Eigen::MatrixXd> const& B);
 
} // namespace Modeling
} // namespace muq
 
 
#endif // #ifndef KRONECKERPRODUCTOPERATOR_H