MsrsbLevelBuilderCoupled.hpp

Go to the documentation of this file.

/*
 * ------------------------------------------------------------------------------------------------------------
 * SPDX-License-Identifier: LGPL-2.1-only
 *
 * Copyright (c) 2018-2019 Lawrence Livermore National Security LLC
 * Copyright (c) 2018-2019 The Board of Trustees of the Leland Stanford Junior University
 * Copyright (c) 2018-2019 Total, S.A
 * Copyright (c) 2019-     GEOS/GEOSX Contributors
 * All right reserved
 *
 * See top level LICENSE, COPYRIGHT, CONTRIBUTORS, NOTICE, and ACKNOWLEDGEMENTS files for details.
 * ------------------------------------------------------------------------------------------------------------
 */
 
#ifndef GEOS_LINEARALGEBRA_MULTISCALE_MSRSBLEVELBUILDERCOUPLED_HPP
#define GEOS_LINEARALGEBRA_MULTISCALE_MSRSBLEVELBUILDERCOUPLED_HPP
 
#include "MsrsbLevelBuilder.hpp"
#include "linearAlgebra/solvers/BlockPreconditioner.hpp"
 
namespace geos
{
namespace multiscale
{
 
template< typename LAI >
class MsrsbLevelBuilderCoupled : public MsrsbLevelBuilderBase< LAI >
{
public:
 
  using Base = MsrsbLevelBuilderBase< LAI >;
 
  using Vector = typename Base::Vector;
 
  using Matrix = typename Base::Matrix;
 
  explicit MsrsbLevelBuilderCoupled( string name, LinearSolverParameters params );
 
  virtual void initializeFineLevel( DomainPartition & domain,
                                    geos::DofManager const & dofManager,
                                    MPI_Comm const & comm ) override;
 
  virtual void initializeCoarseLevel( LevelBuilderBase< LAI > & fineLevel,
                                      Matrix const & fineMatrix ) override;
 
  virtual bool updateProlongation( Matrix const & fineMatrix ) override;
 
  virtual std::unique_ptr< PreconditionerBase< LAI > > makeCoarseSolver() const override;
 
private:
 
  using Base::m_params;
  using Base::m_name;
  using Base::m_prolongation;
  using Base::m_restriction;
  using Base::m_matrix;
  using Base::m_dofManager;
  using Base::m_preSmoother;
  using Base::m_postSmoother;
  using Base::m_fineLevel;
 
  void initializeCommon( DomainPartition & domain, MPI_Comm const & comm );
 
  void createSmoothers();
 
  void buildProlongationStructure( DofManager const & fineDofManager );
 
  std::vector< string > m_fields;
 
  std::vector< Matrix > m_selectors;
 
  std::vector< std::unique_ptr< MsrsbLevelBuilder< LAI > > > m_builders;
 
  std::vector< CRSMatrix< real64, globalIndex > > m_prolongationBlocks;
 
  CRSMatrix< real64, globalIndex > m_localProlongation;
};
 
} // namespace multiscale
 
} // namespace geos
 
#endif //GEOS_LINEARALGEBRA_MULTISCALE_MSRSBLEVELBUILDERCOUPLED_HPP