/* see reporting.h */
struct outbuf *outbuf_new(outbuf_flush_cb_t *flush_cb, void *cb_out)
{
struct outbuf *ob;
ob = xa_malloc(sizeof(*ob));
ob->strbuf = xa_malloc(OUTBUF_MIN_FREE);
ob->strbuf[0] = '\0';
ob->pos = ob->strbuf;
ob->end = ob->strbuf + OUTBUF_MIN_FREE;
ob->flush = flush_cb;
ob->cb_out = cb_out;
return ob;
}
// Generate the texture for error display
void Error3DViewerWidget::genErrorTextures()
{
static const GLint internalformat = GL_R3_G3_B2;
GLubyte *texture;
GLint tw,max_n;
QString tmps;
int i;
// Only in error mapping mode and if not disabled
if (!texture_enabled)
return;
makeCurrent(); // make sure we use the correct GL context
// Allocate texture names (IDs) if not present
if (etex_id == NULL) {
etex_id = (GLuint*) xa_malloc(sizeof(*etex_id)*model->mesh->num_faces);
etex_sz = (int*) xa_malloc(sizeof(*etex_sz)*model->mesh->num_faces);
glGenTextures(model->mesh->num_faces,etex_id);
}
// Get maximum texture size
max_n = 0;
for (i=0; i<model->mesh->num_faces; i++) {
if (max_n < model->fe[i].sample_freq) max_n = model->fe[i].sample_freq;
}
max_n = 1<<ceilLog2(max_n); // round (towards infinity) to power of two
// Test if OpenGL implementation can deal with maximum texture size
// Unfortunately GL_PROXY_TEXTURE_2D fails on IRIX 6.2 for some SGI
// machines, so use older GL_MAX_TEXTURE_SIZE method.
glGetIntegerv(GL_MAX_TEXTURE_SIZE,&tw);
checkGLErrors("error texture size check");
if (tw < max_n) {
tmps.sprintf("The OpenGL implementation does not support\n"
"the required texture size (%ix%i).\n"
"Using plain white color",max_n,max_n);
QMessageBox::critical(this,"OpenGL texture size exceeded",tmps);
// Displaying another window can change the current GL context
makeCurrent();
for (i=0; i<model->mesh->num_faces; i++) {
etex_sz[i] = 1; // avoid having divide by zero texture coords
}
return;
}
// What follows is a potentially slow operation
QApplication::setOverrideCursor(Qt::waitCursor);
// Allocate temporary texture storage
texture = (GLubyte*) xa_malloc(sizeof(*texture)*3*(max_n+2)*(max_n+2));
glPixelStorei(GL_UNPACK_ALIGNMENT,1); /* pixel rows aligned on bytes only */
for (i=0; i<model->mesh->num_faces; i++) {
glBindTexture(GL_TEXTURE_2D,etex_id[i]);
etex_sz[i] = fillTexture(&(model->fe[i]),texture);
glTexImage2D(GL_TEXTURE_2D,0,internalformat,etex_sz[i]+2,etex_sz[i]+2,1,
GL_RGB,GL_UNSIGNED_BYTE,texture);
glTexParameteri(GL_TEXTURE_2D,GL_TEXTURE_WRAP_S,GL_CLAMP);
glTexParameteri(GL_TEXTURE_2D,GL_TEXTURE_WRAP_T,GL_CLAMP);
// Default GL_TEXTURE_MIN_FILTER requires mipmaps!
glTexParameteri(GL_TEXTURE_2D,GL_TEXTURE_MIN_FILTER,GL_LINEAR);
glTexParameteri(GL_TEXTURE_2D,GL_TEXTURE_MAG_FILTER,GL_LINEAR);
}
checkGLErrors("error texture generation");
free(texture);
QApplication::restoreOverrideCursor();
}
/* 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);
}
}
