int main( int argc, char **argv ) { int i; QString m1,n1,o1; struct args pargs; QApplication *a; InitWidget *b; ScreenWidget *c; TextWidget *textOut; QProgressDialog *qProg; QPixmap *qpxMeshIcon=NULL; struct model_error model1,model2; int rcode; struct outbuf *log; struct prog_reporter pr; /* Initialize application */ a = NULL; b = NULL; c = NULL; qProg = NULL; memset(&model1,0,sizeof(model1)); memset(&model2,0,sizeof(model2)); memset(&pr,0,sizeof(pr)); log = NULL; i = 0; while (i<argc) { if (strcmp(argv[i],"-t") == 0) /* text version requested */ break; if (strcmp(argv[i],"-h") == 0) /* just asked for command line help */ break; i++; } if (i == argc) { /* no text version requested, initialize QT */ a = new QApplication( argc, argv ); /* Load pixmap for icon */ qpxMeshIcon = new QPixmap((const char**)meshIcon); if (a != NULL) a->connect( a, SIGNAL(lastWindowClosed()), a, SLOT(quit()) ); } else { a = NULL; /* No QT app needed */ } /* Parse arguments */ parse_args(argc,argv,&pargs); /* Display starting dialog if insufficient arguments */ if (pargs.m1_fname != NULL || pargs.m2_fname != NULL) { if (pargs.m1_fname == NULL || pargs.m2_fname == NULL) { fprintf(stderr,"ERROR: missing file name(s) in command line\n"); exit(1); } if (!pargs.do_wlog) { log = outbuf_new(stdio_puts,stdout); } else { textOut = new TextWidget(); textOut->setIcon(*qpxMeshIcon); log = outbuf_new(TextWidget_puts,textOut); textOut->show(); } if (pargs.no_gui) { pr.prog = stdio_prog; pr.cb_out = stdout; } else { qProg = new QProgressDialog("Calculating distance",0,100); qProg->setIcon(*qpxMeshIcon); qProg->setMinimumDuration(1500); pr.prog = QT_prog; pr.cb_out = qProg; } mesh_run(&pargs, &model1, &model2, log, &pr); } else { b = new InitWidget(pargs, &model1, &model2); b->setIcon(*qpxMeshIcon); b->show(); } if (a != NULL) { if (pargs.m1_fname != NULL || pargs.m2_fname != NULL) { c = new ScreenWidget(&model1, &model2, &pargs); c->setIcon(*qpxMeshIcon); a->setMainWidget(c); c->show(); } rcode = a->exec(); } else { rcode = 0; } /* Free widgets */ outbuf_delete(log); delete qProg; delete qpxMeshIcon; delete b; delete c; delete a; // QApplication must be last QT thing to delete /* Free model data */ if (model1.mesh != NULL) __free_raw_model(model1.mesh); free(model1.verror); free(model1.info); free_face_error(model1.fe); if (model2.mesh != NULL) __free_raw_model(model2.mesh); free(model2.verror); free(model2.info); /* Return exit code */ return rcode; }
/* see mesh_run.h */ void mesh_run(const struct args *args, struct model_error *model1, struct model_error *model2, struct outbuf *out, struct prog_reporter *progress) { clock_t start_time; struct dist_surf_surf_stats stats; struct dist_surf_surf_stats stats_rev; double bbox1_diag,bbox2_diag; struct model_info *m1info,*m2info; double abs_sampling_step,abs_sampling_dens; int nv_empty,nf_empty; /* Read models from input files */ memset(model1,0,sizeof(*model1)); memset(model2,0,sizeof(*model2)); m1info = (struct model_info*) xa_malloc(sizeof(*m1info)); m2info = (struct model_info*) xa_malloc(sizeof(*m2info)); outbuf_printf(out,"Reading %s ... ",args->m1_fname); outbuf_flush(out); start_time = clock(); model1->mesh = read_model_file(args->m1_fname); outbuf_printf(out,"Done (%.2f secs)\n", (double)(clock()-start_time)/CLOCKS_PER_SEC); outbuf_printf(out,"Reading %s ... ",args->m2_fname); outbuf_flush(out); start_time = clock(); model2->mesh = read_model_file(args->m2_fname); outbuf_printf(out,"Done (%.2f secs)\n", (double)(clock()-start_time)/CLOCKS_PER_SEC); outbuf_flush(out); /* Analyze models (we don't need normals for model 1, so we don't request * for it to be oriented). */ start_time = clock(); bbox1_diag = dist_v(&model1->mesh->bBox[0], &model1->mesh->bBox[1]); bbox2_diag = dist_v(&model2->mesh->bBox[0], &model2->mesh->bBox[1]); analyze_model(model1->mesh,m1info,0,args->verb_analysis,out,"model 1"); model1->info = m1info; analyze_model(model2->mesh,m2info,1,args->verb_analysis,out,"model 2"); model2->info = m2info; /* Adjust sampling step size */ abs_sampling_step = args->sampling_step*bbox2_diag; abs_sampling_dens = 1/(abs_sampling_step*abs_sampling_step); /* Print available model information */ outbuf_printf(out,"\n Model information\n" " (degenerate faces ignored for manifold/closed info)\n\n"); outbuf_printf(out,"Number of vertices: \t%11d\t%11d\n", model1->mesh->num_vert,model2->mesh->num_vert); outbuf_printf(out,"Number of triangles: \t%11d\t%11d\n", model1->mesh->num_faces,model2->mesh->num_faces); outbuf_printf(out,"Degenerate triangles: \t%11d\t%11d\n", m1info->n_degenerate,m2info->n_degenerate); outbuf_printf(out,"BoundingBox diagonal: \t%11g\t%11g\n", bbox1_diag,bbox2_diag); outbuf_printf(out,"Number of disjoint parts:\t%11d\t%11d\n", m1info->n_disjoint_parts,m2info->n_disjoint_parts); outbuf_printf(out,"Manifold: \t%11s\t%11s\n", (m1info->manifold ? "yes" : "no"), (m2info->manifold ? "yes" : "no")); outbuf_printf(out,"Originally oriented: \t%11s\t%11s\n", (m1info->orig_oriented ? "yes" : "no"), (m2info->orig_oriented ? "yes" : "no")); outbuf_printf(out,"Orientable: \t%11s\t%11s\n", (m1info->orientable ? "yes" : "no"), (m2info->orientable ? "yes" : "no")); outbuf_printf(out,"Closed: \t%11s\t%11s\n", (m1info->closed ? "yes" : "no"), (m2info->closed ? "yes" : "no")); outbuf_flush(out); /* Compute the distance from one model to the other */ dist_surf_surf(model1,model2->mesh,abs_sampling_dens,args->min_sample_freq, &stats,!args->no_gui,(args->quiet ? NULL : progress)); /* Print results */ outbuf_printf(out,"Surface area: \t%11g\t%11g\n", stats.m1_area,stats.m2_area); outbuf_printf(out,"\n Distance from model 1 to model 2\n\n"); outbuf_printf(out," \t Absolute\t%% BBox diag\n"); outbuf_printf(out," \t \t (Model 2)\n"); outbuf_printf(out,"Min: \t%11g\t%11g\n", stats.min_dist,stats.min_dist/bbox2_diag*100); outbuf_printf(out,"Max: \t%11g\t%11g\n", stats.max_dist,stats.max_dist/bbox2_diag*100); outbuf_printf(out,"Mean: \t%11g\t%11g\n", stats.mean_dist,stats.mean_dist/bbox2_diag*100); outbuf_printf(out,"RMS: \t%11g\t%11g\n", stats.rms_dist,stats.rms_dist/bbox2_diag*100); outbuf_printf(out,"\n"); outbuf_flush(out); if (args->do_symmetric) { /* Invert models and recompute distance */ outbuf_printf(out," Distance from model 2 to model 1\n\n"); dist_surf_surf(model2,model1->mesh,abs_sampling_dens,args->min_sample_freq, &stats_rev,0,(args->quiet ? NULL : progress)); free_face_error(model2->fe); model2->fe = NULL; outbuf_printf(out," \t Absolute\t%% BBox diag\n"); outbuf_printf(out," \t \t (Model 2)\n"); outbuf_printf(out,"Min: \t%11g\t%11g\n", stats_rev.min_dist,stats_rev.min_dist/bbox2_diag*100); outbuf_printf(out,"Max: \t%11g\t%11g\n", stats_rev.max_dist,stats_rev.max_dist/bbox2_diag*100); outbuf_printf(out,"Mean: \t%11g\t%11g\n", stats_rev.mean_dist,stats_rev.mean_dist/bbox2_diag*100); outbuf_printf(out,"RMS: \t%11g\t%11g\n", stats_rev.rms_dist,stats_rev.rms_dist/bbox2_diag*100); outbuf_printf(out,"\n"); /* Print symmetric distance measures */ outbuf_printf(out, " Symmetric distance between model 1 and model 2\n\n"); outbuf_printf(out," \t Absolute\t%% BBox diag\n"); outbuf_printf(out," \t \t (Model 2)\n"); outbuf_printf(out,"Min: \t%11g\t%11g\n", max(stats.min_dist,stats_rev.min_dist), max(stats.min_dist,stats_rev.min_dist)/bbox2_diag*100); outbuf_printf(out,"Max: \t%11g\t%11g\n", max(stats.max_dist,stats_rev.max_dist), max(stats.max_dist,stats_rev.max_dist)/bbox2_diag*100); outbuf_printf(out,"Mean: \t%11g\t%11g\n", max(stats.mean_dist,stats_rev.mean_dist), max(stats.mean_dist,stats_rev.mean_dist)/bbox2_diag*100); outbuf_printf(out,"RMS: \t%11g\t%11g\n", max(stats.rms_dist,stats_rev.rms_dist), max(stats.rms_dist,stats_rev.rms_dist)/bbox2_diag*100); outbuf_printf(out,"\n"); } outbuf_printf(out," \tAbsolute\t %% BBox diag\t " "Expected samples\n" " \t \t model 2 \t " "model 1\tmodel 2\n"); if (!args->do_symmetric) { outbuf_printf(out,"Sampling step: \t%8g\t %7g \t %7d\t%7d\n", abs_sampling_step,abs_sampling_step/bbox2_diag*100, (int)(stats.m1_area*abs_sampling_dens),0); outbuf_printf(out,"\n"); outbuf_printf(out," \t Total\t Avg. / triangle\t\t" "Tot (%%) area of\n" " \t \tmodel 1\tmodel 2 \t\t" "sampled triang.\n"); outbuf_printf(out,"Samples:\t%9d\t%7.2g\t%7.2g\t\t%15.2f\n",stats.m1_samples, ((double)stats.m1_samples)/model1->mesh->num_faces, ((double)stats.m1_samples)/model2->mesh->num_faces, stats.st_m1_area/stats.m1_area*100.0); } else { outbuf_printf(out,"Sampling step: \t%8g\t %7g \t %7d\t%7d\n", abs_sampling_step,abs_sampling_step/bbox2_diag*100, (int)(stats.m1_area*abs_sampling_dens), (int)(stats.m2_area*abs_sampling_dens)); outbuf_printf(out,"\n"); outbuf_printf(out," \t Total\t Avg. / triangle\t\t" "Tot (%%) area of\n" " \t \tmodel 1 \tmodel 2 \t" "sampled triang.\n"); outbuf_printf(out,"Samples (1->2):\t%9d\t%7.2g\t%15.2g\t%18.2f\n", stats.m1_samples, ((double)stats.m1_samples)/model1->mesh->num_faces, ((double)stats.m1_samples)/model2->mesh->num_faces, stats.st_m1_area/stats.m1_area*100.0); outbuf_printf(out,"Samples (2->1):\t%9d\t%7.2g\t%15.2g\t%18.2f\n", stats_rev.m1_samples, ((double)stats_rev.m1_samples)/model1->mesh->num_faces, ((double)stats_rev.m1_samples)/model2->mesh->num_faces, stats_rev.st_m1_area/stats_rev.m1_area*100.0); } outbuf_printf(out,"\n"); if (!args->do_symmetric) { outbuf_printf(out, " \t X\t Y\t Z\t Total\n"); outbuf_printf(out,"Partitioning grid size:\t%6d\t%5d\t%4d\t%8d\n", stats.grid_sz.x,stats.grid_sz.y,stats.grid_sz.z, stats.grid_sz.x*stats.grid_sz.y*stats.grid_sz.z); outbuf_printf(out,"\nAvg. number of triangles per non-empty cell:\t%.2f\n", stats.n_t_p_nec); outbuf_printf(out,"Proportion of non-empty cells: \t%.2f%%\n", (double)stats.n_ne_cells/(stats.grid_sz.x*stats.grid_sz.y* stats.grid_sz.z)*100.0); } else { outbuf_printf(out," \t " "X\t Y\t Z\t Total\n"); outbuf_printf(out,"Partitioning grid size (1 to 2):\t%6d\t%5d\t%4d\t%8d\n", stats.grid_sz.x,stats.grid_sz.y,stats.grid_sz.z, stats.grid_sz.x*stats.grid_sz.y*stats.grid_sz.z); outbuf_printf(out,"Partitioning grid size (2 to 1):\t%6d\t%5d\t%4d\t%8d\n", stats_rev.grid_sz.x,stats_rev.grid_sz.y,stats_rev.grid_sz.z, stats_rev.grid_sz.x*stats_rev.grid_sz.y*stats_rev.grid_sz.z); outbuf_printf(out,"\nAvg. number of triangles per non-empty cell (1 to 2):" "\t%.2f\n",stats.n_t_p_nec); outbuf_printf(out,"Avg. number of triangles per non-empty cell (2 to 1):" "\t%.2f\n",stats_rev.n_t_p_nec); outbuf_printf(out, "Proportion of non-empty cells (1 to 2): \t%.2f%%\n", (double)stats.n_ne_cells/(stats.grid_sz.x*stats.grid_sz.y* stats.grid_sz.z)*100.0); outbuf_printf(out, "Proportion of non-empty cells (2 to 1): \t%.2f%%\n", (double)stats_rev.n_ne_cells/ (stats_rev.grid_sz.x*stats_rev.grid_sz.y* stats_rev.grid_sz.z)*100.0); } outbuf_printf(out,"\n"); outbuf_printf(out,"Analysis and measuring time (secs.):\t%.2f\n", (double)(clock()-start_time)/CLOCKS_PER_SEC); outbuf_flush(out); if(!args->no_gui){ /* Get the per vertex error metric */ nv_empty = nf_empty = 0; /* keep compiler happy */ calc_vertex_error(model1,&nv_empty,&nf_empty); if (nv_empty>0) { outbuf_printf(out, "WARNING: %.2f%% of vertices (%i out of %i) have no error " "samples\n",100.0*nv_empty/model1->mesh->num_vert, nv_empty,model1->mesh->num_vert); } if (nf_empty>0) { outbuf_printf(out, "WARNING: %.2f%% of faces (%i out of %i) have no error " "samples\n",100.0*nf_empty/model1->mesh->num_faces, nf_empty,model1->mesh->num_faces); } outbuf_flush(out); } }