int KMeansClusterizer::getMinimumClusterID(ClusterizerData* train_data_ti)
{
vd dist_v(getCenterNum());
for (size_t ci = 0; ci < dist_v.sz; ci++) {
auto center_ci = m_center.at(ci);
dist_v.at(ci) = train_data_ti->calcDistance(center_ci);
}
return (int)(min_element(all(dist_v)) - begin(dist_v));
}
/* 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);
}
}
