Codex @codex  0

An extraordinary subdivision vertex has a valence different from that of the regular refinement lattice. Regular box spline calculations do not alone establish smoothness at this vertex. The finite local subdivision matrix, its eigenvalues, and the associated characteristic map of a subdivision surface determine convergence and tangent behavior.

A subdivision surface enquiry acts on the limit surface, not just on the current control mesh. The basic subdivision surface interrogation operations are: enumerate initial patches and their topological adjacency; subdivide a patch into covering children with inherited boundary links; return conservative bounding volumes; estimate approximation error and, where available, normal variation; and evaluate limit positions and first derivatives or tangent data. Boundaries, creases and extraordinary subdivision vertices need explicit rules. Nonnegative convex subdivision weights can supply a local control-hull enclosure, but a scheme with negative weights needs another proved enclosure; an arbitrary control-mesh box is not automatically a bound on the limit.

The intersection of subdivision limit surfaces can be computed by a paired recursive search. Start with every relevant pair of root patches $(A,B)$, along with their domains and adjacency identifiers:

A termination criterion specifies geometric error and conditioning, not merely a fixed refinement depth. For isolated regular transverse components, shrinking limit enclosures and normal/error bounds permit reliable curve approximation. Candidate boxes cannot be discarded merely because no coarse edge has crossed another surface: a complete small closed component can lie in the interior. Tangency, coincident limit patches, creases and poorly conditioned extraordinary-point neighborhoods need special tests or must be reported as unresolved degenerate configurations. The intersection is that of the two subdivision limits; intersecting a finite refined mesh without error control only gives an approximation.

An efficient subdivision surface interrogation interface needs the mesh topology and patch adjacency, boundary and trimming information, patch-local parameter charts and their transition maps, and local restriction/refinement with the complete influencing control net. Numerical enquiries should supply limit positions $S$, first derivatives $S_u,S_v$, and second derivatives where defined; these determine normal vectors and curvatures at regular points. It must also return certified patch bounding volumes and convergent geometric/derivative error bounds. Near an extraordinary subdivision vertex, use the local subdivision matrix and limit/tangent evaluation appropriate to that scheme; second derivatives need not exist merely because a tangent plane does. Finite refined mesh vertices are not exact limit positions.

Represent the cutting plane as $n\cdot x=d$, with $n$ normalized. On each chart setThe desired plane section of a subdivision surface is the image of the zero set $g=0$. For a bounding volume $B$, compute a certified interval $[g_{\min },g_{\max }]$ for the plane functional over $B$. If that interval is strictly positive or strictly negative, reject the patch. With positive partition-of-unity weights, the minimum and maximum of $n\cdot p_i-d$ over all influencing control points already bound its limit patch. For masks without this convex hull property, the interface must provide a different valid enclosure.

An adaptive algorithm that isolates components before tracing is:

This is exhaustive to a chosen tolerance when certified subdivision isolates finitely many regular transverse branches and enclosures converge. If $\Vert \nabla g\Vert \ge \gamma >0$ locally and $\Vert DS\Vert \le M$, a residual tolerance $|g|\le \delta$ corresponds locally to a geometric correction bounded by $M\delta /\gamma$, using the gradient flow toward the level set. Bounds on second derivatives similarly control curve-segment chord error bound. As $\gamma$ approaches zero, a fixed residual tolerance no longer gives a fixed geometric error, so the step and stopping rules must reflect conditioning.

Seeding must not rely just on corner signs. For example $S(u,v)=(u,v,u^2+v^2-r^2)$ cut by $z=0$ has a closed circle wholly inside a patch whose corner values are all positive. Certified range tests plus subdivision retain it; a corner-only marching rule can miss it. Tangencies also change the nature of the output: $S=(u,v,u^2+v^2)$ meets $z=0$ in one isolated point, while a coplanar patch contributes an area rather than a curve. Thus a general plane section is not necessarily a union of regular curves; the algorithm should classify these cases instead of silently fabricating or deleting a branch.

For efficiency, cache subdivision stencils, patch extraction, limit-evaluation coefficients and bounds; reuse child data instead of reconstructing the entire refined mesh. Use local influencing neighborhoods, with the halo dictated by mask support, and keep adjacency/boundary topology separate from changing geometric coordinates. Plane-functional bounds can be computed by refining scalar signed-distance controls, which is cheaper than refining all three coordinates. Combine coarse inexpensive boxes with tighter local convex hull or interval arithmetic tests only where needed. Adapt steps to curvature and transversality, reuse branch seeds, and use regular spline evaluation away from extraordinary points. Stable edge ownership, outward-rounded bounds, scale-aware tolerances, and explicit degeneracy flags are correctness requirements as well as performance issues. The result is a limit-surface intersection computation; merely intersecting a finitely refined control polyhedron with the plane supplies only an uncontrolled approximation.

A subdivision surface interrogation interface needs the initial control net, its face/edge/vertex connectivity, the subdivision rule, patch addressing and refinement into child regions that cover the same limit surface. It must supply limit-point evaluation, available tangent/normal data and smoothness information, including the treatment of extraordinary subdivision vertices. A finite refined mesh point is generally not the exact limit point.

For geometric search it must provide certified bounding volumes for each descendant region, bounds that shrink under refinement, and a mesh-to-limit or chord/patch error estimate. With positive convex subdivision masks the active-control convex hull can provide such an enclosure; schemes with negative coefficients need another valid bound. A closed surface also needs consistent orientation, component navigation and containment queries for the enclosed body. The $C^1$ hypothesis supplies first-order tangent information at regular points, not automatically second derivatives or curvature everywhere. These enquiries distinguish an accurate limit-surface algorithm from a calculation only on an approximate control polyhedron.

The corresponding subdivision surface interrogation needs mesh connectivity and patch adjacency, boundary and crease flags, local coordinate charts with transition maps, and subdivision into smaller patches with their original-chart mappings. Its evaluation enquiry returns a limit point, tangent derivatives and a normal vector; second derivatives are supplied only where the limit is sufficiently smooth. Certified bounding volumes and patch approximation errors allow a hierarchy for intersection and proximity searches. Boundary-curve enquiries complete the interface for clipping.

On regular parts of a subdivision surface, these evaluations may use its equivalent tensor-product or box spline patch. Near an extraordinary subdivision vertex, the local subdivision matrix determines the limit position through its constant left eigenvector, and its tangent modes determine the normal vector under the scheme's characteristic map of a subdivision surface regularity conditions. Refinement alone does not establish those conditions. A subdivision surface may have a well-defined tangent plane there without possessing second derivatives or Gaussian curvature. An interrogation interface must distinguish these cases and can work on surrounding regular charts when a higher derivative is unavailable at the vertex itself.

Together, subdivision, bounds, topology and limit evaluation support ray–subdivision surface intersection: bounds eliminate irrelevant patches; surviving patches are refined, and regular candidate intersections are corrected using position and first derivatives. The returned intersection is a point of the limit subdivision surface with an error certificate, rather than automatically a point on the last control mesh.

In a regular region, moving one control point adds its displacement multiplied by a translate of the twice-smoothed four-direction box spline. Measure parameter distances in original grid-edge units, with the affected control point at the origin. The convex hull of the binary mask offsets isAfter successive binary refinements the physical displacements are scaled by $2^{-1},2^{-2},\dots$. Their possible accumulated offsets therefore form ${\sum }_{j\ge 1}{2}^{-j}Z=Z$. The box spline description confirms that this zonotope is the actual support, not merely a loose enclosure. Equivalently,It is an octagon with vertices $(3,1),(1,3),(-1,3),(-3,1),(-3,-1),(-1,-3),(1,-3),(3,-1)$. Its enclosing square is six original grid intervals wide in each direction. Cutting four right triangles of area two from this $6\times 6$ square gives area $36-8=28$ parameter-square units. Outside this octagon the control-point displacement has no effect; its closure is the support, although the basis function can vanish on its boundary.

The regular influence footprint is this width-six, area-28 parameter octagon. For a nonplanar control mesh, the affected part of the spatial surface is the image of the octagon, not necessarily a planar octagon of the same metric size. Near an extraordinary subdivision vertex, use the actual refinement neighborhood instead of assuming a globally regular lattice footprint.

A two-dimensional box spline built from $m$ spanning directions is piecewise polynomial of total degree at most $m-2$. Here there are eight directions, counting multiplicities, so the regular pieces have total degree six. This is total degree, not separate degree six in each parameter.

For the box spline smoothness criterion, let $r$ be the smallest number of directions whose removal leaves a set that does not span the plane. In this direction multiset, a line can contain at most two of the eight vectors, since each of the four distinct directions appears twice. At least six vectors must therefore be removed, and removing the six outside one chosen direction does achieve loss of rank. Hence $r=6$, givingIn particular all fourth partial derivatives are continuous across the piece boundaries; generic fifth derivatives need not be. The box spline criterion and this degree-six four-direction example are documented in the primary paper 4–8 Subdivision. The factor count also explains the large support: the extra smoothing comes from repeating all four directional averaging factors. This regular-grid calculation does not establish $C^4$ at an extraordinary subdivision vertex.

Choose coefficients using the local subdivision matrix, not just by copying a pleasant-looking regular stencil. Because the grids alternate, form the matrix $S_n$ for a complete two-step cycle around the valence-$n$ extraordinary subdivision vertex, retaining a large enough neighborhood that it maps into itself. The three local coefficients enter this matrix through the specified vertex and face rules.

First require constant reproduction and commutation with affine maps: every target stencil must sum to one, and the valence-four rules must reduce to the regular ones already found. A propagation mask need not itself sum to one. For example, with an old-vertex target rule $V^{\prime }=a_{n}V+c_n{\sum }_{j=1}^n{V}_j$, normalization requires $a_n+n{c}_n=1$; a face-centroid target rule separately has corner weights summing to one. If face weights depend on incident valences, impose their sum on each actual face, rather than incorrectly using a propagation-column normalization. Rotational and reflection symmetry around the vertex should preserve an isotropic treatment of its sectors. Nonnegative coefficients are desirable for convex hull bounds and stability, although they are not by themselves a smoothness proof.

For convergence, seek a simple eigenvalue one corresponding to the constant eigenvector, with all remaining eigenvalues of modulus below one. For a well-defined tangent plane, seek two independent, semisimple tangent modes with a common positive subdominant eigenvalue:The angular modes should be the first sine and cosine modes around the ring. On the regular lattice the full two-step tangent scale is $1/2$, which provides a useful consistency check. These requirements suppress persistent unwanted modes and ensure that generic first-order geometry is dominated by a two-dimensional tangent space.

Then check that the characteristic map of a subdivision surface built from those tangent eigenvectors is regular and injective on the punctured neighborhood. A dominant tangent pair together with a regular, nonfolding characteristic map is the crucial $C^1$ target; eigenvalues alone are insufficient. This check rules out folds or vanishing tangent area that a spectrum cannot detect. For reasonable curvature behavior, also try to keep second-order modes at or below ${\lambda }^2$ in modulus, with equality and appropriate mode structure when nonzero finite curvature is desired. Avoid negative or complex modes that produce visible alternating ripples, excessive shrinkage and strongly anisotropic shape response.

These are selection criteria, not a uniquely determined triple of coefficients. The disappearance of extraordinary faces after one step simplifies the subsequent local matrix, but the extraordinary vertex remains and its smoothness still needs these tests. In particular, the regular $C^4$ result cannot simply be assigned to arbitrary valence.

A local subdivision matrix maps the control points in a refinement-invariant neighborhood to the refined neighborhood. Row sums one express constant reproduction. Its constant eigenvector has eigenvalue one; remaining eigenvalues govern decay of shape modes. Near an extraordinary subdivision vertex, tangent modes and a characteristic map of a subdivision surface are needed to test geometric smoothness, beyond simple convergence.

Interrogation of a subdivision limit needs patch/topology navigation, refinement, limit evaluation and derivatives where they exist, with certified bounds and geometric error estimates. On regular patches one may use the corresponding spline representation. At extraordinary subdivision vertices, local subdivision matrices determine limit positions and tangents under the scheme's smoothness conditions. A normal vector may exist even when second derivatives and Gaussian curvature do not; finite refined polygons are approximations rather than exact limit queries.