the present invention can be described as a single parameterized scheme , which is a blend of an approximating and an interpolating scheme . this blend does a better job at preserving a silhouette and a volume of the meshes . when a surface is conventionally subdivided , the subdivided surface is then typically smoothed to maintain the shape of the base mesh . this smoothing nevertheless typically results in a distortion of the original shape . the present invention performs a push back of displaced vertices in order to compensate for this distortion . the present invention performs the scheme in two operations : a linear subdivision , which increase the number of faces , vertices , edges ; and then a smoothing operation , where the vertex positions are modified , regardless on how vertices were created . the prior art typically did everything in one operation . fig2 illustrates the subdividing process . a base mesh is comprised of a triangle 200 and a square 202 . the triangle 200 is divided into a subdivided triangle 204 with four pieces . the square 202 is divided into a subdivided square 206 with four pieces . a smoothing algorithm is then applied to the subdivided triangle 204 and the subdivided square 206 to produce a subdivided smooth surface 208 . this process can be repeated a preset number of times until the desired effect is achieved . note that the present invention allows for both triangles and squares to be subdivided in the same surface . this can be accomplished using methods described later on . note that the iterative process described above may result in a surface smaller than the original mesh , as illustrated in fig1 . to compensate for this effect , a “ push back ” operation can be implemented which moves each original vertex back towards the vertex &# 39 ; s original position ( before smoothing ) by an amount controlled by the user . the push back operation will now be described . fig3 illustrates three segments comprising a polyline made of a b - spline curve . polyline 300 comprises vertices 302 , 304 , 306 , and 308 , which define segments 310 , 312 , and 314 , respectively , fig4 illustrates the segments in fig3 linearly divided into a 6 - segment polyline . this division is typically performed in a subdivision surface in order to increase the number of pieces of a segment ( or surface ) in order to produce smoother results . the polyline 400 comprises vertex 402 ( corresponding to vertex 302 ), new vertex 404 , vertex 406 ( corresponds to vertex 304 ), new vertex 408 , vertex 410 ( corresponds to vertex 306 ), new vertex 412 , and vertex 414 ( corresponds to vertex 308 ). the division can be performed in any conventional manner . fig5 illustrates smoothing of the polyline . the smoothing is typically performed so that the subdivided polyline is smoother and rounder . the original polyline 500 comprises points 502 , 504 , 506 , 508 , 510 , 512 and 514 . the original polyline 514 can be smoothed into a smoothed polyline 516 , comprising points 518 , 520 , 522 , 524 , 526 , 528 , and 530 . the smoothing can be performed using any conventional technique , such as smoothing with a [ 1 2 1 ] mask as illustrated . fig6 illustrates a determination of displacement from the original points . the displacement is determined so that it can be used to subsequently calculate an interpolated displacement , which is ultimately used to compensate for any distortion resulting from the smoothing . subdivided original polyline 600 ( before smoothing ) comprises points 602 , 604 , 606 , 608 , 610 , 612 , and 614 . the smoothed polyline 616 comprises points 618 , 620 , 622 , 624 , 626 , 628 , and 630 . displacement 632 is determined which is the distance between subdivided point 602 and its corresponding smoothed point 618 . displacement 634 is determined which is the distance between subdivided point 606 and its corresponding smoothed point 622 . displacement 632 is determined which is the distance between subdivided point 610 and its corresponding smoothed point 626 . displacement 634 is determined which is the distance between subdivided point 614 and its corresponding smoothed point 630 . displacements are only computed between original points ( 302 , 304 , 306 , 308 ) and the corresponding ones in the subdivided and smoothed surface . the other displacements are not meaningful . missing displacement vectors are interpolated as described in the paragraph below . fig7 illustrates linear interpolation of the displacement . the linear interpolation of the displacement is calculated so it can be subsequently used to compensate for any distortion resulting from the smoothing . interpolated displacement vectors 700 , 702 , and 704 are computed based on the displacements illustrated in fig6 . the interpolation can be performed by any conventional method , typically a linear method . fig8 illustrates moving points ( or “ push back ”) to compensate for distortion during the smoothing . subdivided polyline 800 comprises points 802 , 804 , 806 , 808 , 810 , 812 , and 814 . the smoothed polyline 816 comprises points 818 , 820 , 822 , 824 , 826 , 828 , and 830 . displacements 832 , 834 , 836 , and 838 are determined based on distances between points on the original polyline and their corresponding smoothed subdivided points . interpolated displacement vectors 838 , 840 , and 842 are interpolated displacements based on the determined displacements . points 844 , 846 , and 848 correspond to points 804 , 808 , and 812 after being moved by the corresponding interpolated displacement vector ( 838 , 840 , 842 ). thus , the corrected polyline 850 after the push back operation comprises the points 802 , 844 , 806 , 846 , 810 , 848 , and 814 . fig9 illustrates how the push back effect can be weighted by a volume parameter , herein defined as α . the volume parameter allows a user to specify how much “ push - back effect ” is desired . generally , the greater the value used for α , the larger the volume of the resulting polyline . polyline 900 ( corresponding to polyline 816 from fig8 ) illustrates the case where α = 0 . in this case , no interpolation is reflected , as this is a pure approximation scheme . polyline 902 ( corresponding to polyline 850 from fig8 ) illustrates the case where α = 1 , whereas this is a pure interpolation scheme ( in which the original vertices do not move ). polyline 904 illustrates the case where α = ½ , which comprises a “ blended ” interpolation and approximation scheme . generally speaking , setting α = ½ typically produces the best effects , although of course the user will want to adjust this parameter to suit his preferences . fig1 a , 10 b , and 10 c illustrate the effects of changing α . fig1 a , 10 b , and 10 c illustrate refinement operations 1 , 2 and the limit curve , using degree three interpolation and a division value d - 2 . in fig1 a , α = 0 . note how the limit curve 1002 is smaller than the base mesh 1004 . in fig1 b , α = ½ . note how the limit curve 1006 more closely follows the base mesh 1008 . in fig1 c , α = 1 . note how the limit curve 1010 has a larger volume than the base mesh 1010 . note that the present invention is equally beneficial for three - dimensional ( or higher ) surfaces as well . fig1 illustrates how the push back effect can also be weighted by a rounding parameter , herein defined as β . β is a parameter used to bulge flat areas around segment centers . note that displacements 1100 and 1102 have a sharp angle between them , and an interpolated vector between them may not produce ideal results . an interpolated vector between two vectors with a very sharp angle between them may approach zero . because of this , areas toward the ends of spherical surfaces may look “ boxy ” ( see fig1 a for an example ). thus we introduce a rounding parameter β which weights ( normalizes ) the interpolation results when around segment centers . 1104 illustrates the interpolation result when β = 0 , while 1106 illustrates the interpolation result when β = 1 . more on the rounding parameter will be described below . fig1 illustrates a flowchart of one possible implementation of the methods of the present invention , according to an embodiment of the present invention . a computer implementing the method first subdivides 1200 . after subdividing , the computer smoothes 1202 the subdivided surface . the smoothing operation can be repeated n times , as set by the user . after the smoothing is performed , the computer computes 1204 displacement of the subdivided points and the smoothed points . the computer then interpolates 1206 the displacement between the original subdivided points and the smoothed points . the computer then pushes back 1208 points based on the interpolated displacement . repeating n times can include either 1202 only as shown or all 1202 – 1208 operations . in the later case , the push back will be smoother , as it &# 39 ; s done as many small adjustments , as opposed to one big adjustment at the end . this is generally preferable . repeating only 1202 will be slightly faster to compute , but the surface typically will have a slightly lower quality . optionally , the user can set a push back weight ( or volume parameter ). the user can also optionally set a rounding ( or bulging parameter ). the prior art limits a number of faces into which a surface can be subdivided . for example , the prior art allows a square to be subdivided into a grid of ( n × n ). thus , a square can be divided into 2 × 2 = 4 faces , 4 × 4 = 16 faces , etc . the prior art always subdivided a face into 4 pieces , whether for triangle or quadrilateral schemes . the present invention allows dividing into p × p pieces at each step , thus giving a total of ( p × p ) m . continuous subdividing into faces of a power of 4 results in a growth which some users may consider too high . in another embodiment of the present invention , a line ( or surface ) can be subdivided into an arbitrary number of pieces ( or faces ). a surface can be subdivided into an arbitrary number of pieces not limited to a power of four like the prior art . fig1 illustrates linearly subdividing an original polyline by an odd number , particularly into d = 5 pieces . the top segment of polyline 1500 is divided into pieces 1502 , 1504 , 1506 , 1508 , and 1510 . the center segment is needed for the subdivision operation . the smoothing needs the first neighbors too . fig1 illustrates using a smoothing mask on the subdivided pieces in the odd number subdivisions . a [ 1 2 . . . d . . . 2 1 ] smoothing mask is applied . points along the polyline are numbered 1 - 1400 , 2 - 1402 , 3 - 1404 , 4 - 1406 , 5 - 1408 , 4 - 1410 , 3 - 1412 , 2 - 1414 , 1 - 1416 . the mask requires the position of some of the vertices from the next segments . there are more , but only the ones overlapping with the mask are drawn . fig1 illustrates linear interpolation of push back vectors for the odd number subdivisions . the linear interpolation can be performed by any conventional method as described herein or elsewhere . points 1500 , 1502 , 1504 , 1506 , 1508 , and 1510 are the smoothed points corresponding to points 1404 , 1406 , 1408 , 1410 , 1412 , and 1414 . thus , we have now subdivided the original polyline into 5 pieces ( 1512 , 1514 , 1516 , 1518 , 1520 ). the prior art only divided in 2 , and explicitly defined the resulting position . we will now elaborate on the above methods and describe their implementation in more detail . one type of smoothing method is described in the article titled “ on subdivision schemes generalizing uniform b - spline surfaces of arbitrary degree ,” by jos stam . b - splines of odd degree p can be subdivided by first linearly subdividing the control mesh and then performing smoothing operations . each operation involves averaging a vertex with its immediate neighbors using the for example , from the last line we obtain the subdivision masks for b - spline curves of degree three . every old vertex is updated using the masks . the crucial observation is that these two masks are obtained by simply applying the mask to the second row , which corresponds to the subdivision rules of linear subdivision . in fact , this construction is easily generalized to any number of subdivisions d . in this case we generalize the pascal triangle to take the advantage of the d elements in the row directly above it : this allows us to compute the corresponding masks for these subdivision schemes . first , linearly subdivide each segment into d pieces , then smooth each vertex using the mask 1 d 2 ⁡ [ 1 , 2 , … ⁢ , d , … ⁢ , 2 , 1 ] . in the limit this process generates b - spline curves of degree 2m + 1 if the smoothing is applied m times . note that any conventional smoothing method can be used to present invention we add the push back operation that updates the position of the vertices after smoothing , order to limit the amount of the shrinking characteristic of approximating subdivision schemes . each original vertex is moved back towards its original position by an amount controlled by the user . newly introduced vertices are also adjusted by linear interpolation of the adjusted original vertices . we denote by p j the new vertices obtained by subdividing the original vertices throughout the following description , “:=” denotes assignment , while “=” denotes a true equality of two quantities . in these notations , the first operation is ( note that d = number of pieces , and k is a temporary variable representing a distance from the closest original vertex on the left . this operation is followed by a smoothing operation that modifies the vertices p i : finally , the smoothing operation is followed by a push - back of these new vertices : evaluations are done in parallel in a conventional “ jacobi manner ” to avoid any side effects . in practice this requires the use of an intermediate array to store the vertices &# 39 ; positions . the volume parameter a controls the transition from approximation to interpolation . when α = 0 there is no push - back operation and the subdivision scheme produces uniform b - splines in the limit . on the other hand , when α = 1 the scheme described herein is interpolatory . fig4 shows the influence of the parameter a on the subdivided control vertices after several refinements . when the degree p is greater than three , the smoothing and push - back operations are repeated times . in particular , when p = 3 and d = 2 , one smoothing and one push - back operation is performed . in this case we can explicitly write down the subdivision matrix applied to five consecutive control vertices : in particular , when α = 1 , the p 2i are moved back exactly to their original position p i , and we obtain the well known four point interpolation scheme , with the surface case is similar to the curve one : we perform one bilinear subdivision operation followed by a smoothing operation . we first define rules for binary subdivision schemes when d = 2 . we introduce the following notations . the number of elements in a set a is denoted by a #. the vertices of the mesh before a subdivision operation are denoted by v 1 i . during a subdivision operation these vertices are transformed into new vertices v i . at the same time new vertices e i are introduced by splitting each edge , and new vertices f i are introduced for each face as in fig5 . let p be a vertex of the mesh , then ε ( p ) is the set containing all the vertices sharing an edge with p . the set c ( p ) contains the “ corner vertices :” the vertices sharing a face with p not in ε ( p ). fig1 a and 16b illustrate subdivisions of quadrilaterals and triangles . fig1 a illustrates an initial mesh 1600 . fig1 b illustrates a subdivided mesh 1602 , that is the mesh of fig1 a after one subdivision . each subdivision produces v , e , and f vertex types . faces with 5 and more vertices use the quadrilateral subdivision rule . to illustrate the above definitions refer to fig1 a and 16b , wherein ε ( v 2 )={ e 2 , e 3 , e 5 }, and c ( e 3 )={ e 1 , e 2 }. we will now focus entirely on quadrilateral schemes . however , triangular schemes can be treated in a similar way , with the exception that there are no face vertices f i , and c ( v ) is always empty . the subdivision operation should distinguish between these two types of faces . the stam article previously mentioned provides different smoothing rules of the vertices that result in uniform b - spline surfaces in the limit on the regular part of the mesh . the simplest smoothing algorithm which corresponds to “ repeated averaging ” replaces each vertex by a weighted average of its direct neighbors : catmull - clark surfaces are obtained with a different choice for the weights : we observe that formulae 9 and 10 are identical when n i = 4 . this comes as no surprise since both of these schemes produce uniform b - spline surfaces on regular meshes ( n i = 4 everywhere ). we further observe that , when n i ≠ 4 , the catmull - clark rule can be obtained by following formula 9 with an adjustment of all the extraordinary vertices : the parameter γ allows us to interpolate between the two schemes . not only does this adjustment unify these two schemes , but it also simplifies the implementation of the catmull - clark subdivision : a simple smoothing followed by a vertex update operation . following the curve case , the simplest push - back operation is to compute the differences δ i between v i 1 , and v i , followed by a ( bi -) linear interpolation of these differences for the new e i and f i vertices . however , close to very sharp corners the scheme of the present invention tends to create flat areas around the face centers . the reason is that a bilinear interpolation of vectors of the same length with different angles produces smaller vectors at the center of the faces . this is a well known artifact of certain renderers , which do not renormalize vertex normals after interpolation , and consequently produce darker areas in the face centers . we can fix this problem by introducing a renormalization operation for the interpolation of the δ vectors . this is achieved by interpolating the length and direction of the δ vectors separately . to smooth the transition between these new rules and the ones without the normalization , we introduce a rounding factor parameter β . when β = 0 there is no renormalization , while when β = 1 the lengths of the δ are exactly interpolated . in the case β ≠ 0 these subdivision rules do not reproduce uniform b - splines on the regular part of the mesh in the limit . this doesn &# 39 ; t matter since we do not use these rules to generate limit surfaces . fig1 a , 17 b , and 17 c illustrate how meshes are deformed when β is increased . fig1 a illustrates a case where β = 0 . note how the sphere 1700 has a square - like shape . fig1 b illustrates a case where β = ½ . note how the sphere 1702 has a rounder shape . fig1 c illustrates a case where β = 1 . note how the shape of the sphere 1704 is affected . for a cube , produces the most “ rounded ” meshes . note that in this example we have γ = 0 ( no catmull - clark ) to emphasize the flattening problem . the present invention allows for subdivision into an arbitrary number of pieces . for regular meshes the corresponding limited surfaces σ ( s , t ) are equal to a tensor product of uniform b - spline curves . therefore , the subdivision scheme for these surfaces is simply a linear subdivision operation followed by a smoothing operation with a mask equal to the tensor product of the mask these rules can be naturally extended to extraordinary regions . ordinary regions are made either of triangles only , where all vertices have 6 neighbors ; or made of quadrilaterals only , where all vertices have exactly 4 neighbors . in practice , it turns out that it is easier to decompose the smoothing operation into two simple averaging operations . the averaging operation is different depending on whether d is odd or even . in the odd case we replace each vertex by a simple average of its k - ring neighborhood , where , where when d is even , each averaging operation replaces each face with a vertex that is the average of the k - ring of vertices surrounding it , where the new vertices after this operation form the dual of the initial mesh . in practice , however , the dual is never explicitly computed since the averaging operation is always performed twice ( an even amount in general ). after two dualizations the vertices are again “ in place .” more formally , let v k ( v i ) be the set of all vertices which can be reached from v i by traversing at most k faces and let f k ( v i ) denote the corresponding set of the faces traversed . fig1 a , 18 b , and 18 c illustrate examples of vertex ( or face ) neighborhoods . vertex ( or face ) neighborhoods v k [ f k ] are defined by adding one more ring to the previous set . for some examples . we also define a set of face neighborhoods by v k ( f i )= u vεfi v k ( v ). using these definitions we can explicitly state the smoothing operations . when d is odd , we apply the following rule p − 1 times , where when d is even the procedure typically only works for odd degrees p . we set the neighborhood to in practice , we prefer odd degrees so that no constraint is necessary on the number of subdivisions d . the catmull - clark correction operation defined by formula 11 was introduced for the case d = 2 and is typically applied to the extraordinary vertices of the mesh . for arbitrary divisions d we observe that this correction influences a small neighborhood around each extraordinary vertex . fig1 a and 19b illustrate catmull - clark correction weights of the first two operations , along the vector . weights in the right image should be divided by 64 . more precisely , this correction never propagates further than two rings of faces around the extraordinary vertex as shown in fig1 a and 19b . in addition , the corrections are only noticeable in the first couple of subdivision operations . the first subdivision operation produces the most visible change , which from formula 11 is equal to c i = n i - 4 n i ⁢ ( v i ′ - v i ) . subsequent subdivisions produce changes , w i c i , which are proportional to the first one by a weight w i . it is possible to compute these weights exactly for the first couple of subdivision operations the catmull - clark result can be computed on a plane , with all points assigned a z value of 0 , except for one . then the rules described herein can be applied to the same plane , and the difference can be computed . this gives the contribution of that point . this operation can be repeated for each point in the plane . as all calculations are linear , the same weights apply to any objects . this is what was done in fig1 a and 19b . these sampled weights then define a piecewise bilinear function on the unit square that can be used to compute the corresponding weight values when d is a power of two . for more general d values the weights can be interpolated from this function . in practice , however , we found that a similar behavior can be achieved using the push - back operation described in the next section . the effect of the catmull - clark correction can be emulated by using a higher α value and by adjusting the β parameter . this is apparent in fig1 b , where a value the push - back is similar to the d = 2 case described above : we first compute the δ values for the original vertices and then update the newly introduced vertices using bilinear interpolation . in a similar fashion we can use the normalized interpolation of the δ values to keep the lengths equal . for an even d , the push - back operation should be applied after equations ( 23 ) and ( 24 ) have been applied . this is because it doesn &# 39 ; t make sense to apply the push - back to the “ intermediate ” vertices f i , which are used temporarily to compute the new vertex positions . to make the algorithm consistent for every number of divisions d , we prefer the algorithm to perform the push - back for odd d when equation ( 22 ) is applied twice . after some experimentation with higher order interpolation schemes we concluded that the differences were likely too small to prefer a more expensive interplant . the present invention provides users with a simple smoothing tool for polygonal meshes . the smoothing operation allows users to create refined versions of their models . crucial to the success of such a model is that the transitions between the different resolutions of the meshes are almost imperceptible . in practice , we found that the new subdivision scheme of the present invention works best when we used a push - back operation with α = 1 2 , β = 1 2 ⁢ ⁢ and ⁢ ⁢ γ = 0 . of course , these parameters should typically be set by the artist , who can freely explore the effect of varying the parameters to meet particular needs . although this might be tricky , it is a huge improvement over current practice , where artists sometimes have to adjust individual vertices at each level of refinement . with the present invention , on the other hand , artists only have to worry about setting a few parameters at each level . in another embodiment , the present methods described herein can also be applied to not just quadrilateral surfaces but also polygon surfaces with five sides or more . fig2 a illustrates the subdividing process of an n - sided polygon . this method works with n ≧ 4 . first , the process subdivides 2000 edges of the polygon into d pieces . this creates new vertices along each edge . this can be accomplished by using any method , including ones described herein . after the subdivide operation 2000 , then for each edge , the process joins the first vertices on the adjacent edges , creating new edges . after the join operation 2002 , the process then connects 2004 the remaining vertices ( not the first vertices ) to the new edges . fig2 b illustrates an example of the subdividing process of an n - sided polygon using the process illustrated in fig2 a . original polygon 2006 is a four sided polygon to which we will apply the subdividing process of fig2 a . after performing the subdividing operation 2000 ( from fig2 a ), the result is subdivided edge polygon 2008 . note here we have divided each ( or edge ) into three pieces , as new vertices are created . next the joining operation 2002 ( from fig2 a ) is performed , which results in the joined edge polygon 2010 . note for simplicity , only one edge is joined in fig2 b , although typically all edges of the polygon will be joined . new edge 2011 is a newly created edge . next , the connecting operation 2004 ( from fig2 a ) is performed , which results in the connected edge polygon 2012 . note that for simplicity , the connecting is only applied to one edge in fig2 b , although typically the connecting operation 2004 applies to all edges . note that newly created points a 2014 and b 2016 are created , which are then connected to the subdivided points 2018 2020 on the outside edge . newly created points a 2014 and b 2016 can be calculated in any manner . one way would be to maintain a one to one correspondence with the original points 2018 2020 as shown here . if the new edge 2011 is not the same length as the original edge , then points a and b may be computed to maintain the same ratio along the new edge 2011 as original points 2018 2020 maintain along the original edge . numerous other connecting methods exist which could determine points a and b , and any of these methods can be applied , after the above - described operations are applied to all sides ( edges ) of a polygon , this leaves in the center an n - sided polygon ( same n as original polygon ) that should now be subdivided into d - 2 pieces . operations 2000 , 2002 , and 2004 are repeated until d = 1 or d = 0 . in the case of the quadrilateral shown in fig2 b , since d = 3 pieces , and d - 2 = 1 , the method is considered complete and no repeating is necessary . fig2 a illustrates the subdividing process of an n - sided polygon as applied to a square , with d = 6 . an outside part 2100 , a middle part 2102 , and a center part 2104 represent each different iteration of the process being applied . fig2 b illustrates the subdividing process of an n - sides polygon as applied to a pentagon , with d = 4 . an outside part 2106 and a center part 2108 represent each different iteration of the process being applied . fig2 c illustrates the subdividing process of an n - sides polygon as applied to a pentagon , with d = 5 . an outside part 2110 and a center part 2112 represent each different iteration of the process being applied . fig2 d illustrates the subdividing process of an n - sides polygon as applied to a hexagon , with d = 6 . an outside part 2114 , a middle part 2116 , and a center part 2118 represent each different iteration of the process being applied . in a further embodiment , the present methods can by applied to a collection of polytopes of arbitrary dimension . a polytope is a known math term and is a generalization of a polygonal mesh to arbitrary dimensions . for example in three dimensions , a polytope is a collection of closed polyhedra , where a polyhedron is a solid whose boundary is a closed polygonal mesh . in general , a k - dimensional polytope is defined recursively in terms of a collection of simple polytopes whose boundaries are ( k - 1 )- dimensional closed polytopes . in this embodiment , every polytope is first linearly subdivided into smaller polytopes and then smoothed as follows . each vertex of the polytope is replaced by the centroid of the polytope defined by the centroids of the polytopes adjacent to the vertex . alternatively , every polytope can be subdivided and then the vertices are smoothed by replacing them by the average of all centroids of the neighboring polytopes . the positions of the vertices of the polytopes can also be adjusted as described herein . each of these individual operations applied to polytopes may be accomplished by conventional methods and / or using the methods described herein which can simply be applied to k - dimensional polytopes . fig2 a illustrates a standard linear subdivision operation for cubes . original cube 2200 is subdivided into a subdivided cube 2202 , which comprises 8 smaller cubes . fig2 b illustrates a standard linear subdivision operation for tetrahedra . original tetrahedron 2204 is subdivided into a subdivided tetrahedron 2206 , which comprises 8 smaller tetrahedra . fig2 illustrates an example of linear subdivision of a 3d polytope . original polytope 2300 is subdivided into a subdivided polytope 2302 . fig2 is a block diagram illustrating one example of a configuration of hardware used to implement the present invention . a display monitor 2400 is connected to a computer 2402 . the computer performs the operational processes described herein based upon input from a keyboard 2406 and / or a mouse 2408 . a drawing tablet 2404 can also be connected to the computer 2402 . the computer 2402 has connected a rom 2410 , a ram 2412 , and a disk drive 2414 . in addition , a drawing pen 2416 and / or a puck 2418 can also be used as input devices for the tablet . of course , any applicable configuration of hardware can be used to implement the present invention . the system can also include any type of conventional peripherals , including permanent or removable storage , such as magnetic and optical discs , etc . further , any storage used with the computer ( disk drive , ram , etc .) can store the process and data structures of the present invention . the processes can also be distributed via , for example , downloading over a network such as the internet . the present invention has been described with respect to a general set of subdivision rules , which provide users with more control over the subdivision process . this application incorporates by reference the article entitled , “ a unified subdivision scheme for polygonal modeling ,” by jerome maillot and jos stam , eurographics 2001 , volume 20 ( 2001 ), number 3 , published 9 / 2001 . this application also incorporates by reference the article entitled , “ on subdivision schemes generalizing uniform b - spline surfaces of arbitrary degree ,” by jos stam , computer aided geometric design 18 ( 5 ), published june 2001 . the many features and advantages of the invention are apparent from the detailed specification and , thus , it is intended by the appended claims to cover all such features and advantages of the invention that fall within the true spirit and scope of the invention . further , since numerous modifications and changes will readily occur to those skilled in the art , it is not desired to limit the invention to the exact construction and operation illustrated and described , and accordingly all suitable modifications and equivalents may be resorted to , falling within the