Ifpack2Preconditioner.h

Go to the documentation of this file.

//
// Class Ifpack2Preconditioner
//   Interface to Ifpack2 preconditioners.
//
// Copyright (c) 2017 - 2020, Matthias Frey, Paul Scherrer Institut, Villigen PSI, Switzerland
// All rights reserved
//
// Implemented as part of the PhD thesis
// "Precise Simulations of Multibunches in High Intensity Cyclotrons"
//
// This file is part of OPAL.
//
// OPAL 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.
//
// You should have received a copy of the GNU General Public License
// along with OPAL. If not, see <https://www.gnu.org/licenses/>.
//
#ifndef IFPACK_2_PRECONDITIONER_H
#define IFPACK_2_PRECONDITIONER_H
 
#include "AmrPreconditioner.h"
 
#include <Ifpack2_Factory.hpp>
 
template <class Level>
class Ifpack2Preconditioner : public AmrPreconditioner<amr::matrix_t, Level>
{
public:
    typedef amr::Preconditioner Preconditioner;
    typedef Ifpack2::Preconditioner<
        amr::scalar_t,
        amr::local_ordinal_t,
        amr::global_ordinal_t,
        amr::node_t
    > precond_t;
    
    typedef std::map<std::string, Preconditioner> map_t;
    
public:
    
    Ifpack2Preconditioner(Preconditioner prec);
    
    void create(const Teuchos::RCP<amr::matrix_t>& A, Level* level_p = nullptr);
    
    Teuchos::RCP<amr::operator_t> get();
    
    static void fillMap(map_t& map);
    
private:
    void init_m(Preconditioner prec);
    
private:
    std::string prectype_m;
    
    Teuchos::RCP<Teuchos::ParameterList> params_mp;
    
    Teuchos::RCP<precond_t> prec_mp;
};
 
#include "Ifpack2Preconditioner.hpp"
 
#endif