AmrSmoother.h

Go to the documentation of this file.

#ifndef AMR_SMOOTHER_H
#define AMR_SMOOTHER_H

#include <string>

#include "AmrMultiGridDefs.h"

#include "Ifpack2_Factory.hpp"

class AmrSmoother {
    
public:
    typedef amr::global_ordinal_t   go_t;
    typedef amr::local_ordinal_t    lo_t;
    typedef amr::scalar_t           scalar_t;
    typedef amr::node_t             node_t;
    typedef amr::matrix_t           matrix_t;
    typedef amr::vector_t           vector_t;
    typedef Ifpack2::Preconditioner<scalar_t,
                                    lo_t,
                                    go_t,
                                    node_t
            > preconditioner_t;
    
    enum Smoother {
        GAUSS_SEIDEL = 0,
        SGS,    // symmetric Gauss-Seidel
        JACOBI //,
//         SOR
    };
    
public:
    AmrSmoother(const Teuchos::RCP<const matrix_t>& A,
                const Smoother& smoother,
                lo_t nSweeps);
    
    ~AmrSmoother();
    
    void smooth(const Teuchos::RCP<vector_t>& x,
                const Teuchos::RCP<matrix_t>& A,
                const Teuchos::RCP<vector_t>& b);
    
    static Smoother convertToEnumSmoother(const std::string& smoother);
    
private:
    void initParameter_m(const Smoother& smoother,
                         lo_t nSweeps);
    
    
private:
    Teuchos::RCP<preconditioner_t> prec_mp;
    
    Teuchos::RCP<Teuchos::ParameterList> params_mp;
};

#endif