int main(int argc, char *argv[]) {
if(argc < 2) {
std::cerr << "ERROR: input file required.\nusage: ./geonum_TP9 [dataname] [numsteps=1]\n";
return -1;
}
const std::string dataname(argv[1]);
const std::string in_dir = "../data/";
const std::string in_filename = in_dir+dataname+".off";
const std::string out_dir = "../plots/";
const std::string out_filename = out_dir+dataname+".off";
// mesh data : vertices and faces
MatdX3 V0, V1;
MatiX3 F0, F1;
// read mesh from OFF file
if( !readOFF( in_filename, V0, F0 )) {
std::cerr << "ERROR: cannot read file " << in_filename << std::endl;
return -1;
}
unsigned int num_steps = 1;
if( argc > 2 )
num_steps = atoi(argv[2]);
for( int s=1; s<=num_steps; s++ ){
std::cout << "\nLoop step " << s << "...";
assert(V0.rows() < 20000); // limit
begin = std::clock();
LoopStep(V0,F0,V1,F1);
end = std::clock();
V0 = V1;
F0 = F1;
std::cout <<
" [" << double(end - begin) / CLOCKS_PER_SEC << "s]\n"
" #V" << s << " = " << V0.rows() << std::endl <<
" #F" << s << " = " << F0.rows() << std::endl;
}
std::cout << "\n" << std::endl;
// write mesh to OFF file
if( !writeOFF( out_filename, V0, F0 )) {
std::cerr << "ERROR: cannot write file " << out_filename << std::endl;
return -1;
}
//*******************//
//** SimpleViewer **//
//*****************//
SimpleViewer viewer;
viewer.set_mesh(V0,F0);
viewer.set_facecolor( 255, 163, 0 ); // [255,163, 0] = dark orange
viewer.set_edgecolor( 55, 55, 55 ); // [ 55, 55, 55] = gray
return viewer.show();
}
/**
* writes a surface as a given format
*/
int writeSurface(surface *surf, char *filename, enum SURFACE_FORMAT format) {
switch(format) {
case OFF:
return writeOFF(surf,filename);
break;
case OBJ:
return writeOBJ(surf,filename);
break;
case MGHSURF:
return writeMGHsurf(surf,filename);
break;
default:
fprintf(stderr,"writeSurface(): unsupported format\n");
break;
}
return SR_FAILURE;
}
/**
* main
*/
int main(int argc, char **args) {
list *slices;
surface *surf;
/* turn on memory tracing */
putenv("MALLOC_TRACE=memtest_memtrace.txt");
mtrace();
/* read the slices from the file */
slices = readSliceContours(args[1]);
if(slices == NULL) {
printf("couldn't read slices from %s!\n",args[1]);
return 0;
}
preprocessSliceContours(slices);
/* builds a guess at contour correspondence */
// correspondenceMethod = HISTOGRAM;
// correspondenceScope = LOCAL;
buildCorrespondenceGuess(slices);
/* tile the slices */
// minSliceInd = 0;
// maxSliceInd = 40;
surf = tileSlices(slices);
deleteSliceContours(slices);
fillBranchedHoles(surf);
writeOFF(surf,"testmem.off");
deleteSurface(surf);
return 0;
}
