DenseMatrixOperator.h

Go to the documentation of this file.

/*
 * This program is free software; you can redistribute it and/or modify
 * it under the terms of the GNU General Public License as published by
 * the Free Software Foundation; either version 3 of the License, or
 * (at your option) any later version.
 * 
 * Written (W) 2013 Soumyajit De
 */

#ifndef DENSE_MATRIX_OPERATOR_H_
#define DENSE_MATRIX_OPERATOR_H_

#include <shogun/lib/config.h>

#ifdef HAVE_EIGEN3
#include <shogun/mathematics/linalg/linop/MatrixOperator.h>

namespace shogun
{
template<class T> class SGVector;
template<class T> class SGMatrix;

template<class T> class CDenseMatrixOperator : public CMatrixOperator<T>
{
typedef bool supports_complex128_t;

public:
    CDenseMatrixOperator();

    explicit CDenseMatrixOperator(SGMatrix<T> op);

    CDenseMatrixOperator(const CDenseMatrixOperator<T>& orig);

    ~CDenseMatrixOperator();

    virtual SGVector<T> apply(SGVector<T> b) const;

    virtual void set_diagonal(SGVector<T> diag);

    virtual SGVector<T> get_diagonal() const;

    SGMatrix<T> get_matrix_operator() const;

    template<class Scalar>
    inline operator CDenseMatrixOperator<Scalar>*() const
    {
        REQUIRE(m_operator.matrix, "Matrix is not initialized!\n");

        SGMatrix<Scalar> casted_m(m_operator.num_rows, m_operator.num_cols);
        for (index_t i=0; i<m_operator.num_cols; ++i)
        {
            for (index_t j=0; j<m_operator.num_rows; ++j)
                casted_m(j,i)=static_cast<Scalar>(m_operator(j,i));
        }
        SG_SDEBUG("DenseMatrixOperator::static_cast(): Creating casted operator!\n");

        return new CDenseMatrixOperator<Scalar>(casted_m);
    }

    virtual const char* get_name() const
    {
        return "DenseMatrixOperator";
    }

private:
    SGMatrix<T> m_operator;

    void init();

};

}

#endif // HAVE_EIGEN3
#endif // DENSE_MATRIX_OPERATOR_H_