Revision: 42127
http://brlcad.svn.sourceforge.net/brlcad/?rev=42127&view=rev
Author: brlcad
Date: 2011-01-12 03:32:55 +0000 (Wed, 12 Jan 2011)
Log Message:
-----------
massive restructure in order to fix a bug assuming that variables are
accessible after a longjmp. quite a chore to quell the warning due to a bug in
pre 4.3 gcc, but seems to be possible to quiet the warning if we start a new
frame (and don't access the var/arg after that call). some testing, seems to
work.
Modified Paths:
--------------
brlcad/trunk/src/conv/g-acad.c
Modified: brlcad/trunk/src/conv/g-acad.c
===================================================================
--- brlcad/trunk/src/conv/g-acad.c 2011-01-12 03:28:53 UTC (rev 42126)
+++ brlcad/trunk/src/conv/g-acad.c 2011-01-12 03:32:55 UTC (rev 42127)
@@ -46,227 +46,37 @@
#define V3ARGSIN(a) (a)[X]/25.4, (a)[Y]/25.4, (a)[Z]/25.4
-BU_EXTERN(union tree *do_region_end, (struct db_tree_state *tsp, const struct
db_full_path *pathp, union tree *curtree, genptr_t client_data));
-static char usage[] = "\
-Usage: %s [-v][-i][-xX lvl][-a abs_tess_tol][-r rel_tess_tol][-n
norm_tess_tol]\n\
-[-e error_file ][-D dist_calc_tol] -o output_file_name brlcad_db.g
object(s)\n";
+static int NMG_debug; /* saved arg of -X, for longjmp handling */
+static int verbose;
+static int ncpu = 1; /* Number of processors */
+static char *output_file = NULL; /* output filename */
+static char *error_file = NULL; /* error filename */
+static FILE *fp; /* Output file pointer */
+static FILE *fpe; /* Error file pointer */
+static struct db_i *dbip;
+static struct rt_tess_tol ttol;
+static struct bn_tol tol;
+static struct model *the_model;
-static int NMG_debug; /* saved arg of -X, for longjmp handling */
-static int verbose;
-static int ncpu = 1; /* Number of processors */
-static char *output_file = NULL; /* output filename */
-static char *error_file = NULL; /* error filename */
-static FILE *fp; /* Output file pointer */
-static FILE *fpe; /* Error file pointer */
-static struct db_i *dbip;
-static struct rt_tess_tol ttol;
-static struct bn_tol tol;
-static struct model *the_model;
+static struct db_tree_state tree_state; /* includes tol & model */
-static struct db_tree_state tree_state; /* includes tol & model */
+static int regions_tried = 0;
+static int regions_converted = 0;
+static int regions_written = 0;
+static int inches = 0;
+static long tot_polygons = 0;
-static int regions_tried = 0;
-static int regions_converted = 0;
-static int regions_written = 0;
-static int inches = 0;
-static long tot_polygons = 0;
-/*
- * M A I N
- */
-int
-main(int argc, char **argv)
+static void
+usage(const char *argv0)
{
- int c;
- double percent;
- int i;
+ bu_log("Usage: %s [-v][-i][-xX lvl][-a abs_tess_tol][-r rel_tess_tol][-n
norm_tess_tol]\n", argv0);
+ bu_log("[-e error_file ][-D dist_calc_tol] -o output_file_name brlcad_db.g
object(s)\n");
+ bu_exit(1, NULL);
+}
- bu_setlinebuf( stderr );
- tree_state = rt_initial_tree_state; /* struct copy */
- tree_state.ts_tol = &tol;
- tree_state.ts_ttol = &ttol;
- tree_state.ts_m = &the_model;
-
- ttol.magic = RT_TESS_TOL_MAGIC;
- /* Defaults, updated by command line options. */
- ttol.abs = 0.0;
- ttol.rel = 0.01;
- ttol.norm = 0.0;
-
- /* FIXME: These need to be improved */
- tol.magic = BN_TOL_MAGIC;
- tol.dist = 0.0005;
- tol.dist_sq = tol.dist * tol.dist;
- tol.perp = 1e-5;
- tol.para = 1 - tol.perp;
-
- rt_init_resource( &rt_uniresource, 0, NULL );
-
- the_model = nmg_mm();
- BU_LIST_INIT( &rt_g.rtg_vlfree ); /* for vlist macros */
-
- /* Get command line arguments. */
- while ((c = bu_getopt(argc, argv, "a:n:o:r:vx:D:P:X:e:i")) != EOF) {
- switch (c) {
- case 'a': /* Absolute tolerance. */
- ttol.abs = atof(bu_optarg);
- ttol.rel = 0.0;
- break;
- case 'n': /* Surface normal tolerance. */
- ttol.norm = atof(bu_optarg);
- ttol.rel = 0.0;
- break;
- case 'o': /* Output file name. */
- output_file = bu_optarg;
- break;
- case 'r': /* Relative tolerance. */
- ttol.rel = atof(bu_optarg);
- break;
- case 'v':
- verbose++;
- break;
- case 'P':
- ncpu = atoi( bu_optarg );
- rt_g.debug = 1; /* NOTE: enabling DEBUG_ALLRAYS to get core
dumps */
- break;
- case 'x':
- sscanf( bu_optarg, "%x", (unsigned int *)&rt_g.debug );
- break;
- case 'D':
- tol.dist = atof(bu_optarg);
- tol.dist_sq = tol.dist * tol.dist;
- rt_pr_tol( &tol );
- break;
- case 'X':
- sscanf( bu_optarg, "%x", (unsigned int *)&rt_g.NMG_debug );
- NMG_debug = rt_g.NMG_debug;
- break;
- case 'e': /* Error file name. */
- error_file = bu_optarg;
- break;
- case 'i':
- inches = 1;
- break;
- default:
- bu_exit(1, usage, argv[0]);
- break;
- }
- }
-
- if (bu_optind+1 >= argc) {
- bu_exit(1, usage, argv[0]);
- }
-
- if (!output_file) {
- fp = stdout;
-#if defined(_WIN32) && !defined(__CYGWIN__)
- setmode(fileno(fp), O_BINARY);
-#endif
- } else {
- /* Open output file */
- if ((fp=fopen( output_file, "wb+")) == NULL) {
- perror(argv[0]);
- bu_exit(1, "Cannot open output file (%s) for writing\n",
output_file);
- }
- }
-
- /* Open g-acad error log file */
- if (!error_file) {
- fpe = stderr;
-#if defined(_WIN32) && !defined(__CYGWIN__)
- setmode(fileno(fpe), O_BINARY);
-#endif
- } else if ((fpe=fopen( error_file, "wb")) == NULL) {
- perror(argv[0]);
- bu_exit(1, "Cannot open output file (%s) for writing\n", error_file);
- }
-
-
- /* Open BRL-CAD database */
- argc -= bu_optind;
- argv += bu_optind;
- if ((dbip = db_open(argv[0], "r")) == DBI_NULL) {
- perror(argv[0]);
- bu_exit(1, "Cannot open geometry file (%s) for reading\n", argv[0]);
- }
- if ( db_dirbuild( dbip ) ) {
- bu_exit(1, "db_dirbuild failed\n" );
- }
-
- BN_CK_TOL(tree_state.ts_tol);
- RT_CK_TESS_TOL(tree_state.ts_ttol);
-
- fprintf( fpe, "Model: %s\n", argv[0] );
- fprintf( fpe, "Objects:" );
- for ( i=1; i<argc; i++ )
- fprintf( fpe, " %s", argv[i] );
- fprintf( fpe, "\nTesselation tolerances:\n\tabs = %g mm\n\trel =
%g\n\tnorm = %g\n",
- tree_state.ts_ttol->abs, tree_state.ts_ttol->rel,
tree_state.ts_ttol->norm );
- fprintf( fpe, "Calculational tolerances:\n\tdist = %g mm perp = %g\n",
- tree_state.ts_tol->dist, tree_state.ts_tol->perp );
-
- bu_log( "Model: %s\n", argv[0] );
- bu_log( "Objects:" );
- for ( i=1; i<argc; i++ )
- bu_log( " %s", argv[i] );
- bu_log( "\nTesselation tolerances:\n\tabs = %g mm\n\trel = %g\n\tnorm =
%g\n",
- tree_state.ts_ttol->abs, tree_state.ts_ttol->rel,
tree_state.ts_ttol->norm );
- bu_log( "Calculational tolerances:\n\tdist = %g mm perp = %g\n",
- tree_state.ts_tol->dist, tree_state.ts_tol->perp );
-
-/* Write out ACAD facet header */
- if (inches)
- fprintf(fp, "BRL-CAD generated ACAD FACET FILE (Units in)\n");
- else
- fprintf(fp, "BRL-CAD generated ACAD FACET FILE (Units mm)\n");
-
-/* Generate space for number of facet entities, will write over later */
-
- fprintf(fp, " ");
-
- /* Walk indicated tree(s). Each region will be output separately */
- (void) db_walk_tree(dbip, argc-1, (const char **)(argv+1),
- 1, /* ncpu */
- &tree_state,
- 0, /* take all regions */
- do_region_end,
- nmg_booltree_leaf_tess,
- (genptr_t)NULL); /* in librt/nmg_bool.c */
-
- percent = 0;
- if (regions_tried>0) {
- percent = ((double)regions_converted * 100) / regions_tried;
- printf("Tried %d regions, %d converted to NMG's successfully. %g%%\n",
- regions_tried, regions_converted, percent);
- }
- percent = 0;
-
- if ( regions_tried > 0 ) {
- percent = ((double)regions_written * 100) / regions_tried;
- printf( " %d triangulated successfully. %g%%\n",
- regions_written, percent );
- }
-
- bu_log( "%ld triangles written\n", tot_polygons );
- fprintf( fpe, "%ld triangles written\n", tot_polygons );
-
- /* Write out number of facet entities to .facet file */
-
- rewind(fp);
- fseek(fp, 46, 0); /* Re-position pointer to 2nd line */
- fprintf(fp, "%d\n", regions_written); /* Write out number of regions */
- fclose(fp);
-
- /* Release dynamic storage */
- nmg_km(the_model);
- rt_vlist_cleanup();
- db_close(dbip);
-
- return 0;
-}
-
static void
nmg_to_acad(struct nmgregion *r, const struct db_full_path *pathp, int
region_id)
{
@@ -280,20 +90,20 @@
int tricount = 0; /* Triangle number */
int i;
- NMG_CK_REGION( r );
+ NMG_CK_REGION(r);
RT_CK_FULL_PATH(pathp);
- region_name = db_path_to_string( pathp );
+ region_name = db_path_to_string(pathp);
m = r->m_p;
- NMG_CK_MODEL( m );
+ NMG_CK_MODEL(m);
/* triangulate model */
- nmg_triangulate_model( m, &tol );
+ nmg_triangulate_model(m, &tol);
/* list all vertices in result */
- nmg_vertex_tabulate( &verts, &r->l.magic );
+ nmg_vertex_tabulate(&verts, &r->l.magic);
/* Get number of vertices */
@@ -303,64 +113,56 @@
/* Check vertices */
- for ( i=0; i<numverts; i++ )
- {
- v = (struct vertex *)BU_PTBL_GET( &verts, i );
- NMG_CK_VERTEX( v );
+ for (i=0; i<numverts; i++) {
+ v = (struct vertex *)BU_PTBL_GET(&verts, i);
+ NMG_CK_VERTEX(v);
}
/* Check triangles */
- for ( BU_LIST_FOR( s, shell, &r->s_hd ) )
- {
+ for (BU_LIST_FOR(s, shell, &r->s_hd)) {
struct faceuse *fu;
- NMG_CK_SHELL( s );
+ NMG_CK_SHELL(s);
- for ( BU_LIST_FOR( fu, faceuse, &s->fu_hd ) )
- {
+ for (BU_LIST_FOR(fu, faceuse, &s->fu_hd)) {
struct loopuse *lu;
- NMG_CK_FACEUSE( fu );
+ NMG_CK_FACEUSE(fu);
- if ( fu->orientation != OT_SAME )
+ if (fu->orientation != OT_SAME)
continue;
- for ( BU_LIST_FOR( lu, loopuse, &fu->lu_hd ) )
- {
+ for (BU_LIST_FOR(lu, loopuse, &fu->lu_hd)) {
struct edgeuse *eu;
int vert_count=0;
- NMG_CK_LOOPUSE( lu );
+ NMG_CK_LOOPUSE(lu);
- if ( BU_LIST_FIRST_MAGIC( &lu->down_hd ) != NMG_EDGEUSE_MAGIC )
+ if (BU_LIST_FIRST_MAGIC(&lu->down_hd) != NMG_EDGEUSE_MAGIC)
continue;
/* check vertex numbers for each triangle */
- for ( BU_LIST_FOR( eu, edgeuse, &lu->down_hd ) )
- {
- NMG_CK_EDGEUSE( eu );
+ for (BU_LIST_FOR(eu, edgeuse, &lu->down_hd)) {
+ NMG_CK_EDGEUSE(eu);
v = eu->vu_p->v_p;
- NMG_CK_VERTEX( v );
+ NMG_CK_VERTEX(v);
vert_count++;
- i = bu_ptbl_locate( &verts, (long *)v );
- if ( i < 0 )
- {
- bu_ptbl_free( &verts);
- bu_free( region_name, "region name" );
+ i = bu_ptbl_locate(&verts, (long *)v);
+ if (i < 0) {
+ bu_ptbl_free(&verts);
+ bu_free(region_name, "region name");
bu_log("Vertex from eu x%x is not in nmgregion x%x\n",
eu, r);
bu_exit(1, "ERROR: Triangle vertex was not located\n");
}
}
- if ( vert_count > 3 )
- {
- bu_ptbl_free( &verts);
- bu_free( region_name, "region name" );
- bu_log( "lu x%x has too many (%d) vertices!\n", lu,
vert_count );
+ if (vert_count > 3) {
+ bu_ptbl_free(&verts);
+ bu_free(region_name, "region name");
+ bu_log("lu x%x has too many (%d) vertices!\n", lu,
vert_count);
bu_exit(1, "ERROR: LU is not a triangle\n");
- }
- else if ( vert_count < 3 )
+ } else if (vert_count < 3)
continue;
numtri++;
}
@@ -370,120 +172,170 @@
/* END CHECK SECTION */
/* Write pertinent info for this region */
- fprintf( fp, "%s\n", (region_name+1));
+ fprintf(fp, "%s\n", (region_name+1));
/* No mirror plane */
- fprintf( fp, "%d\n", 0);
+ fprintf(fp, "%d\n", 0);
/* Number of vertices */
- fprintf( fp, "%d\n", numverts);
+ fprintf(fp, "%d\n", numverts);
/* Write numverts, then vertices */
- for ( i=0; i<numverts; i++ )
- {
- v = (struct vertex *)BU_PTBL_GET( &verts, i );
- NMG_CK_VERTEX( v );
+ for (i=0; i<numverts; i++) {
+ v = (struct vertex *)BU_PTBL_GET(&verts, i);
+ NMG_CK_VERTEX(v);
if (inches)
- fprintf( fp, "%f %f %f\n", V3ARGSIN( v->vg_p->coord ));
+ fprintf(fp, "%f %f %f\n", V3ARGSIN(v->vg_p->coord));
else
- fprintf( fp, "%f %f %f\n", V3ARGS( v->vg_p->coord ));
+ fprintf(fp, "%f %f %f\n", V3ARGS(v->vg_p->coord));
}
/* Number of sub-parts (always 1 with BRL-CAD) */
- fprintf( fp, "%d\n", 1);
+ fprintf(fp, "%d\n", 1);
/* Write out name again */
- fprintf( fp, "%s\n", (region_name+1));
+ fprintf(fp, "%s\n", (region_name+1));
/* Number of triangles, number of vert/tri (3) */
- fprintf( fp, "%d %d\n", numtri, 3);
+ fprintf(fp, "%d %d\n", numtri, 3);
/* output triangles */
- for ( BU_LIST_FOR( s, shell, &r->s_hd ) )
- {
+ for (BU_LIST_FOR(s, shell, &r->s_hd)) {
struct faceuse *fu;
- NMG_CK_SHELL( s );
+ NMG_CK_SHELL(s);
- for ( BU_LIST_FOR( fu, faceuse, &s->fu_hd ) )
- {
+ for (BU_LIST_FOR(fu, faceuse, &s->fu_hd)) {
struct loopuse *lu;
- NMG_CK_FACEUSE( fu );
+ NMG_CK_FACEUSE(fu);
- if ( fu->orientation != OT_SAME )
+ if (fu->orientation != OT_SAME)
continue;
- for ( BU_LIST_FOR( lu, loopuse, &fu->lu_hd ) )
- {
+ for (BU_LIST_FOR(lu, loopuse, &fu->lu_hd)) {
struct edgeuse *eu;
int vert_count=0;
- NMG_CK_LOOPUSE( lu );
+ NMG_CK_LOOPUSE(lu);
- if ( BU_LIST_FIRST_MAGIC( &lu->down_hd ) != NMG_EDGEUSE_MAGIC )
+ if (BU_LIST_FIRST_MAGIC(&lu->down_hd) != NMG_EDGEUSE_MAGIC)
continue;
/* list vertex numbers for each triangle */
- for ( BU_LIST_FOR( eu, edgeuse, &lu->down_hd ) )
- {
- NMG_CK_EDGEUSE( eu );
+ for (BU_LIST_FOR(eu, edgeuse, &lu->down_hd)) {
+ NMG_CK_EDGEUSE(eu);
v = eu->vu_p->v_p;
- NMG_CK_VERTEX( v );
+ NMG_CK_VERTEX(v);
vert_count++;
- i = bu_ptbl_locate( &verts, (long *)v );
- if ( i < 0 )
- {
- bu_ptbl_free( &verts);
- bu_log( "Vertex from eu x%x is not in nmgregion x%x\n",
eu, r );
- bu_free( region_name, "region name" );
+ i = bu_ptbl_locate(&verts, (long *)v);
+ if (i < 0) {
+ bu_ptbl_free(&verts);
+ bu_log("Vertex from eu x%x is not in nmgregion x%x\n",
eu, r);
+ bu_free(region_name, "region name");
bu_exit(1, "ERROR: Can't find vertex in list!\n");
}
- fprintf( fp, " %d", i+1 );
+ fprintf(fp, " %d", i+1);
}
/* Output other info. for triangle ICOAT, component#, facet# */
/* Map Icoat from material table later */
- /* fprintf( fp, "%s icomp=%d material=%d:\n", (region_name+1),
region_id );*/
+ /* fprintf(fp, "%s icomp=%d material=%d:\n", (region_name+1),
region_id);*/
- fprintf( fp, " %d %d %d\n", 0, region_id, ++tricount);
+ fprintf(fp, " %d %d %d\n", 0, region_id, ++tricount);
- if ( vert_count > 3 )
- {
- bu_ptbl_free( &verts);
- bu_free( region_name, "region name" );
- bu_log( "lu x%x has %d vertices!\n", lu, vert_count );
+ if (vert_count > 3) {
+ bu_ptbl_free(&verts);
+ bu_free(region_name, "region name");
+ bu_log("lu x%x has %d vertices!\n", lu, vert_count);
bu_exit(1, "ERROR: LU is not a triangle\n");
- }
- else if ( vert_count < 3 )
+ } else if (vert_count < 3)
continue;
tot_polygons++;
}
}
}
-/* regions_converted++;
- printf("Processed region %s\n", region_name);
- printf("Regions attempted = %d Regions done = %d\n", regions_tried,
regions_converted);
- fflush(stdout);
+/* regions_converted++;
+ printf("Processed region %s\n", region_name);
+ printf("Regions attempted = %d Regions done = %d\n", regions_tried,
regions_converted);
+ fflush(stdout);
*/
- bu_ptbl_free( &verts);
- bu_free( region_name, "region name" );
+ bu_ptbl_free(&verts);
+ bu_free(region_name, "region name");
}
+
+static union tree *
+process_region(const struct db_full_path *pathp, union tree *curtree, struct
db_tree_state *tsp)
+{
+ /* Begin bomb protection */
+ if (!BU_SETJUMP) {
+ /* try */
+
+ union tree *ret_tree;
+
+ printf("Attempting to process region %s\n", db_path_to_string(pathp));
+ fflush(stdout);
+ ret_tree = nmg_booltree_evaluate(curtree, tsp->ts_tol, &rt_uniresource);
+ db_free_tree(curtree, &rt_uniresource); /* Does an nmg_kr()
*/
+
+ return ret_tree;
+ } else {
+ /* catch */
+
+ /* Error, bail out */
+ char *sofar;
+ BU_UNSETJUMP; /* Relinquish the protection */
+
+ sofar = db_path_to_string(pathp);
+ bu_log("FAILED in Boolean evaluation: %s\n", sofar);
+ fprintf(fpe, "Failed Bool. Eval.: %s\n", sofar);
+ fflush(fpe);
+ bu_free((char *)sofar, "sofar");
+
+ /* Sometimes the NMG library adds debugging bits when
+ * it detects an internal error, before bombing out.
+ */
+ rt_g.NMG_debug = NMG_debug; /* restore mode */
+
+ /* Release any intersector 2d tables */
+ nmg_isect2d_final_cleanup();
+
+ /* Release the tree memory & input regions */
+
+ /* FIXME: memory leak? */
+ /* db_free_tree(curtree);*/ /* Does an nmg_kr() */
+
+ /* Get rid of (m)any other intermediate structures */
+ if ((*tsp->ts_m)->magic == NMG_MODEL_MAGIC) {
+ nmg_km(*tsp->ts_m);
+ } else {
+ bu_log("WARNING: tsp->ts_m pointer corrupted, ignoring it.\n");
+ }
+
+ /* Now, make a new, clean model structure for next pass. */
+ *tsp->ts_m = nmg_mm();
+
+ return TREE_NULL;
+ } BU_UNSETJUMP; /* Relinquish the protection */
+}
+
+
/*
- * D O _ R E G I O N _ E N D
+ * D O _ R E G I O N _ E N D
*
- * Called from db_walk_tree().
+ * Called from db_walk_tree().
*
- * This routine must be prepared to run in parallel.
+ * This routine must be prepared to run in parallel.
*/
union tree *
do_region_end(struct db_tree_state *tsp, const struct db_full_path *pathp,
union tree *curtree, genptr_t UNUSED(client_data))
{
- union tree *ret_tree;
- struct bu_list vhead;
- struct nmgregion *r;
+ union tree *ret_tree = NULL;
+ struct bu_list vhead;
+ /* static due to longjmp */
+ static struct nmgregion *r = NULL;
RT_CK_FULL_PATH(pathp);
RT_CK_TREE(curtree);
@@ -494,7 +346,7 @@
BU_LIST_INIT(&vhead);
if (RT_G_DEBUG&DEBUG_TREEWALK || verbose) {
- char *sofar = db_path_to_string(pathp);
+ char *sofar = db_path_to_string(pathp);
bu_log("\ndo_region_end(%d %d%%) %s\n",
regions_tried,
regions_tried>0 ? (regions_converted * 100) / regions_tried : 0,
@@ -507,82 +359,37 @@
regions_tried++;
- /* Begin bomb protection */
- if ( ncpu == 1 ) {
- if ( BU_SETJUMP ) {
- /* Error, bail out */
- char *sofar;
- BU_UNSETJUMP; /* Relinquish the protection */
+ /* do the deed */
+ ret_tree = process_region(pathp, curtree, tsp);
- sofar = db_path_to_string(pathp);
- bu_log( "FAILED in Boolean evaluation: %s\n", sofar );
- fprintf(fpe, "Failed Bool. Eval.: %s\n", sofar);
- fflush(fpe);
- bu_free( (char *)sofar, "sofar" );
-
- /* Sometimes the NMG library adds debugging bits when
- * it detects an internal error, before bombing out.
- */
- rt_g.NMG_debug = NMG_debug; /* restore mode */
-
- /* Release any intersector 2d tables */
- nmg_isect2d_final_cleanup();
-
- /* Release the tree memory & input regions */
-
- /* FIXME: memory leak? */
- /* db_free_tree(curtree);*/ /* Does an nmg_kr() */
-
- /* Get rid of (m)any other intermediate structures */
- if ( (*tsp->ts_m)->magic == NMG_MODEL_MAGIC ) {
- nmg_km(*tsp->ts_m);
- } else {
- bu_log("WARNING: tsp->ts_m pointer corrupted, ignoring it.\n");
- }
-
- /* Now, make a new, clean model structure for next pass. */
- *tsp->ts_m = nmg_mm();
- goto out;
- }
- }
- printf("Attempting to process region %s\n", db_path_to_string( pathp ));
- fflush(stdout);
- ret_tree = nmg_booltree_evaluate( curtree, tsp->ts_tol, &rt_uniresource );
/* librt/nmg_bool.c */
-
- if ( ret_tree )
+ if (ret_tree)
r = ret_tree->tr_d.td_r;
- else
- {
+ else {
if (verbose) {
- printf( "\tNothing left of this region after Boolean evaluation\n"
);
- fprintf( fpe, "WARNING: Nothing left after Boolean evaluation:
%s\n",
- db_path_to_string( pathp ));
+ printf("\tNothing left of this region after Boolean evaluation\n");
+ fprintf(fpe, "WARNING: Nothing left after Boolean evaluation: %s\n",
+ db_path_to_string(pathp));
fflush(fpe);
}
regions_written++; /* don't count as a failure */
r = (struct nmgregion *)NULL;
}
- BU_UNSETJUMP; /* Relinquish the protection */
regions_converted++;
- if (r != (struct nmgregion *)NULL)
- {
+ if (r != (struct nmgregion *)NULL) {
struct shell *s;
int empty_region=0;
int empty_model=0;
/* Kill cracks */
- s = BU_LIST_FIRST( shell, &r->s_hd );
- while ( BU_LIST_NOT_HEAD( &s->l, &r->s_hd ) )
- {
+ s = BU_LIST_FIRST(shell, &r->s_hd);
+ while (BU_LIST_NOT_HEAD(&s->l, &r->s_hd)) {
struct shell *next_s;
- next_s = BU_LIST_PNEXT( shell, &s->l );
- if ( nmg_kill_cracks( s ) )
- {
- if ( nmg_ks( s ) )
- {
+ next_s = BU_LIST_PNEXT(shell, &s->l);
+ if (nmg_kill_cracks(s)) {
+ if (nmg_ks(s)) {
empty_region = 1;
break;
}
@@ -591,24 +398,30 @@
}
/* kill zero length edgeuses */
- if ( !empty_region )
- {
- empty_model = nmg_kill_zero_length_edgeuses( *tsp->ts_m );
+ if (!empty_region) {
+ empty_model = nmg_kill_zero_length_edgeuses(*tsp->ts_m);
}
- if ( !empty_region && !empty_model )
- {
- if ( BU_SETJUMP )
- {
+ if (!empty_region && !empty_model) {
+ if (!BU_SETJUMP) {
+ /* try */
+
+ /* Write the region to the TANKILL file */
+ nmg_to_acad(r, pathp, tsp->ts_regionid);
+ regions_written++;
+
+ } else {
+ /* catch */
+
char *sofar;
BU_UNSETJUMP;
sofar = db_path_to_string(pathp);
- bu_log( "FAILED in triangulator: %s\n", sofar );
+ bu_log("FAILED in triangulator: %s\n", sofar);
fprintf(fpe, "Failed in triangulator: %s\n", sofar);
fflush(fpe);
- bu_free( (char *)sofar, "sofar" );
+ bu_free((char *)sofar, "sofar");
/* Sometimes the NMG library adds debugging bits when
* it detects an internal error, before bombing out.
@@ -619,41 +432,33 @@
nmg_isect2d_final_cleanup();
/* Get rid of (m)any other intermediate structures */
- if ( (*tsp->ts_m)->magic == NMG_MODEL_MAGIC )
- {
+ if ((*tsp->ts_m)->magic == NMG_MODEL_MAGIC) {
nmg_km(*tsp->ts_m);
- }
- else
- {
+ } else {
bu_log("WARNING: tsp->ts_m pointer corrupted, ignoring
it.\n");
}
/* Now, make a new, clean model structure for next pass. */
*tsp->ts_m = nmg_mm();
- goto out;
- }
- /* Write the region to the TANKILL file */
- nmg_to_acad( r, pathp, tsp->ts_regionid );
- regions_written++;
+ /* FIXME: leaking memory with curtree */
- BU_UNSETJUMP;
+ return TREE_NULL;
+ } BU_UNSETJUMP;
}
- if ( !empty_model )
- nmg_kr( r );
+ if (!empty_model)
+ nmg_kr(r);
}
- out:
/*
- * Dispose of original tree, so that all associated dynamic
- * memory is released now, not at the end of all regions.
- * A return of TREE_NULL from this routine signals an error,
- * and there is no point to adding _another_ message to our output,
- * so we need to cons up an OP_NOP node to return.
+ * Dispose of original tree, so that all associated dynamic
+ * memory is released now, not at the end of all regions.
+ * A return of TREE_NULL from this routine signals an error,
+ * and there is no point to adding _another_ message to our output,
+ * so we need to cons up an OP_NOP node to return.
*/
-
if (regions_tried>0) {
float npercent, tpercent;
@@ -663,15 +468,211 @@
regions_tried, regions_converted, regions_written, npercent,
tpercent);
}
- db_free_tree(curtree, &rt_uniresource); /* Does an nmg_kr() */
-
BU_GETUNION(curtree, tree);
curtree->magic = RT_TREE_MAGIC;
curtree->tr_op = OP_NOP;
+
return curtree;
}
+
/*
+ * M A I N
+ */
+int
+main(int argc, char **argv)
+{
+ int c;
+ double percent;
+ int i;
+
+ bu_setlinebuf(stderr);
+
+ tree_state = rt_initial_tree_state; /* struct copy */
+ tree_state.ts_tol = &tol;
+ tree_state.ts_ttol = &ttol;
+ tree_state.ts_m = &the_model;
+
+ ttol.magic = RT_TESS_TOL_MAGIC;
+ /* Defaults, updated by command line options. */
+ ttol.abs = 0.0;
+ ttol.rel = 0.01;
+ ttol.norm = 0.0;
+
+ /* FIXME: These need to be improved */
+ tol.magic = BN_TOL_MAGIC;
+ tol.dist = 0.0005;
+ tol.dist_sq = tol.dist * tol.dist;
+ tol.perp = 1e-5;
+ tol.para = 1 - tol.perp;
+
+ rt_init_resource(&rt_uniresource, 0, NULL);
+
+ the_model = nmg_mm();
+ BU_LIST_INIT(&rt_g.rtg_vlfree); /* for vlist macros */
+
+ /* Get command line arguments. */
+ while ((c = bu_getopt(argc, argv, "a:n:o:r:vx:D:P:X:e:i")) != EOF) {
+ switch (c) {
+ case 'a': /* Absolute tolerance. */
+ ttol.abs = atof(bu_optarg);
+ ttol.rel = 0.0;
+ break;
+ case 'n': /* Surface normal tolerance. */
+ ttol.norm = atof(bu_optarg);
+ ttol.rel = 0.0;
+ break;
+ case 'o': /* Output file name. */
+ output_file = bu_optarg;
+ break;
+ case 'r': /* Relative tolerance. */
+ ttol.rel = atof(bu_optarg);
+ break;
+ case 'v':
+ verbose++;
+ break;
+ case 'P':
+ ncpu = atoi(bu_optarg);
+ rt_g.debug = 1; /* NOTE: enabling DEBUG_ALLRAYS to get core
dumps */
+ break;
+ case 'x':
+ sscanf(bu_optarg, "%x", (unsigned int *)&rt_g.debug);
+ break;
+ case 'D':
+ tol.dist = atof(bu_optarg);
+ tol.dist_sq = tol.dist * tol.dist;
+ rt_pr_tol(&tol);
+ break;
+ case 'X':
+ sscanf(bu_optarg, "%x", (unsigned int *)&rt_g.NMG_debug);
+ NMG_debug = rt_g.NMG_debug;
+ break;
+ case 'e': /* Error file name. */
+ error_file = bu_optarg;
+ break;
+ case 'i':
+ inches = 1;
+ break;
+ default:
+ usage(argv[0]);
+ break;
+ }
+ }
+
+ if (bu_optind+1 >= argc) {
+ usage(argv[0]);
+ }
+
+ if (!output_file) {
+ fp = stdout;
+#if defined(_WIN32) && !defined(__CYGWIN__)
+ setmode(fileno(fp), O_BINARY);
+#endif
+ } else {
+ /* Open output file */
+ if ((fp=fopen(output_file, "wb+")) == NULL) {
+ perror(argv[0]);
+ bu_exit(1, "Cannot open output file (%s) for writing\n",
output_file);
+ }
+ }
+
+ /* Open g-acad error log file */
+ if (!error_file) {
+ fpe = stderr;
+#if defined(_WIN32) && !defined(__CYGWIN__)
+ setmode(fileno(fpe), O_BINARY);
+#endif
+ } else if ((fpe=fopen(error_file, "wb")) == NULL) {
+ perror(argv[0]);
+ bu_exit(1, "Cannot open output file (%s) for writing\n", error_file);
+ }
+
+
+ /* Open BRL-CAD database */
+ argc -= bu_optind;
+ argv += bu_optind;
+ if ((dbip = db_open(argv[0], "r")) == DBI_NULL) {
+ perror(argv[0]);
+ bu_exit(1, "Cannot open geometry file (%s) for reading\n", argv[0]);
+ }
+ if (db_dirbuild(dbip)) {
+ bu_exit(1, "db_dirbuild failed\n");
+ }
+
+ BN_CK_TOL(tree_state.ts_tol);
+ RT_CK_TESS_TOL(tree_state.ts_ttol);
+
+ fprintf(fpe, "Model: %s\n", argv[0]);
+ fprintf(fpe, "Objects:");
+ for (i=1; i<argc; i++)
+ fprintf(fpe, " %s", argv[i]);
+ fprintf(fpe, "\nTesselation tolerances:\n\tabs = %g mm\n\trel = %g\n\tnorm
= %g\n",
+ tree_state.ts_ttol->abs, tree_state.ts_ttol->rel,
tree_state.ts_ttol->norm);
+ fprintf(fpe, "Calculational tolerances:\n\tdist = %g mm perp = %g\n",
+ tree_state.ts_tol->dist, tree_state.ts_tol->perp);
+
+ bu_log("Model: %s\n", argv[0]);
+ bu_log("Objects:");
+ for (i=1; i<argc; i++)
+ bu_log(" %s", argv[i]);
+ bu_log("\nTesselation tolerances:\n\tabs = %g mm\n\trel = %g\n\tnorm =
%g\n",
+ tree_state.ts_ttol->abs, tree_state.ts_ttol->rel,
tree_state.ts_ttol->norm);
+ bu_log("Calculational tolerances:\n\tdist = %g mm perp = %g\n",
+ tree_state.ts_tol->dist, tree_state.ts_tol->perp);
+
+/* Write out ACAD facet header */
+ if (inches)
+ fprintf(fp, "BRL-CAD generated ACAD FACET FILE (Units in)\n");
+ else
+ fprintf(fp, "BRL-CAD generated ACAD FACET FILE (Units mm)\n");
+
+/* Generate space for number of facet entities, will write over later */
+
+ fprintf(fp, " ");
+
+ /* Walk indicated tree(s). Each region will be output separately */
+ (void) db_walk_tree(dbip, argc-1, (const char **)(argv+1),
+ 1, /* ncpu */
+ &tree_state,
+ 0, /* take all regions */
+ do_region_end,
+ nmg_booltree_leaf_tess,
+ (genptr_t)NULL); /* in librt/nmg_bool.c */
+
+ percent = 0;
+ if (regions_tried>0) {
+ percent = ((double)regions_converted * 100) / regions_tried;
+ printf("Tried %d regions, %d converted to NMG's successfully. %g%%\n",
+ regions_tried, regions_converted, percent);
+ }
+ percent = 0;
+
+ if (regions_tried > 0) {
+ percent = ((double)regions_written * 100) / regions_tried;
+ printf(" %d triangulated successfully. %g%%\n",
+ regions_written, percent);
+ }
+
+ bu_log("%ld triangles written\n", tot_polygons);
+ fprintf(fpe, "%ld triangles written\n", tot_polygons);
+
+ /* Write out number of facet entities to .facet file */
+
+ rewind(fp);
+ fseek(fp, 46, 0); /* Re-position pointer to 2nd line */
+ fprintf(fp, "%d\n", regions_written); /* Write out number of regions */
+ fclose(fp);
+
+ /* Release dynamic storage */
+ nmg_km(the_model);
+ rt_vlist_cleanup();
+ db_close(dbip);
+
+ return 0;
+}
+
+
+/*
* Local Variables:
* mode: C
* tab-width: 8
This was sent by the SourceForge.net collaborative development platform, the
world's largest Open Source development site.
------------------------------------------------------------------------------
Protect Your Site and Customers from Malware Attacks
Learn about various malware tactics and how to avoid them. Understand
malware threats, the impact they can have on your business, and how you
can protect your company and customers by using code signing.
http://p.sf.net/sfu/oracle-sfdevnl
_______________________________________________
BRL-CAD Source Commits mailing list
[email protected]
https://lists.sourceforge.net/lists/listinfo/brlcad-commits
