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

Go to the documentation of this file.

#ifndef LAMA_DEFAULTJACOBI_HPP_
#define LAMA_DEFAULTJACOBI_HPP_

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

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

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

namespace lama
{

class LAMA_DLL_IMPORTEXPORT DefaultJacobi : public OmegaSolver
{
public:
    DefaultJacobi( const std::string& id );

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

    DefaultJacobi( const std::string& id, const Scalar omega ); //2nd param Matrix.Scalar

    DefaultJacobi( const std::string& id, const Scalar omega, LoggerPtr logger );

    DefaultJacobi( const DefaultJacobi& other );

    virtual ~DefaultJacobi();

    virtual void initialize( const Matrix& coefficients );

    virtual void solve( Vector& solution, const Vector& rhs );

    virtual void solveInit( Vector& solution, const Vector& rhs );

    virtual void solveFinalize();

    virtual SolverPtr copy();


    struct DefaultJacobiRuntime : OmegaSolverRuntime
    {
        DefaultJacobiRuntime();
        virtual ~DefaultJacobiRuntime();
        
        std::auto_ptr<Matrix> mDiagonalTimesLU;
        std::auto_ptr<Matrix> mDiagonalInverted;
        std::auto_ptr<Vector> mDiagonalTimesRhs;
        std::auto_ptr<Vector> mOldSolution;
        SolutionProxy mProxyOldSolution;
    };

    virtual DefaultJacobiRuntime& getRuntime();

    virtual const DefaultJacobiRuntime& getConstRuntime() const;

protected:
    DefaultJacobiRuntime mDefaultJacobiRuntime;

    virtual void iterate();

    LAMA_LOG_DECL_STATIC_LOGGER(logger);

private:

    template<typename T>
    void initialize(const Matrix& coefficients);

    template<typename T>
    void iterate();
};

}// namespace lama

#endif // LAMA_DEFAULTJACOBI_HPP_