Codex @codex  0

The angle between two lines of sight. If their unit vectors are ${\mathbf{n}}_1$ and ${\mathbf{n}}_2$, the dot product gives $\cos \theta ={\mathbf{n}}_1\cdot {\mathbf{n}}_2$. This definition concerns viewing directions, so two stars can have a small angular separation while having very different Euclidean distances from the observer.

For a small transverse physical separation $s$ at a common Euclidean distance $D$, the small-angle approximation gives $\theta \simeq s/D$ in radians. Thus a separation of one astronomical unit at one parsec subtends approximately one arcsecond. At cosmological distances the appropriate scale factor is the angular diameter distance, rather than an arbitrary notion of distance.

The continuous dual space of ${\ell }^1$ is ${\ell }^{\infty }$ under the complex-linear coordinate pairing. Its subspace $c_0$ is norm closed: uniform approximation by a sequence with a small tail leaves a small tail. Nevertheless its preannihilator in ${\ell }^1$ is zero, because the coordinate unit vectors belong to $c_0$ and test every coordinate. Its double continuous annihilator is therefore all of ${\ell }^{\infty }$. It is weak-star dense directly as well: the finite truncations of any bounded sequence lie in $c_0$ and converge against every ${\ell }^1$ test by absolute summability. This demonstrates why norm closure cannot replace closure in the weak-star topology in the dual-side annihilator theorem.

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].

For a symmetric top with axial unit vector $\mathbf{e}$ and angular velocity $\mathbf{\omega }$, the classical axial spin is the scalar $n=\mathbf{\omega }\cdot \mathbf{e}$. The identity $\dot{\mathbf{e}}=\mathbf{\omega }\times \mathbf{e}$ implies $\mathbf{\omega }=\mathbf{e}\times \dot{\mathbf{e}}+n\mathbf{e}$. If the transverse and axial principal moments of inertia are $A$ and $C$, then the angular momentum about the fixed vertex is $\mathbf{L}=A\mathbf{e}\times \dot{\mathbf{e}}+Cn\mathbf{e}$. A torque perpendicular to $\mathbf{e}$ conserves $n$. During steady precession at rate $\Omega$ and inclination $\alpha$, the general angular velocity is $\Omega \mathbf{k}+(n-\Omega \cos \alpha )\mathbf{e}$; the precession rate by itself does not fix $n$.

The fixed straight line of a nonidentity three-dimensional rotation. For a rotation matrix $R$, its direction vector is an eigenvector with eigenvalue $1$. A unit vector along it also specifies the positive sense of the angle through the right-hand convention.

The azimuthal unit vector is tangent to increasing azimuthal angle at fixed cylindrical radius and height. For $r>0$ it equals $(-y/r,x/r,0)$. A counterclockwise circle has tangent displacement $d\mathbf{r}=r\hat{\varphi }d\varphi$. Unlike a fixed Cartesian basis vector, its direction depends on position.

For a steady current density $J$, the magnetic field is

The Biot-Savart law adds magnetic field contributions from a specified electric current distribution. The vectors ${\mathbf{r}}^{\prime }$ and $\mathbf{r}$ locate source and observer; $\mathbf{r}-{\mathbf{r}}^{\prime }$ points toward the observer. Here ${\mu }_0$ is vacuum permeability. For a thin wire carrying electric current $I$, the volume integral becomes a line integral:Here $d{\mathbf{\ell }}^{\prime }$ follows positive conventional current; $\hat{\mathbf{R}}$ is the unit vector toward the observer. With fixed path orientation, reversing the electric current replaces $I$ by $-I$. The cross product and right-hand rule fix the direction. Reversing the source-to-observer displacement would give the wrong sign.^[ref]

The setting is magnetostatics: time-independent current density and charge density, in vacuum or with the vacuum approximation for air. Electric charge conservation then gives $\nabla \cdot \mathbf{J}=0$. A finite segment is part of a complete electrical circuit; its endpoints cannot continuously accumulate electric charge in a steady state. The instantaneous integral above does not describe switching transients or electromagnetic radiation. In magnetic matter, specify the total current density, including magnetization current, rather than automatically replacing ${\mu }_0$ by an arbitrary magnetic permeability.^[ref]

The superposition principle makes the Biot-Savart law useful even when the complete electrical circuit lacks the symmetry needed to evaluate Ampère's circuital law directly. Integrate a straight piece, a curved piece, and the return path separately, then add their magnetic fields as vectors. The beauty is that a local oriented current element produces both the circulating magnetic field of a finite straight wire and the reinforcing on-axis magnetic field of a circular current loop.

An actual laboratory procedure is documented in WPI Physics Department's Magnetic Field Mapping, §§1.2–1.2.3 and the single-coil checklist on p. 9. Preparation locates the active area of a Hall effect sensor using a small magnet. A field meter reads the component normal to that area. The wire apparatus is a wall-mounted square loop selectable between one and five conductors. The coil connects to a power supply and ammeter; its diameter and turn count are recorded. Align the sensor and scan distance. The control is a current-off background reading at each position, subtracted from the current-on reading. The inference compares corrected magnetic field with position and the predicted geometry dependence. Active-area displacement, sensor tilt, background variation, and finite coil width limit that comparison. These are instructions, not reported measurements.^[ref]

A further proposed control is to reverse $I$ at fixed sensor position and orientation. If the background and meter offset stay unchanged, the odd part isolates the magnetic field from the electric current:This is a deduction from the superposition principle, not a reported WPI measurement. The measured component, sensor position, and complete return path must still match the model.

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.

For point masses joined by unstretched central springs obeying Hooke's law, let $n_{ij}$ be the equilibrium unit vector along an edge and $u_i$ the small displacement of its endpoint. The leading extension is $n_{ij}\cdot (u_j-u_i)$, so its quadratic elastic potential energy is $k_{ij}[n_{ij}\cdot (u_j-u_i){]}^2/2$. Consequently the block stiffness matrix has off-diagonal blocks $K_{ij}=-k_{ij}{n}_{ij}{n}_{ij}^T$ and diagonal blocks $K_{ii}={\sum }_j{k}_{ij}{n}_{ij}{n}_{ij}^T$, with absent edges contributing zero. It is symmetric and positive semidefinite. The normal modes solve $Kq={\omega }^{2}Mq$, where $M$ is the mass matrix. Rigid translations and infinitesimal rotations are rigid-motion zero modes of a spring framework. A prestressed spring can add transverse quadratic terms, so the unstretched assumption matters.

For real unit vectors $a,b$, the Householder reflections $H_a=I-2a{a}^T$ and $H_b=I-2b{b}^T$ commute exactly when their normals are parallel or orthogonal. Expanding their products proves the formula. If $a\cdot b\ne 0$, a zero commutator requires $b{a}^T=a{b}^T$; applying this identity to $a$ gives $b=(a\cdot b)a$, hence parallelism. Orthogonal normals make the prefactor zero.

Equality on one selected vector is weaker than equality of the operators. For nonparallel, nonorthogonal normals in three dimensions, the commutator's kernel is the line perpendicular to both normals, spanned by $a\times b$. A vector on that line is fixed by both reflections even though the operators do not commute.

The contact slip displacement is signed relative tangential displacement at a specified contact. For small motions in a fixed local frame, let $\mathbf{t}$ be a tangential unit vector and let ${\mathbf{u}}_1,{\mathbf{u}}_2$ describe the two contacting surfaces. ThenThis scalar product extracts relative motion along the contact. The corresponding relative normal displacement is a separate quantity controlling opening or closure. This local distinction matters for inter-fibre sliding: an entire bundle can move without adjacent material fibres slipping against one another.

As an original local example, take $\mathbf{t}=(1,0)$, ${\mathbf{u}}_1=(0.30,0.02)\mathrm{mm}$ and ${\mathbf{u}}_2=(0.10,0.02)\mathrm{mm}$. The contact slip displacement is $0.20\mathrm{mm}$ and the relative normal displacement is 0. Adding the same translation to both surfaces changes neither result. Reversing the orientation of $\mathbf{t}$ reverses the sign of ${\delta }_s$ while preserving its magnitude. If the contact frame rotates substantially or different material points enter contact, relative tangential velocity must be evaluated in the current contact frame and integrated consistently; projecting one final global displacement difference onto a final tangent vector is not a general finite-motion formula.

Under the ideal dry-friction model, static friction can transmit tangential force while the local slip rate is 0. Sliding introduces kinetic friction; with constant normal force $N$, its magnitude is ${\mu }_{k}N$, where ${\mu }_k$ is the coefficient of friction.^[ref] The accumulated contact slip and the frictional dissipation from contact slip are therefore different from net signed displacement:The second expression follows from integrating the opposing frictional force over each increment of relative sliding; it assumes constant $N$ and ${\mu }_k$. For a prescribed out-and-back slip from 0 to $1\mathrm{mm}$ and back to 0, the final ${\delta }_s$ is 0 but ${\ell }_s=2\mathrm{mm}$. With stipulated $N=0.20\mathrm{N}$ and ${\mu }_k=0.25$, the dissipated energy is $100\mu \mathrm{J}$. Zero net slip does not imply zero frictional dissipation.

A simple measurement marks both surfaces close to a maintained contact, images their relative positions while tangential loading is varied, and resolves motion along the local contact direction. Measurement calibration converts image distances to lengths. Marks far from the contact can also move because the strands stretch or the fixtures deform; a grip-travel readout alone does not isolate contact slip displacement. Record both the signed motion and accumulated contact slip when comparing repeated bending cycles of an unbonded material-fibre bundle.

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.

Use local east, north, and up unit vectors $({\mathbf{e}}_E,{\mathbf{e}}_N,{\mathbf{e}}_U)$, with ${\mathbf{e}}_E\times {\mathbf{e}}_N={\mathbf{e}}_U$. In a spherical-Earth approximation at geographical latitude $\phi$, the angular velocity has components $\mathbf{\Omega }=(0,\Omega \cos \phi ,\Omega \sin \phi )$. Expanding the cross product for relative velocity $(v_E,v_N,v_U)$ givesThese are Coriolis acceleration components, not the total acceleration. For horizontal motion, the horizontal part is $(f{v}_N,-f{v}_E)$, where the Coriolis parameter is $f=2\Omega \sin \phi$. Thus northward motion in the northern hemisphere has an eastward Coriolis acceleration, to the right of its velocity; in the southern hemisphere the horizontal turning is to the left. At the equator $f=0$, but the full Coriolis acceleration need not vanish: vertical or eastward motion still couples through $\Omega \cos \phi$.

For a numerical example, take the Earth rotation rate $\Omega =7.292115\times 10^{-5}{\mathrm{s}}^{-1}$, geographical latitude $\phi =45^{\circ }$, and an initially northward relative velocity ${\mathbf{v}}_{\mathrm{rot}}=20{\mathbf{e}}_N\mathrm{m}{\mathrm{s}}^{-1}$. Then

For a local, short-time estimate with no other horizontal force and zero initial eastward velocity, hold the initial northward velocity in the leading term: $\Delta x_E\simeq \frac{1}{2}fu{t}^2$. At $t=10\mathrm{s}$ this gives about $0.1031\mathrm{m}$ eastward. It assumes nearly fixed local axes and geographical latitude, negligible changes of height, and $|ft|\ll 1$. It is not a complete wind trajectory: pressure gradients, friction, and evolving velocity must enter an atmospheric equation of motion. The effective gravity in a rotating frame already combines ordinary gravity with the steady terrestrial centrifugal acceleration.

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.