testMetrics-benchmark.cpp

Go to the documentation of this file.

 
#include <iostream>
#include "DGtal/base/Common.h"
#include "DGtal/helpers/StdDefs.h"
#include "DGtal/geometry/volumes/distance/ExactPredicateLpSeparableMetric.h"
#include "DGtal/geometry/volumes/distance/InexactPredicateLpSeparableMetric.h"
#include <boost/lexical_cast.hpp>
 
using namespace std;
using namespace DGtal;
 
// Functions for testing class Metrics-benchmark.
 
template<int p>
bool runATest( unsigned int maxTest)
{
  ExactPredicateLpSeparableMetric<Z3i::Space, p> exactMetric;
  InexactPredicateLpSeparableMetric<Z3i::Space> approxMetric(p);
  double sum=0.0, d=0,d2=0, MSE=0;
  std::string txt = "Testing Exponent" +   boost::lexical_cast<string>(p);
  
    
  DGtal::Dimension dim = 1000;
  trace.beginBlock(txt);
  trace.beginBlock("Exact Predicate run");
  for(unsigned int i=0; i< maxTest; ++i)
    {
      Z3i::Point P( rand() % dim,rand() % dim ,rand() % dim );
      Z3i::Point Q( rand() % dim, rand() % dim, rand() % dim);
      sum += exactMetric(P,Q);
    }
 
  trace.endBlock();
  
  trace.beginBlock("Inexact Predicate run");
 for(unsigned int i=0; i< maxTest; ++i)
    {
      Z3i::Point P( rand() % dim,rand() % dim ,rand() % dim );
      Z3i::Point Q( rand() % dim, rand() % dim, rand() % dim);
      sum += approxMetric(P,Q);
    }
  trace.endBlock();
  
 
  trace.beginBlock("Both with comparison run");
  for(unsigned int i=0; i< maxTest; ++i)
    {
      Z3i::Point P( rand() % dim,rand() % dim ,rand() % dim );
      Z3i::Point Q( rand() % dim, rand() % dim, rand() % dim);
      d = exactMetric(P,Q);
      d2 = approxMetric(P,Q);
      MSE += (d-d2)*(d-d2);
    }
  trace.endBlock();
  trace.info() << "Mean Square Error = "<<MSE<<std::endl;
  trace.endBlock();
  return true;
}
 
 
 
// Standard services - public :
 
int main( int argc, char** argv )
{
  trace.beginBlock ( "Testing class Metrics-benchmark" );
  trace.info() << "Args:";
  for ( int i = 0; i < argc; ++i )
    trace.info() << " " << argv[ i ];
  trace.info() << endl;
 
  bool res = runATest<2>(1000000) &&
    runATest<4>(10000000);
    runATest<6>(1000000);
  trace.emphase() << ( res ? "Passed." : "Error." ) << endl;
  trace.endBlock();
  return res ? 0 : 1;
}
//                                                                           //