void set_diameter(node* root)
{
if (!root)
return ;
set_diameter(root->left);
int l = find_height(root->left);
int r = find_height(root->right);
if (l+r+1 > max_diameter)
max_diameter = l+r+1+1;// 1 is added to count the parent node in that path, and another 1 is added because when we are calculating height of LST and RST then both of them are starting from 0, but while counting the diameter we need to count the number of nodes, so we add one.
set_diameter(root->right);
}
Surface& Surface::operator=(const Surface& rSurf)
{
_pMaterialNext=0;
_pMaterialPrev=0;
set_type(rSurf.type());
set_conic(rSurf.conic());
set_radius_curvature(rSurf.radius_curvature());
set_R4(rSurf.R4());
set_R6(rSurf.R6());
set_R8(rSurf.R8());
set_R10(rSurf.R10());
set_comment(rSurf.comment());
set_x(rSurf.x());
set_y(rSurf.y());
set_z(rSurf.z());
set_diameter(rSurf.diameter());
set_auto_diameter(rSurf.get_auto_diameter());
set_inner_diameter(rSurf.inner_diameter());
set_auto_inner_diameter(rSurf.get_auto_inner_diameter());
return *this;
}
int main()
{
node* root = NULL;
/*
root = insert(root, 100);
root = insert(root, 70);
root = insert(root, 110);
root = insert(root, 50);
root = insert(root, 75);
root = insert(root, 78);
root = insert(root, 80);
root = insert(root, 82);
root = insert(root, 45);
root = insert(root, 40);
root = insert(root, 57);
*/
/*
root = insert(root, 10);
root = insert(root, 5);
root = insert(root, 15);
root = insert(root, 16);
root = insert(root, 4);
*/
/*
root = insert(root, 2);
root = insert(root, 1);
root = insert(root, 3);
*/
root = insert(root, 10);
root = insert(root, 5);
root = insert(root, 30);
root = insert(root, 4);
root = insert(root, 6);
//inorder(root);
std::cout<<std::endl;
set_diameter(root);
std::cout<<"Diameter = "<<max_diameter;
std::cout<<std::endl;
return 0;
}
void
use_air(char *buffer, com_table *ctp, struct rt_i *UNUSED(rtip))
{
int new_use = 0; /* current position on the *buffer */
char response[128];
char *rp = response;
char db_title[TITLE_LEN+1]; /* title from MGED database */
struct rt_i *my_rtip;
extern char *db_name; /* Name of MGED database file */
while (isspace((int)*buffer))
++buffer;
if (*buffer == '\0') {
/* display current value of use_of_air */
bu_log("use_air = %d\n", ap.a_rt_i->useair);
return;
}
if (!isdigit((int)*buffer)) {
com_usage(ctp);
return;
}
while (isdigit((int)*buffer)) {
new_use *= 10;
new_use += *buffer++ - '0';
}
if (new_use && (new_use != 1)) {
bu_log("Warning: useair=%d specified, will set to 1\n",
new_use);
new_use = 1;
}
if (rti_tab[new_use] == RTI_NULL) {
bu_log(" Air %s in the current directory of database objects.\n",
new_use ? "is not included" : "is included");
bu_log(
" To set useair=%d requires building/prepping another directory.\n",
new_use);
bu_log(" Do you want to do that now (y|n)[n]? ");
bu_fgets(response, sizeof(response), stdin);
while ((*rp == ' ') || (*rp == '\t'))
++rp;
if ((*rp != 'y') && (*rp != 'Y')) {
bu_log("useair remains %d\n", ap.a_rt_i->useair);
return;
}
bu_log("Building the directory...");
if ((my_rtip = rt_dirbuild(db_name, db_title, TITLE_LEN)) == RTI_NULL) {
bu_log("Could not load file %s\n", db_name);
printusage();
bu_exit(1, NULL);
}
rti_tab[new_use] = my_rtip;
my_rtip->useair = new_use;
my_rtip->rti_save_overlaps = (overlap_claims > 0);
bu_log("Prepping the geometry...");
do_rt_gettrees(my_rtip, NULL, 0, &need_prep);
}
ap.a_rt_i = rti_tab[new_use];
ap.a_resource = &res_tab[new_use];
set_diameter(ap.a_rt_i);
}
void
shoot(char *UNUSED(buffer), com_table *UNUSED(ctp), struct rt_i *rtip)
{
int i;
double bov = 0.0; /* back out value */
extern void init_ovlp();
if (!rtip)
return;
if (need_prep) {
rt_clean(rtip);
do_rt_gettrees(rtip, NULL, 0, &need_prep);
}
if (do_backout) {
point_t ray_point;
vect_t ray_dir;
vect_t center_bsphere;
fastf_t dist_to_target;
vect_t dvec;
fastf_t delta;
for (i = 0; i < 3; ++i) {
ray_point[i] = target(i);
ray_dir[i] = direct(i);
}
if (bsphere_diameter < 0)
set_diameter(rtip);
/*
* calculate the distance from a plane normal to the ray direction through the center of
* the bounding sphere and a plane normal to the ray direction through the aim point.
*/
VADD2SCALE(center_bsphere, rtip->mdl_max, rtip->mdl_min, 0.5);
dist_to_target = DIST_PT_PT(center_bsphere, ray_point);
VSUB2(dvec, ray_point, center_bsphere);
VUNITIZE(dvec);
delta = dist_to_target*VDOT(ray_dir, dvec);
/*
* this should put us about a bounding sphere radius in front of the bounding sphere
*/
bov = bsphere_diameter + delta;
}
for (i = 0; i < 3; ++i) {
target(i) = target(i) + (bov * -direct(i));
ap.a_ray.r_pt[i] = target(i);
ap.a_ray.r_dir[i] = direct(i);
}
if (nirt_debug & DEBUG_BACKOUT) {
bu_log("Backing out %g units to (%g %g %g), shooting dir is (%g %g %g)\n",
bov * base2local,
ap.a_ray.r_pt[0] * base2local,
ap.a_ray.r_pt[1] * base2local,
ap.a_ray.r_pt[2] * base2local,
V3ARGS(ap.a_ray.r_dir));
}
init_ovlp();
(void) rt_shootray(&ap);
/* Restore pre-backout target values */
if (do_backout) {
for (i = 0; i < 3; ++i) {
target(i) = target(i) - (bov * -direct(i));
}
}
}
int
main(int argc, char *argv[])
{
struct rt_i *rtip = NULL;
char db_title[TITLE_LEN+1];/* title from MGED file */
const char *tmp_str;
extern char local_u_name[65];
extern double base2local;
extern double local2base;
FILE *fPtr;
int Ch; /* Option name */
int mat_flag = 0; /* Read matrix from stdin? */
int use_of_air = 0;
int print_ident_flag = 1;
char ocastring[1024] = {0};
struct bu_list script_list; /* For -e and -f options */
struct script_rec *srp;
extern outval ValTab[];
/* from if.c, callback functions for overlap, hit, and miss shots */
int if_overlap(struct application *, struct partition *, struct region *, struct region *, struct partition *);
int if_hit(struct application *, struct partition *, struct seg *);
int if_miss(struct application *);
BU_LIST_INIT(&script_list);
bu_setprogname(argv[0]);
ocname[OVLP_RESOLVE] = "resolve";
ocname[OVLP_REBUILD_FASTGEN] = "rebuild_fastgen";
ocname[OVLP_REBUILD_ALL] = "rebuild_all";
ocname[OVLP_RETAIN] = "retain";
*ocastring = '\0';
bu_optind = 1; /* restart */
/* Handle command-line options */
while ((Ch = bu_getopt(argc, argv, OPT_STRING)) != -1) {
if (bu_optopt == '?') Ch='h';
switch (Ch) {
case 'A':
attrib_add(bu_optarg, &need_prep);
break;
case 'B':
rt_bot_minpieces = atoi(bu_optarg);
break;
case 'T':
setenv("LIBRT_BOT_MINTIE", bu_optarg, 1);
break;
case 'b':
do_backout = 1;
break;
case 'E':
if (nirt_debug & DEBUG_SCRIPTS)
show_scripts(&script_list, "before erasure");
while (BU_LIST_WHILE(srp, script_rec, &script_list)) {
BU_LIST_DEQUEUE(&(srp->l));
free_script(srp);
}
if (nirt_debug & DEBUG_SCRIPTS)
show_scripts(&script_list, "after erasure");
break;
case 'e':
enqueue_script(&script_list, READING_STRING, bu_optarg);
if (nirt_debug & DEBUG_SCRIPTS)
show_scripts(&script_list, "after enqueueing a literal");
break;
case 'f':
enqueue_script(&script_list, READING_FILE, bu_optarg);
if (nirt_debug & DEBUG_SCRIPTS)
show_scripts(&script_list, "after enqueueing a file name");
break;
case 'L':
listformats();
bu_exit(EXIT_SUCCESS, NULL);
case 'M':
mat_flag = 1;
break;
case 'O':
sscanf(bu_optarg, "%1024s", ocastring);
break;
case 's':
silent_flag = SILENT_YES; /* Positively yes */
break;
case 'v':
silent_flag = SILENT_NO; /* Positively no */
break;
case 'x':
sscanf(bu_optarg, "%x", (unsigned int *)&RTG.debug);
break;
case 'X':
sscanf(bu_optarg, "%x", (unsigned int *)&nirt_debug);
break;
case 'u':
if (sscanf(bu_optarg, "%d", &use_of_air) != 1) {
(void) fprintf(stderr,
"Illegal use-air specification: '%s'\n", bu_optarg);
return 1;
}
break;
case 'H':
if (sscanf(bu_optarg, "%d", &print_ident_flag) != 1) {
(void) fprintf(stderr,
"Illegal header output option specified: '%s'\n", bu_optarg);
return 1;
}
break;
default:
printusage();
bu_exit (Ch != 'h', NULL);
}
} /* end while getopt */
if (argc - bu_optind < 2) {
printusage();
return 1;
}
if (isatty(0)) {
if (silent_flag != SILENT_YES)
silent_flag = SILENT_NO;
} else {
/* stdin is not a TTY */
if (silent_flag != SILENT_NO)
silent_flag = SILENT_YES;
}
if (silent_flag != SILENT_YES && print_ident_flag)
(void) fputs(brlcad_ident("Natalie's Interactive Ray Tracer"), stdout);
if (use_of_air && (use_of_air != 1)) {
fprintf(stderr,
"Warning: useair=%d specified, will set to 1\n", use_of_air);
use_of_air = 1;
}
switch (*ocastring) {
case '\0':
overlap_claims = OVLP_RESOLVE;
break;
case '0':
case '1':
case '2':
case '3':
if (ocastring[1] == '\0') {
sscanf(ocastring, "%d", &overlap_claims);
} else {
fprintf(stderr,
"Illegal overlap_claims specification: '%s'\n", ocastring);
return 1;
}
break;
case 'r':
if (BU_STR_EQUAL(ocastring, "resolve"))
overlap_claims = OVLP_RESOLVE;
else if (BU_STR_EQUAL(ocastring, "rebuild_fastgen"))
overlap_claims = OVLP_REBUILD_FASTGEN;
else if (BU_STR_EQUAL(ocastring, "rebuild_all"))
overlap_claims = OVLP_REBUILD_ALL;
else if (BU_STR_EQUAL(ocastring, "retain"))
overlap_claims = OVLP_RETAIN;
else {
fprintf(stderr,
"Illegal overlap_claims specification: '%s'\n", ocastring);
return 1;
}
break;
default:
fprintf(stderr,
"Illegal overlap_claims specification: '%s'\n", ocastring);
return 1;
}
db_name = argv[bu_optind];
/* build directory for target object */
if (silent_flag != SILENT_YES) {
printf("Database file: '%s'\n", db_name);
printf("Building the directory...");
}
if ((rtip = rt_dirbuild(db_name, db_title, TITLE_LEN)) == RTI_NULL) {
fflush(stdout);
fprintf(stderr, "Could not load file %s\n", db_name);
return 1;
}
rti_tab[use_of_air] = rtip;
rti_tab[1 - use_of_air] = RTI_NULL;
rtip->useair = use_of_air;
rtip->rti_save_overlaps = (overlap_claims > 0);
++bu_optind;
do_rt_gettrees(rtip, argv + bu_optind, argc - bu_optind, &need_prep);
/* Initialize the table of resource structures */
rt_init_resource(&res_tab, 0, rtip);
/* initialization of the application structure */
RT_APPLICATION_INIT(&ap);
ap.a_hit = if_hit; /* branch to if_hit routine */
ap.a_miss = if_miss; /* branch to if_miss routine */
ap.a_overlap = if_overlap;/* branch to if_overlap routine */
ap.a_logoverlap = rt_silent_logoverlap;
ap.a_onehit = 0; /* continue through shotline after hit */
ap.a_resource = &res_tab;
ap.a_purpose = "NIRT ray";
ap.a_rt_i = rtip; /* rt_i pointer */
ap.a_zero1 = 0; /* sanity check, sayth raytrace.h */
ap.a_zero2 = 0; /* sanity check, sayth raytrace.h */
ap.a_uptr = (void *)a_tab.attrib;
/* initialize variables */
azimuth() = 0.0;
elevation() = 0.0;
direct(X) = -1.0;
direct(Y) = 0.0;
direct(Z) = 0.0;
grid(HORZ) = 0.0;
grid(VERT) = 0.0;
grid(DIST) = 0.0;
grid2targ();
set_diameter(rtip);
/* initialize the output specification */
default_ospec();
/* initialize NIRT's local units */
base2local = rtip->rti_dbip->dbi_base2local;
local2base = rtip->rti_dbip->dbi_local2base;
tmp_str = bu_units_string(local2base);
if (tmp_str) {
bu_strlcpy(local_u_name, tmp_str, sizeof(local_u_name));
} else {
bu_strlcpy(local_u_name, "Unknown units", sizeof(local_u_name));
}
if (silent_flag != SILENT_YES) {
printf("Database title: '%s'\n", db_title);
printf("Database units: '%s'\n", local_u_name);
printf("model_min = (%g, %g, %g) model_max = (%g, %g, %g)\n",
rtip->mdl_min[X] * base2local,
rtip->mdl_min[Y] * base2local,
rtip->mdl_min[Z] * base2local,
rtip->mdl_max[X] * base2local,
rtip->mdl_max[Y] * base2local,
rtip->mdl_max[Z] * base2local);
}
/* Run the run-time configuration file, if it exists */
if ((fPtr = fopenrc()) != NULL) {
interact(READING_FILE, fPtr, rtip);
fclose(fPtr);
}
/* Run all scripts specified on the command line */
run_scripts(&script_list, rtip);
/* Perform the user interface */
if (mat_flag) {
read_mat(rtip);
return 0;
} else {
interact(READING_FILE, stdin, rtip);
}
return 0;
}
