/home/brandes/workspace/LAMA/src/lama/solver/InverseSolver.hpp

Go to the documentation of this file.

#ifndef LAMA_INVERSESOLVER_HPP_
#define LAMA_INVERSESOLVER_HPP_

// for dll_import
#include <lama/config.hpp>

// base classes
#include <lama/solver/Solver.hpp>

// others
#include <lama/matrix/DenseMatrix.hpp>

// logging
#include <logging/Logger.hpp>

namespace lama
{

class LAMA_DLL_IMPORTEXPORT InverseSolver : public Solver
{
public:
    InverseSolver(const std::string& id);

    InverseSolver(const std::string& id, LoggerPtr logger);

    InverseSolver( const InverseSolver& other );

    virtual ~InverseSolver();

    virtual void initialize(const Matrix& coefficients);

    virtual void solveImpl();

    virtual void setContext( ContextPtr context );

    void computeInverse( Matrix& matrix ) const;

    struct InverseSolverRuntime : SolverRuntime
    {
        InverseSolverRuntime();
        virtual ~InverseSolverRuntime();
        
        boost::shared_ptr<Matrix> mInverse;
    };

    virtual SolverPtr copy();

    virtual InverseSolverRuntime& getRuntime();

    virtual const InverseSolverRuntime& getConstRuntime() const;

protected:
    InverseSolverRuntime mInverseSolverRuntime;

private:

    template<typename T>
    void invert( DenseMatrix<T>& matrix, IndexType* const permutation ) const;

    template<typename T>
    void decompose( DenseMatrix<T>& matrix,IndexType* const permutation ) const;

    void logStartSolve();
    void logEndSolve();

    LAMA_LOG_DECL_STATIC_LOGGER(logger);
};

}

#endif // LAMA_INVERSESOLVER_HPP_