3#ifndef DUNE_ISTL_TWOLEVELMETHODCPR_HH
4#define DUNE_ISTL_TWOLEVELMETHODCPR_HH
11#include<dune/istl/operators.hh>
14#include<dune/istl/paamg/amg.hh>
15#include<dune/istl/paamg/galerkin.hh>
16#include<dune/istl/solver.hh>
18#include<dune/common/unused.hh>
19#include<dune/common/version.hh>
42template<
class FO,
class CO>
151template<
class O,
class C>
155 typedef Dune::Amg::AggregatesMap<typename O::matrix_type::size_type> AggregatesMap;
159 typedef SequentialInformation ParallelInformation;
167 prolongDamp_ = criterion_.getProlongationDampingFactor();
168 GalerkinProduct<ParallelInformation> productBuilder;
169 typedef typename Dune::Amg::MatrixGraph<const typename O::matrix_type> MatrixGraph;
170 typedef typename Dune::Amg::PropertiesGraph<MatrixGraph,Dune::Amg::VertexProperties,
171 Dune::Amg::EdgeProperties,Dune::IdentityMap,Dune::IdentityMap> PropertiesGraph;
172 MatrixGraph mg(fineOperator.getmat());
173 PropertiesGraph pg(mg,Dune::IdentityMap(),Dune::IdentityMap());
174 typedef NegateSet<typename ParallelInformation::OwnerSet> OverlapFlags;
176 aggregatesMap_.reset(
new AggregatesMap(pg.maxVertex()+1));
178 int noAggregates, isoAggregates, oneAggregates, skippedAggregates;
180 std::tie(noAggregates, isoAggregates, oneAggregates, skippedAggregates) =
181 aggregatesMap_->buildAggregates(fineOperator.getmat(), pg, criterion_,
true);
182 std::cout<<
"no aggregates="<<noAggregates<<
" iso="<<isoAggregates<<
" one="<<oneAggregates<<
" skipped="<<skippedAggregates<<std::endl;
184 Dune::Amg::IndicesCoarsener<Dune::Amg::SequentialInformation,int> renumberer;
185 typedef std::vector<bool>::iterator Iterator;
186 typedef Dune::IteratorPropertyMap<Iterator, Dune::IdentityMap> VisitedMap;
187 std::vector<bool> excluded(fineOperator.getmat().N(),
false);
188 VisitedMap vm(excluded.begin(), Dune::IdentityMap());
189 ParallelInformation pinfo;
191 std::size_t aggregates = coarsen(renumberer, pinfo, pg, vm,*aggregatesMap_, noAggregates,
true);
192 std::vector<bool>& visited=excluded;
194 typedef std::vector<bool>::iterator Iterator;
196 for(Iterator iter= visited.begin(), end=visited.end();
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,
476 smoother_ = smoother;
479 coarsePolicy.setCoarseOperator(*policy_);
480 coarseSolver_->updatePreconditioner();
484 coarseSolver_ = coarsePolicy.createCoarseLevelSolver(*policy_);
501 LevelContext context;
502 SequentialInformation info;
506 context.smoother=smoother_;
508 context.matrix=operator_;
510 presmooth(context, preSteps_);
513 InverseOperatorResult res;
517 *context.update += *context.lhs;
519 postsmooth(context, postSteps_);
523 virtual SolverCategory::Category category()
const
525 return SolverCategory::sequential;
535 typedef S SmootherType;
537 std::shared_ptr<SmootherType> smoother;
553 SequentialInformation* pinfo;
565 std::shared_ptr<S> smoother_;
567 LevelTransferPolicyCpr<FO,typename CSP::Operator>* policy_;
569 std::size_t preSteps_;
571 std::size_t postSteps_;
A LeveTransferPolicy that used aggregation to construct the coarse level system.
Definition twolevelmethodcpr.hh:154
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:165
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:44
CoarseOperatorType::range_type CoarseRangeType
The type of the range of the coarse level operator.
Definition twolevelmethodcpr.hh:67
CoarseDomainType & getCoarseLevelLhs()
Get the coarse level left hand side.
Definition twolevelmethodcpr.hh:93
virtual LevelTransferPolicyCpr * clone() const =0
Clone the current object.
FO FineOperatorType
The linear operator of the finel level system.
Definition twolevelmethodcpr.hh:50
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:71
FineOperatorType::range_type FineRangeType
The type of the range of the fine level operator.
Definition twolevelmethodcpr.hh:54
CoarseRangeType & getCoarseLevelRhs()
Get the coarse level right hand side.
Definition twolevelmethodcpr.hh:84
virtual ~LevelTransferPolicyCpr()
Destructor.
Definition twolevelmethodcpr.hh:135
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:141
CoarseRangeType rhs_
The coarse level rhs.
Definition twolevelmethodcpr.hh:139
FineOperatorType::domain_type FineDomainType
The type of the domain of the fine level operator.
Definition twolevelmethodcpr.hh:58
CO CoarseOperatorType
The linear operator of the finel level system.
Definition twolevelmethodcpr.hh:63
std::shared_ptr< CoarseOperatorType > & getCoarseLevelOperator()
Get the coarse level operator.
Definition twolevelmethodcpr.hh:76
std::shared_ptr< CoarseOperatorType > operator_
the coarse level linear operator.
Definition twolevelmethodcpr.hh:143
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