generateSimplicityTables3D.cpp

#include <iostream>
#include <vector>
#include <bitset>
#include "DGtal/topology/Object.h"
#include "DGtal/helpers/StdDefs.h"


using namespace std;
using namespace DGtal;


template <typename Object, typename Map>
void 
generateSimplicityTable( const typename Object::DigitalTopology & dt,
                         Map & map )
{
  typedef typename Object::DigitalTopology DigitalTopology;
  typedef typename Object::DigitalSet DigitalSet;
  typedef typename Object::Point Point;
  typedef typename DigitalSet::Domain Domain;
  typedef typename Domain::ConstIterator DomainConstIterator;

  Point p1 = Point::diagonal( -1 );
  Point p2 = Point::diagonal(  1 );
  Point c = Point::diagonal( 0 );
  Domain domain( p1, p2 );
  DigitalSet shapeSet( domain );
  Object shape( dt, shapeSet );
  unsigned int k = 0;
  for ( DomainConstIterator it = domain.begin(); it != domain.end(); ++it )
    if ( *it != c ) ++k;
  ASSERT( ( k < 32 )
          && "[generateSimplicityTable] number of configurations is too high." );
  unsigned int nbCfg = 1 << k;
  for ( unsigned int cfg = 0; cfg < nbCfg; ++cfg )
    {
      if ( ( cfg % 1000 ) == 0 )
        {
          trace.progressBar( (double) cfg, (double) nbCfg );
        }
      shape.pointSet().clear();
      shape.pointSet().insert( c );
      unsigned int mask = 1;
      for ( DomainConstIterator it = domain.begin(); it != domain.end(); ++it )
        {
          if ( *it != c )
            {
              if ( cfg & mask ) shape.pointSet().insert( *it );
              mask <<= 1;
            }
        }
      bool simple = shape.isSimple( c );
      map[ cfg ] = simple;
    }
}



int main( int argc, char** argv )
{
  typedef std::bitset<67108864> ConfigMap; // 2^26

  using namespace Z3i;
  trace.beginBlock ( "Generate 3d table for 26_6 topology" );
  // Too big for stack. Use heap instead.
  ConfigMap* table26_6 = new ConfigMap;
  generateSimplicityTable< Object26_6 >( dt26_6, *table26_6 );
  trace.endBlock();

  ofstream file26_6( "simplicity_table26_6.txt" );
  file26_6 << *table26_6;
  file26_6.close();

  delete table26_6;
  return 0;
}
//                                                                           //