3 #ifndef DUNE_GLOBALAGGREGATES_HH
4 #define DUNE_GLOBALAGGREGATES_HH
18 #include <dune/common/parallel/indexset.hh>
19 #include <dune/common/unused.hh>
26 template<
typename T,
typename TI>
36 typedef typename ParallelIndexSet::LocalIndex
LocalIndex;
41 const GlobalLookupIndexSet<ParallelIndexSet>& indexset)
42 : aggregates_(aggregates), indexset_(indexset)
47 const Vertex& aggregate = aggregates_[index];
50 return isolatedMarker;
52 const Dune::IndexPair<GlobalIndex,LocalIndex >* pair = indexset_.pair(aggregate);
54 return pair->global();
61 const Vertex& aggregate = aggregates_[index];
63 const Dune::IndexPair<GlobalIndex,LocalIndex >* pair = indexset_.pair(aggregate);
71 Proxy(
const GlobalLookupIndexSet<ParallelIndexSet>& indexset,
Vertex& aggregate)
72 : indexset_(&indexset), aggregate_(&aggregate)
77 if(global==isolatedMarker)
81 *aggregate_ = indexset_->operator[](global).local();
86 const GlobalLookupIndexSet<ParallelIndexSet>* indexset_;
92 return Proxy(indexset_, aggregates_[index]);
97 aggregates_[i]=indexset_[global].local();
103 const GlobalLookupIndexSet<ParallelIndexSet>& indexset_;
107 template<
typename T,
typename TI>
108 const typename TI::GlobalIndex GlobalAggregatesMap<T,TI>::isolatedMarker =
109 std::numeric_limits<typename TI::GlobalIndex>::max();
111 template<
typename T,
typename TI>
128 template<
typename T,
typename O,
typename I>
140 template<
class T1,
class T2>
157 template<
typename T,
typename O,
typename T1,
typename T2>
171 GlobalMap gmap(aggregates, globalLookup);
178 typedef typename ParallelInformation::RemoteIndices::const_iterator Lists;
187 for(Iter i=aggregates.
begin(), end=aggregates.
end(); i!=end; ++i)
188 maxAggregate = std::max(maxAggregate, *i);
191 std::map<Vertex,Vertex> newMapping;
194 typedef typename ParallelInformation::RemoteIndices::RemoteIndexList
195 ::const_iterator RIter;
196 for(RIter ri=lists->second.first->begin(), rend = lists->second.first->end();
198 if(O::contains(ri->localIndexPair().local().attribute()))
199 newMapping.insert(std::make_pair(aggregates[ri->localIndexPair().local()],
202 typedef typename std::map<Vertex,Vertex>::iterator MIter;
203 for(MIter mi=newMapping.begin(), mend=newMapping.end();
205 mi->second=++maxAggregate;
208 for(RIter ri=lists->second.first->begin(), rend = lists->second.first->end();
210 if(O::contains(ri->localIndexPair().local().attribute()))
211 aggregates[ri->localIndexPair().local()] =
212 newMapping[aggregates[ri->localIndexPair().local()]];
218 template<
typename T,
typename O>
229 DUNE_UNUSED_PARAMETER(aggregates);
230 DUNE_UNUSED_PARAMETER(pinfo);
231 DUNE_UNUSED_PARAMETER(globalLookup);
239 template<
typename T,
typename TI>
240 struct CommPolicy<Amg::GlobalAggregatesMap<T,TI> >
Provides classes for the Coloring process of AMG.
static const GlobalIndex & gather(const GlobalAggregatesMap< T, TI > &ga, size_t i)
Definition: globalaggregates.hh:117
ParallelInformation::GlobalLookupIndexSet GlobalLookupIndexSet
Definition: globalaggregates.hh:163
O OverlapFlags
Definition: globalaggregates.hh:161
static int getSize(const Type &, int)
Definition: globalaggregates.hh:245
GlobalIndex & get(std::size_t index)
Definition: globalaggregates.hh:59
Amg::GlobalAggregatesMap< T, TI >::IndexedType IndexedType
Definition: globalaggregates.hh:243
T Vertex
Definition: globalaggregates.hh:160
static void publish(AggregatesMap< Vertex > &aggregates, ParallelInformation &pinfo, const GlobalLookupIndexSet &globalLookup)
Definition: globalaggregates.hh:225
SequentialInformation ParallelInformation
Definition: globalaggregates.hh:222
ParallelIndexSet::GlobalIndex GlobalIndex
Definition: globalaggregates.hh:32
ParallelInformation::ParallelIndexSet IndexSet
Definition: globalaggregates.hh:164
void put(const GlobalIndex &global, size_t i)
Definition: globalaggregates.hh:95
T Vertex
Definition: globalaggregates.hh:38
GlobalAggregatesMap(AggregatesMap< Vertex > &aggregates, const GlobalLookupIndexSet< ParallelIndexSet > &indexset)
Definition: globalaggregates.hh:40
TI ParallelIndexSet
Definition: globalaggregates.hh:30
Amg::AggregatesMap< T > Type
Definition: globalaggregates.hh:242
static void scatter(GlobalAggregatesMap< T, TI > &ga, GlobalIndex global, size_t i)
Definition: globalaggregates.hh:122
OwnerOverlapCopyCommunication< T1, T2 > ParallelInformation
Definition: globalaggregates.hh:162
const_iterator begin() const
Definition: aggregates.hh:723
SizeOne IndexedTypeFlag
Definition: globalaggregates.hh:244
Proxy(const GlobalLookupIndexSet< ParallelIndexSet > &indexset, Vertex &aggregate)
Definition: globalaggregates.hh:71
const_iterator end() const
Definition: aggregates.hh:728
static void publish(AggregatesMap< Vertex > &aggregates, ParallelInformation &pinfo, const GlobalLookupIndexSet &globalLookup)
Definition: globalaggregates.hh:166
ParallelIndexSet::GlobalIndex GlobalIndex
Definition: globalaggregates.hh:115
ParallelIndexSet::LocalIndex LocalIndex
Definition: globalaggregates.hh:36
TI ParallelIndexSet
Definition: globalaggregates.hh:114
T Vertex
Definition: globalaggregates.hh:221
const GlobalIndex & operator[](std::size_t index) const
Definition: globalaggregates.hh:45
ParallelInformation::GlobalLookupIndexSet GlobalLookupIndexSet
Definition: globalaggregates.hh:223
Proxy & operator=(const GlobalIndex &global)
Definition: globalaggregates.hh:75
ParallelIndexSet::GlobalIndex IndexedType
Definition: globalaggregates.hh:34
Proxy operator[](std::size_t index)
Definition: globalaggregates.hh:90
Definition: allocator.hh:7
A class setting up standard communication for a two-valued attribute set with owner/overlap/copy sema...
Definition: owneroverlapcopy.hh:172
const CollectiveCommunication< MPI_Comm > & communicator() const
Definition: owneroverlapcopy.hh:311
void copyCopyToAll(const T &source, T &dest) const
Communicate values from copy data points to all other data points.
Definition: owneroverlapcopy.hh:340
Dune::GlobalLookupIndexSet< ParallelIndexSet > GlobalLookupIndexSet
The type of the reverse lookup of indices.
Definition: owneroverlapcopy.hh:466
const RemoteIndices & remoteIndices() const
Get the underlying remote indices.
Definition: owneroverlapcopy.hh:481
void copyOwnerToAll(const T &source, T &dest) const
Communicate values from owner data points to all other data points.
Definition: owneroverlapcopy.hh:323
Dune::ParallelIndexSet< GlobalIdType, LI, 512 > ParallelIndexSet
The type of the parallel index set.
Definition: owneroverlapcopy.hh:459
Definition: globalaggregates.hh:28
Definition: globalaggregates.hh:69
Definition: globalaggregates.hh:113
Definition: globalaggregates.hh:130
int GlobalLookupIndexSet
Definition: pinfo.hh:58
@ nonoverlapping
Category for non-overlapping solvers.
Definition: solvercategory.hh:25
static Category category(const OP &op, decltype(op.category()) *=nullptr)
Helperfunction to extract the solver category either from an enum, or from the newly introduced virtu...
Definition: solvercategory.hh:32