amgcl_mpi_solver.h

Go to the documentation of this file.

//    |  /           |
//    ' /   __| _` | __|  _ \   __|
//    . \  |   (   | |   (   |\__ `
//   _|\_\_|  \__,_|\__|\___/ ____/
//                   Multi-Physics
//
//  License:         BSD License
//                     Kratos default license: kratos/license.txt
//
//  Main authors:    Denis Demidov
//                   Riccardo Rossi
//
 
#if !defined (KRATOS_AMGCL_MPI_SOLVER_H_INCLUDED)
#define KRATOS_AMGCL_MPI_SOLVER_H_INCLUDED
 
#ifndef AMGCL_PARAM_UNKNOWN
#include "input_output/logger.h"
#  define AMGCL_PARAM_UNKNOWN(name)                                            \
    Kratos::Logger("AMGCL") << KRATOS_CODE_LOCATION << Kratos::Logger::Severity::WARNING << "Unknown parameter " << name << std::endl
#endif
 
// System includes
#include <iostream>
#include <fstream>
#include <utility>
 
// Project includes
#include "includes/define.h"
#include "includes/kratos_parameters.h"
#include "external_includes/amgcl_mpi_solve_functions.h"
#include "linear_solvers/amgcl_solver.h"
 
namespace Kratos
{
 
template< class TSparseSpaceType, class TDenseSpaceType,
          class TReordererType = Reorderer<TSparseSpaceType, TDenseSpaceType> >
class AmgclMPISolver : public AMGCLSolver< TSparseSpaceType,TDenseSpaceType, TReordererType >
{
public:
 
    KRATOS_CLASS_POINTER_DEFINITION( AmgclMPISolver );
 
    typedef typename TSparseSpaceType::MatrixType SparseMatrixType;
 
    typedef typename TSparseSpaceType::VectorType VectorType;
 
    typedef typename TDenseSpaceType::MatrixType DenseMatrixType;
 
    typedef ModelPart::DofsArrayType DofsArrayType;
 
    typedef std::size_t IndexType;
 
    typedef std::size_t SizeType;
 
 
    AmgclMPISolver(Parameters ThisParameters = Parameters(R"({})"))
        :
        AMGCLSolver< TSparseSpaceType,TDenseSpaceType, TReordererType>(ThisParameters) { }
 
    AmgclMPISolver(const AmgclMPISolver& Other) = delete;
 
    ~AmgclMPISolver() override = default;
 
 
    AmgclMPISolver& operator=(const AmgclMPISolver& Other) = delete;
 
 
    void  PrintInfo(std::ostream& rOStream) const override
    {
        rOStream << "AMGCL-MPI-Solver";
    }
 
 
}; // Class AmgclMPISolver
 
template<class TSparseSpaceType, class TDenseSpaceType, class TReordererType>
inline std::ostream& operator << (std::ostream& rOStream,
                                  const AmgclMPISolver<TSparseSpaceType,
                                  TDenseSpaceType, TReordererType>& rThis)
{
    rThis.PrintInfo(rOStream);
    rOStream << std::endl;
    rThis.PrintData(rOStream);
 
    return rOStream;
}
 
}  // namespace Kratos.
 
#endif // KRATOS_AMGCL_MPI_SOLVER_H_INCLUDED defined