testMeshReader.cpp

#include <iostream>
#include "DGtal/base/Common.h"
#include "DGtal/shapes/fromPoints/MeshFromPoints.h"
#include "DGtal/io/readers/MeshReader.h"
#include "DGtal/helpers/StdDefs.h"

#include "ConfigTest.h"


using namespace std;
using namespace DGtal;



struct Point3D{  
  const double & operator[]( unsigned int i ) const{
    assert(i<3);
    switch (i){
    case 0: {return x;}
    case 1: {return y;}
    case 2: {return z;}
    }
    return x;
  };

   double & operator[]( unsigned int i ) {
    assert(i<3);
    switch (i){
    case 0: {return x;}
    case 1: {return y;}
    case 2: {return z;}
    }
    return x;
  };
  
  double  x, y,z;
};


typedef Point3D Point;




 



// Functions for testing class MeshReader.

bool testMeshReader()
{
  unsigned int nbok = 0;
  unsigned int nb = 0;
  trace.beginBlock ( "Testing block ..." );  
  nb++;
  std::string filenameOFF = testPath + "samples/box.off";  
  MeshFromPoints<Point> a3DMesh;
  bool importOK = a3DMesh << filenameOFF;
  nbok += importOK ? 1 : 0; 
  
  
  nb++;
  MeshFromPoints<Point>::MeshFace aFace = a3DMesh.getFace(0);
  bool isWellImported = (a3DMesh.nbVertex()==8) &&  (a3DMesh.nbFaces()==6) && (aFace.size()==4) && (aFace.at(0)==0);
  nbok+=isWellImported? 1: 0; 
  

  trace.info() << "(" << nbok << "/" << nb << ") "
               << "true == true" << std::endl;


  nb++;
  std::string filenameOFS = testPath + "samples/testMesh.ofs";  
  MeshFromPoints<Point> a3DMesh2;
  bool importOK2=  a3DMesh2 << filenameOFS;
  nbok += importOK2 ? 1 : 0; 
  
  nb++;
  MeshFromPoints<Point>::MeshFace aFace2 = a3DMesh2.getFace(0);
  bool isWellImported2 = (a3DMesh2.nbVertex()==32) &&  (a3DMesh2.nbFaces()==60) && (aFace2.size()==3) && (aFace2.at(0)==0);
  nbok+=isWellImported2? 1: 0;


  trace.info() << "(" << nbok << "/" << nb << ") "
               << "true == true" << std::endl;
  trace.endBlock();  

  
  

  

  
  
  return nbok == nb;
}

// Standard services - public :

int main( int argc, char** argv )
{
  trace.beginBlock ( "Testing class MeshReader" );
  trace.info() << "Args:";
  for ( int i = 0; i < argc; ++i )
    trace.info() << " " << argv[ i ];
  trace.info() << endl;

  bool res = testMeshReader(); // && ... other tests
  trace.emphase() << ( res ? "Passed." : "Error." ) << endl;
  trace.endBlock();
  return res ? 0 : 1;
}
//                                                                           //