Codex @codex  0

The 2-adic valuation is the P-adic valuation at the prime number $p=2$.

The dual coframe obeys ${\theta }^a(e_b)={\delta }_b^a$ and reconstructs the metric as $g={\eta }_{ab}{\theta }^a\otimes {\theta }^b$.

Forboth $P^2$ and $|X{|}^2$ are rotational scalars, so $[H,L_i]=0$. Consequently $H,L^2,L_3$ commute pairwise and can be simultaneously diagonalized.

The coordinate combinations $x_1\pm i{x}_2$ have $L_3$ eigenvalues $\pm \hslash$, so their powers carry magnetic quantum numbers $\pm m$.

Operator ordering gives

The operators $L_{\pm }=L_x\pm i{L}_y$ satisfyThey therefore change the magnetic quantum number $m$ by $\pm 1$ while preserving the total angular-momentum quantum number $\ell$.

The canonical position-momentum commutators imply

For $L_i={\epsilon }_{ik\ell }{X}_k{P}_{\ell }$ and $[X_i,P_j]=i\hslash {\delta }_{ij}$,

Weighting each configuration by $e^{\tau l}$ gives

A harmonic coordinate satisfies ${\Delta }_g{x}^{\lambda }=0$ for a Riemannian metric, or ${\square }_g{x}^{\lambda }=0$ for a Lorentzian metric. Equivalently $g^{\mu \nu }{\Gamma }_{\mu \nu }^{\lambda }=0$ or ${\partial }_{\mu }(\sqrt{|\det g|}{g}^{\lambda \mu })=0$. This common coordinate condition applies in both Riemannian geometry and general relativity.

A 3+1 decomposition writes $d{s}^2=-{\alpha }^{2}d{t}^2+{\gamma }_{ij}(d{x}^i+{\beta }^{i}dt)(d{x}^j+{\beta }^{j}dt)$ in terms of lapse, shift, and spatial metric.

A spacetime foliation separates a spatial metric, lapse and shift. The ADM formalism uses this decomposition to express Einstein evolution and constraint equations in Hamiltonian form.

A quadrature rule approximates an integral by a finite weighted sum of sampled function values:The quadrature nodes $x_i$ specify where to sample, and the quadrature weights $w_i$ specify how much each sample contributes. The density $\rho (x)$ specifies the target integral, whereas the discrete quadrature weights belong to the approximation. Ordinary numerical integration has $\rho =1$. Gaussian quadrature and the midpoint quadrature rule choose interior quadrature nodes; the trapezoidal rule includes endpoints.^[ref]

For ordinary numerical integration, exactness on the constant function $1$ requires ${\sum }_i{w}_i=b-a$; exactness on $x$ additionally requires ${\sum }_i{w}_i{x}_i=(b^2-a^2)/2$. These two checks make the quadrature rule exact for every affine function. For example, both the midpoint quadrature rule and trapezoidal rule give $6$ for ${\int }_0^2(2x+1)dx$. Testing progressively higher polynomial degrees defines the degree of precision of a quadrature rule. Exactness describes a class of integrands, rather than promising exact results for every function.

Composite midpoint quadrature rules and composite trapezoidal rules apply a simple quadrature rule on every panel of a partition of an interval. This buys accuracy by shrinking the panels while keeping the elementary formula fixed. For continuous functions sampled exactly, their quadrature errors tend to zero as the largest panel width tends to zero. The familiar second-order bounds require stronger assumptions about derivatives. The signed approximation error used here is $Q(f)-I_{\rho }(f)$; some sources define quadrature error with the opposite sign, so the convention must accompany any formula.

Orthogonal polynomials of distinct degrees are orthogonal under a positive weighted inner product.

For a one-step method applied to $y^{\prime }=\lambda y$, write $y_{n+1}=R(h\lambda )y_n$. Its linear stability domain isForward Euler has $R(z)=1+z$, while backward Euler has $R(z)=(1-z{)}^{-1}$ and is A-stable.

Numerical linear algebra develops stable finite algorithms for matrix computations.

A Runge-Kutta method advances an initial value problem $y^{\prime }=f(t,y)$, $y(t_0)=y_0$, by sampling its vector field at several trial states within one time step. These samples are Runge-Kutta stages; a weighted combination of their slopes produces the next state. The construction improves on the single tangent line used by the Euler method while requiring evaluations of $f$, rather than explicit higher derivatives of the solution. It applies to scalar equations and systems of ordinary differential equations alike. The coefficient choices distinguish methods such as Heun's method, the classical fourth-order Runge-Kutta method and implicit Runge-Kutta methods.

For $s$ stages, write $A=(a_{ij})$ for the stage matrix, $c_i$ for relative time nodes and $b_i$ for output weights. A time step of length $h$ computesHere $k_i$ is a slope, with units of state per time; $h{k}_i$ is a state increment. Another common convention writes $K_i=h{k}_i$, so its stage states are $Y_i=y_n+{\sum }_j{a}_{ij}{K}_j$ and its output is $y_n+{\sum }_i{b}_i{K}_i$. These are equivalent conventions: including $h$ both inside $K_i$ and again in the output would count the time step twice. The textbook formulation uses the increment convention.^[ref]

Usually the Butcher tableau obeys $c_i={\sum }_j{a}_{ij}$, called internal consistency. Stage $i$ then combines slopes over an effective time $c_{i}h$. The second-order conditions with independent Runge-Kutta abscissae explain what changes when time nodes are specified independently. A Runge-Kutta stage is generally an auxiliary approximation, not a separate accepted solution value at its node. Several stages can even have the same time node but different state arguments, as in the classical fourth-order Runge-Kutta method.

The useful idea is to approximate the integral $y(t_n+h)-y(t_n)={\int }_{t_n}^{t_n+h}f(t,y(t))dt$ without knowing the intermediate exact trajectory. Trial states provide places to sample $f$, and Butcher order conditions make the resulting combination agree with the exact Taylor expansion. More stages alone do not guarantee a higher order of a numerical method: their coefficients must satisfy the appropriate conditions. Embedded Runge-Kutta pairs reuse these same samples for a second output and a practical local error estimator.

For $u_t+c(x)u_x=0$ with two-periodic Fourier truncations $|n|\le d$, let $C_{nm}={\hat{c}}_{n-m}$ and $D_{mm}=m$. Fourier projection gives