1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
|
// -*- tab-width: 4; indent-tabs-mode: nil; c-basic-offset: 2 -*-
// vi: set et ts=4 sw=2 sts=2:
// SPDX-FileCopyrightInfo: Copyright © DUNE Project contributors, see file LICENSE.md in module root
// SPDX-License-Identifier: LicenseRef-GPL-2.0-only-with-DUNE-exception
#include <iostream>
#include <dune/common/enumset.hh>
#include <dune/common/timer.hh>
#include <dune/common/parallel/indexset.hh>
#include <dune/common/parallel/plocalindex.hh>
#include <dune/common/parallel/selection.hh>
enum GridFlags {
owner, overlap, border
};
template<class T>
int measure(const T& selection)
{
/*
return measure<1>(selection);
}
template<int LOOPS, class T>
int measure(const T& selection)
{*/
typedef typename T::const_iterator iterator;
const iterator end = selection.end();
int count=0;
Dune::Timer timer;
timer.reset();
for(int i=0; i<10; i++)
for(iterator iter = selection.begin(); iter != end; ++iter)
count+=*iter;
std::cout<<" took "<< timer.elapsed()<<" seconds"<<std::endl;
return count;
}
template<int SIZE>
void test()
{
const int Nx = SIZE;
const int Ny = SIZE;
// Process configuration
const int ALSIZE=55;
Dune::ParallelIndexSet<int,Dune::ParallelLocalIndex<GridFlags>,ALSIZE> distIndexSet;
distIndexSet.beginResize();
for(int y=0, i=0; y < Ny; y++)
for(int x=0; x < Nx; x++, i++) {
GridFlags flag = owner;
if(x==0 || x == Nx-1 || y ==0 || y==Ny-1)
flag = overlap;
distIndexSet.add(i, Dune::ParallelLocalIndex<GridFlags> (i, flag, true));
}
distIndexSet.endResize();
Dune::UncachedSelection<Dune::EnumItem<GridFlags,owner>,int,Dune::ParallelLocalIndex<GridFlags>,ALSIZE>
ownerUncached(distIndexSet);
Dune::Selection<Dune::EnumItem<GridFlags,owner>,int,Dune::ParallelLocalIndex<GridFlags>,ALSIZE>
ownerCached(distIndexSet);
Dune::UncachedSelection<Dune::EnumItem<GridFlags,overlap>,int,Dune::ParallelLocalIndex<GridFlags>,ALSIZE>
overlapUncached(distIndexSet);
Dune::Selection<Dune::EnumItem<GridFlags,overlap>,int,Dune::ParallelLocalIndex<GridFlags>,ALSIZE>
overlapCached(distIndexSet);
int count=0;
std::cout<<" Owner selection uncached:";
count+=measure(ownerUncached);
std::cout<<" Owner selection cached:";
count+=measure(ownerCached);
std::cout<<" Overlap selection uncached:";
count+=measure(overlapUncached);
std::cout<<" Overlap selection cached:";
count+=measure(overlapCached);
std::cout<<count<<std::endl;
}
int main()
{
test<1000>();
}
|
