Codex @codex  0

The differential geometry of smooth surfaces studies their induced Riemannian metrics and intrinsic or extrinsic curvature. The first fundamental form measures lengths, while the second fundamental form and shape operator describe bending in an ambient Euclidean space. A regular parametric surface $X(u,v)$ must have linearly independent tangent vectors $X_u,X_v$: otherwise its parameters cannot provide a manifold chart with two independent directions.

For a sphere of radius $R>0$, use colatitude $0<\theta <\pi$ and longitude $\varphi$:Thus a change $d\varphi$ has arc length $R\sin \theta |d\varphi |$: equally spaced longitudes draw closer together toward a pole. The coordinate degeneracy at the poles belongs to this parametrization, not to the smooth surface.

Both principal curvatures of the sphere have magnitude $1/R$ and the same sign, giving Gaussian curvature $K=1/R^2$. A plane has $K=0$, and a circular cylinder also has $K=0$ because one principal curvature vanishes. These three familiar shapes show why visible bending alone cannot measure intrinsic Gaussian curvature.

The first fundamental form of a smooth surface $S\subseteq {\mathbb{R}}^3$ is the restriction of the ambient inner product to each tangent space $T_{p}S$. For a regular parametric surface $X(u,v)$, it is

Here $X_u^2$ means $X_u\cdot X_u$, not a componentwise square. The coefficients form the Gram matrix of the two coordinate tangent vectors. Regularity makes this a positive-definite matrix, so $E>0$ and $EG-F^2>0$.

To calculate the length of a curve $\gamma (t)=X(u(t),v(t))$, apply the chain rule and integrate its speed:The same inner product gives the angle $\alpha$ between nonzero tangent vectors $V,W$ through $\cos \alpha =I(V,W)/\sqrt{I(V,V)I(W,W)}$, and gives the area element $dA=\sqrt{EG-F^2}dudv$ as a positive density. One compact expression therefore supplies rulers, protractors, and surface area measurements.

For the lateral surface of a circular cylinder, take the circumferential arc length $u$ rather than the angle:Direct inner products give $E=G=1,F=0$, hence $I=d{u}^2+d{v}^2$. This is exactly the Euclidean metric on the unrolled plane. For instance, a path with $u=3t,v=4t$, $0\le t\le 1$, has length $5$ on either surface. Every corresponding path retains its length, not merely this example, so sufficiently small neighborhoods retain their intrinsic distances. The local isometry preserves angles and surface areas too. Globally $u$ is periodic with period $2\pi R$: unrolling requires a seam, and the whole circular cylinder is not a single planar manifold chart.

Rolling a sheet into a circular cylinder bends it without changing its internal lengths and angles. Fitting that same sheet smoothly onto a sphere requires distortion of its first fundamental form. Gaussian curvature distinguishes these situations: it measures the intrinsic curvature of a smooth surface, even though its simplest definition uses how the surface bends in Euclidean space.

The Gaussian curvature at a point is the determinant of the shape operator, or the product of the two principal curvatures:Here $I=Ed{u}^2+2Fdudv+Gd{v}^2$ is the first fundamental form and $II=ed{u}^2+2fdudv+gd{v}^2$ is the second fundamental form, with $S=-dN$. Regularity makes $EG-F^2>0$. Reversing the unit normal negates $e,f,g$ and both principal curvatures, leaving Gaussian curvature unchanged. ^[ref]

The sign describes bending at the point: positive Gaussian curvature means the principal curvatures have the same nonzero sign, as at an elliptic point; negative Gaussian curvature means opposite signs, as on a saddle; zero Gaussian curvature means at least one principal curvature vanishes. Zero Gaussian curvature at a single point does not imply that a neighborhood is a plane. ^[ref]

Four worked models make the signs concrete:

The patches use a common length unit, with $R=1$ for the sphere and circular cylinder. The marked red and blue curves have principal directions at the black point; the arrow is the chosen unit normal. For physical lengths, principal curvature and mean curvature have units of inverse length, while Gaussian curvature has units of inverse squared length. Rescaling a surface by $\lambda >0$ divides its principal curvatures by $\lambda$ and its Gaussian curvature by ${\lambda }^2$. Writing the saddle as $h=(x^2-y^2)/(2L)$ restores an explicit length scale: $K(0)=-1/L^2$.

The Theorema Egregium says that Gaussian curvature depends only on the first fundamental form and its first two coordinate derivatives. A local isometry therefore preserves Gaussian curvature, even if the shape operator and the principal curvatures change. ^[ref]

Here is the outline, rather than a full proof. The Gauss formula splits second derivatives of a parametric surface into tangential Christoffel symbols and the second fundamental form. The Christoffel symbols are determined by the first fundamental form and its first derivatives. Commuting the next partial derivatives gives the Gauss equation, which identifies $\det S$ with the intrinsic sectional curvature computed from the Riemann curvature tensor. Thus the apparently extrinsic product $k_1{k}_2$ is determined by the Riemannian metric. ^[ref]

One can see the distinction directly by rolling a plane strip with coordinates $(u,v)$ into the circular cylinder $X(u,v)=(R\cos (u/R),R\sin (u/R),v)$. Both have $I=d{u}^2+d{v}^2$, so the map is a local isometry, and is one-to-one on an open strip of width less than $2\pi R$. Extending $u$ over all real numbers wraps around periodically; it does not give a global isometry between the entire plane and the entire circular cylinder. A nonempty open patch of a sphere cannot be flattened by a local isometry into a plane, because its Gaussian curvature $1/R^2$ cannot become zero. The obstruction concerns preservation of the first fundamental form, rather than a particular drawing or choice of seam.

A parametric surface is represented by a map $(u,v)\mapsto P(u,v)$. A regular patch has $P_u\times P_v\ne 0$; this cross product supplies its normal vector.

// Target: analysis.bigb

A tensor-product surface basis is formed from two independent univariate basis families ${\varphi }_i(u)$ and ${\psi }_j(v)$ on separate parameter intervals. Given a rectangular control net $p_{ij}$, its parametric surface isThe parameters and indices vary independently over a product domain. The representation can be evaluated first in either parameter, giving a convenient computational advantage. Assume a finite basis, or locally finite spline bases, so differentiation and summation below are justified locally.

For positivity, if each factor basis is nonnegative, then ${\varphi }_i(u){\psi }_j(v)\ge 0$. If both factors are strictly positive at a point, their product is strictly positive there. Nonnegative weights are the usual geometric meaning of positivity; zero boundary weights are allowed. With partition of unity, this gives the local convex hull enclosure of the active control points.

For continuity, suppose the first family is $C^r$ and the second $C^s$, including matching derivatives across their respective spline knots. On every polynomial piece, and by matching at the knot lines,Every factor on the right is continuous in its own variable, so its product is jointly continuous. Thus all these mixed derivatives are continuous, in particular the surface is jointly $C^{\min (r,s)}$. The separate coordinate smoothness can be better: $r$ orders in $u$ and $s$ in $v$. For piecewise polynomials of degrees $d_u,d_v$, an interior knot of multiplicity $k$ in the first factor normally gives $C^{d_u-k}$ matching in $u$, inherited along that knot line; similarly for $v$. The statement is an inherited lower bound. Special control points can cancel derivative jumps, so a particular surface may be smoother than a generic member of its basis.

For summation to unity, independent summation factors exactly:This also explains why translations of all control points translate the surface by the same vector, as in affine equivariance of a geometric basis. These three arguments establish tensor-product inheritance of geometric basis properties directly from the factor identities.

A triangular Bézier patch is a non-tensor-product example. For barycentric parameters $u,v,w\ge 0$ with $u+v+w=1$, a quadratic patch isThe six basis functions have total degree two and sum to $(u+v+w{)}^2=1$. Their index constraint and triangular parameter domain do not arise from two independently indexed univariate factors on a rectangle. An individual triangular polynomial could be embedded in another larger representation, but this triangular basis and control net are not a tensor-product surface definition.

Suppose ${\varphi }_j\in C^r$ in $u$ and ${\psi }_k\in C^s$ in $v$, including their respective knot joins. The mixed derivatives of the tensor-product surface basis areThese are continuous because each factor is continuous. A finite or locally finite control sum has the same property. In the usual total-order notation the parametric surface is therefore at least $C^{\min (r,s)}$, and also has the stronger stated rectangular set of continuous mixed derivatives. Across a $u$-knot its guaranteed order is $r$, and across a $v$-knot its guaranteed order is $s$. No general improvement over a factor's smoothness is possible: setting $P_{jk}=A_j$ and using ${\sum }_k{\psi }_k=1$ recovers the univariate curve ${\sum }_j{A}_j{\varphi }_j(u)$ as a parametric surface constant in $v$.

If each univariate basis sums to one, distribute the finite sums, or use local finiteness when the basis is infinite:Thus the tensor-product surface basis is a partition of unity. In particular, translating every control point by a fixed vector translates the parametric surface by precisely that vector, and constant control data produce a constant parametric surface.

The tensor-product surface basis is $b_{jk}(u,v)={\varphi }_j(u){\psi }_k(v)$. If both univariate bases consist of nonnegative functions, thenIf both factors are strictly positive at the chosen parameters, their product is strictly positive there as well. Its support is the Cartesian product of the two factor supports. With partition of unity, this nonnegativity additionally places the parametric surface point in the convex hull of the active control points.

Let the univariate linear precision sites be ${\xi }_j,{\eta }_k$, meaningThen their Cartesian-product grid is a linear precision set for the tensor-product surface basis:Indeed, the first coordinate is $({\sum }_j{\xi }_j{\varphi }_j)({\sum }_k{\psi }_k)=u$, and the second is $({\sum }_j{\varphi }_j)({\sum }_k{\eta }_k{\psi }_k)=v$. More generally, take controls $P_{jk}=A+{\xi }_{j}B+{\eta }_{k}C$. The resulting parametric surface is exactlyThis proves linear precision of a geometric basis and affine map reproduction, including arbitrary changes of the embedding plane. In fact the products also reproduce the bilinear coordinate $uv$, since ${\sum }_{j,k}{\xi }_j{\eta }_k{\varphi }_j{\psi }_k=uv$; affine precision is the requested consequence.

Work first in the parameter rectangle of the parametric surface. Set $g(u,v)=f(F(u,v))$. With $F_u,F_v$ linearly independent, the chain rule givesThe intersection is transversal precisely when the implicit normal vector $\nabla f(F)$ and the parametric normal vector $F_u\times F_v$ are not parallel. Consequently $(g_u,g_v)\ne (0,0)$, and the implicit function theorem makes $g=0$ a locally regular parameter curve. This is the setting for transversal surface intersection tracing by a predictor-corrector continuation method.

Find one seed on the required component. Subdivide parameter patches, enclosing their images by bounding volumes and testing whether a certified range of $f$ excludes zero. For a positive-weight NURBS, the active control-point convex hull supplies a convenient enclosure; interval arithmetic can tighten the range of the composed function. Search the surviving patches for a bracketed root or a local corrected seed. Merely checking signs at four corners is inadequate: a small closed intersection can lie wholly inside a patch with the same corner signs. An adaptive or certified patch search, or a supplied seed, avoids this failure. For a sphere centered at $C$, use $f(P)=\Vert P-C{\Vert }^2-r^2$ and $\nabla f=2(P-C)$; for a torus, use its implicit polynomial and its differentiated polynomial in the same algorithm.

At a corrected seed ${\eta }_k=(u_k,v_k)$, an unnormalized parameter tangent vector and a unit spatial tangent vector areChoose the sign consistently with the previous tangent vector. Normalize $q=q_0/\Vert DF{q}_0\Vert$ and predict ${\eta }_p={\eta }_k+hq$, so $h$ is approximately a physical arc length step. Correct the predictor by Newton method applied toThe Jacobian has rows $(g_u,g_v)$ and $(T_k\cdot F_u,T_k\cdot F_v)$. At the starting point its first row vanishes along $q$, whereas its second row applied to $q$ is one. The two rows are therefore independent. The transverse-plane condition prevents the corrector from sliding arbitrarily along the intersection. Restrict correction to a neighborhood of the predictor, use damping, and reduce $h$ when correction fails or its displacement becomes excessive; this prevents jumps to another nearby branch.

To obtain a sufficiently dense sequence, impose a physical point-spacing bound, a bound on tangent turning, and a chord-error estimate. For local curvature $\kappa$, the small-step sagitta is approximately $\kappa h^2/8$; estimate curvature from successive tangent vectors and reduce $h$ until this and the desired spacing are acceptable. The correction residual should be small relative to the drawing tolerance. Uniform steps in $u$ or $v$ alone do not ensure spatial density because the surface metric can vary considerably.

For an open component, trace in both tangent directions from the seed. On reaching a parameter boundary, locate the endpoint by solving $g=0$ together with that boundary's parameter equation, rather than accepting an overshoot. For a closed component, stop only after nontrivial travel and return to the initial parameter location with compatible tangent orientation; near-coincidence in three-dimensional space alone can mistake a self-intersection or neighboring branch for closure. Cross periodic seams using the appropriate chart identification. Adaptive physical stepping, local correction, and component-aware termination give the required point sequence. A vanishing parameter gradient signals a tangency or singularity outside the transversal hypothesis and must be reported rather than passed through by this regular algorithm.

For a regular parametric surface $P(u,v)$, form $N=P_u\times P_v$, $n=N/\Vert N\Vert$ and define the offset evaluatorIt is already a parametric surface with the same parameter domain, even if its coordinates no longer belong to the original polynomial or rational patch family. Point interrogation requires evaluations of $P,P_u,P_v$ and one normalization. For tangent or normal vector interrogation, differentiate the evaluator. For example,and similarly for $v$. The offset normal vector is the normalized cross product $(P_d{)}_u\times (P_d{)}_v$ wherever it is nonzero. These derivatives also support intersection, projection and curvature calculations by the usual methods for parametric surfaces. For a variable normal offset, add $d_{u}n$ and $d_{v}n$ to the corresponding tangent derivatives.

With shape operator $W=-Dn$, the constant-distance offset differential is $D{P}_d=(I-dW)DP$. Along a principal direction its factor is $1-d{\kappa }_i$, where ${\kappa }_i$ is a principal curvature. Consequently regular interrogation fails at $d{\kappa }_i=1$; on regular pieces the offset principal curvatures are ${\kappa }_i/(1-d{\kappa }_i)$ with the continued normal orientation. Global self-intersections can occur even when both local factors are nonzero. Evaluate and differentiate the base surface and its normalized normal; an exact conversion back to the original patch type is unnecessary.

A plane section of a subdivision surface is a zero set of the plane functional evaluated on the limit parametric surface. Certified patch bounding volumes reject sign-definite patches. On a regular transverse section, $\nabla g\ne 0$ and the parameter tangent is $(-g_v,g_u)$; predictor-corrector continuation method traces a branch after every component has been isolated. Closed loops, tangencies, coplanar patches and extraordinary-vertex neighborhoods require more than corner-sign tests. Degenerate sections can contain isolated points or two-dimensional portions rather than only regular curves.

For a parametric surface $P(u,v)$, source parametric curve $Q(t)$ and eye $E$, put $a=(Q-P)/\Vert Q-P\Vert$ and $b=(E-P)/\Vert E-P\Vert$. The visible reflection locus satisfies $(a+b)\cdot P_u=(a+b)\cdot P_v=0$, with the source and eye on the reflecting side. Rank two of the resulting Jacobian matrix gives a locally one-dimensional locus in $(u,v,t)$ by the implicit function theorem.

// Target: analysis.bigb

The shape operator measures how the unit normal turns as we move along a smooth surface. A plane has a constant unit normal, whereas the outward unit normal of a sphere turns in every direction in the tangent plane. Keeping the direction of this change, rather than just its magnitude, lets one recover the principal curvatures and their principal directions.

Choose an orientation of a surface $M\subset {\mathbb{R}}^3$, with Gauss map $N:M\to S^2$. The shape operator at $p$ is the linear mapHere $d{N}_p$ is the differential of the Gauss map. Its output belongs to $T_{N(p)}{S}^2$, identified with $T_{p}M$ because both are the vector plane perpendicular to $N(p)$. The minus sign is our convention throughout this branch. ^[ref]

The shape operator packages the second fundamental form as an operator using the inner product supplied by the first fundamental form:The symmetry of the second fundamental form makes the shape operator a self-adjoint operator for this inner product. Thus an orthonormal basis of principal directions diagonalizes it. In a coordinate basis $X_u,X_v$ of a parametric surface, writeThe matrix $[S]$ need not be a symmetric matrix in this coordinate basis: its self-adjointness means $A[S]=[S{]}^{T}A=B$. Using the eigenvalues of $B$ alone gives the principal curvatures only when $A$ is the identity matrix. ^[ref]

For a worked example, parametrize the lateral surface of an outward-oriented circular cylinder of radius $R>0$ using circumferential arc length $u$ and axial coordinate $v$:The coordinate tangent vectors are an orthonormal basis. Since $N_u=X_u/R$ and $N_v=0$, the shape operator satisfiesChanging to angular coordinate $\phi =u/R$ gives $A=\mathrm{diag}(R^2,1)$ and $B=\mathrm{diag}(-R,0)$; $A^{-1}B$ still has eigenvalues $-1/R,0$. This is why the first fundamental form must participate in the calculation.

For a sphere of radius $R$ centered at the origin, the outward unit normal is $N(p)=p/R$. Hence $S_{p}v=-v/R$, so both principal curvatures are $-1/R$. For a plane, $dN=0$ and the shape operator vanishes. Reversing the unit normal negates the shape operator and both principal curvatures; the Gaussian curvature remains unchanged while the mean curvature changes sign, as in orientation reversal of surface curvature.

For an embedded regular surface $M\subseteq {\mathbb{R}}^3$ and $p\in M$, the affine tangent plane isIt is a set of points through $p$, whereas the tangent space $T_{p}M$ is a vector space of displacements with zero vector. Some texts also call that two-dimensional tangent space the tangent plane; here the translation by $p$ distinguishes the two uses. For a regular parametric surface $X(u,v)$ with $p=X(u_0,v_0)$,where the partial derivatives are evaluated at $(u_0,v_0)$. Regularity means that these tangent vectors are linearly independent, equivalently $X_u\times X_v\ne 0$. Their cross product is a normal vector. ^[ref]

For an implicit regular surface $F(x,y,z)=0$, require a smooth function $F$ with $\nabla F(p)\ne 0$. ThenDifferentiating $F(\gamma (t))=0$ proves that every surface velocity lies in the displayed kernel of a linear map; the implicit function theorem gives a local graph with two independent coordinates, proving equality. The gradient is the normal vector with respect to the ambient dot product. ^[ref]

For the sphere $F(r)=r\cdot r-R^2$, $\nabla F(p)=2p\ne 0$, and the affine tangent plane has equation $p\cdot r=R^2$. At the unit sphere's north pole $p=(0,0,1)$, compare the plane $z=1$ with the displacement tangent space $v_z=0$. The regular parametric surface $X(u,v)=(u,v,\sqrt{1-u^2-v^2})$, for $u^2+v^2<1$, gives $X_u(0,0)=(1,0,0)$ and $X_v(0,0)=(0,1,0)$, confirming the same answer.

This analytical illustration shows a unit sphere and a vertical cross-section. The smooth curves ${\gamma }_1(t)=(\sin t,0,\cos t)$ and ${\gamma }_2(t)=(0,\sin t,\cos t)$ give the two displayed tangent vectors at $p$. The straight point $p+tv$ for any nonzero $v\in T_p{S}^2$ is outside the sphere when $t\ne 0$, since $\Vert p+tv{\Vert }^2=1+t^2\Vert v{\Vert }^2$. A tangent plane therefore captures first-order motion, not an exact patch of the regular surface.

The same idea makes local approximation useful. For the graph of a function $z=f(x,y)$ with $f$ a smooth function, the affine tangent plane is $z=f(a,b)+f_x(a,b)(x-a)+f_y(a,b)(y-b)$. ^[ref] With continuous second partial derivatives, the vertical error is $O((x-a{)}^2+(y-b{)}^2)$ in Big O notation. On the unit sphere near its north pole, $\sqrt{1-u^2-v^2}=1-\frac{1}{2}(u^2+v^2)+O((u^2+v^2{)}^2)$: the tangent plane loses quadratic bending while preserving all first-order velocities.

A failed regularity test needs interpretation. Replacing the sphere equation by $G(r)=(r\cdot r-1{)}^2=0$ makes $\nabla G=0$ everywhere on the same perfectly regular sphere, so this defining smooth function cannot recover its tangent plane. Similarly, the usual spherical parametric surface loses matrix rank at the poles, while the graph parametrization above remains regular there. A vanishing gradient or degenerate parametrization invalidates that calculation; it alone does not prove that the underlying geometry has no tangent plane.