| author | marci |
| Sat, 24 Jul 2004 14:01:31 +0000 | |
| changeset 738 | 56e60e9eb2da |
| parent 711 | b6c56353832c |
| child 742 | 235fd36336b7 |
| permissions | -rw-r--r-- |
| alpar@711 | 1 |
// -*- mode:C++ -*- |
| alpar@711 | 2 |
#ifndef HUGO_BENCH_TEST_H |
| alpar@711 | 3 |
#define HUGO_BENCH_TEST_H |
| alpar@711 | 4 |
|
| alpar@711 | 5 |
#include<vector> |
| alpar@718 | 6 |
#include<iostream> |
| alpar@718 | 7 |
|
| alpar@718 | 8 |
#include<hugo/time_measure.h> |
| alpar@711 | 9 |
|
| alpar@711 | 10 |
///An experimental typedef factory |
| alpar@711 | 11 |
#define GRAPH_TYPEDEF_FACTORY(Graph) \ |
| alpar@711 | 12 |
typedef typename Graph:: Node Node;\ |
| alpar@711 | 13 |
typedef typename Graph:: NodeIt NodeIn;\ |
| alpar@711 | 14 |
typedef typename Graph:: Edge Edge;\ |
| alpar@711 | 15 |
typedef typename Graph:: EdgeIt EdgeIt;\ |
| alpar@711 | 16 |
typedef typename Graph:: InEdgeIt InEdgeIt;\ |
| alpar@711 | 17 |
typedef typename Graph::OutEdgeIt OutEdgeIt; |
| alpar@711 | 18 |
|
| alpar@711 | 19 |
#define GRAPH_TYPEDEF_FACTORY_NOTYPENAME(Graph) \ |
| alpar@711 | 20 |
typedef Graph:: Node Node;\ |
| alpar@711 | 21 |
typedef Graph:: NodeIt NodeIn;\ |
| alpar@711 | 22 |
typedef Graph:: Edge Edge;\ |
| alpar@711 | 23 |
typedef Graph:: EdgeIt EdgeIt;\ |
| alpar@711 | 24 |
typedef Graph:: InEdgeIt InEdgeIt;\ |
| alpar@711 | 25 |
typedef Graph::OutEdgeIt OutEdgeIt; |
| alpar@711 | 26 |
|
| alpar@711 | 27 |
|
| alpar@711 | 28 |
///A primitive primtest |
| alpar@718 | 29 |
|
| alpar@718 | 30 |
///\bug 2 is not a prime according to this function! |
| alpar@711 | 31 |
bool isPrim(int n) |
| alpar@711 | 32 |
{
|
| alpar@711 | 33 |
if(n%2) {
|
| alpar@711 | 34 |
for(int k=3;n/k>=k;k+=2) |
| alpar@711 | 35 |
if(!(n%k)) return false; |
| alpar@711 | 36 |
return true; |
| alpar@711 | 37 |
} |
| alpar@711 | 38 |
return false; |
| alpar@711 | 39 |
} |
| alpar@711 | 40 |
|
| alpar@711 | 41 |
///Finds the smallest prime not less then \c n. |
| alpar@711 | 42 |
int nextPrim(int n) |
| alpar@711 | 43 |
{
|
| alpar@711 | 44 |
for(n+=!(n%2);!isPrim(n);n+=2) ; |
| alpar@711 | 45 |
return n; |
| alpar@711 | 46 |
} |
| alpar@711 | 47 |
|
| alpar@711 | 48 |
|
| alpar@711 | 49 |
/// Class to generate consecutive primes |
| alpar@711 | 50 |
class Primes |
| alpar@711 | 51 |
{
|
| alpar@711 | 52 |
std::vector<int> primes; |
| alpar@711 | 53 |
int n; |
| alpar@711 | 54 |
|
| alpar@711 | 55 |
bool isPrime(int m) |
| alpar@711 | 56 |
{
|
| alpar@711 | 57 |
for(int i=0;m<primes[i]*primes[i];i++) if(!(m%primes[i])) return false; |
| alpar@711 | 58 |
return true; |
| alpar@711 | 59 |
} |
| alpar@711 | 60 |
public: |
| alpar@711 | 61 |
Primes() : n(1) {}
|
| alpar@711 | 62 |
|
| alpar@711 | 63 |
int operator() () |
| alpar@711 | 64 |
{
|
| alpar@711 | 65 |
if(primes.size()==0) {
|
| alpar@711 | 66 |
primes.push_back(2); |
| alpar@711 | 67 |
return 2; |
| alpar@711 | 68 |
} |
| alpar@711 | 69 |
else {
|
| alpar@711 | 70 |
do n+=2; while(!isPrime(n)); |
| alpar@711 | 71 |
primes.push_back(n); |
| alpar@711 | 72 |
return n; |
| alpar@711 | 73 |
} |
| alpar@711 | 74 |
} |
| alpar@711 | 75 |
}; |
| alpar@711 | 76 |
|
| alpar@718 | 77 |
inline void PrintTime(char *ID,hugo::Timer &T) |
| alpar@718 | 78 |
{
|
| alpar@718 | 79 |
hugo::TimeStamp S(T); |
| alpar@718 | 80 |
std::cout << ID << ' ' << S.getUserTime() << ' ' |
| alpar@718 | 81 |
<< S.getSystemTime() << ' ' << S.getRealTime() << std::endl; |
| alpar@718 | 82 |
} |
| alpar@718 | 83 |
|
| alpar@718 | 84 |
|
| alpar@711 | 85 |
|
| alpar@711 | 86 |
#endif |