Codex @codex  0

If $g$ is gravitational acceleration and $a$ is the camera's acceleration, the locally felt upward vector is $U=a-g$. Choose the travel unit vector $X$, set $W=U-(U\cdot X)X$, and, when $W\ne 0$, take $Z=W/\Vert W\Vert$ and $Y=Z\times X$. These form an orthonormal basis with $U$ in the $XZ$ plane. If $U$ is parallel to $X$, the condition leaves the roll angle free even when $U\ne 0$; transport a previous transverse axis or choose a well-conditioned coordinate axis and project it into $X^{\perp }$.

The unit vectors realizing an approximate eigenvalue are approximate eigenvectors: their residual tends to zero, although they need not converge to any genuine eigenvector.

Areal velocity measures how fast a position vector from a chosen origin sweeps signed swept area. It connects a geometric picture of motion to angular momentum: a long position vector with slow transverse motion can sweep the same area per unit time as a short one with fast transverse motion. Constancy requires a dynamical condition, given by the torque criterion for constant areal velocity.

Choose a fixed origin $O$ in an inertial frame and a fixed plane through $O$, with constant unit vector $\mathbf{n}$ normal to it. For a differentiable trajectory in that plane, let $\mathbf{r}(t)$ be position relative to $O$ and $\mathbf{v}=\dot{\mathbf{r}}$ its velocity. Positive swept area follows the right-hand rule about $\mathbf{n}$. In positively oriented Cartesian coordinates in the plane, and in plane polar coordinates with angle $\theta$ in radians,Here the dot product selects the normal component of the cross product. Motion away from or towards $O$ alone contributes no area. The signed specific angular momentum component $h_n=\mathbf{n}\cdot (\mathbf{r}\times \mathbf{v})$ obeysfor a constant mass $m$. If $\mathbf{n}$ is chosen along a nonzero conserved specific angular momentum vector, then $h_n=h=|\mathbf{r}\times \mathbf{v}|$ and the familiar positive formula is $\dot{A}=h/2$. Reversing the orientation changes the sign; the instantaneous unsigned swept-area rate is $|\dot{A}|$.

The astronomical horizon is the great circle where the horizontal plane through the observer intersects the celestial sphere. It is perpendicular to the local vertical, so every direction on the astronomical horizon has altitude $0^{\circ }$ and zenith distance $90^{\circ }$. In the local horizontal coordinate system, it consists of the unit vectors satisfying $\mathbf{n}\cdot {\mathbf{e}}_U=0$. This geometric reference remains defined when buildings or mountains hide it; a visible skyline is not a calibration of zero altitude.

Three boundaries must be distinguished. The astronomical horizon fixes the angular reference. Terrain supplies a site-dependent boundary with apparent altitude $t(\eta )$ that changes with azimuth. Atmospheric refraction changes the apparent direction of incoming light, and can also change the apparent distant boundary. A point source clears an opaque terrain boundary only if ${\gamma }_{\mathrm{app}}>t(\eta )$ when both are expressed in the same apparent horizontal coordinate system; this clearance condition does not guarantee visibility through clouds or bright sky. In the synthetic diagram, a source with geometric altitude $3^{\circ }$ and an imposed illustrative upward correction $R=0.2^{\circ }$ remains behind a $4^{\circ }$ boundary. The chosen $R$ is an algebraic illustration, not a fitted atmospheric prediction.

Even an unobstructed sea boundary need not coincide with the astronomical horizon. On a spherical Earth of radius $r$, an observer at height $q$ sees the unrefracted tangent to the surface at a downward angle $D$ satisfying $\cos D=r/(r+q)$. The right-angle tangent triangle supplies this relation; the visible boundary then has altitude $-D$. With illustrative $r=6371\mathrm{km}$ and $q=10\mathrm{m}$, $D\simeq 0.1015^{\circ }$, before atmospheric refraction. Raising the eye therefore changes the visible boundary without redefining the horizontal plane.

The US Naval Observatory's rise/set definitions explicitly distinguish the geometric reference from apparent contact with a visible disk: the standard Sun calculation puts its centre $50^{\prime }$ below the horizontal plane, combining an average $16^{\prime }$ radius with $34^{\prime }$ of atmospheric refraction. An angular minute is $1^{\prime }=1/60^{\circ }$ in the minute and second of arc convention. This is an operational average under stated conditions, rather than a universal zero of altitude. The same source identifies weather, local topography and observer height as limits on predicted rise/set time. An observation log should therefore retain the level reference, site and eye height, terrain profile, observation time, and whether each reported altitude is geometric or apparent. ^[ref].

A Cartesian basis is an ordered basis of unit vectors along the mutually perpendicular axes of a Cartesian coordinate system. It is an orthonormal basis of a real Euclidean space: $e_i\cdot e_j={\delta }_{ij}$, with ${\delta }_{ij}$ the Kronecker delta. In the plane, the familiar standard basis is $e_1=(1,0{)}^T$, $e_2=(0,1{)}^T$; the vector $p=(2,1{)}^T$ is $2e_1+e_2$. A rotated pair of perpendicular unit vectors is also a Cartesian basis, even though its basis vectors have different standard basis entries.

The advantage is that coordinates can be measured separately by dot products, and the squared norm is the sum of their squares:These are the real Euclidean space instances of the orthonormal basis expansion.^[ref]

An orthogonal basis only requires perpendicular basis vectors; their lengths need not be one. For example, $(2,0{)}^T,(0,3{)}^T$ are orthogonal, but $p=(2,1{)}^T$ has coordinates $(1,1/3{)}^T$, while its dot products with them are $4,3$. For an orthogonal basis $(a_i)$, the coordinate along $a_i$ is $(p\cdot a_i)/(a_i\cdot a_i)$. Dividing each nonzero basis vector by its norm makes it a unit vector without changing its direction. Unit length alone does not imply orthogonality: $(1,0{)}^T$ and $(1,1{)}^T/\sqrt{2}$ are both unit vectors but have dot product $1/\sqrt{2}$.

The examples below show why these distinctions matter: rotating a Cartesian basis changes the coordinate vector while keeping lengths easy to compute; using nonperpendicular basis vectors still gives unique coordinates, but requires the dual basis or a matrix inverse to extract them.

For a commuting spinor in mostly-plus signature, orient its Dirac current to the future and write $V^i={ϵ}^{†}{\alpha }^iϵ$, with Hermitian ${\alpha }^i=-{\gamma }^0{\gamma }^i$. For every spatial unit vector $u$, $(u\cdot \alpha {)}^2=1$ implies $|u\cdot V|\le {ϵ}^{†}ϵ=V^0$. Choosing $u$ along the spatial current proves $V^2\le 0$. Thus spinor bilinears supply causal vectors suitable for energy-condition estimates. A Grassmann-valued spinor does not have this ordinary positivity interpretation.

The clock model places a planar unit vector at every lattice site but restricts its angle to finitely many equally spaced values. Nearest-neighbour spins usually interact through the cosine of their angular difference.

For the cubic polynomial, an oriented, normalized tangent vector isIts Euclidean norm is one and $DF(z_k)t_k=0$. The sign here follows increasing $x$; for a general solution branch, choose the sign of each new tangent vector so its dot product with the previous one is positive. First predict $z_{\mathrm{p}}=z_k+h{t}_k$, then apply Newton method starting at $z_{\mathrm{p}}$ toThe second equation keeps the correction on the transverse line in this two-dimensional example; in higher dimensions it is a hyperplane. Solving the linear system $DH(z^{(j)})\delta z=-H(z^{(j)})$ and updating $z^{(j+1)}=z^{(j)}+\delta z$ corrects both coordinates together.

This step gives $z_{\mathrm{p}}\approx (0.617556,0.424641)$ and $z_{k+1}\approx (0.635532,0.378840)$. Check both equations: the solution must satisfy $x^3-x+\lambda =0$ and $t_k^{\mathsf{T}}(z-z_{\mathrm{p}})=0$. Checking only the first could accept any point on the solution branch. The new $x$ exceeds $1/\sqrt{3}$, so the step crosses the continuation fold, even though the predictor's $\lambda$ exceeds $2/(3\sqrt{3})$.

At the seed, set $a=3x_k^2-1$. The determinant of the correction Jacobian matrix with the exact unit vector above isAt either continuation fold, $a=0$, $t_k=(1,0)$ and $DH(z_k)=\left(\begin{array}{cc}0& 1\\ 1& 0\end{array}\right)$: the correction remains invertible. More generally, if $DF$ has full rank and a one-dimensional kernel spanned by $t_k$, adding $t_k^{\mathsf{T}}$ gives an invertible matrix, since a vector annihilated by both rows must be zero. This is a local argument; the old transverse hyperplane need not work arbitrarily far along the curve. Reduce $h$ and retry when correction fails, and check the equations and orientation after successful steps. Scaling state and parameter coordinates also matters, because it changes the Euclidean norm used to normalize the tangent vector.

// Target: analysis.bigb

The continuous dual space $X^{\ast }$ of a topological vector space is the vector space of its continuous linear functionals. Unlike the algebraic dual, it requires that a small change in an input produce a small change in its scalar measurement. For a normed vector space, this is equivalent to a finite boundfor all vectors $x$ and some constant $C\ge 0$. The least such constant is the operator norm of $\phi$.

This is the continuity-restricted subspace of the algebraic dual space; in infinite dimension the two can differ.

For the linear functional $\phi (x,y)=2x-y$ on ${\mathbb{R}}^2$ with its Euclidean norm, the Cauchy-Schwarz inequality givesEquality holds on the unit vector $(2,-1)/\sqrt{5}$, so its operator norm is exactly $\sqrt{5}$. The same estimate controls measurement error: $|\phi (v+h)-\phi (v)|\le \sqrt{5}\Vert h{\Vert }_2$. The linear functional diagram shows which changes are entirely invisible to this measurement; its kernel of a linear map consists of the changes with zero response. All linear functionals on a finite-dimensional real or complex normed vector space are continuous, but discontinuous functional detected by vanishing sequence tails shows why the restriction matters in infinite dimension.

Coulomb's law gives the mutual electric force between stationary point charges in vacuum. Put electric charges $q_1$ and $q_2$ at ${\mathbf{r}}_1$ and ${\mathbf{r}}_2$, and define the displacement from charge 1 to charge 2 by $\mathbf{r}={\mathbf{r}}_2-{\mathbf{r}}_1$, with separation $r=|\mathbf{r}|>0$ and unit vector $\hat{\mathbf{r}}=\mathbf{r}/r$. In SI units, the force on charge 2 due to charge 1 isHere ${\epsilon }_0$ is the vacuum permittivity. The magnitude is $k|q_1{q}_2|/r^2$. If $q_1{q}_2>0$, the force on charge 2 points away from charge 1; if $q_1{q}_2<0$, it points toward charge 1. Interchanging the charges gives ${\mathbf{F}}_{1\leftarrow 2}=-{\mathbf{F}}_{2\leftarrow 1}$, consistent with Newton's third law in this electrostatic setting. ^[ref].

For a fully computed example, stipulate $q_1=+2\mathrm{nC}$, $q_2=-3\mathrm{nC}$, ${\mathbf{r}}_1=0$, ${\mathbf{r}}_2=0.20\mathrm{m}\hat{\mathbf{x}}$ and the approximate value $k=8.99\times 10^9\mathrm{N}{\mathrm{m}}^2/{\mathrm{C}}^2$. ThenThe force is attractive, of magnitude $1.3485\times 10^{-6}\mathrm{N}$. The source electric field at charge 2 is ${\mathbf{E}}_1({\mathbf{r}}_2)=449.5N/C\hat{\mathbf{x}}$, pointing away from the positive source. Multiplication by the negative $q_2$ reverses the direction to give the attractive force. This field excludes charge 2's own self-field.

Holding the electric charges fixed while changing the separation from $R$ to $2R$ or $3R$ gives the normalized inverse-square law valuesFor the stipulated example, separations $0.20$, $0.40$ and $0.60\mathrm{m}$ give magnitudes $1.3485\times 10^{-6}$, $3.37125\times 10^{-7}$ and $1.498333\dots \times 10^{-7}\mathrm{N}$. The diagram computes these ratios and uses one common arrow-length scale for both sign choices.

The inverse-square law has a geometric explanation through Gauss's law: a source point charge gives fixed total electric flux through a concentric sphere, whose surface area is $4\pi r^2$. Dividing that flux by the area gives a radial electric field proportional to $1/r^2$. This connects a pairwise force calculation to a global conservation-like accounting of field sources. For several stationary point charges, the superposition principle adds the individual electric fields as vectors, and ${\mathbf{F}}_i=q_i{\mathbf{E}}_{\mathrm{other}}({\mathbf{r}}_i)$ uses all sources except charge $i$ itself. The Coulomb potential energy gives an alternative route when the desired quantity is work rather than instantaneous force.

The convex hull of the positive and negative coordinate unit vectors. Its boundary triangulates a sphere, and each face contains at most one member of each antipodal pair of coordinate unit vectors. This makes its barycentric subdivision useful for triangulating Real projective space.

Declination $\delta$ is the signed angle from the reference celestial equator to a direction, measured along its hour circle: $-90^{\circ }\le \delta \le +90^{\circ }$. Positive declination means celestial north, rather than necessarily north of the observer's zenith. A direction with $\delta =0$ lies on the celestial equator; $\delta =+90^{\circ }$ is the north reference celestial pole. Declination is an angle, not a height or a measure of Euclidean distance. ^[ref].

The connection to an observation follows from a dot product. In the elementary spherical-Earth, unrefracted model, an observer at geographical latitude $\varphi$ has zenith direction $\mathbf{z}=(\cos \varphi \cos \Theta ,\cos \varphi \sin \Theta ,\sin \varphi )$ in the same reference axes. Since $\mathbf{n}\cdot \mathbf{z}=\sin a$ for altitude $a$, the cosine addition formula givesAt upper meridian transit, $H=0$, so $a_{\mathrm{upper}}=90^{\circ }-|\varphi -\delta |$. For the synthetic site $\varphi =52^{\circ }$ and synthetic direction $\delta =20^{\circ }$, this yields $a_{\mathrm{upper}}=58^{\circ }$. A direction with $\delta =70^{\circ }$ instead reaches $72^{\circ }$, north of the zenith; positive declination does not imply a crossing on the southern half of the celestial meridian. A negative computed upper altitude means the direction never rises above the geometric astronomical horizon at that site.

For an angular readout, define signed zenith distance $z_s$ on the upper celestial meridian, positive toward the south and negative toward the north. Then $\delta =\varphi -z_s$. A synthetic southward readout $z_s=32^{\circ }$ at $\varphi =52^{\circ }$ recovers $\delta =20^{\circ }$. An unsigned distance $|z_s|$ alone would leave two possibilities, $\varphi \pm |z_s|$, when both are in range. A calibrated optical telescope must therefore record which side of the zenith it points toward. In this linear relation, a $1^{\prime }$ error in the assumed geographical latitude produces a $1^{\prime }$ error in inferred declination; averaging repeated readings cannot remove that shared offset. Near the astronomical horizon, uncorrected atmospheric refraction also invalidates the geometric readout.

The unit vector construction explains a subtle feature of sky maps. At fixed declination, changing right ascension traces a parallel whose radius on the unit sphere is $\cos \delta$. Its arc length is $\cos \delta |d\alpha |$, with angles in radians; at $\delta =60^{\circ }$, a $15^{\circ }$ right ascension increment spans a $7.5^{\circ }$ arc along that parallel. This arc is not generally the shortest great circle separation. At either celestial pole, all right ascension values represent the same direction: the coordinate singularity belongs to the coordinate grid. The unit vector remains regular there.

This collects the gradient, divergence, curl and Laplacian formulas in cylindrical and spherical coordinates. Coordinate unit vectors vary with position, so taking their derivatives is essential when translating vector calculus operators from Cartesian coordinates. Curl in spherical coordinates is one of these formulas.

If $F:{\mathbb{R}}^n\to \mathbb{R}$ is differentiable at the point $\gamma (t)$ of a differentiable parametric curve, the multivariable chain rule givesThe gradient supplies the local spatial rates; the velocity vector ${\gamma }^{\prime }(t)$ supplies the direction and rate of movement. Their dot product is the change in the sampled scalar field per unit parameter. If the parameter is time, this is a rate per unit time; a directional derivative per unit distance instead uses a unit vector.

Take $F(x,y)=x^2+2y^2$ and the unit circle path $\gamma (t)=(\cos t,\sin t)$. Along the path,At $t=\pi /4$, the gradient is $(\sqrt{2},2\sqrt{2})$, the velocity vector is $(-1/\sqrt{2},1/\sqrt{2})$, and their dot product is $-1+2=1$. The sampled value is $3/2$ and increases at rate one, even though one coordinate is decreasing. This is why a moving observation must combine every coordinate's contribution.

For a level curve, $F(\gamma (t))$ is constant, so the derivative along a curve is zero and its velocity vector is orthogonal to the gradient. Along other paths, the same formula combines a field's geometry with the observer's motion. ^[ref].

The directional derivative tests a function along a chosen straight line, including lines that follow neither coordinate axis. For a fixed vector $v$ and a point $a$ in an open subset of ${\mathbb{R}}^n$, defineif the limit exists. With a unit vector $v$, the parameter $t$ measures signed distance in the input Euclidean space, so the directional derivative is change per unit distance. For a general $v$, it is change per unit parameter; the norm of $v$ also sets the rate of travel.^[ref]

For a differentiable real-valued function, the multivariable chain rule givesThus one gradient predicts the local response in every direction. This formula requires differentiability, rather than merely separate partial derivatives.^[ref]

Return to $f(x,y)=x^2+xy+2y^2$ at $a=(1,1)$. Its gradient is $(3,5)$. In the unit direction $u=(3/5,4/5)$,One can check this without invoking the gradient formula: substitution givesThe difference quotient is $29/5+(53/25)t$, whose limit is $29/5$. This calculation separates the exact change along a finite path from its first-order rate at the starting point.

Changing the travel speed changes the parameter rate. Taking $v=(3,4)=5u$ gives $D_{v}f(a)=29$, because the input travels five distance units per parameter unit. For a differentiable function, $D_{cv}f(a)=c{D}_{v}f(a)$ for any real $c$: negative $c$ reverses the orientation as well as scaling the rate. Normalizing $v$ before comparing different slopes avoids confusing direction with speed.

The Cauchy-Schwarz inequality bounds the rate in any unit direction by $\Vert \nabla f(a)\Vert$. Here the largest rate is $\sqrt{34}\approx 5.831$, attained in the direction $(3,5)/\sqrt{34}$; the chosen $(3/5,4/5)$ gives the slightly smaller rate $5.8$. The perpendicular unit direction $(-5,3)/\sqrt{34}$ gives zero. That means no first-order height change, rather than no change at all: along this direction the example has $f(a+tu)=4+(28/34)t^2$.

This is useful for following a height or temperature scalar field, choosing a local search direction in mathematical optimization, and distinguishing a coordinate sensitivity from a response along an actual path. For a curved path, derivative along a curve uses its instantaneous tangent vector. Even existence of all directional derivatives does not guarantee one common linear approximation; directional derivatives do not imply differentiability gives a concrete counterexample.