/*
* Call parse_args to handle command line arguments first, then
* process input.
*/
int
main(int ac, char *av[])
{
struct bn_tol tol;
FILE *fdplot;
struct rt_wdb *fdmodel;
parse_args(ac, av);
if (ac==1) usage((char *)NULL,0);
if (!manifold[0] && !manifold[1] && !manifold[2] && !manifold[3])
usage("No manifolds selected\n",1);
m = nmg_mm();
r = nmg_mrsv(m);
s = BU_LIST_FIRST(shell, &r->s_hd);
tol.magic = BN_TOL_MAGIC;
tol.dist = 0.01;
tol.dist_sq = tol.dist * tol.dist;
tol.perp = 0.001;
tol.para = 0.999;
if (manifold[3]) make_3manifold_bits(&tol);
if (manifold[2]) make_2manifold_bits(&tol);
if (manifold[1]) make_1manifold_bits();
if (manifold[0]) make_0manifold_bits();
NMG_CK_MODEL(m);
/* write a plot file */
if ((fdplot = fopen(plotfilename, "w")) == (FILE *)NULL)
perror(plotfilename);
else {
nmg_pl_m(fdplot, m);
fclose(fdplot);
}
/* write the database */
if ((fdmodel = wdb_fopen(mfilename)) == NULL)
perror(mfilename);
else {
mk_id(fdmodel, "hairy NMG");
mk_nmg(fdmodel, "s.NMG", m); /* releases m, boo */
/* build a database region mentioning the solid */
mk_comb1(fdmodel, "r.NMG", "s.NMG", 1);
wdb_close(fdmodel);
}
return 0;
}
int off2nmg(FILE *fpin, struct rt_wdb *fpout)
{
char title[64], geom_fname[64];
char rname[67], sname[67];
char buf[200], buf2[200];
FILE *fgeom;
struct model *m;
title[0] = geom_fname[0] = '\0';
bu_fgets(buf, sizeof(buf), fpin);
while (!feof(fpin)) {
/* Retrieve the important data */
if (sscanf(buf, "name %[^\n]s", buf2) > 0)
bu_strlcpy(title, buf2, sizeof(title));
if (sscanf(buf, "geometry %200[^\n]s", buf2) > 0) {
char dtype[40], format[40];
if (sscanf(buf2, "%40s %40s %64s", dtype, format, geom_fname) != 3)
bu_exit(1, "Incomplete geometry field in input file.");
if (!BU_STR_EQUAL(dtype, "indexed_poly"))
bu_exit(1, "Unknown geometry data type. Must be \"indexed_poly\".");
}
bu_fgets(buf, sizeof(buf), fpin);
}
if (strlen(title) < (unsigned)1)
fprintf(stderr, "Warning: no title\n");
if (strlen(geom_fname) < (unsigned)1)
bu_exit(1, "ERROR: no geometry filename given");
if ((fgeom = fopen(geom_fname, "rb")) == NULL) {
bu_exit(1, "off2nmg: cannot open %s (geometry description) for reading\n",
geom_fname);
}
m = nmg_mm();
read_faces(m, fgeom);
fclose(fgeom);
snprintf(sname, 67, "s.%s", title);
snprintf(rname, 67, "r.%s", title);
mk_id(fpout, title);
mk_nmg(fpout, sname, m);
mk_comb1(fpout, rname, sname, 1);
nmg_km(m);
return 0;
}
/*
* Write the nmg to a BRL-CAD style data base.
*/
void
create_brlcad_db(struct rt_wdb *fpout, struct model *m, char *reg_name, char *grp_name)
{
char *rname, *sname;
int size = sizeof(reg_name) + 3;
mk_id(fpout, "Ascii NMG");
rname = (char *)bu_malloc(size, "rname"); /* Region name. */
sname = (char *)bu_malloc(size, "sname"); /* Solid name. */
snprintf(sname, size, "s.%s", reg_name);
mk_nmg(fpout, sname, m); /* Make nmg object. */
snprintf(rname, size, "r.%s", reg_name);
mk_comb1(fpout, rname, sname, 1); /* Put object in a region. */
if (grp_name) {
mk_comb1(fpout, grp_name, rname, 1); /* Region in group. */
}
}
int
main(int argc, char **argv) /* really has no arguments */
{
struct nmgregion *r;
char * id_name = "BRL-CAD t-NURBS NMG Example";
char * tea_name = "UtahTeapot";
char * uplot_name = "teapot.pl";
struct bu_list vhead;
FILE *fp;
int i;
tol.magic = BN_TOL_MAGIC;
tol.dist = 0.005;
tol.dist_sq = tol.dist * tol.dist;
tol.perp = 1e-6;
tol.para = 1 - tol.perp;
BU_LIST_INIT( &rt_g.rtg_vlfree );
outfp = wdb_fopen( "tea_nmg.g" );
rt_g.debug |= DEBUG_ALLRAYS; /* Cause core dumps on bu_bomb(), but no extra messages */
while ((i=bu_getopt(argc, argv, "d")) != EOF) {
switch (i) {
case 'd' :
rt_g.debug |= DEBUG_MEM | DEBUG_MEM_FULL;
break;
default :
(void)fprintf(stderr,
"Usage: %s [-d] > database.g\n", *argv);
return(-1);
}
}
mk_id( outfp, id_name);
m = nmg_mm();
NMG_CK_MODEL( m );
r = nmg_mrsv( m );
NMG_CK_REGION( r );
s = BU_LIST_FIRST( shell, &r->s_hd );
NMG_CK_SHELL( s );
/* Step through each patch and create a NMG TNURB face
* representing the patch then dump them out.
*/
for ( i = 0; i < PATCH_COUNT; i++)
{
dump_patch( patches[i] );
}
/* Connect up the coincident vertexuses and edges */
(void)nmg_model_fuse( m, &tol );
/* write NMG to output file */
(void)mk_nmg( outfp, tea_name, m );
wdb_close(outfp);
/* Make a vlist drawing of the model */
BU_LIST_INIT( &vhead );
nmg_m_to_vlist( &vhead, m, 0 );
/* Make a UNIX plot file from this vlist */
if ( (fp=fopen( uplot_name, "w" )) == NULL )
{
bu_log( "Cannot open plot file: %s\n", uplot_name );
perror( "teapot_nmg" );
}
else
rt_vlist_to_uplot( fp, &vhead );
return(0);
}
int
main(int argc, char **argv)
{
int c;
int i;
struct pshell *psh;
struct pbar *pbp;
struct wmember head;
struct wmember all_head;
char *nastran_file = "Converted from NASTRAN file (stdin)";
bu_setprogname(argv[0]);
fpin = stdin;
units = INCHES;
/* 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-6;
tol.para = 1 - tol.perp;
while ((c=bu_getopt(argc, argv, "x:X:t:ni:o:mh?")) != -1) {
switch (c) {
case 'x':
sscanf(bu_optarg, "%x", (unsigned int *)&RTG.debug);
bu_printb("librt RT_G_DEBUG", RT_G_DEBUG, DEBUG_FORMAT);
bu_log("\n");
break;
case 'X':
sscanf(bu_optarg, "%x", (unsigned int *)&RTG.NMG_debug);
bu_printb("librt RTG.NMG_debug", RTG.NMG_debug, NMG_DEBUG_FORMAT);
bu_log("\n");
break;
case 't': /* calculational tolerance */
tol.dist = atof(bu_optarg);
tol.dist_sq = tol.dist * tol.dist;
break;
case 'n':
polysolids = 0;
break;
case 'm':
units = MM;
break;
case 'i':
fpin = fopen(bu_optarg, "rb");
if (fpin == (FILE *)NULL) {
bu_log("Cannot open NASTRAN file (%s) for reading!\n", bu_optarg);
bu_exit(1, Usage, argv[0]);
}
nastran_file = bu_optarg;
break;
case 'o':
output_file = bu_optarg;
break;
default:
bu_exit(1, Usage, argv[0]);
}
}
fpout = wdb_fopen(output_file);
if (fpout == NULL) {
bu_log("Cannot open BRL-CAD file (%s) for writing!\n", output_file);
bu_exit(1, Usage, argv[0]);
}
if (!fpin || !fpout) {
bu_exit(1, Usage, argv[0]);
}
line = (char *)bu_malloc(MAX_LINE_SIZE, "line");
next_line = (char *)bu_malloc(MAX_LINE_SIZE, "next_line");
prev_line = (char *)bu_malloc(MAX_LINE_SIZE, "prev_line");
curr_rec = (char **)bu_calloc(NO_OF_FIELDS, sizeof(char *), "curr_rec");
for (i=0; i<NO_OF_FIELDS; i++)
curr_rec[i] = (char *)bu_malloc(sizeof(char)*FIELD_LENGTH, "curr_rec[i]");
prev_rec = (char **)bu_calloc(NO_OF_FIELDS, sizeof(char *), "prev_rec");
for (i=0; i<NO_OF_FIELDS; i++)
prev_rec[i] = (char *)bu_malloc(sizeof(char)*FIELD_LENGTH, "prev_rec[i]");
/* first pass, find start of NASTRAN "bulk data" */
start_off = (-1);
bulk_data_start_line = 0;
while (bu_fgets(line, MAX_LINE_SIZE, fpin)) {
bulk_data_start_line++;
if (bu_strncmp(line, "BEGIN BULK", 10))
continue;
start_off = bu_ftell(fpin);
break;
}
if (start_off < 0) {
bu_log("Cannot find start of bulk data in NASTRAN file!\n");
bu_exit(1, Usage, argv[0]);
}
/* convert BULK data deck into something reasonable */
fptmp = bu_temp_file(NULL, 0);
if (fptmp == NULL) {
perror(argv[0]);
bu_exit(1, "Cannot open temporary file\n");
}
convert_input();
/* initialize some lists */
BU_LIST_INIT(&coord_head.l);
BU_LIST_INIT(&pbar_head.l);
BU_LIST_INIT(&pshell_head.l);
BU_LIST_INIT(&all_head.l);
nmg_model = (struct model *)NULL;
/* count grid points */
bu_fseek(fptmp, 0, SEEK_SET);
while (bu_fgets(line, MAX_LINE_SIZE, fptmp)) {
if (!bu_strncmp(line, "GRID", 4))
grid_count++;
}
if (!grid_count) {
bu_exit(1, "No geometry in this NASTRAN file!\n");
}
/* get default values and properties */
bu_fseek(fptmp, 0, SEEK_SET);
while (get_next_record(fptmp, 1, 0)) {
if (!bu_strncmp(curr_rec[0], "BAROR", 5)) {
/* get BAR defaults */
bar_def_pid = atoi(curr_rec[2]);
} else if (!bu_strncmp(curr_rec[0], "PBAR", 4)) {
struct pbar *pb;
BU_ALLOC(pb, struct pbar);
pb->pid = atoi(curr_rec[1]);
pb->mid = atoi(curr_rec[2]);
pb->area = atof(curr_rec[3]);
BU_LIST_INIT(&pb->head.l);
BU_LIST_INSERT(&pbar_head.l, &pb->l);
} else if (!bu_strncmp(curr_rec[0], "PSHELL", 6)) {
BU_ALLOC(psh, struct pshell);
psh->s = (struct shell *)NULL;
psh->pid = atoi(curr_rec[1]);
psh->mid = atoi(curr_rec[2]);
psh->thick = atof(curr_rec[3]);
BU_LIST_INSERT(&pshell_head.l, &psh->l);
pshell_count++;
}
}
/* allocate storage for grid points */
g_pts = (struct grid_point *)bu_calloc(grid_count, sizeof(struct grid_point), "grid points");
/* get all grid points */
bu_fseek(fptmp, 0, SEEK_SET);
while (get_next_record(fptmp, 1, 0)) {
int gid;
int cid;
double tmp[3];
if (bu_strncmp(curr_rec[0], "GRID", 4))
continue;
gid = atoi(curr_rec[1]);
cid = atoi(curr_rec[2]);
for (i=0; i<3; i++) {
tmp[i] = atof(curr_rec[i+3]);
}
g_pts[grid_used].gid = gid;
g_pts[grid_used].cid = cid;
g_pts[grid_used].v = (struct vertex **)bu_calloc(pshell_count + 1, sizeof(struct vertex *), "g_pts vertex array");
VMOVE(g_pts[grid_used].pt, tmp);
grid_used++;
}
/* find coordinate systems */
bu_fseek(fptmp, 0, SEEK_SET);
while (get_next_record(fptmp, 1, 0)) {
if (bu_strncmp(curr_rec[0], "CORD", 4))
continue;
get_coord_sys();
}
/* convert everything to BRL-CAD coordinate system */
i = 0;
while (convert_all_cs() || convert_all_pts()) {
i++;
if (i > 10) {
bu_exit(1, "Cannot convert to default coordinate system, check for circular definition\n");
}
}
mk_id(fpout, nastran_file);
/* get elements */
bu_fseek(fptmp, 0, SEEK_SET);
while (get_next_record(fptmp, 1, 0)) {
if (!bu_strncmp(curr_rec[0], "CBAR", 4))
get_cbar();
else if (!bu_strncmp(curr_rec[0], "CROD", 4))
get_cbar();
else if (!bu_strncmp(curr_rec[0], "CTRIA3", 6))
get_ctria3();
else if (!bu_strncmp(curr_rec[0], "CQUAD4", 6))
get_cquad4();
}
if (nmg_model) {
nmg_rebound(nmg_model, &tol);
if (polysolids)
mk_bot_from_nmg(fpout, "pshell.0", nmg_shell);
else
mk_nmg(fpout, "pshell.0", nmg_model);
}
BU_LIST_INIT(&head.l);
for (BU_LIST_FOR(psh, pshell, &pshell_head.l)) {
struct model *m;
char name[NAMESIZE+1];
if (!psh->s)
continue;
m = nmg_find_model(&psh->s->l.magic);
nmg_rebound(m, &tol);
nmg_fix_normals(psh->s, &tol);
if (psh->thick > tol.dist) {
nmg_model_face_fuse(m, &tol);
nmg_hollow_shell(psh->s, psh->thick*conv[units], 1, &tol);
}
sprintf(name, "pshell.%d", psh->pid);
if (polysolids)
mk_bot_from_nmg(fpout, name, psh->s);
else
mk_nmg(fpout, name, m);
mk_addmember(name, &head.l, NULL, WMOP_UNION);
}
if (BU_LIST_NON_EMPTY(&head.l)) {
mk_lfcomb(fpout, "shells", &head, 0);
mk_addmember("shells", &all_head.l, NULL, WMOP_UNION);
}
BU_LIST_INIT(&head.l);
for (BU_LIST_FOR(pbp, pbar, &pbar_head.l)) {
char name[NAMESIZE+1];
if (BU_LIST_IS_EMPTY(&pbp->head.l))
continue;
sprintf(name, "pbar_group.%d", pbp->pid);
mk_lfcomb(fpout, name, &pbp->head, 0);
mk_addmember(name, &head.l, NULL, WMOP_UNION);
}
if (BU_LIST_NON_EMPTY(&head.l)) {
mk_lfcomb(fpout, "pbars", &head, 0);
mk_addmember("pbars", &all_head.l, NULL, WMOP_UNION);
}
if (BU_LIST_NON_EMPTY(&all_head.l)) {
mk_lfcomb(fpout, "all", &all_head, 0);
}
wdb_close(fpout);
return 0;
}
int
brep(int entityno)
{
int sol_num; /* IGES solid type number */
int shell_de; /* Directory sequence number for a shell */
int orient; /* Orientation of shell */
int *void_shell_de; /* Directory sequence number for an void shell */
int *void_orient; /* Orientation of void shell */
int num_of_voids; /* Number of inner void shells */
struct model *m; /* NMG model */
struct nmgregion *r; /* NMG region */
struct shell **void_shells; /* List of void shells */
struct shell *s_outer; /* Outer shell */
struct iges_vertex_list *v_list;
struct iges_edge_list *e_list;
int i;
/* Acquiring Data */
if (dir[entityno]->param <= pstart) {
bu_log("Illegal parameter pointer for entity D%07d (%s)\n" ,
dir[entityno]->direct , dir[entityno]->name);
return 0;
}
Readrec(dir[entityno]->param);
Readint(&sol_num , "");
Readint(&shell_de , "");
Readint(&orient , "");
Readint(&num_of_voids , "");
if (num_of_voids) {
void_shell_de = (int *)bu_calloc(num_of_voids , sizeof(int) , "BREP: void shell DE's");
void_orient = (int *)bu_calloc(num_of_voids , sizeof(int) , "BREP: void shell orients");
void_shells = (struct shell **)bu_calloc(num_of_voids , sizeof(struct shell *) , "BREP: void shell pointers");
for (i = 0; i < num_of_voids; i++) {
Readint(&void_shell_de[i] , "");
Readint(&void_orient[i] , "");
}
} else {
void_shell_de = NULL;
void_orient = NULL;
void_shells = NULL;
}
/* start building */
m = nmg_mmr();
r = BU_LIST_FIRST(nmgregion, &m->r_hd);
/* Put outer shell in region */
if ((s_outer = Get_outer_shell(r , (shell_de - 1)/2)) == (struct shell *)NULL)
goto err;
ON_Brep* outer = ON_Brep::New();
if (Get_outer_brep(outer, (shell_de - 1)/2, orient))
goto err;
/* Put voids in */
for (i = 0; i < num_of_voids; i++) {
if ((void_shells[i] = Add_inner_shell(r, (void_shell_de[i] - 1)/2))
== (struct shell *)NULL)
goto err;
}
/* orient loops */
Orient_loops(r);
/* orient shells */
nmg_fix_normals(s_outer , &tol);
for (i = 0; i < num_of_voids; i++) {
nmg_fix_normals(void_shells[i] , &tol);
nmg_invert_shell(void_shells[i]);
}
if (do_bots) {
/* Merge all shells into one */
for (i = 0; i < num_of_voids; i++)
nmg_js(s_outer, void_shells[i], &tol);
/* write out BOT */
if (mk_bot_from_nmg(fdout, dir[entityno]->name, s_outer))
goto err;
} else {
/* Compute "geometry" for region and shell */
nmg_region_a(r , &tol);
/* Write NMG solid */
if (mk_nmg(fdout , dir[entityno]->name , m))
goto err;
}
if (num_of_voids) {
bu_free((char *)void_shell_de , "BREP: void shell DE's");
bu_free((char *)void_orient , "BREP: void shell orients");
bu_free((char *)void_shells , "brep: void shell list");
}
v_list = vertex_root;
while (v_list != NULL) {
bu_free((char *)v_list->i_verts , "brep: iges_vertex");
bu_free((char *)v_list , "brep: vertex list");
v_list = v_list->next;
}
vertex_root = NULL;
e_list = edge_root;
while (e_list != NULL) {
bu_free((char *)e_list->i_edge , "brep:iges_edge");
bu_free((char *)e_list , "brep: edge list");
e_list = e_list->next;
}
edge_root = NULL;
return 1;
err :
if (num_of_voids) {
bu_free((char *)void_shell_de , "BREP: void shell DE's");
bu_free((char *)void_orient , "BREP: void shell orients");
bu_free((char *)void_shells , "brep: void shell list");
}
nmg_km(m);
return 0;
}
