3#ifndef DUNE_ISTL_TWOLEVELMETHODCPR_HH
4#define DUNE_ISTL_TWOLEVELMETHODCPR_HH
9#include <dune/istl/operators.hh>
12#include <dune/istl/paamg/amg.hh>
13#include <dune/istl/paamg/galerkin.hh>
14#include <dune/istl/solver.hh>
16#include <dune/common/unused.hh>
17#include <dune/common/version.hh>
47template<
class FO,
class CO>
156template<
class O,
class C>
160 typedef Dune::Amg::AggregatesMap<typename O::matrix_type::size_type> AggregatesMap;
164 typedef SequentialInformation ParallelInformation;
172 prolongDamp_ = criterion_.getProlongationDampingFactor();
173 GalerkinProduct<ParallelInformation> productBuilder;
174 typedef typename Dune::Amg::MatrixGraph<const typename O::matrix_type> MatrixGraph;
175 typedef typename Dune::Amg::PropertiesGraph<MatrixGraph,Dune::Amg::VertexProperties,
176 Dune::Amg::EdgeProperties,Dune::IdentityMap,Dune::IdentityMap> PropertiesGraph;
177 MatrixGraph mg(fineOperator.getmat());
178 PropertiesGraph pg(mg,Dune::IdentityMap(),Dune::IdentityMap());
179 typedef NegateSet<typename ParallelInformation::OwnerSet> OverlapFlags;
181 aggregatesMap_.reset(
new AggregatesMap(pg.maxVertex()+1));
183 int noAggregates, isoAggregates, oneAggregates, skippedAggregates;
185 std::tie(noAggregates, isoAggregates, oneAggregates, skippedAggregates) =
186 aggregatesMap_->buildAggregates(fineOperator.getmat(), pg, criterion_,
true);
187 std::cout<<
"no aggregates="<<noAggregates<<
" iso="<<isoAggregates<<
" one="<<oneAggregates<<
" skipped="<<skippedAggregates<<std::endl;
189 Dune::Amg::IndicesCoarsener<Dune::Amg::SequentialInformation,int> renumberer;
190 typedef std::vector<bool>::iterator Iterator;
191 typedef Dune::IteratorPropertyMap<Iterator, Dune::IdentityMap> VisitedMap;
192 std::vector<bool> excluded(fineOperator.getmat().N(),
false);
193 VisitedMap vm(excluded.begin(), Dune::IdentityMap());
194 ParallelInformation pinfo;
196 std::size_t aggregates = coarsen(renumberer, pinfo, pg, vm,*aggregatesMap_, noAggregates,
true);
198 std::fill(excluded.begin(), excluded.end(),
false);
199 matrix_.reset(productBuilder.build(mg, vm,
200 SequentialInformation(),
204 productBuilder.calculate(fineOperator.getmat(), *aggregatesMap_, *matrix_, pinfo, OverlapFlags());
205 this->
lhs_.resize(this->matrix_->M());
206 this->
rhs_.resize(this->matrix_->N());
212 Transfer<std::size_t,typename FatherType::FineRangeType,ParallelInformation>
213 ::restrictVector(*aggregatesMap_, this->
rhs_, fineRhs, ParallelInformation());
219 Transfer<std::size_t,typename FatherType::FineRangeType,ParallelInformation>
220 ::prolongateVector(*aggregatesMap_, this->
lhs_, fineLhs,
221 prolongDamp_, ParallelInformation());
231 template <
typename RN,
typename PI,
typename PG,
typename VM,
typename AM>
232 auto coarsen(RN& renumberer,
238 bool useFixedOrder) ->
239 decltype(renumberer.coarsen(pinfo, pg, vm, aggregatesMap, pinfo, noAggregates, useFixedOrder))
241 return renumberer.coarsen(pinfo, pg, vm, aggregatesMap, pinfo, noAggregates, useFixedOrder);
244 template <
typename RN,
typename PI,
typename PG,
typename VM,
typename AM>
245 auto coarsen(RN& renumberer,
250 int noAggregates, ...)
252 return renumberer.coarsen(pinfo, pg, vm, aggregatesMap, pinfo, noAggregates);
255 typename O::matrix_type::field_type prolongDamp_;
256 std::shared_ptr<AggregatesMap> aggregatesMap_;
257 Criterion criterion_;
258 std::shared_ptr<typename O::matrix_type> matrix_;
267template<
class O,
class S,
class C>
274 typedef typename O::range_type
X;
289 : smootherArgs_(args), criterion_(c)
293 : coarseOperator_(other.coarseOperator_), smootherArgs_(other.smootherArgs_),
294 criterion_(other.criterion_)
303 struct AMGInverseOperator :
public InverseOperator<X,X>
305 AMGInverseOperator(
const typename AMGType::Operator& op,
307 const typename AMGType::SmootherArgs& args)
308 : amg_(op, crit,args), first_(true)
311 void apply(
X& x,
X& b, [[maybe_unused]]
double reduction, [[maybe_unused]] InverseOperatorResult& res)
322 void apply(
X& x,
X& b, InverseOperatorResult& res)
324 return apply(x,b,1e-8,res);
327 virtual SolverCategory::Category category()
const
329 return amg_.category();
332 ~AMGInverseOperator()
337 AMGInverseOperator(
const AMGInverseOperator& other)
338 : x_(other.x_), amg_(other.amg_), first_(other.first_)
361 coarseOperator_=transferPolicy.getCoarseLevelOperator();
362 AMGInverseOperator* inv =
new AMGInverseOperator(*coarseOperator_,
372 std::shared_ptr<Operator> coarseOperator_;
384template<
class FO,
class CSP,
class S>
386 public Preconditioner<typename FO::domain_type, typename FO::range_type>
439 std::shared_ptr<SmootherType> smoother,
443 std::size_t preSteps=1, std::size_t postSteps=1)
444 : operator_(&op), smoother_(smoother),
445 preSteps_(preSteps), postSteps_(postSteps)
447 policy_ = policy.
clone();
449 coarseSolver_=coarsePolicy.createCoarseLevelSolver(*policy_);
453 : operator_(other.operator_), coarseSolver_(new
CoarseLevelSolver(*other.coarseSolver_)),
454 smoother_(other.smoother_), policy_(other.policy_->clone()),
455 preSteps_(other.preSteps_), postSteps_(other.postSteps_)
462 delete coarseSolver_;
466 std::shared_ptr<SmootherType> smoother,
469 updatePreconditioner(smoother, coarsePolicy);
472 void updatePreconditioner(std::shared_ptr<SmootherType> smoother,
475 smoother_ = smoother;
480 coarsePolicy.setCoarseOperator(*policy_);
481 coarseSolver_->updatePreconditioner();
485 coarseSolver_ = coarsePolicy.createCoarseLevelSolver(*policy_);
498 bool hasPerfectUpdate()
const
501 return smoother_->hasPerfectUpdate() && coarseSolver_->hasPerfectUpdate();
508 LevelContext context;
509 SequentialInformation info;
513 context.smoother=smoother_;
515 context.matrix=operator_;
517 presmooth(context, preSteps_);
520 InverseOperatorResult res;
524 *context.update += *context.lhs;
526 postsmooth(context, postSteps_);
530 virtual SolverCategory::Category category()
const
532 return SolverCategory::sequential;
542 typedef S SmootherType;
544 std::shared_ptr<SmootherType> smoother;
560 SequentialInformation* pinfo;
572 std::shared_ptr<S> smoother_;
574 LevelTransferPolicyCpr<FO,typename CSP::Operator>* policy_;
576 std::size_t preSteps_;
578 std::size_t postSteps_;
A LeveTransferPolicy that used aggregation to construct the coarse level system.
Definition twolevelmethodcpr.hh:159
void moveToFineLevel(typename FatherType::FineDomainType &fineLhs)
Updates the fine level linear system after the correction of the coarse levels system.
Definition twolevelmethodcpr.hh:217
void createCoarseLevelSystem(const O &fineOperator)
Algebraically creates the coarse level system.
Definition twolevelmethodcpr.hh:170
AggregationLevelTransferPolicyCpr * clone() const
Clone the current object.
Definition twolevelmethodcpr.hh:224
Abstract base class for transfer between levels and creation of the coarse level system.
Definition twolevelmethodcpr.hh:49
CoarseOperatorType::range_type CoarseRangeType
The type of the range of the coarse level operator.
Definition twolevelmethodcpr.hh:72
CoarseDomainType & getCoarseLevelLhs()
Get the coarse level left hand side.
Definition twolevelmethodcpr.hh:98
virtual LevelTransferPolicyCpr * clone() const =0
Clone the current object.
FO FineOperatorType
The linear operator of the finel level system.
Definition twolevelmethodcpr.hh:55
virtual void createCoarseLevelSystem(const FineOperatorType &fineOperator)=0
Algebraically creates the coarse level system.
virtual void moveToCoarseLevel(const FineRangeType &fineRhs)=0
Transfers the data to the coarse level.
CoarseOperatorType::domain_type CoarseDomainType
The type of the domain of the coarse level operator.
Definition twolevelmethodcpr.hh:76
FineOperatorType::range_type FineRangeType
The type of the range of the fine level operator.
Definition twolevelmethodcpr.hh:59
CoarseRangeType & getCoarseLevelRhs()
Get the coarse level right hand side.
Definition twolevelmethodcpr.hh:89
virtual ~LevelTransferPolicyCpr()
Destructor.
Definition twolevelmethodcpr.hh:140
virtual void moveToFineLevel(FineDomainType &fineLhs)=0
Updates the fine level linear system after the correction of the coarse levels system.
virtual void calculateCoarseEntries(const FineOperatorType &fineOperator)=0
???.
CoarseDomainType lhs_
The coarse level lhs.
Definition twolevelmethodcpr.hh:146
CoarseRangeType rhs_
The coarse level rhs.
Definition twolevelmethodcpr.hh:144
FineOperatorType::domain_type FineDomainType
The type of the domain of the fine level operator.
Definition twolevelmethodcpr.hh:63
CO CoarseOperatorType
The linear operator of the finel level system.
Definition twolevelmethodcpr.hh:68
std::shared_ptr< CoarseOperatorType > & getCoarseLevelOperator()
Get the coarse level operator.
Definition twolevelmethodcpr.hh:81
std::shared_ptr< CoarseOperatorType > operator_
the coarse level linear operator.
Definition twolevelmethodcpr.hh:148
A policy class for solving the coarse level system using one step of AMG.
Definition twolevelmethodcpr.hh:269
AMG< Operator, X, Smoother > AMGType
The type of the AMG construct on the coarse level.
Definition twolevelmethodcpr.hh:282
OneStepAMGCoarseSolverPolicyCpr(const OneStepAMGCoarseSolverPolicyCpr &other)
Copy constructor.
Definition twolevelmethodcpr.hh:292
AMGInverseOperator CoarseLevelSolver
The type of solver constructed for the coarse level.
Definition twolevelmethodcpr.hh:349
S Smoother
The type of the smoother used in AMG.
Definition twolevelmethodcpr.hh:278
C Criterion
The type of the crition used for the aggregation within AMG.
Definition twolevelmethodcpr.hh:276
Dune::Amg::SmootherTraits< S >::Arguments SmootherArgs
The type of the arguments used for constructing the smoother.
Definition twolevelmethodcpr.hh:280
CoarseLevelSolver * createCoarseLevelSolver(P &transferPolicy)
Constructs a coarse level solver.
Definition twolevelmethodcpr.hh:359
O::range_type X
The type of the range and domain of the operator.
Definition twolevelmethodcpr.hh:274
OneStepAMGCoarseSolverPolicyCpr(const SmootherArgs &args, const Criterion &c)
Constructs the coarse solver policy.
Definition twolevelmethodcpr.hh:288
O Operator
The type of the linear operator used.
Definition twolevelmethodcpr.hh:272
Definition twolevelmethodcpr.hh:387
S SmootherType
The type of the fine level smoother.
Definition twolevelmethodcpr.hh:422
FineOperatorType::range_type FineRangeType
The type of the range of the fine level operator.
Definition twolevelmethodcpr.hh:401
TwoLevelMethodCpr(const FineOperatorType &op, std::shared_ptr< SmootherType > smoother, const LevelTransferPolicyCpr< FineOperatorType, CoarseOperatorType > &policy, CoarseLevelSolverPolicy &coarsePolicy, std::size_t preSteps=1, std::size_t postSteps=1)
Constructs a two level method.
Definition twolevelmethodcpr.hh:438
FineOperatorType::domain_type FineDomainType
The type of the domain of the fine level operator.
Definition twolevelmethodcpr.hh:405
CoarseOperatorType::domain_type CoarseDomainType
The type of the domain of the coarse level operator.
Definition twolevelmethodcpr.hh:418
CSP::Operator CoarseOperatorType
The linear operator of the finel level system.
Definition twolevelmethodcpr.hh:410
CSP CoarseLevelSolverPolicy
The type of the policy for constructing the coarse level solver.
Definition twolevelmethodcpr.hh:390
CoarseLevelSolverPolicy::CoarseLevelSolver CoarseLevelSolver
The type of the coarse level solver.
Definition twolevelmethodcpr.hh:392
CoarseOperatorType::range_type CoarseRangeType
The type of the range of the coarse level operator.
Definition twolevelmethodcpr.hh:414
FO FineOperatorType
The linear operator of the finel level system.
Definition twolevelmethodcpr.hh:397