DGCGCoupling.h

Go to the documentation of this file.

// This file is part of the Peano's MGHyPE extension. For conditions of
// distribution and use, please see the copyright notice at
// www.peano-framework.org
#pragma once
 
 
#include "tarch/la/Matrix.h"
#include "tarch/la/Vector.h"
#include "peano4/utils/Globals.h"
 
 
namespace mghype {
  namespace matrixfree {
    namespace solvers {
      namespace dgcgcoupling {
        void rollOverAndPrepareCGVertex(
          auto&  vertex,
          bool   restrictRhsAndInjectSolution,
          bool   fas
        );
 
        extern tarch::logging::Log _log;
 
        void updateRestrictionCounters(
          auto&  vertices
        );
 
        template <class DGSolver, class CGSolver>
        void injectSolution(
          const tarch::la::Matrix< TwoPowerD * CGSolver::VertexUnknowns, DGSolver::NodesPerCell * DGSolver::UnknownsPerCellNode, double >&  injectionMatrix,
          const tarch::la::Vector< DGSolver::NodesPerCell * DGSolver::UnknownsPerCellNode, double >&                                        dgSolution,
          auto&           cgVertices
        );
 
        template <class DGSolver, class CGSolver>
        tarch::la::Vector< DGSolver::NodesPerCell * DGSolver::UnknownsPerCellNode, double > prolongate(
          const tarch::la::Matrix< DGSolver::NodesPerCell * DGSolver::UnknownsPerCellNode, TwoPowerD * CGSolver::VertexUnknowns, double >&  prolongationMatrix,
          const tarch::la::Vector< DGSolver::NodesPerCell * DGSolver::UnknownsPerCellNode, double >&                                        dgSolution,
          const auto&           cgVertices
        );
 
        template <class DGSolver, class CGSolver>
        void computeAndRestrictResidual(
          const tarch::la::Matrix< TwoPowerD * CGSolver::VertexUnknowns, DGSolver::NodesPerCell * DGSolver::UnknownsPerCellNode, double >&  restrictionMatrix,
          const tarch::la::Matrix< DGSolver::NodesPerCell * DGSolver::UnknownsPerCellNode, DGSolver::NodesPerCell * DGSolver::UnknownsPerCellNode, double >&  massMatrix,
          const tarch::la::Matrix< DGSolver::NodesPerCell * DGSolver::UnknownsPerCellNode, DGSolver::NodesPerCell * DGSolver::UnknownsPerCellNode, double >&  systemMatrix,
          const tarch::la::Vector< DGSolver::NodesPerCell * DGSolver::UnknownsPerCellNode, double >&                                        dgSolution,
          const tarch::la::Vector< DGSolver::NodesPerCell * DGSolver::UnknownsPerCellNode, double >&                                        dgRhs,
          const tarch::la::Vector< DGSolver::NodesPerCell * DGSolver::UnknownsPerCellNode, double >&                                        faceSolutions,
          auto&           cgVertices
        );
 
        template <class DGSolver, class CGSolver>
        void computeAndRestrictHierarchicalResidual(
          const tarch::la::Matrix< DGSolver::NodesPerCell * DGSolver::UnknownsPerCellNode, TwoPowerD * CGSolver::VertexUnknowns, double >&  prolongationMatrix,
          const tarch::la::Matrix< TwoPowerD * CGSolver::VertexUnknowns, DGSolver::NodesPerCell * DGSolver::UnknownsPerCellNode, double >&  restrictionMatrix,
          const tarch::la::Matrix< DGSolver::NodesPerCell * DGSolver::UnknownsPerCellNode, DGSolver::NodesPerCell * DGSolver::UnknownsPerCellNode, double >&  massMatrix,
          const tarch::la::Matrix< DGSolver::NodesPerCell * DGSolver::UnknownsPerCellNode, DGSolver::NodesPerCell * DGSolver::UnknownsPerCellNode, double >&  systemMatrix,
          const tarch::la::Vector< DGSolver::NodesPerCell * DGSolver::UnknownsPerCellNode, double >&                                        dgSolution,
          const tarch::la::Vector< DGSolver::NodesPerCell * DGSolver::UnknownsPerCellNode, double >&                                        dgRhs,
          auto&           cgVertices
        );
 
        template <class CGSolver>
        tarch::la::Vector< TwoPowerD * CGSolver::VertexUnknowns, double > getCoarseGridSolution(
          auto&           cgVertices
        );
      }
    }
  }
}
 
 
#include "DGCGCoupling.cpph"