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

Go to the documentation of this file.

#ifndef LAMA_GMRES_HPP_
#define LAMA_GMRES_HPP_

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

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

// boost
#include <boost/scoped_array.hpp>

namespace lama
{

class LAMA_DLL_IMPORTEXPORT GMRES : public IterativeSolver
{
public:

    GMRES(const std::string& id);

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

    GMRES( const GMRES& other );

    virtual ~GMRES();

    virtual void initialize(const Matrix&  coefficients);

    void setKrylovDim( unsigned int krylovDim );

    virtual SolverPtr copy();

    struct GMRESRuntime : IterativeSolverRuntime
    {
        GMRESRuntime();
        virtual ~GMRESRuntime();
        
        // arrays to store rotations
        boost::scoped_array<double>         mCC;
        boost::scoped_array<double>         mSS;

        // array for Hessenberg equation
        // H*y=g
        boost::scoped_array<double>         mG;
        boost::scoped_array<double>         mY;

        // Hessenberg matrix
        // mH:  Upper triangular (columnwise)
        // mHd: diagonal band h(i+1,i)
        boost::scoped_array<double>         mH;
        boost::scoped_array<double>         mHd;

        // krylov space
        std::vector<Vector*>                *mV;

        // temp-arrays
        Vector                              *mW;
        Vector                              *mT;

        // remember starting solution
        // only needed if x is modified within krylov loop
        Vector                              *mX0;
    };

    virtual GMRESRuntime& getRuntime();

    virtual const GMRESRuntime& getConstRuntime() const;

protected:

    virtual void iterate();

    GMRESRuntime mGMRESRuntime;

    LAMA_LOG_DECL_STATIC_LOGGER(logger);

private:

    void updateX( unsigned int i );

    // krylov dimension
    unsigned int                        mKrylovDim;
};

}  // namespace lama

#endif // LAMA_GMRES_HPP_