Codex @codex  0

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.

Rotate the standard basis through $30^{\circ }$ using the rotation matrix whose columns are the new basis vectors:Here $u\cdot v=0$ and $u\cdot u=v\cdot v=1$, so $Q$ is an orthogonal matrix with $Q^{T}Q=I$, where $T$ denotes transpose and $I$ is the identity matrix. For the same vector $p=(2,1{)}^T$, its coordinate vector in this Cartesian basis isThus $\alpha \approx 2.232$ and $\beta \approx -0.134$: travel forwards along $u$ and a little backwards along $v$ to reconstruct $p=\alpha u+\beta v$. The orthogonal projections onto the two axes give these signed coordinates because the measuring vectors are both perpendicular and of unit length.

The squared norm remainsThis is useful when choosing axes aligned with a geometric feature: the coordinates change but distances and dot products retain the same sum-of-products formula. The two directions of the change of basis must be kept straight: $Q$ reconstructs the vector in old coordinates, $p=Q[p{]}_{\mathcal{U}}$, and $Q^T$ extracts its new coordinates. Applying $Q$ to the old coordinate vector $p$ instead would actively rotate the vector.

A short proof explains both formulas. If $p=Qc$, multiplication by $Q^T$ gives $c=Q^{T}p$. Then $\Vert p{\Vert }^2=c^T{Q}^{T}Qc=c^{T}c$. Perpendicular unit-length columns make the matrix inverse equal the transpose; this is the computational convenience supplied by a Cartesian basis.

A Cartesian frame combines a chosen origin point $O$ with a Cartesian basis $\mathcal{E}$, including its ordering as an ordered basis. A point $P$ has coordinate vector $c$ when $\overrightarrow{OP}={\sum }_i{c}_i{e}_i$. The basis supplies measuring directions and unit lengths; the coordinate vector is the resulting tuple of numbers; the Cartesian frame additionally fixes where measurement starts. A vector already belongs to a vector space with a zero vector, so expressing it in a basis requires no additional choice of spatial origin.

For example, keeping the standard basis but moving the origin from $O=(0,0)$ to $O^{\prime }=(1,0)$ changes the coordinates of the physical point $P=(2,1)$ from $(2,1{)}^T$ to $(1,1{)}^T$. This happens because $\overrightarrow{O^{\prime }P}=\overrightarrow{OP}-\overrightarrow{O{O}^{\prime }}$ is a different displacement vector. By contrast, the two diagrams above keep $O$ fixed and change only the basis. In a fixed reference Cartesian frame, a new frame with origin column $o$ and basis-column matrix $Q$ describes points by $p=o+Qc$; the new coordinates are $c=Q^T(p-o)$ when $Q$ is an orthogonal matrix.

In exact classical identification of an unknown linear Boolean function $f_a(x)=a\cdot x$, a function oracle call returns one bit $f_a(x)$. The query complexity is exactly $n$ for $n\ge 1$. Queries on the standard basis inputs $e_j$ reveal $a_j$, proving sufficiency. For necessity, a depth-$k$ decision tree has at most $2^k$ answer histories and must distinguish $2^n$ hidden binary strings; hence $k\ge n$, even with adaptive inputs. A reversible circuit that preserves $x$ and updates a known target by exclusive or reveals the same one bit per call.

There is also a useful linear algebra proof. Along a fixed history, the chosen inputs form the rows of a matrix $M$ over the two-element finite field, and the replies specify $Ma$. If there are $k<n$ queries, the rank-nullity theorem gives a nonzero vector $v$ with $Mv=0$. Then $a$ and $a+v$ give identical replies throughout this history. This also explains why adaptive choices cannot evade the bound: equal earlier replies force the same later choices.

One of the straight lines through the origin parallel to a standard basis vector in a Cartesian coordinate system. The positive direction corresponds to positive values of that coordinate.

A coordinate functional extracts one coefficient in a specified finite basis. It is an element of that basis's dual basis, rather than a coordinate rule independent of the basis. For the standard basis of $F^n$, the $i$th coordinate functional sends $(a_1,\dots ,a_n)$ to $a_i$.

The determinant assigns a scalar to a square matrix. For a linear map on the real plane, it measures signed area scaling: its absolute value tells how much area grows, and its sign records whether the orientation of a vector space is preserved or reversed. In three dimensions, the same idea measures signed volume. Over a field, a square matrix is an invertible matrix exactly when its determinant is nonzero. These are two views of one phenomenon: losing an independent direction collapses full-dimensional volume and makes the original input impossible to recover uniquely.^[ref]

Start with an explicit matrix and its action on a vector:The columns are the images of the standard basis vectors: $A{e}_1=(2,1{)}^T$ and $A{e}_2=(1,3{)}^T$. Consequently the unit square with corners $(0,0),(1,0),(1,1),(0,1)$ becomes the parallelogram with corners $(0,0),(2,1),(3,4),(1,3)$. Its area is five. Every small patch of this square undergoes the same linear map, so the factor describes the whole deformation, rather than just one specially chosen boundary. The matrix mixes horizontal and vertical components, but its determinant compresses that combined effect into one area factor.^[ref]

Now exchange the columns:The parallelogram occupies exactly the same region, but its generating vectors arrive in the opposite order. Moving from the first column to the second now turns clockwise rather than counterclockwise. Ordinary area remains five; the oriented volume in two dimensions changes sign. A negative determinant therefore signals reversed orientation of a vector space, not a negative amount of material. Keeping the sign makes successive reversals calculable: two reversals restore the original orientation of a vector space.

Compare a singular matrix:Its columns are linearly dependent: the second is twice the first. The unit square collapses onto the segment from $(0,0)$ to $(3,6)$, which has zero area despite nonzero length. Since $C(-2,1{)}^T=0$, the kernel of a linear map contains a nonzero vector; inputs differing by that vector have identical outputs. The collapsed picture explains why an inverse matrix cannot exist.

For $A$, recovery is possible, and the complete matrix multiplication check isHere $I$ is the identity matrix. The inverse matrix restores the original square, shrinking its area by $1/5$. Generally, the adjugate identity gives $A^{-1}=\mathrm{adj}(A)/\det A$ whenever the determinant is nonzero over a field; this is why the denominator matters.^[ref]

The multiplicativity of the determinant describes two successive deformations. Matrix multiplication means that $AB$ applies $B$ first, then $A$. TakeThe first linear map doubles area; the second multiplies the resulting area by five. Thus the unit square passes through areas $1,2,10$, with the last parallelogram bounded by $(0,0),(2,1),(4,7),(2,6)$. Direct calculation agrees:This explains physically why area factors multiply instead of adding. The general identity holds for square matrices of the same size, including singular matrices. If either stage loses a direction, later processing cannot recover the discarded information.^[ref]

The same interpretation works in three dimensions without drawing a new picture. The diagonal matrixsends the unit cube to a rectangular box with side lengths two, three, and one. Its ordinary volume is six. The minus sign comes from reversing the third direction, so the oriented volume is minus six. Replacing the last diagonal entry by zero flattens the cube into a planar region: its area may be positive, but its three-dimensional volume is zero. This distinction matters whenever a linear map changes shape: the determinant measures size in the full input dimension, rather than the size of whichever lower-dimensional object remains.

The determinant measures neither trace nor perimeter. For $A$, the trace happens also to be five, but the identity matrix in two dimensions has trace two and determinant one. The first parallelogram above has perimeter $2(\sqrt{5}+\sqrt{10})$, compared with four for the unit square; boundary length follows a different rule. Nor does determinant one require every length to stay fixed: the diagonal matrix $\mathrm{diag}(2,1/2)$ preserves area while stretching one direction and compressing the other.

These properties make the determinant useful when coordinates obscure what a deformation does. A nonzero determinant certifies that linearly independent input directions remain distinguishable; its magnitude counts the area or volume represented by a transformed cell. That second role reappears in the change of variables formula, where the Jacobian determinant corrects the sizes of small cells. Signed volume also separates a deformation that reverses the orientation of a vector space from one that merely stretches it, even when both give identical ordinary volume factors.^[ref]

For general $n$, regard the columns $v_1,\dots ,v_n$ as ordered arguments. The determinant is the unique alternating multilinear map to the scalar field normalized by $\det I=1$. The multilinearity of the determinant means linearity in each column separately, with all others held fixed:Alternation means equal columns give zero; exchanging columns negates the result. Normalization sets the oriented volume of the standard basis to one. These properties also hold in rows. For real three-dimensional columns, their determinant is the scalar triple product, and its absolute value is the enclosed volume. Over other fields, the algebraic definition still works even without a geometric notion of positive volume.^[ref]

Separate linearity does not mean the determinant is a linear map on the whole matrix. Scaling every column of an $n$-by-$n$ matrix by $t$ gives $\det (tA)=t^n\det A$. For our two-dimensional example, $\det (2A)=4\det A=20$: doubling both directions quadruples area. The power records how many independent directions are being scaled together.

The properties suggest practical computations. Adding a multiple of one column to another preserves the determinant: the extra term has repeated columns and vanishes. This is the algebra behind an area-preserving shear mapping. Alternatively, expansion along a row uses cofactors, each a matrix minor multiplied by a sign:where $A^{(ij)}$ deletes row $i$ and column $j$. The Leibniz formula for determinants supplies the general entry formula. Both calculations encode the same geometric invariant; the defining properties explain their structure without requiring memorization of a separate recipe for every size.^[ref]

A determinantal function on $M_n(F)$ is an alternating multilinear map of the $n$ columns, normalized to take value one at the identity matrix. It is uniquely the determinant. Indeed, expand each column in the standard basis. Terms with a repeated basis vector vanish, and all remaining terms differ from the value on the standard basis only by the sign of a permutation. Thus an alternating multilinear map of $n$ columns is determined by its value on the identity matrix.

The exterior algebra of a vector space $V$ combines scalars, vectors, and higher-dimensional oriented elements in one graded algebra. In a finite-dimensional vector space it iswith multiplication given by the exterior product. Degree zero is the scalar field, degree one is $V$, and degree two is the space of bivectors. Multiplication adds degrees: combining two vectors produces a bivector, rather than another ordinary vector.^[ref]

The essential rule is $v\wedge v=0$. Repeating the same direction cannot produce a nonzero oriented area. Bilinearity then gives $u\wedge v=-v\wedge u$: expand $(u+v)\wedge (u+v)=0$ and remove the two zero self-products. The alternating rule remains part of the definition in characteristic two, where sign reversal alone would not force a self-product to vanish.

If $V={\mathbb{R}}^2$ with standard basis $e_1,e_2$, every element of its exterior algebra has the formThere are four independent coefficients: one scalar, two vector components, and one area component. In ${\mathbb{R}}^3$, the respective degree dimensions are $1,3,3,1$, giving eight components in total. More generally, the exterior power basis gives total dimension ${\sum }_{k=0}^n\left(\genfrac{}{}{0ex}{}{n}{k}\right)=2^n$.

This language is useful when direction and orientation matter together. Determinants measure the action on top-degree exterior powers; differential forms use exterior powers of a dual space to describe quantities integrated along paths, surfaces, and higher-dimensional regions. The Grassmann algebra viewpoint also makes the sign changes in anticommuting generators systematic. A single multiplication rule connects geometry and algebra.

The exterior product, also called the wedge product, combines oriented elements through the multiplicationFor two vectors, it is bilinear and alternating:These rules imply $v\wedge u=-u\wedge v$. They distinguish the exterior product from both the symmetric dot product and the unconstrained tensor product.^[ref]

In an oriented orthonormal plane, write $u=(a,b)$ and $v=(c,d)$. Expanding in its standard basis givesThe determinant $ad-bc$ is the signed area of the parallelogram formed by the two inputs. Its absolute value is the ordinary area. Reversing the input order reverses the sign; parallel inputs produce zero.

For the original example $u=(3,1)$ and $v=(1,2)$,Replacing $v$ by $v+u=(4,3)$ changes the shape but preserves the product: $3\cdot 3-1\cdot 4=5$. Both parallelograms have the same oriented area element. If $v=2u$, the exterior product is zero, since $u\wedge (2u)=2(u\wedge u)=0$.

The same rule adds directions in higher dimensions. With $w=(2,-1,1)$ and the three-dimensional inputs under bivector,Its coefficient is the oriented volume of the three-dimensional parallelepiped. Repeating or linearly depending on earlier inputs collapses the volume, giving zero.

For homogeneous elements of a graded algebra of degrees $p,q$, the general commutation rule isSwapping a bivector with another bivector therefore produces a positive sign, since $pq=4$. It is the number of individual vector swaps that controls the sign, not a rule that every pair of objects must anticommute.

Here is an original reproducible finite set with 40 points and 512 coordinates. Let $e_1,\dots ,e_4$ be the first four vectors of the standard basis of ${\mathbb{R}}^{512}$. With point indices $i=0,\dots ,39$, generateThe four groups have different centers, and the standard normal random variables perturb all 512 coordinates. Generate the coordinates in row order using NumPy's PCG64 pseudorandom number generator with seed 20261013. Then keep this finite set fixed. A second PCG64 stream with seed 424242 generates a $32\times 512$ matrix of independent random variables with the standard normal distribution; divide it by $\sqrt{32}$ to obtain the Gaussian Johnson–Lindenstrauss transform.

For all 780 unordered pairs, compute the distance ratio $r_{ij}=\Vert T{x}_i-T{x}_j{\Vert }_2/\Vert x_i-x_j{\Vert }_2$. In this draw, its minimum is approximately $0.733010$ and its maximum $1.478785$, giving a maximum relative error in Euclidean distance of about $0.478785$. The squared-distance ratios $r_{ij}^2$ instead range from approximately $0.537304$ to $2.186806$: their maximum relative error is about $1.186806$. Thus this particular draw fails a 20% distance tolerance. The Johnson–Lindenstrauss lemma does not assert that every draw into 32 dimensions must succeed.

The two bands of original Euclidean distances distinguish pairs within a group from pairs between groups. Their spread after the linear map shows why checking an average distance is insufficient: the Johnson–Lindenstrauss lemma controls every pair. The dashed 20% distance band is an illustrative tolerance, not a confidence interval or a squared-distance band. The affine hull of 40 points has dimension at most 39; 512 is the ambient coordinate count, and no full-dimensionality assumption is required. For this draw the bi-Lipschitz distortion is approximately $1.478785/0.733010=2.017414$.

Now keep the same points and draw 48 fresh Gaussian Johnson–Lindenstrauss transforms for each target dimension $k\in \{8,16,32,64,128\}$. NumPy's SeedSequence(314159).spawn(240) supplies distinct PCG64 streams, ordered first by target dimension, then by draw. Each stream generates a new $k\times 512$ random matrix, normalized by $1/\sqrt{k}$; no map is a row prefix of another map. Each plotted value is

The observed median maximum relative errors in Euclidean distance are approximately $0.735$, $0.520$, $0.381$, $0.267$, and $0.184$, respectively. All pairs fall inside the 20% distance tolerance in $0$, $0$, $0$, $0$, and $35$ of the 48 draws, respectively. These are computations on one synthetic finite set, rather than optimal dimension bounds or measured physical data. The displayed quantiles describe the sampled maps; they do not certify a success probability for other sets. Even at 128 dimensions, 13 of the sampled maps fail this tolerance.

The practical benefit is that each point acquires fewer coordinates while its pairwise Euclidean distances can remain useful. The choice of target dimension balances compression, accuracy, and failure probability. For example, a conservative sufficient squared-distance tolerance $\epsilon =0.2$ with $n=40$ and $\delta =0.05$ gives $k\ge 1195$ by the explicit bound above. Here that exceeds the ambient dimension, so the identity map is better for a certified guarantee. The much smaller dimensions in the plots illustrate observed behavior without claiming that this sufficient bound is sharp.

On $K^4$, use the standard basis $\mathcal{U}=(u_1,u_2,u_3,u_4)$ and the alternating bilinear form with matrix $J$ above for $n=2$. A second symplectic basis isFor instance, $\omega (f_1,f_2)=1-1=0$, while $\omega (e_1,f_1)=\omega (e_2,f_2)=1$ and the cross-pairings vanish. The change of basis matrix, whose columns are these new basis vectors in $\mathcal{U}$, isThe symmetric matrix $C$ makes $S^{T}JS=J$, so $S$ is a symplectic matrix. Explicitly, multiplying the blocks gives a lower-right block $C^T-C=0$; the other blocks agree with $J$. For a fixed vector, change of basis means $[v{]}_{\mathcal{U}}=S[v{]}_{\mathcal{B}}$, henceTheir symplectic form gives $3\cdot 0+2\cdot 1-1\cdot 1-1\cdot 0=1$ in $\mathcal{U}$, and $1\cdot 0+1\cdot 1-1\cdot 0-1\cdot 0=1$ in $\mathcal{B}$. The separate coordinate-plane contributions change, but their sum does not.

The same matrix also represents an active symplectic linear map: keeping one basis fixed, it sends $(q,p)$ to $(q+Cp,p)$. The following diagram shows this active interpretation over the real numbers. Its preserved total pairing is precisely the test used for linear canonical transformations.

A linear Boolean function is a linear map from ${\mathbb{F}}_2^n$ to ${\mathbb{F}}_2$. Here ${\mathbb{F}}_2$ is the two-element finite field, whose addition is exclusive or. Every such function has a unique coefficient binary string $a$ and takes the formThe dot is the binary inner product. To prove the representation, expand $x={\sum }_j{x}_j{e}_j$ in the standard basis and use linearity, obtaining $f_a(x)={\sum }_j{x}_j{f}_a(e_j)$; thus $a_j=f_a(e_j)$. For example, $a=101$ gives $f_a(x)=x_1\oplus x_3$. Adding a nonzero constant produces an affine map, not a linear map.

For an invertible matrix $A\in {\mathbb{R}}^{n\times n}$ and real vectors $u,v\in {\mathbb{R}}^n$, the matrix determinant lemma statesThe expression $v^T{A}^{-1}u$ is a scalar, so a rank-one matrix update changes the determinant by a single multiplicative factor. One need not expand the updated determinant entry by entry.

For the preceding example with $A=I$, the dot product is $v^{T}u=3-2=1$. Hence the updated determinant is $1+1=2$. Direct evaluation confirmsMore generally, let $B(t)=I+tu{v}^T$ for a real scalar $t$. The same calculation gives $\det B(t)=1+t$. The determinant is one at $t=0$, two at $t=1$, and zero at $t=-1$.

A compact proof explains the scalar reduction. Put $w=A^{-1}u$ and factor $A+u{v}^T=A(I+w{v}^T)$. If $w\ne 0$, choose a basis starting with $w$. In that basis, $(I+w{v}^T)w=(1+v^{T}w)w$, while every other basis vector changes only by a multiple of $w$. Its matrix is triangular with diagonal $1+v^{T}w,1,\dots ,1$, giving the required determinant. If $w=0$, the correction vanishes. This uses multiplicativity of the determinant to restore the factor $\det A$.

The denominator in the Sherman–Morrison formula is exactly this factor. The rank-one update singularity criterion therefore identifies when the updated matrix inverse ceases to exist. Batson, Spielman and Srivastava use both identities for the symmetric update $A+v{v}^T$ in Twice-Ramanujan Sparsifiers, section 2.3, Lemmas 2.1–2.2.

For an $n$-by-$n$ square matrix $A$ over a field, its matrix inverse is an $n$-by-$n$ matrix $A^{-1}$ satisfying both identitieswhere $I_n$ is the identity matrix. It exists exactly when $A$ is an invertible matrix, equivalently when its matrix rank is $n$ or its determinant is nonzero. It reverses the linear map: applying $A$ and then $A^{-1}$ recovers the input, while applying $A^{-1}$ and then $A$ recovers the output.^[ref]

The matrix inverse is unique. If $B$ and $C$ both satisfy the identities, then $B=B(AC)=(BA)C=C$. Although the definition is two-sided, for finite square matrices either identity alone suffices: $BA=I$ forces a trivial kernel, whereas $AB=I$ forces surjectivity; either makes $A$ an invertible matrix by the rank-nullity theorem. Rectangular matrices may have only one of these properties.

For a concrete example, takeDirect matrix multiplication givesConsequently the vector $v=(1,2{)}^T$ follows a reversible journey:The columns $(2,1{)}^T$ and $(1,1{)}^T$ are the images of the standard basis vectors. They turn the unit square into a parallelogram; the matrix inverse restores both the square and its surrounding grid. Since $\det A=1$, this example changes shape while preserving area, illustrating that invertibility does not require preserving lengths or angles.

For comparison, the singular matrixcannot have a matrix inverse: the output $(1,0{)}^T$ cannot specify which of the two inputs occurred. Their difference $(1,-1{)}^T$ belongs to its kernel. Indeed $S(x,y{)}^T=(x+y,0{)}^T$ collapses the plane to a line, its matrix rank is one, and its determinant is zero. No subsequent function of the output alone can reconstruct every input.

An important use of the matrix inverse is describing the unique solution of a linear system:For $b=(4,3{)}^T$, the preceding computation solves $2x_1+x_2=4$ and $x_1+x_2=3$, giving $x=(1,2{)}^T$. In computation, usually solve the linear system using Gaussian elimination and a factorization such as an LU decomposition, rather than forming every entry of the matrix inverse. The inverse formula specifies the answer; a factorization provides a reusable way to compute it for each right-hand side.^[ref]

Here the LU decomposition is especially simple:Solve $Ly=b$ by forward substitution in a triangular system: $y_1=4$ and $y_2=3-\frac{1}{2}\cdot 4=1$. Then solve $Ux=y$ by backward substitution in a triangular system: $x_2=2$ and $x_1=(4-2)/2=1$. The same factors work when $b$ changes. If the whole matrix inverse is actually needed, its $j$th column is the solution of $Ax=e_j$, where $e_j$ is the $j$th standard basis vector.

Reversing several operations means undoing the last one first. For invertible matrices of the same size,because $(AB)(B^{-1}{A}^{-1})=A(B{B}^{-1})A^{-1}=I$, and the product in the other order is also $I$. This is why a change-of-basis matrix uses its matrix inverse to reverse the direction of coordinate conversion.

Two square matrices $A,B$ over the same field are similar if $B=S^{-1}AS$ for an invertible matrix $S$. They can describe one linear operator using two different ordered bases, with the same basis used for input and output in each description. Fix the convention $x=Sy$: $x$ is the old coordinate vector, $y$ the new coordinate vector, and the columns of $S$ are the new basis vectors in old coordinates. This is the convention of the change-of-basis matrix.^[ref]

For a concrete linear operator on the physical plane, let $T$ double horizontal components and multiply vertical components by five. In the old standard basis $\mathcal{E}=(e_1,e_2)$, its matrix of a linear map is $A=\mathrm{diag}(2,5)$. Choose the new ordered basis $\mathcal{F}=(b_1,b_2)$ with $b_1=e_1$, $b_2=e_1+e_2$. ThenThe second column of $B$ says $T(b_2)=-3b_1+5b_2=2e_1+5e_2$. Its off-diagonal entry records how the chosen basis vectors describe the output; the physical linear operator still stretches along the old horizontal and vertical directions.

Take the physical vector $p=2e_1+e_2=b_1+b_2$. Its two coordinate vectors, and the two descriptions of $T(p)$, areReconstruction gives $T(p)=4e_1+5e_2=-b_1+5b_2$. The output coordinate vectors differ, but the endpoint in the physical plane is the same.

This is why matrix similarity is useful: a convenient basis can simplify calculations without changing the linear operator. A diagonalizable matrix can be replaced by a diagonal matrix, making repeated application a calculation with powers of its eigenvalues. In the example, $A$ makes the two independent stretch factors visible, while $B$ expresses the same action in directions chosen for another purpose. Matrix similarity invariants separate properties of the linear operator from features of its address system.

A partial derivative measures the local response of a function to one input while its other inputs stay fixed. For a function $f$ on an open subset of ${\mathbb{R}}^2$, its partial derivatives at $(a,b)$ are the limitswhen the respective limits exist. Each is an ordinary derivative of a one-variable slice: $x\mapsto f(x,b)$ or $y\mapsto f(a,y)$. This is why the other variable is treated as a constant during partial differentiation.^[ref]

For an original example, take the polynomial functionHolding $y$ fixed gives $f_x=2x+y$, while holding $x$ fixed gives $f_y=x+4y$. At $(1,1)$, the height is $f=4$ and the two partial derivatives are $f_x=3$ and $f_y=5$. The values answer different experiments on the same mathematical model: change $x$ alone, or change $y$ alone.

The two exact slices through that point areDividing either height change by $t\ne 0$ gives $3+t$ or $5+2t$; taking the limit $t\to 0$ recovers the two slopes. At $t=0.1$, the finite height changes are $0.31$ and $0.52$, although their linear approximations are $0.3$ and $0.5$. A derivative measures the limiting response, rather than an exact finite change.

For a function of $n$ variables, ${\partial }_{j}f(a)$ varies only its $j$th coordinate: it is the same difference quotient along $e_j$, the $j$th standard basis vector. This gives a practical sensitivity analysis: one can compare how a mathematical model responds to changes in its individual inputs, provided the chosen system of units and the quantities held fixed are stated. Comparing the numerical sizes of two partial derivatives in different units does not by itself identify a physically preferable input.

Existence of the partial derivatives at one point is weaker than differentiability. The example under partial derivatives do not imply continuity shows how both coordinate tests can miss a bad diagonal approach. When the first partial derivatives are continuous in a neighborhood, they do give differentiability there, allowing the total differential to combine all small input changes.

Write $a=a_1{a}_2$ and $x=x_1{x}_2$, with the top wire carrying $x_1$, the middle wire $x_2$, and the bottom wire the ancilla bit $y$. The binary inner product is $a\cdot x=a_1{x}_1\oplus a_2{x}_2$, where $\oplus$ is exclusive or. Let $C_{j\to y}$ denote a CNOT gate controlled by wire $j$ and targeting $y$. The four linear Boolean functions and their Boolean quantum oracles are:

The CNOT gates in the last row commute: their common target acquires the exclusive or of the two controls. Each input wire is preserved, and applying any of these reversible circuits twice restores $y$.

For classical reversible computation, call the unknown function oracle on $|10⟩|0⟩$ and $|01⟩|0⟩$. The resulting ancilla bits are respectively $a_1$ and $a_2$:Use a fresh zero ancilla bit for each call, or preserve the first result in workspace before restoring the target to zero. Thus two classical calls determine the two-bit string.

For $n$ input bits, query the standard basis strings $e_j$, since $a\cdot e_j=a_j$. This gives an upper bound of $n$ calls. It is also optimal for exact identification: every classical query reveals at most one new bit, since $x$ is unchanged and $y\oplus f_a(x)$ reveals just $f_a(x)$ when $y$ is known. A depth-$k$ decision tree has at most $2^k$ answer histories, but must distinguish $2^n$ possible hidden strings. Hence $k\ge n$, including adaptive strategies. The classical query complexity of linear Boolean function identification is therefore

The increments ${\Delta }_i=d_i-d_{i-1}$ are positive integers by part (a)(i), are nonincreasing by part (a)(ii), and sum to $d_m=4$. Thus they form an integer partition of $4$. The complete list is $(4)$, $(3,1)$, $(2,2)$, $(2,1,1)$, and $(1,1,1,1)$.

To prove all five occur over every field, let $J_s(0)$ be a Nilpotent Jordan block of size $s$. It sends each standard basis vector to its predecessor, with the first sent to zero, soFor a block diagonal matrix, the nullities add. The following explicit nilpotent matrices therefore realize every possible sequence:In each row the largest Jordan block size is exactly the nilpotency index $m$, so the displayed sequence stops at the correct place.

Take positive angles counterclockwise. The rotation sends the standard basis vectors to $(\cos \alpha ,\sin \alpha )$ and $(-\sin \alpha ,\cos \alpha )$. These images form the columns of its matrix, soThe columns are an orthonormal basis and the determinant is one, as required for a plane rotation.