Type at least 2 letters. Enter opens the first result, Esc closes.

    Glossary

    The technical terms of this track, with short definitions (and the Persian term where available). Each term links to the lesson where it is introduced.

    AA-conjugate
    Vectors with pi⊤Apj=0p_i^\top A p_j = 0 for i≠ji \ne j: orthogonal in the inner product defined by an SPD matrix AA. Introduced in: Lesson 14
    affine map نگاشت آفین
    A linear map plus a fixed shift, x↦Wx+bx \mapsto Wx + b. It is not linear when b≠0b \ne 0. Introduced in: Lesson 1
    algebraic multiplicity
    The number of times λ\lambda is a root of the characteristic polynomial. Introduced in: Lesson 5
    Arnoldi iteration
    Builds an orthonormal basis QkQ_k of the Krylov subspace with Gram–Schmidt, and the upper Hessenberg matrix Hk=Qk⊤AQkH_k = Q_k^\top AQ_k. Introduced in: Lesson 14
    backward error
    The smallest relative change of the data that makes the computed answer exact. It can be computed from the residual. Introduced in: Lesson 13
    backward stable
    An algorithm whose backward error is O(u)O(u) for every input. Its forward error is then at most about κ⋅u\kappa \cdot u. Introduced in: Lesson 13
    basis پایه
    An independent list of vectors that spans a subspace. Every vector of the subspace has unique coordinates in it. Introduced in: Lesson 1
    Bauer–Fike theorem
    For diagonalizable A=XΛX−1A = X\Lambda X^{-1}, every eigenvalue of A+EA + E is within κ2(X)∥E∥2\kappa_2(X)\|E\|_2 of an eigenvalue of AA. Introduced in: Lesson 17
    BBP phase transition
    In the spiked covariance model, the top sample eigenvalue separates from the noise bulk only when the spike strength θ\theta is above γ\sqrt\gamma; below it the spike cannot be detected. Introduced in: Lesson 18
    Bessel's correction
    Dividing by N−1N - 1 instead of NN in the sample covariance. It makes the estimate unbiased, because the mean was fitted to the same data. Introduced in: Lesson 9
    catastrophic cancellation
    The loss of correct digits when two nearly equal rounded numbers are subtracted. The fix is to rewrite the formula. Introduced in: Lesson 13
    Cauchy interlacing
    If B=V⊤SVB = V^\top S V with V∈Rn×mV \in \mathbb{R}^{n \times m} with orthonormal columns (for example a principal submatrix), then λk(S)≥λk(B)≥λk+n−m(S)\lambda_k(S) \ge \lambda_k(B) \ge \lambda_{k+n-m}(S). Introduced in: Lesson 6
    Cauchy–Schwarz inequality
    ∣u⊤v∣≤∥u∥ ∥v∥|u^\top v| \le \|u\|\,\|v\|, with equality only when one vector is a multiple of the other. It makes the angle between two vectors well defined. Introduced in: Lesson 3
    centering
    Subtracting the mean of each feature from the data: Xc=X−1μ⊤X_c = X - \mathbf{1}\mu^\top. PCA needs centered data. Introduced in: Lesson 9
    change-of-basis matrix
    The matrix that turns coordinates in one basis into coordinates in another. If the new basis vectors are the columns of MM, then x=Mcx = Mc and c=M−1xc = M^{-1}x. Introduced in: Lesson 4
    characteristic polynomial
    p(λ)=det⁡(A−λI)p(\lambda) = \det(A - \lambda I), a polynomial of degree nn for an n×nn \times n matrix. Its roots are the eigenvalues. Introduced in: Lesson 5
    Cholesky factorization
    S=LL⊤S = LL^\top with LL lower triangular with a positive diagonal. It exists, and is unique, exactly when SS is positive definite. Introduced in: Lesson 7
    circulant matrix
    A matrix CC with Cjl=cj−lC_{jl} = c_{j-l} (indices mod nn): each column is the previous one shifted down by one. CxCx is the circular convolution c⊛xc \circledast x. Introduced in: Lesson 19
    circular convolution
    (c⊛x)j=∑lcj−l xl(c \circledast x)_j = \sum_l c_{j-l}\,x_l with indices mod nn. The DFT turns it into an entrywise product: c⊛x=F−1(c^⊙Fx)c \circledast x = F^{-1}(\hat c \odot Fx). Introduced in: Lesson 19
    column space فضای ستونی
    C(A)C(A), the span of the columns of AA. It is the set of all outputs AxAx. Introduced in: Lesson 1
    composition ترکیب نگاشت‌ها
    Applying one map after another. ABAB is the matrix of "first BB, then AA": (AB)x=A(Bx)(AB)x = A(Bx). Introduced in: Lesson 1
    condition number
    The worst amplification of a small relative input change. For an invertible matrix, κ(A)=∥A∥ ∥A−1∥\kappa(A) = \|A\|\,\|A^{-1}\|, and κ2(A)=σmax⁡/σmin⁡\kappa_2(A) = \sigma_{\max}/\sigma_{\min}. Introduced in: Lesson 13
    conjugate gradient (CG)
    A Krylov method for Ax=bAx = b with AA symmetric positive definite. Step kk minimizes the AA-norm of the error over x0+Kk(A,r0)x_0 + \mathcal{K}_k(A, r_0); convergence depends on κ\sqrt{\kappa}. Introduced in: Lesson 14
    conjugate transpose
    A∗A^*: the transpose of AA with every entry replaced by its complex conjugate. For a complex vector, ∥x∥2=x∗x\|x\|^2 = x^* x. Introduced in: Lesson 4
    coordinates مختصات
    The numbers cjc_j in x=∑jcjvjx = \sum_j c_j v_j for a basis v1,…,vkv_1, \dots, v_k. Introduced in: Lesson 1
    Courant–Fischer theorem
    The min–max form of the eigenvalues of a symmetric matrix, sorted downward: λk=max⁡dim⁡U=k min⁡x∈U, x≠0R(x)\lambda_k = \max_{\dim U = k}\, \min_{x \in U,\, x \ne 0} R(x). Introduced in: Lesson 6
    CP decomposition
    A tensor written as a sum of RR rank-1 tensors, X≈∑rar∘br∘cr\mathcal{X} \approx \sum_r a_r \circ b_r \circ c_r. The smallest exact RR is the CP rank. Introduced in: Lesson 15
    Davis–Kahan theorem
    Bounds the rotation of an eigenvector or eigenspace of a symmetric matrix under noise EE: sin⁡θ≲∥E∥2/gap\sin\theta \lesssim \|E\|_2 / \text{gap}, for example sin⁡θ(q^k,qk)≤2∥E∥2/δk\sin\theta(\hat q_k, q_k) \le 2\|E\|_2/\delta_k. Introduced in: Lesson 17
    defective matrix
    A square matrix with fewer independent eigenvectors than its size, such as the shear (1101)\begin{pmatrix} 1 & 1 \\ 0 & 1 \end{pmatrix}. It is not diagonalizable. Introduced in: Lesson 5
    deflation
    Setting a tiny subdiagonal entry to 0 in the QR algorithm, which splits off a converged eigenvalue so the work continues on a smaller block. Introduced in: Lesson 14
    determinant
    The factor det⁡A\det A by which a square matrix changes volumes, with a minus sign when the map mirrors space. det⁡A=0\det A = 0 exactly when AA is not invertible. Introduced in: Lesson 3
    diagonalization
    Writing A=XΛX−1A = X\Lambda X^{-1} with eigenvectors as the columns of XX and eigenvalues on the diagonal of Λ\Lambda. It is possible exactly when AA has nn independent eigenvectors. Introduced in: Lesson 5
    differential
    The linear part dfdf of the change f(X+dX)−f(X)f(X + dX) - f(X). If df=tr⁡(G⊤dX)df = \operatorname{tr}(G^\top dX) for every dXdX, then ∇Xf=G\nabla_X f = G. Introduced in: Lesson 16
    dimension بُعد
    The number of vectors in a basis of a subspace. All bases of a subspace have the same size. Introduced in: Lesson 1
    dot product
    u⊤v=u1v1+⋯+unvn=∥u∥ ∥v∥cos⁡θu^\top v = u_1 v_1 + \dots + u_n v_n = \|u\|\,\|v\|\cos\theta. Its sign tells whether the angle between uu and vv is below, at or above 90°90°. Introduced in: Lesson 0
    dual norm
    ∥y∥∗=max⁡∥x∥≤1x⊤y\|y\|_* = \max_{\|x\| \le 1} x^\top y. The dual of ℓp\ell_p is ℓq\ell_q with 1/p+1/q=11/p + 1/q = 1. Introduced in: Lesson 11
    echelon form
    The staircase form that elimination produces. The first nonzero entry (the pivot) of each nonzero row lies to the right of the pivot in the row above. Introduced in: Lesson 2
    Eckart–Young–Mirsky theorem
    The truncated SVD AkA_k is a best approximation of AA of rank at most kk in the spectral and the Frobenius norm. The errors are σk+1\sigma_{k+1} and σk+12+⋯+σr2\sqrt{\sigma_{k+1}^2 + \dots + \sigma_r^2}. Introduced in: Lesson 8
    eigengap
    The distance from an eigenvalue (or a group of eigenvalues) to the rest of the spectrum. A small gap makes the eigenvector sensitive to noise. Introduced in: Lesson 17
    eigenspace
    Eλ=N(A−λI)E_\lambda = N(A - \lambda I): all eigenvectors for the eigenvalue λ\lambda, together with 00. Introduced in: Lesson 5
    eigenvalue
    A number λ\lambda with Av=λvAv = \lambda v for some nonzero vector vv. The eigenvalues are the roots of det⁡(A−λI)\det(A - \lambda I). Introduced in: Lesson 5
    eigenvector
    A nonzero vector vv with Av=λvAv = \lambda v. The matrix only stretches, shrinks or flips it; its line does not change. Introduced in: Lesson 5
    elimination حذف
    Subtracting multiples of rows from other rows to create zeros below each pivot (Gaussian elimination). Introduced in: Lesson 1
    entry درایه
    One number of a matrix; AijA_{ij} is the entry in row ii and column jj. Introduced in: Lesson 1
    equivariant map
    A map ff with f(ρ(g)x)=ρ′(g)f(x)f(\rho(g)x) = \rho'(g)f(x) for every group element gg: transforming the input transforms the output in a matching way. It is invariant if ρ′(g)=I\rho'(g) = I. Introduced in: Lesson 19
    Euclidean norm
    The length of a vector, ∥v∥=∥v∥2=v⊤v=v12+⋯+vn2\|v\| = \|v\|_2 = \sqrt{v^\top v} = \sqrt{v_1^2 + \dots + v_n^2}. Introduced in: Lesson 0
    explained variance ratio
    The share of the total variance along principal direction ii: λi/∑jλj=σi2/∥Xc∥F2\lambda_i / \sum_j \lambda_j = \sigma_i^2 / \|X_c\|_F^2. Introduced in: Lesson 9
    Fiedler vector
    The eigenvector of the second smallest Laplacian eigenvalue λ2\lambda_2. Its signs split the graph in two (spectral bisection). Introduced in: Lesson 19
    forward error
    ∥x^−x∥/∥x∥\|\hat x - x\| / \|x\|: how wrong the computed answer is. It is usually unknown, because xx is unknown. Introduced in: Lesson 13
    forward mode
    Automatic differentiation that multiplies the chain of Jacobians from the input side (JVPs). One pass per input direction. Introduced in: Lesson 15
    four fundamental subspaces
    The column space C(A)C(A), null space N(A)N(A), row space C(A⊤)C(A^\top) and left null space N(A⊤)N(A^\top) of a matrix AA. Introduced in: Lesson 2
    Fourier neural operator (FNO)
    A network whose layers apply a learned filter to the lowest Fourier modes plus a pointwise linear map: a global convolution that can be used on grids of different size. Introduced in: Lesson 19
    free variable
    A variable of Ax=0Ax = 0 whose column has no pivot. You choose it freely; back substitution then gives the pivot variables. Introduced in: Lesson 2
    frequency response
    The eigenvalues c^=Fc\hat c = Fc of a circulant (the DFT of its first column). They tell how the filter scales each Fourier mode. Introduced in: Lesson 19
    Frobenius norm
    ∥A∥F=∑ijaij2=∑iσi2\|A\|_F = \sqrt{\sum_{ij} a_{ij}^2} = \sqrt{\sum_i \sigma_i^2}: the Euclidean norm of AA seen as a long vector. Introduced in: Lesson 11
    Gavish–Donoho threshold
    A cut-off for singular values of a noisy matrix that minimizes the mean squared error; for a square matrix it is (4/3)n σ(4/\sqrt3)\sqrt n\,\sigma, a little above the bulk edge 2n σ2\sqrt n\,\sigma. Introduced in: Lesson 18
    generalized eigenvector
    A nonzero vector vv with (A−λI)kv=0(A - \lambda I)^k v = 0 for some k≥1k \ge 1. A defective matrix needs them to complete a basis. Introduced in: Lesson 5
    generalized Rayleigh quotient
    x⊤Sx/x⊤Bxx^\top S x / x^\top B x with SS symmetric and BB symmetric positive definite. Its critical values are the eigenvalues of Sx=λBxSx = \lambda Bx. Introduced in: Lesson 6
    geometric multiplicity
    dim⁡N(A−λI)\dim N(A - \lambda I), the number of independent eigenvectors for λ\lambda. It is at least 1 and at most the algebraic multiplicity. Introduced in: Lesson 5
    Gershgorin disc theorem
    Every eigenvalue of AA lies in at least one disc with center aiia_{ii} and radius ∑j≠i∣aij∣\sum_{j \ne i} |a_{ij}|. Introduced in: Lesson 5
    Givens rotation
    A rotation that changes only two coordinates. It is chosen to set one entry to zero. Introduced in: Lesson 12
    GMRES
    A Krylov method for any nonsingular AA: it picks xk∈x0+Kk(A,r0)x_k \in x_0 + \mathcal{K}_k(A, r_0) with the smallest residual ∥b−Axk∥2\|b - Ax_k\|_2. Introduced in: Lesson 14
    Golub–Kahan bidiagonalization
    A=UBV⊤A = UBV^\top with BB bidiagonal, computed with Householder reflections from both sides. It is the first phase of computing the SVD. Introduced in: Lesson 14
    gradient check
    A test of a derivative: compare f(x+h d)−f(x)f(x + h\,d) - f(x) with hh times the claimed directional derivative for decreasing hh. The error must fall like h2h^2. Introduced in: Lesson 16
    Gram–Schmidt
    Builds orthonormal vectors from the columns of AA by subtracting projections onto the earlier ones. The modified version (MGS) is more stable than the classical one (CGS). Introduced in: Lesson 12
    graph Fourier transform
    x^=U⊤x\hat x = U^\top x, the coordinates of a graph signal in the eigenvectors of the Laplacian L=UΛU⊤L = U\Lambda U^\top. Small eigenvalues belong to smooth signals. Introduced in: Lesson 19
    graph Laplacian
    L=D−AL = D - A for a graph with adjacency matrix AA and degree matrix DD. It is PSD, x⊤Lx=∑edgesAij(xi−xj)2x^\top Lx = \sum_{\text{edges}} A_{ij}(x_i - x_j)^2, and its number of zero eigenvalues is the number of connected components. Introduced in: Lesson 19
    Grassmann manifold
    Gr⁡(n,p)\operatorname{Gr}(n, p), the set of pp-dimensional subspaces of Rn\mathbb{R}^n. One point can be stored as the projector XX⊤XX^\top for any XX in St⁡(n,p)\operatorname{St}(n, p) that spans it. Introduced in: Lesson 19
    group representation
    A rule ρ\rho that gives each element gg of a group an invertible matrix ρ(g)\rho(g), with ρ(gh)=ρ(g)ρ(h)\rho(gh) = \rho(g)\rho(h). Introduced in: Lesson 19
    Hessenberg form
    A matrix with zeros below the first subdiagonal. An orthogonal similarity A=QHQ⊤A = QHQ^\top (Householder reflections) brings any square matrix to this form, and a symmetric matrix to tridiagonal form. Introduced in: Lesson 14
    Hoffman–Wielandt inequality
    For symmetric AA and A+EA + E with sorted eigenvalues: ∑k(λk(A+E)−λk(A))2≤∥E∥F2\sum_k (\lambda_k(A + E) - \lambda_k(A))^2 \le \|E\|_F^2. Introduced in: Lesson 17
    HOSVD
    Higher-order SVD: in each mode, take the top left singular vectors of the unfolding, then project the tensor onto them. The truncated HOSVD is quasi-optimal, not optimal. Introduced in: Lesson 15
    Householder reflection
    H=I−2vv⊤/(v⊤v)H = I - 2vv^\top / (v^\top v), the reflection across the hyperplane orthogonal to vv. It is symmetric and orthogonal and can map a vector to a multiple of e1e_1. Introduced in: Lesson 12
    image (of a vector) تصویر (یک بردار)
    The output T(x)T(x) of a map for the input xx. The columns of AA are the images of the standard basis vectors. Introduced in: Lesson 1
    indefinite matrix
    A symmetric matrix with both positive and negative eigenvalues. The graph of its quadratic form is a saddle. Introduced in: Lesson 7
    induced norm
    The matrix norm ∥A∥=max⁡x≠0∥Ax∥/∥x∥\|A\| = \max_{x \ne 0} \|Ax\|/\|x\| given by a vector norm. ∥A∥1\|A\|_1 is the largest absolute column sum and ∥A∥∞\|A\|_\infty the largest absolute row sum. Introduced in: Lesson 11
    inverse iteration
    Power iteration with (A−σI)−1(A - \sigma I)^{-1}: solve (A−σI)y=xk(A - \sigma I)y = x_k and normalize. It converges to the eigenvector whose eigenvalue is closest to the shift σ\sigma. Introduced in: Lesson 14
    invertible matrix
    A square matrix AA with a matrix A−1A^{-1} such that A−1A=AA−1=IA^{-1}A = AA^{-1} = I. Equivalently, N(A)={0}N(A) = \{0\}, or det⁡A≠0\det A \ne 0. Introduced in: Lesson 3
    Jacobian
    The m×nm \times n matrix Jij=∂fi/∂xjJ_{ij} = \partial f_i / \partial x_j of a map f:Rn→Rmf: \mathbb{R}^n \to \mathbb{R}^m. Near xx, f(x+δ)≈f(x)+Jδf(x + \delta) \approx f(x) + J\delta. Introduced in: Lesson 15
    Johnson–Lindenstrauss lemma
    A random projection to k=O(ε−2log⁡n)k = O(\varepsilon^{-2}\log n) dimensions keeps all pairwise distances of nn points within a factor 1±ε1 \pm \varepsilon, with high probability. kk does not depend on the original dimension. Introduced in: Lesson 18
    Jordan form
    A block-diagonal matrix of Jordan blocks λI+N\lambda I + N, where NN has ones just above the diagonal. Every square matrix is similar to a Jordan form over C\mathbb{C}. Introduced in: Lesson 5
    JVP (Jacobian-vector product)
    The product JvJv of the Jacobian with a tangent vector vv. Forward mode computes it in one pass. Introduced in: Lesson 15
    Kronecker product
    The block matrix A⊗BA \otimes B with blocks aijBa_{ij}B. Its eigenvalues are the products λi(A)μj(B)\lambda_i(A)\mu_j(B). Introduced in: Lesson 16
    Kronecker sum
    A⊗Im+In⊗BA \otimes I_m + I_n \otimes B for AA of size n×nn \times n and BB of size m×mm \times m. Its eigenvalues are the sums λi(A)+μj(B)\lambda_i(A) + \mu_j(B). Introduced in: Lesson 16
    Krylov subspace
    Kk(A,b)=span⁡{b,Ab,…,Ak−1b}={p(A)b:deg⁡p≤k−1}\mathcal{K}_k(A, b) = \operatorname{span}\{b, Ab, \dots, A^{k-1}b\} = \{p(A)b : \deg p \le k - 1\}. It needs only matrix-vector products. Introduced in: Lesson 14
    Lanczos iteration
    Arnoldi for a symmetric matrix. HkH_k becomes tridiagonal, and each new vector needs only the two previous ones. Introduced in: Lesson 14
    least squares
    Find x^\hat x that minimizes ∥Ax−b∥2\|Ax - b\|_2. Then Ax^A\hat x is the orthogonal projection of bb onto C(A)C(A). Introduced in: Lesson 12
    left null space
    N(A⊤)={y:A⊤y=0}N(A^\top) = \{y : A^\top y = 0\}, a subspace of Rm\mathbb{R}^m for an m×nm \times n matrix AA. It is the orthogonal complement of C(A)C(A). Introduced in: Lesson 2
    linear combination ترکیب خطی
    A sum c1v1+⋯+ckvk=Vcc_1 v_1 + \dots + c_k v_k = Vc of vectors multiplied by numbers. Introduced in: Lesson 1
    linear map نگاشت خطی
    A map TT with T(αx+βy)=αT(x)+βT(y)T(\alpha x + \beta y) = \alpha T(x) + \beta T(y) for all vectors x,yx, y and numbers α,β\alpha, \beta. With a fixed basis, it is the same as a matrix. Introduced in: Lesson 1
    linearly dependent وابسته خطی
    Not independent; some nonzero combination gives 0, so one vector is a combination of the others. Introduced in: Lesson 1
    linearly independent مستقل خطی
    Vectors v1,…,vkv_1, \dots, v_k with c1v1+⋯+ckvk=0c_1 v_1 + \dots + c_k v_k = 0 only for c=0c = 0. Each one adds a new direction. Introduced in: Lesson 1
    Lipschitz constant
    The smallest LL with ∥f(x)−f(y)∥≤L∥x−y∥\|f(x) - f(y)\| \le L\|x - y\| for all x,yx, y. For a linear map x↦Wxx \mapsto Wx in the 2-norm it is σ1(W)\sigma_1(W). Introduced in: Lesson 11
    Loewner order
    A⪯BA \preceq B means that B−AB - A is positive semidefinite: x⊤Ax≤x⊤Bxx^\top A x \le x^\top B x for every xx. It is a partial order; some pairs cannot be compared. Introduced in: Lesson 7
    LoRA
    Low-rank adaptation: fine-tune W0+αrBAW_0 + \tfrac{\alpha}{r}BA with B∈Rd×rB \in \mathbb{R}^{d \times r}, A∈Rr×kA \in \mathbb{R}^{r \times k} trainable and W0W_0 frozen, so only r(d+k)r(d + k) numbers are trained. Introduced in: Lesson 15
    LU factorization
    A=LUA = LU with LL unit lower triangular (it stores the elimination multipliers) and UU upper triangular (it holds the pivots). With row exchanges, PA=LUPA = LU. Introduced in: Lesson 10
    machine epsilon
    The gap ε=21−p\varepsilon = 2^{1-p} between 1 and the next larger float, for pp bits of precision. In double precision ε=2−52≈2.2×10−16\varepsilon = 2^{-52} \approx 2.2 \times 10^{-16}. Introduced in: Lesson 13
    Marchenko–Pastur law
    The limit of the eigenvalue histogram of a sample covariance of pure noise when p/n→γp/n \to \gamma. For γ≤1\gamma \le 1 it fills [σ2(1−γ)2,σ2(1+γ)2][\sigma^2(1 - \sqrt\gamma)^2, \sigma^2(1 + \sqrt\gamma)^2]. Introduced in: Lesson 18
    matrix exponential
    eA=∑k≥0Ak/k!e^A = \sum_{k \ge 0} A^k / k!, defined for every square AA. x(t)=etAx0x(t) = e^{tA}x_0 solves x′=Axx' = Ax. Introduced in: Lesson 16
    matrix factorization
    Writing a matrix as a product of simple matrices (triangular, orthogonal or diagonal), for example LULU, QRQR, LL⊤LL^\top or UΣV⊤U\Sigma V^\top. Introduced in: Lesson 10
    matrix function
    For symmetric S=QΛQ⊤S = Q\Lambda Q^\top, f(S)=Qf(Λ)Q⊤f(S) = Qf(\Lambda)Q^\top, with ff applied to each eigenvalue (not to each entry). Examples are eSe^S, S1/2S^{1/2} and log⁡S\log S. Introduced in: Lesson 16
    mode-nn unfolding
    The matrix X(n)X_{(n)} whose columns are the mode-nn fibers of a tensor (all indices fixed except the nn-th). For a third-order tensor, X(1)X_{(1)} is I×JKI \times JK. Introduced in: Lesson 15
    multilinear rank
    The tuple (rank⁡X(1),rank⁡X(2),rank⁡X(3))(\operatorname{rank} X_{(1)}, \operatorname{rank} X_{(2)}, \operatorname{rank} X_{(3)}) of a tensor. Its entries can differ. Introduced in: Lesson 15
    norm
    A function ∥⋅∥\|\cdot\| with ∥x∥≥0\|x\| \ge 0, ∥x∥=0\|x\| = 0 only for x=0x = 0, ∥cx∥=∣c∣ ∥x∥\|cx\| = |c|\,\|x\| and ∥x+y∥≤∥x∥+∥y∥\|x + y\| \le \|x\| + \|y\|. Introduced in: Lesson 11
    normal equations
    A⊤Ax^=A⊤bA^\top A\hat x = A^\top b. They say that the error b−Ax^b - A\hat x is orthogonal to every column of AA, and they give the least-squares solution. Introduced in: Lesson 3
    normalized Laplacian
    Lsym=I−D−1/2AD−1/2L_{\text{sym}} = I - D^{-1/2}AD^{-1/2}. Its eigenvalues lie in [0,2][0, 2]. Introduced in: Lesson 19
    nuclear norm
    ∥A∥∗=∑iσi\|A\|_* = \sum_i \sigma_i, the sum of the singular values. It is the convex surrogate for rank used in low-rank models. Introduced in: Lesson 11
    null space فضای پوچ
    N(A)={x:Ax=0}N(A) = \{x : Ax = 0\}. The columns of AA are independent exactly when N(A)={0}N(A) = \{0\}. Introduced in: Lesson 1
    numerical rank
    The number of singular values above a cut-off, such as rcond⋅σ1\text{rcond} \cdot \sigma_1. Smaller singular values are treated as zero. Introduced in: Lesson 8
    orthogonal complement
    The set V⊥V^\perp of all vectors that are orthogonal to every vector of a subspace VV. In Rn\mathbb{R}^n, N(A)N(A) and C(A⊤)C(A^\top) are orthogonal complements. Introduced in: Lesson 2
    orthogonal matrix
    A real square matrix QQ with Q⊤Q=IQ^\top Q = I. Its columns are orthonormal, Q−1=Q⊤Q^{-1} = Q^\top, and it keeps lengths and angles. It is a rotation (det⁡Q=1\det Q = 1) or a rotation combined with a reflection (det⁡Q=−1\det Q = -1). Introduced in: Lesson 4
    orthogonal projection
    The point pp of a subspace that is closest to a vector bb. The error b−pb - p is orthogonal to the subspace. Introduced in: Lesson 3
    orthogonal projector
    A matrix PP with P2=PP^2 = P and P⊤=PP^\top = P. For AA with independent columns, P=A(A⊤A)−1A⊤P = A(A^\top A)^{-1}A^\top projects onto C(A)C(A). Introduced in: Lesson 3
    orthogonal vectors
    Two vectors with dot product zero, u⊤v=0u^\top v = 0. In the plane and in space they are perpendicular. Introduced in: Lesson 0
    orthonormal
    Vectors that are pairwise orthogonal and have length 1: qi⊤qj=0q_i^\top q_j = 0 for i≠ji \ne j and qi⊤qi=1q_i^\top q_i = 1. Introduced in: Lesson 4
    outer product ضرب بیرونی
    The matrix uv⊤uv^\top for column vectors u,vu, v. Its column jj is vjuv_j u, so its rank is 1 when u,v≠0u, v \ne 0. Introduced in: Lesson 1
    over-smoothing
    In a deep GCN without residual connections, repeated propagation acts like power iteration, so the node features become almost the same up to a degree-dependent factor. Introduced in: Lesson 19
    oversampling
    Taking a few more random vectors than the target rank in a randomized SVD, so that one bad draw does not matter. Introduced in: Lesson 9
    partial pivoting
    At each elimination step, swap in the row with the largest ∣aik∣|a_{ik}| in the current column. Then all multipliers satisfy ∣ℓij∣≤1|\ell_{ij}| \le 1. Introduced in: Lesson 10
    Perron–Frobenius theorem
    A square matrix with all entries positive has a real, positive, simple eigenvalue that is larger in absolute value than all others. Its eigenvector can be chosen with all entries positive. Introduced in: Lesson 5
    pivot محور
    The first nonzero entry of a row after elimination. A column without a pivot is a free column. Introduced in: Lesson 1
    polar decomposition
    For square AA: A=QPA = QP with QQ orthogonal and PP symmetric positive semidefinite. From the SVD, Q=UV⊤Q = UV^\top and P=VΣV⊤P = V\Sigma V^\top. Introduced in: Lesson 8
    positive definite
    S≻0S \succ 0: a symmetric matrix with x⊤Sx>0x^\top S x > 0 for every x≠0x \ne 0. Equivalently, all its eigenvalues are positive. Introduced in: Lesson 7
    positive semidefinite
    S⪰0S \succeq 0: a symmetric matrix with x⊤Sx≥0x^\top S x \ge 0 for every xx. Equivalently, all its eigenvalues are ≥0\ge 0. Introduced in: Lesson 7
    power iteration
    Repeat xk+1=Axk/∥Axk∥x_{k+1} = Ax_k / \|Ax_k\|. If ∣λ1∣>∣λ2∣|\lambda_1| > |\lambda_2|, almost every start vector turns toward the top eigenvector, with error of order ∣λ2/λ1∣k|\lambda_2/\lambda_1|^k. Introduced in: Lesson 5
    preconditioning
    Solving an equivalent system with a cheap M≈AM \approx A, so that the iteration sees M−1AM^{-1}A, which has a smaller condition number or clustered eigenvalues. Introduced in: Lesson 14
    principal angles
    The angles 0≤θ1≤⋯≤θd≤90°0 \le \theta_1 \le \dots \le \theta_d \le 90° between two dd-dimensional subspaces with orthonormal bases VV and V^\hat V: cos⁡θi=σi(V^⊤V)\cos\theta_i = \sigma_i(\hat V^\top V). Introduced in: Lesson 17
    principal component analysis (PCA)
    Finds the orthogonal directions of largest variance of a data set: the eigenvectors of the sample covariance matrix, or the right singular vectors of the centered data matrix. Introduced in: Lesson 9
    principal directions
    The eigenvectors qiq_i of the sample covariance matrix, equal to the right singular vectors of XcX_c. They are sorted by the variance λi=σi2/(N−1)\lambda_i = \sigma_i^2/(N - 1) along them. Introduced in: Lesson 9
    principal minor
    The determinant of a submatrix with the same set of row and column indices. A leading principal minor uses the indices 1,…,k1, \dots, k. Introduced in: Lesson 7
    pseudoinverse
    The Moore–Penrose pseudoinverse A+=VΣ+U⊤A^+ = V\Sigma^+ U^\top, where Σ+\Sigma^+ inverts each nonzero σi\sigma_i and keeps the zeros. A+bA^+ b is the minimum-norm least-squares solution. Introduced in: Lesson 8
    push-through identity
    X⊤(XX⊤+λI)−1=(X⊤X+λI)−1X⊤X^\top(XX^\top + \lambda I)^{-1} = (X^\top X + \lambda I)^{-1}X^\top. It lets ridge regression invert the smaller of the two matrices. Introduced in: Lesson 10
    QK and OV circuits
    The two low-rank d×dd \times d matrices of an attention head: WQKW_{QK} decides where a token looks, and WOVW_{OV} decides what is read and written. Both have rank at most dhd_h. Introduced in: Lesson 15
    QR algorithm
    Repeat: factor Ak=QkRkA_k = Q_kR_k, then set Ak+1=RkQkA_{k+1} = R_kQ_k. With Hessenberg reduction, shifts and deflation it is the standard dense eigenvalue method. Introduced in: Lesson 14
    QR factorization
    A=QRA = QR with QQ having orthonormal columns and RR upper triangular. The reduced (thin) QR has QQ of size m×nm \times n. Introduced in: Lesson 12
    QR with column pivoting
    AP=QRAP = QR, where at each step the remaining column of largest norm goes first. Small diagonal entries of RR reveal the numerical rank. Introduced in: Lesson 12
    quadratic form
    The function x↦x⊤Sxx \mapsto x^\top S x of a symmetric matrix SS, also called the energy. Introduced in: Lesson 7
    random matrix theory
    The study of the eigenvalues and singular values of large random matrices. It shows how pure noise looks, so you can tell which eigenvalues carry signal. Introduced in: Lesson 18
    randomized SVD
    An approximate truncated SVD: multiply AA by a random matrix, orthonormalize the result to get QQ, then take the small SVD of Q⊤AQ^\top A. Introduced in: Lesson 9
    range finder
    The first stage of the randomized SVD: Q=orth⁡(AΩ)Q = \operatorname{orth}(A\Omega) for a random Ω\Omega, an orthonormal basis for (almost) the top part of the range of AA. Introduced in: Lesson 9
    rank رتبه
    The dimension of the column space of a matrix: the number of independent columns. Introduced in: Lesson 1
    rank–nullity theorem
    For A∈Rm×nA \in \mathbb{R}^{m \times n}: rank⁡(A)+dim⁡N(A)=n\operatorname{rank}(A) + \dim N(A) = n. Introduced in: Lesson 2
    Rayleigh quotient
    R(x)=x⊤Sx/x⊤xR(x) = x^\top S x / x^\top x for x≠0x \ne 0. For symmetric SS it lies between the smallest and the largest eigenvalue and reaches them at the eigenvectors. Introduced in: Lesson 6
    Rayleigh quotient iteration
    Inverse iteration with the shift updated to the Rayleigh quotient at every step. For symmetric matrices it converges cubically. Introduced in: Lesson 14
    reduced row echelon form (RREF)
    An echelon form in which every pivot is 1 and every other entry in a pivot column is 0. Each matrix has exactly one RREF. Introduced in: Lesson 2
    residual
    r=b−Ax^r = b - A\hat x for an approximate solution x^\hat x. In least squares it is orthogonal to every column of AA. Introduced in: Lesson 12
    retraction
    A smooth map that brings a point X+VX + V back onto the manifold after a step, for example the Q factor of a QR factorization or the polar factor. Introduced in: Lesson 19
    reverse mode
    Automatic differentiation that multiplies the chain of Jacobians from the output side (VJPs). One pass gives the full gradient of a scalar loss; backpropagation is reverse mode. Introduced in: Lesson 15
    Riemannian gradient
    The Euclidean gradient projected onto the tangent space of the manifold. On the Stiefel manifold it is G−Xsym⁡(X⊤G)G - X\operatorname{sym}(X^\top G). Introduced in: Lesson 19
    Ritz values
    The eigenvalues of Hk=Qk⊤AQkH_k = Q_k^\top AQ_k, the matrix compressed onto a Krylov subspace. They approximate the extreme eigenvalues of AA first. Introduced in: Lesson 14
    row space فضای سطری
    The span of the rows of AA. Every row of ABAB lies in the row space of BB. Introduced in: Lesson 1
    sample covariance matrix
    S=Xc⊤Xc/(N−1)S = X_c^\top X_c / (N - 1) for NN centered samples in the rows of XcX_c. It is symmetric positive semidefinite. Introduced in: Lesson 9
    Schatten pp-norm
    The ℓp\ell_p norm of the vector of singular values. p=1p = 1, 22 and ∞\infty give the nuclear, Frobenius and spectral norms. Introduced in: Lesson 11
    Schur complement
    For M=(ABB⊤C)M = \begin{pmatrix} A & B \\ B^\top & C \end{pmatrix} with AA invertible: M/A=C−B⊤A−1BM/A = C - B^\top A^{-1} B. Then M≻0M \succ 0 exactly when A≻0A \succ 0 and M/A≻0M/A \succ 0. Introduced in: Lesson 7
    scores
    The coordinates of the samples on the principal directions: XcV=UΣX_c V = U\Sigma. Introduced in: Lesson 9
    shear برش
    A map that keeps one axis fixed and slides the points parallel to it, e.g. (x,y)↦(x+y,y)(x, y) \mapsto (x + y, y). Introduced in: Lesson 1
    Sherman–Morrison formula
    The rank-1 case of the Woodbury formula: (A+uv⊤)−1=A−1−A−1uv⊤A−11+v⊤A−1u(A + uv^\top)^{-1} = A^{-1} - \frac{A^{-1}uv^\top A^{-1}}{1 + v^\top A^{-1}u}, if 1+v⊤A−1u≠01 + v^\top A^{-1}u \ne 0. Introduced in: Lesson 10
    significand
    The number (1.f)2(1.f)_2 that a float multiplies by a power of 2. Its leading 1 is not stored (the hidden bit). Introduced in: Lesson 13
    similar matrices
    AA and B=M−1AMB = M^{-1}AM for an invertible MM: one linear map written in two bases. They have the same eigenvalues, trace and determinant, but in general not the same singular values. Introduced in: Lesson 4
    singular value
    A diagonal entry σi≥0\sigma_i \ge 0 of Σ\Sigma in the SVD. The σi2\sigma_i^2 are the eigenvalues of A⊤AA^\top A, and σ1\sigma_1 is the largest stretch of AA. Introduced in: Lesson 8
    singular value decomposition (SVD)
    A=UΣV⊤A = U\Sigma V^\top with UU and VV orthogonal and Σ\Sigma diagonal with σ1≥σ2≥⋯≥0\sigma_1 \ge \sigma_2 \ge \dots \ge 0. Every real matrix has one. Introduced in: Lesson 8
    singular vectors
    The columns uiu_i of UU (left singular vectors) and viv_i of VV (right singular vectors) in the SVD, with Avi=σiuiAv_i = \sigma_i u_i and A⊤ui=σiviA^\top u_i = \sigma_i v_i. Introduced in: Lesson 8
    sketching
    Multiplying the data by a short random matrix SS and working with the small result SASA, for example to solve a huge least-squares problem. Introduced in: Lesson 18
    soft-thresholding
    Sλ(z)=sign⁡(z)max⁡(∣z∣−λ,0)S_\lambda(z) = \operatorname{sign}(z)\max(|z| - \lambda, 0), the minimizer of 12(x−z)2+λ∣x∣\tfrac12(x - z)^2 + \lambda|x|. It sets small values exactly to zero, as in the lasso. Introduced in: Lesson 11
    span فضای پوشش
    The set of all linear combinations of some vectors; for the columns of VV it is the column space C(V)C(V). Introduced in: Lesson 1
    spanning set مجموعه‌ی پوشا
    A set of vectors whose span is the whole subspace. It may contain extra (dependent) vectors. Introduced in: Lesson 1
    special solution
    The solution of Ax=0Ax = 0 with one free variable equal to 1 and the other free variables equal to 0. The special solutions form a basis of N(A)N(A). Introduced in: Lesson 2
    spectral norm
    ∥A∥2=max⁡x≠0∥Ax∥/∥x∥\|A\|_2 = \max_{x \ne 0} \|Ax\|/\|x\|, the largest stretch of AA. It equals the largest singular value σ1\sigma_1; for symmetric AA it is max⁡i∣λi∣\max_i |\lambda_i|. Introduced in: Lesson 6
    spectral normalization
    Dividing each weight matrix by an estimate of its largest singular value, so that every layer is (about) 1-Lipschitz. Introduced in: Lesson 11
    spectral theorem
    Every real symmetric matrix is S=QΛQ⊤S = Q\Lambda Q^\top with QQ orthogonal and Λ\Lambda real and diagonal: real eigenvalues and an orthonormal basis of eigenvectors. Introduced in: Lesson 6
    standard basis vector بردار پایه‌ی استاندارد
    The vector eje_j with a 1 in position jj and 0 elsewhere. AejAe_j is column jj of AA. Introduced in: Lesson 1
    standardization
    Dividing each centered feature by its standard deviation, so that all features have variance 1. PCA of standardized data is PCA of the correlation matrix. Introduced in: Lesson 9
    Stiefel manifold
    St⁡(n,p)={X∈Rn×p:X⊤X=Ip}\operatorname{St}(n, p) = \{X \in \mathbb{R}^{n \times p} : X^\top X = I_p\}, the set of matrices with orthonormal columns. Its dimension is np−p(p+1)/2np - p(p+1)/2. Introduced in: Lesson 19
    subspace زیرفضا
    A set of vectors that contains 0 and is closed under addition and multiplication by numbers. Introduced in: Lesson 1
    subspace embedding
    A matrix SS with (1−ε)∥y∥2≤∥Sy∥2≤(1+ε)∥y∥2(1 - \varepsilon)\|y\|^2 \le \|Sy\|^2 \le (1 + \varepsilon)\|y\|^2 for all yy in a subspace. It is oblivious if it works for any fixed subspace without looking at the data. Introduced in: Lesson 18
    subspace iteration
    Power iteration on a block of vectors, with an orthonormalization after every product. In the randomized SVD it turns the weights σj\sigma_j into σj2q+1\sigma_j^{2q+1}. Introduced in: Lesson 9
    superposition
    A network represents more features than it has dimensions by giving each feature an almost-orthogonal direction. It works when features are sparse. Introduced in: Lesson 18
    swamping
    The loss of a small term tt added to a large sum ss, when ∣t∣|t| is less than half the gap between floats near ss. Introduced in: Lesson 13
    Sylvester equation
    AX+XB=CAX + XB = C, linear in XX. It has a unique solution exactly when no λi(A)+μj(B)\lambda_i(A) + \mu_j(B) is zero. Introduced in: Lesson 16
    symmetric matrix
    A square matrix with S⊤=SS^\top = S, that is, sij=sjis_{ij} = s_{ji}. Introduced in: Lesson 6
    tensor
    An array with three or more indices, such as X∈RI×J×K\mathcal{X} \in \mathbb{R}^{I \times J \times K} with entries xijkx_{ijk}. Introduced in: Lesson 15
    thin SVD
    The reduced SVD of an m×nm \times n matrix: it keeps only p=min⁡(m,n)p = \min(m, n) columns of UU and VV and the p×pp \times p block of Σ\Sigma. Introduced in: Lesson 8
    truncated SVD
    Ak=∑i≤kσiuivi⊤A_k = \sum_{i \le k} \sigma_i u_i v_i^\top: the SVD with only the kk largest terms. Introduced in: Lesson 8
    Tucker decomposition
    X≈G×1A×2B×3C\mathcal{X} \approx \mathcal{G} \times_1 A \times_2 B \times_3 C: a small core tensor G\mathcal{G} multiplied by a factor matrix in each mode. Introduced in: Lesson 15
    underflow
    A result too small for the float format: it loses precision or becomes 0. Small fp16 gradients underflow without loss scaling. Introduced in: Lesson 13
    unit ball
    The set {x:∥x∥≤1}\{x : \|x\| \le 1\} of a norm: a disk for ℓ2\ell_2, a diamond for ℓ1\ell_1 and a square for ℓ∞\ell_\infty in the plane. Introduced in: Lesson 11
    unit lower triangular
    A lower triangular matrix with ones on the diagonal. Introduced in: Lesson 10
    unit roundoff
    u=ε/2u = \varepsilon/2, the largest relative error of rounding to nearest: fl⁡(x)=x(1+δ)\operatorname{fl}(x) = x(1 + \delta) with ∣δ∣≤u|\delta| \le u. Introduced in: Lesson 13
    unitary matrix
    A complex square matrix UU with U∗U=IU^* U = I, the complex version of an orthogonal matrix. Every eigenvalue has ∣λ∣=1|\lambda| = 1. Introduced in: Lesson 4
    vec
    Stacks the columns of a matrix into one long vector. Then vec⁡(AXB)=(B⊤⊗A)vec⁡(X)\operatorname{vec}(AXB) = (B^\top \otimes A)\operatorname{vec}(X). Introduced in: Lesson 16
    VJP (vector-Jacobian product)
    The product J⊤uJ^\top u (the row vector u⊤Ju^\top J). Reverse mode computes it in one pass. Introduced in: Lesson 15
    Wedin's theorem
    The Davis–Kahan bound for singular vectors of any matrix: the top-kk singular subspaces rotate by at most about ∥E∥2\|E\|_2 divided by the gap σk−σk+1\sigma_k - \sigma_{k+1}. Introduced in: Lesson 17
    Weyl's inequality
    For symmetric AA and EE, each sorted eigenvalue moves by at most the spectral norm of the change: ∣λk(A+E)−λk(A)∣≤∥E∥2|\lambda_k(A + E) - \lambda_k(A)| \le \|E\|_2. Introduced in: Lesson 6
    Wilkinson shift
    The eigenvalue of the trailing 2×22 \times 2 block that is closer to the last diagonal entry. The QR algorithm uses it as the shift. Introduced in: Lesson 14
    Woodbury formula
    (A+UCV)−1=A−1−A−1U(C−1+VA−1U)−1VA−1(A + UCV)^{-1} = A^{-1} - A^{-1}U(C^{-1} + VA^{-1}U)^{-1}VA^{-1}. It updates an inverse after a rank-kk change at the cost of a k×kk \times k inverse. Introduced in: Lesson 10
    ℓp\ell_p norm
    ∥x∥p=(∑i∣xi∣p)1/p\|x\|_p = (\sum_i |x_i|^p)^{1/p} for p≥1p \ge 1, and ∥x∥∞=max⁡i∣xi∣\|x\|_\infty = \max_i |x_i|. For p<1p < 1 it is not a norm. Introduced in: Lesson 11