alpar@711: // -*- mode:C++ -*-
alpar@711: #ifndef HUGO_BENCH_TEST_H
alpar@711: #define HUGO_BENCH_TEST_H
alpar@711: 
alpar@711: #include<vector>
alpar@718: #include<iostream>
alpar@718: 
alpar@718: #include<hugo/time_measure.h>
alpar@711: 
alpar@711: ///An experimental typedef factory
alpar@711: #define GRAPH_TYPEDEF_FACTORY(Graph) \
alpar@711:    typedef typename Graph::   Node      Node;\
alpar@711:    typedef typename Graph::   NodeIt    NodeIn;\
alpar@711:    typedef typename Graph::   Edge      Edge;\
alpar@711:    typedef typename Graph::   EdgeIt    EdgeIt;\
alpar@711:    typedef typename Graph:: InEdgeIt  InEdgeIt;\
alpar@711:    typedef typename Graph::OutEdgeIt OutEdgeIt;
alpar@711: 
alpar@711: #define GRAPH_TYPEDEF_FACTORY_NOTYPENAME(Graph) \
alpar@711:    typedef Graph::   Node      Node;\
alpar@711:    typedef Graph::   NodeIt    NodeIn;\
alpar@711:    typedef Graph::   Edge      Edge;\
alpar@711:    typedef Graph::   EdgeIt    EdgeIt;\
alpar@711:    typedef Graph:: InEdgeIt  InEdgeIt;\
alpar@711:    typedef Graph::OutEdgeIt OutEdgeIt;
alpar@711:  
alpar@711: 
alpar@711: ///A primitive primtest
alpar@718: 
alpar@718: ///\bug 2 is not a prime according to this function!
alpar@711: bool isPrim(int n)
alpar@711: {
alpar@711:   if(n%2) {
alpar@711:     for(int k=3;n/k>=k;k+=2)
alpar@711:       if(!(n%k)) return false;
alpar@711:     return true;
alpar@711:   }
alpar@711:   return false;
alpar@711: }
alpar@711: 
alpar@711: ///Finds the smallest prime not less then \c n.
alpar@711: int nextPrim(int n)
alpar@711: {
alpar@711:   for(n+=!(n%2);!isPrim(n);n+=2) ;
alpar@711:   return n;
alpar@711: }
alpar@711: 
alpar@711: 
alpar@711: /// Class to generate consecutive primes
alpar@711: class Primes 
alpar@711: {
alpar@711:   std::vector<int> primes;
alpar@711:   int n;
alpar@711:   
alpar@711:   bool isPrime(int m) 
alpar@711:   {
alpar@711:     for(int i=0;m<primes[i]*primes[i];i++) if(!(m%primes[i])) return false;
alpar@711:     return true;
alpar@711:   }
alpar@711: public:
alpar@711:   Primes() : n(1) {}
alpar@711:   
alpar@711:   int operator() ()
alpar@711:     {
alpar@711:       if(primes.size()==0) {
alpar@711: 	primes.push_back(2);
alpar@711: 	return 2;
alpar@711:       }
alpar@711:       else {
alpar@711: 	do n+=2; while(!isPrime(n));
alpar@711: 	primes.push_back(n);
alpar@711: 	return n;
alpar@711:       }
alpar@711:     }
alpar@711: };
alpar@711: 
alpar@718: inline void PrintTime(char *ID,hugo::Timer &T) 
alpar@718: {
alpar@718:   hugo::TimeStamp S(T);
alpar@718:   std::cout << ID << ' ' << S.getUserTime() << ' '
alpar@718: 	    << S.getSystemTime() << ' ' << S.getRealTime() << std::endl;
alpar@718: }
alpar@718: 
alpar@718:   
alpar@711: 
alpar@711: #endif