BEZIER_NODE *BicubicPatch::bezier_tree_builder(const ControlPoints *Patch, DBL u0, DBL u1, DBL v0, DBL v1, int depth, int& max_depth_reached)
{
ControlPoints Lower_Left, Lower_Right;
ControlPoints Upper_Left, Upper_Right;
BEZIER_CHILDREN *Children;
BEZIER_VERTICES *Vertices;
BEZIER_NODE *Node = create_new_bezier_node();
if (depth > max_depth_reached)
{
max_depth_reached = depth;
}
/* Build the bounding sphere for this subpatch. */
bezier_bounding_sphere(Patch, Node->Center, &(Node->Radius_Squared));
/*
* If the patch is close to being flat, then just perform
* a ray-plane intersection test.
*/
if (flat_enough(Patch))
{
/* The patch is now flat enough to simply store the corners. */
Node->Node_Type = BEZIER_LEAF_NODE;
Vertices = create_bezier_vertex_block();
Vertices->Vertices[0] = (*Patch)[0][0];
Vertices->Vertices[1] = (*Patch)[0][3];
Vertices->Vertices[2] = (*Patch)[3][3];
Vertices->Vertices[3] = (*Patch)[3][0];
Vertices->uvbnds[0] = u0;
Vertices->uvbnds[1] = u1;
Vertices->uvbnds[2] = v0;
Vertices->uvbnds[3] = v1;
Node->Data_Ptr = reinterpret_cast<void *>(Vertices);
}
else
{
if (depth >= U_Steps)
{
if (depth >= V_Steps)
{
/* We are at the max recursion depth. Just store corners. */
Node->Node_Type = BEZIER_LEAF_NODE;
Vertices = create_bezier_vertex_block();
Vertices->Vertices[0] = (*Patch)[0][0];
Vertices->Vertices[1] = (*Patch)[0][3];
Vertices->Vertices[2] = (*Patch)[3][3];
Vertices->Vertices[3] = (*Patch)[3][0];
Vertices->uvbnds[0] = u0;
Vertices->uvbnds[1] = u1;
Vertices->uvbnds[2] = v0;
Vertices->uvbnds[3] = v1;
Node->Count = 0;
Node->Data_Ptr = reinterpret_cast<void *>(Vertices);
}
else
{
bezier_split_up_down(Patch, &Lower_Left, &Upper_Left);
Node->Node_Type = BEZIER_INTERIOR_NODE;
Children = create_bezier_child_block();
Children->Children[0] = bezier_tree_builder(&Lower_Left, u0, u1, v0, (v0 + v1) / 2.0, depth + 1, max_depth_reached);
Children->Children[1] = bezier_tree_builder(&Upper_Left, u0, u1, (v0 + v1) / 2.0, v1, depth + 1, max_depth_reached);
Node->Count = 2;
Node->Data_Ptr = reinterpret_cast<void *>(Children);
}
}
else
{
if (depth >= V_Steps)
{
bezier_split_left_right(Patch, &Lower_Left, &Lower_Right);
Node->Node_Type = BEZIER_INTERIOR_NODE;
Children = create_bezier_child_block();
Children->Children[0] = bezier_tree_builder(&Lower_Left, u0, (u0 + u1) / 2.0, v0, v1, depth + 1, max_depth_reached);
Children->Children[1] = bezier_tree_builder(&Lower_Right, (u0 + u1) / 2.0, u1, v0, v1, depth + 1, max_depth_reached);
Node->Count = 2;
Node->Data_Ptr = reinterpret_cast<void *>(Children);
}
else
{
bezier_split_left_right(Patch, &Lower_Left, &Lower_Right);
bezier_split_up_down(&Lower_Left, &Lower_Left, &Upper_Left);
bezier_split_up_down(&Lower_Right, &Lower_Right, &Upper_Right);
Node->Node_Type = BEZIER_INTERIOR_NODE;
Children = create_bezier_child_block();
Children->Children[0] = bezier_tree_builder(&Lower_Left, u0, (u0 + u1) / 2.0, v0, (v0 + v1) / 2.0, depth + 1, max_depth_reached);
Children->Children[1] = bezier_tree_builder(&Upper_Left, u0, (u0 + u1) / 2.0, (v0 + v1) / 2.0, v1, depth + 1, max_depth_reached);
Children->Children[2] = bezier_tree_builder(&Lower_Right, (u0 + u1) / 2.0, u1, v0, (v0 + v1) / 2.0, depth + 1, max_depth_reached);
Children->Children[3] = bezier_tree_builder(&Upper_Right, (u0 + u1) / 2.0, u1, (v0 + v1) / 2.0, v1, depth + 1, max_depth_reached);
Node->Count = 4;
Node->Data_Ptr = reinterpret_cast<void *>(Children);
}
}
}
return (Node);
}
static BEZIER_NODE *bezier_tree_builder(BICUBIC_PATCH *Object, VECTOR (*Patch)[4][4], DBL u0, DBL u1, DBL v0, DBL v1, int depth, int& max_depth_reached)
{
VECTOR Lower_Left[4][4], Lower_Right[4][4];
VECTOR Upper_Left[4][4], Upper_Right[4][4];
BEZIER_CHILDREN *Children;
BEZIER_VERTICES *Vertices;
BEZIER_NODE *Node = create_new_bezier_node();
if (depth > max_depth_reached)
{
max_depth_reached = depth;
}
/* Build the bounding sphere for this subpatch. */
bezier_bounding_sphere(Patch, Node->Center, &(Node->Radius_Squared));
/*
* If the patch is close to being flat, then just perform
* a ray-plane intersection test.
*/
if (flat_enough(Object, Patch))
{
/* The patch is now flat enough to simply store the corners. */
Node->Node_Type = BEZIER_LEAF_NODE;
Vertices = create_bezier_vertex_block();
Assign_Vector(Vertices->Vertices[0], (*Patch)[0][0]);
Assign_Vector(Vertices->Vertices[1], (*Patch)[0][3]);
Assign_Vector(Vertices->Vertices[2], (*Patch)[3][3]);
Assign_Vector(Vertices->Vertices[3], (*Patch)[3][0]);
Vertices->uvbnds[0] = u0;
Vertices->uvbnds[1] = u1;
Vertices->uvbnds[2] = v0;
Vertices->uvbnds[3] = v1;
Node->Data_Ptr = (void *)Vertices;
}
else
{
if (depth >= Object->U_Steps)
{
if (depth >= Object->V_Steps)
{
/* We are at the max recursion depth. Just store corners. */
Node->Node_Type = BEZIER_LEAF_NODE;
Vertices = create_bezier_vertex_block();
Assign_Vector(Vertices->Vertices[0], (*Patch)[0][0]);
Assign_Vector(Vertices->Vertices[1], (*Patch)[0][3]);
Assign_Vector(Vertices->Vertices[2], (*Patch)[3][3]);
Assign_Vector(Vertices->Vertices[3], (*Patch)[3][0]);
Vertices->uvbnds[0] = u0;
Vertices->uvbnds[1] = u1;
Vertices->uvbnds[2] = v0;
Vertices->uvbnds[3] = v1;
Node->Data_Ptr = (void *)Vertices;
}
else
{
bezier_split_up_down(Patch, (VECTOR(*)[4][4])Lower_Left, (VECTOR(*)[4][4])Upper_Left);
Node->Node_Type = BEZIER_INTERIOR_NODE;
Children = create_bezier_child_block();
Children->Children[0] = bezier_tree_builder(Object, (VECTOR(*)[4][4])Lower_Left, u0, u1, v0, (v0 + v1) / 2.0, depth + 1, max_depth_reached);
Children->Children[1] = bezier_tree_builder(Object, (VECTOR(*)[4][4])Upper_Left, u0, u1, (v0 + v1) / 2.0, v1, depth + 1, max_depth_reached);
Node->Count = 2;
Node->Data_Ptr = (void *)Children;
}
}
else
{
if (depth >= Object->V_Steps)
{
bezier_split_left_right(Patch, (VECTOR(*)[4][4])Lower_Left, (VECTOR(*)[4][4])Lower_Right);
Node->Node_Type = BEZIER_INTERIOR_NODE;
Children = create_bezier_child_block();
Children->Children[0] = bezier_tree_builder(Object, (VECTOR(*)[4][4])Lower_Left, u0, (u0 + u1) / 2.0, v0, v1, depth + 1, max_depth_reached);
Children->Children[1] = bezier_tree_builder(Object, (VECTOR(*)[4][4])Lower_Right, (u0 + u1) / 2.0, u1, v0, v1, depth + 1, max_depth_reached);
Node->Count = 2;
Node->Data_Ptr = (void *)Children;
}
else
{
bezier_split_left_right(Patch, (VECTOR(*)[4][4])Lower_Left, (VECTOR(*)[4][4])Lower_Right);
bezier_split_up_down((VECTOR(*)[4][4])Lower_Left, (VECTOR(*)[4][4])Lower_Left, (VECTOR(*)[4][4])Upper_Left);
bezier_split_up_down((VECTOR(*)[4][4])Lower_Right, (VECTOR(*)[4][4])Lower_Right, (VECTOR(*)[4][4])Upper_Right);
Node->Node_Type = BEZIER_INTERIOR_NODE;
Children = create_bezier_child_block();
Children->Children[0] = bezier_tree_builder(Object, (VECTOR(*)[4][4])Lower_Left, u0, (u0 + u1) / 2.0, v0, (v0 + v1) / 2.0, depth + 1, max_depth_reached);
Children->Children[1] = bezier_tree_builder(Object, (VECTOR(*)[4][4])Upper_Left, u0, (u0 + u1) / 2.0, (v0 + v1) / 2.0, v1, depth + 1, max_depth_reached);
Children->Children[2] = bezier_tree_builder(Object, (VECTOR(*)[4][4])Lower_Right, (u0 + u1) / 2.0, u1, v0, (v0 + v1) / 2.0, depth + 1, max_depth_reached);
Children->Children[3] = bezier_tree_builder(Object, (VECTOR(*)[4][4])Upper_Right, (u0 + u1) / 2.0, u1, (v0 + v1) / 2.0, v1, depth + 1, max_depth_reached);
Node->Count = 4;
Node->Data_Ptr = (void *)Children;
}
}
}
return (Node);
}