Topics (384k) Articles (398k) Users (410) Discussions (242) Comments (383) Files (13k) New article
A network circulation is a flow with zero net supply at every vertex. A feasible residual network circulation changes one feasible flow into another with the same prescribed supplies. It decomposes into directed cycles, giving the negative-cycle optimality test for minimum-cost flow.
Catastrophic cancellation occurs when subtraction of nearly equal approximations removes their matching leading digits and exposes earlier roundoff error. The subtraction can itself be exact on the stored operands: the information may already have been lost while obtaining them. Goldberg distinguishes this from cancellation of accurately known quantities, which can be harmless.[ref]
https://docs.oracle.com/cd/E19957-01/806-3568/ncg_goldberg.html. Goldberg, 1991, §Cancellation
For , consider the single real-valued functionThe rationalization of a square-root difference proves equality in exact arithmetic. The two evaluations use the same input , but the direct version first forms a value near one and then subtracts one. The rearranged version keeps the small input in its numerator and divides by a denominator near two.
Direct and rationalized evaluation of a small square-root difference against a high-precision reference for the same binary64 inputs
. Details. For binary64
x = 1e-16, the rounded sum 1.0 + x is one, so sqrt(1.0 + x) - 1.0 is zero. The rationalized evaluation returns approximately . A 100-digit Decimal calculation using the exact stored input gives relative errors of one and about respectively. At x = 1e-12, the direct result is about , with relative error ; the rearranged result is about , with relative error .The plot uses
Decimal.from_float(x) to preserve the binary input exactly, computes at 100 decimal digits, and verifies the reference at 160 digits. Exact conversion of a stored binary value is documented by Python.[ref]https://docs.python.org/3.14/library/decimal.html#decimal.Decimal.from_float. Python 3.14 documentation, Decimal.from_float
Each error is , with the computed result converted exactly to Decimal before comparison. Hence the curves measure evaluation error for one common input, rather than a mixture of input conversion error and evaluation error. The plotted interval is far above the subnormal number range.The original problem is well conditioned with respect to relative changes in . Using the derivative and rationalization of a square-root difference, its relative condition number isIt tends to one as tends to zero. The direct evaluation instead introduces the intermediate difference with ; an error of order unit roundoff in is large relative to . This is a failure of the chosen algorithm's accuracy despite a small condition number for the same real-valued function. The rationalized evaluation avoids that intermediate catastrophic cancellation, retaining relative accuracy near the unit roundoff scale throughout the plotted interval. Recognizing this structure can improve a calculation far more than merely printing extra digits.
A subnormal number is a nonzero representable value smaller in magnitude than the least positive normalized value . Keep the exponent scale fixed at and allow a zero leading digit: the positive values are for integers . This fills the gap toward zero with constant absolute spacing, at the cost of fewer significant digits. Goldberg describes this gradual approach to zero and why it matters for differences of small numbers.[ref]
https://docs.oracle.com/cd/E19957-01/806-3568/ncg_goldberg.html. Goldberg, 1991, §Denormalized Numbers
For binary64 the least positive normalized value is , and subnormal numbers are spaced by . With nearest rounding and gradual underflow enabled, an exact result with has absolute rounding error at most . A uniform relative bound by unit roundoff fails here: rounds to zero under the even-last-digit tie rule, giving relative error one. Machine epsilon measures resolution near one; it does not measure the distance from zero to the first positive value.
For floating-point arithmetic with precision in radix and rounding to nearest, unit roundoff is the half-gap scale . If a nonzero exact real number lies in the normal finite range, its nearest rounded value satisfiesIndeed, the absolute distance to the nearest value is at most , and gives the bound. For a correctly rounded addition, subtraction, multiplication or division, apply this estimate to the exact result of that one operation on the stored operands. The bound does not include earlier roundoff error or automatically bound a whole algorithm. Results outside the finite range and tiny subnormal numbers require separate treatment. Rounding in a directed mode also changes the bound.[ref]
https://docs.oracle.com/cd/E19957-01/806-3568/ncg_goldberg.html. Goldberg, 1991, §§Relative Error and Ulps, Exactly Rounded Operations
Here machine epsilon means the gap from one to the next larger value in normalized floating-point arithmetic:It is a relative resolution scale near one, not the smallest positive representable value and not a fixed absolute error for all magnitudes. In binary64, giveswhere is unit roundoff for rounding to nearest. Some authors instead call their machine epsilon; MIT 18.335J uses this half-gap convention. Always check which definition an error bound uses.[ref]
https://ocw.mit.edu/courses/18-335j-introduction-to-numerical-methods-spring-2019/2f313023ae3404bc217a81a31b227170_MIT18_335JS19_lec2.pdf. Persson, MIT 18.335J, Lecture 2, slides 3, 8–9: the half-gap convention
The gap below one is smaller: the preceding binary64 value is , whereas the next value above one is . Under nearest-even rounding,
1.0 + 2.0**-53 is exactly 1.0: the exact sum is a tie and one has the even final binary digit. By contrast, 1.0 + 2.0**-52 is the next value. Thus searching for the smallest increment that changes one by trial addition also depends on the rounding convention. At the next larger binary64 value is , illustrating why a small absolute increment can disappear at large magnitude.The predator-prey first integral changes by the displayed amount when prey are removed. At a fixed prey threshold this increment is constant; sufficiently high thresholds reduce oscillation energy until the threshold is no longer crossed. At a fixed predator threshold, the increment depends on the prey abundance and on the crossing direction. A rising-predator balanced-quota cycle has negative return multiplier of magnitude greater than one; its leading small-oscillation approximation is marginal rather than a robust attractor.
Automorphisms and generate a dihedral group of order . The invariant has rational map degree , so . The fixed-field theorem for a finite automorphism group gives the same degree over its fixed field, proving equality. For odd prime , the nontrivial proper subgroups are the rotation group and the reflection groups. Their fixed fields are and for reflections , respectively.
For potential energy of two-layer entrainment with constant , the fixed-energy turbulent mixed layer assumption gives . If and the interfacial power model for entrainment applies with fixed coefficients, integration gives the displayed depth law with . The relation applies only before the layer reaches the tank bottom. Its interfacial Richardson number stays constant for fixed interface thickness, whereas it decreases like in the constant-speed entrainment model.
The four commuting independent stabilizer generators define a two-dimensional stabilizer code. The identity and fifteen single-site Pauli errors have sixteen distinct four-bit error syndromes, so their two-dimensional images form an orthogonal direct sum of the physical space. This proves nondegenerate single-site correction and the perfect quantum error-correcting code property. The operators and commute with the stabilizer and anticommute with each other. Multiplying by gives, up to phase, , a nontrivial weight-three logical operator, proving distance exactly three.
For , an orthonormal-frame curvature calculation gives . Expansion gives the displayed invariant, diverging at zero radius. It is finite at every positive root of V. Thus horizon coordinate singularities are distinct from the central curvature singularity.
When an experiment's configuration is observed, inference conditions on the experiment actually performed rather than averaging over its possible configurations. Conditioning on an ancillary statistic is a central application; freedom from a nuisance parameter must still be checked in the resulting conditional distribution.
For first-order logic, every sentence true in every model of a theory is derivable from that theory. Equivalently, every consistent theory has a model.
The fragment permits at most distinct variable names, but allows bound names to be reused. It can contain formulas of arbitrarily large quantifier depth. On a linear order, alternating two variable names can express an increasing chain of any prescribed finite length.
Let be a commutative algebra over a field , with a homomorphism restricting to the identity on , and let be its kernel. If is generated by finitely many elements, subtracting their images in gives finitely many generators in . Every polynomial in those generators, modulo , is a linear combination of their classes. Thus , the algebraic cotangent space at this point, is finite-dimensional. For and , the classes , , form an infinite basis modulo . This is a subalgebra of a polynomial ring which is not a finitely generated algebra.
Consider the advection equation with constant speed , no sources, and periodic boundary conditions on . Use equal cells of width , time step , and indices modulo . The exact solution is the periodic translation . At a face, the Riemann problem transports the left state to the right, so the two-state upwind numerical flux is , with .[ref]
Riemann Problems and Jupyter Solutions by David I. Ketcheson, Randall J. LeVeque, and Mauricio J. del Razo. Advection, equations adv:advection and adv:solution, §The Riemann problem; Approximate solvers, §High-resolution methods: reduction to upwind advection
The finite volume method with this numerical flux and an explicit Euler method step givesHere represents a cell average; the recurrence also occurs in the upwind finite difference scheme. Impose the specific CFL condition . Each updated cell average is then a convex combination of two old ones. This proves preservation of their initial minimum and maximum, including nonnegativity for nonnegative data, and contraction of the maximum norm of the difference of two solutions on this periodic grid. These statements concern this constant-speed scalar update with this time step.
For an original synthetic example, take , , , , hence Courant number . Choose a function constant on each successive cell, with the following values, so these are exact initial cell averages. One step givesFor example, cell zero receives from cell seven across the periodic face: . All eight face transfers cancel in the total by shared face flux cancellation. More generally this periodic recurrence conserves the total algebraically for any , even though the CFL condition is needed for the convex combination argument. With the same initial data and , one step gives and , despite preserving the total.
Original periodic upwind profile at zero, one and four steps: the integral remains three quarters while jumps spread
. Details. After four steps, , the exact translated profile has moved two cells to the right. Repeated averaging instead gives the explicit binomial coefficient formulawith all indices modulo eight. The diagram compares these computed cell averages with the exact translated cell averages. Their integrals agree, but the numerical jumps spread across cells. The finite volume method is useful for jump discontinuities because it retains the integral balance even at a jump; conservation alone does not guarantee a sharp profile. At this particular update is an exact one-cell shift, whereas for it mixes neighboring states. No accuracy conclusion at a jump follows just from the conserved total.
The Jacobian matrix of a deformation map carries an infinitesimal reference line element to its current image. Its determinant is the local volume ratio. An orientation-preserving deformation has , and incompressibility imposes .
A stretch tensor isolates local changes of material length from the rotation contained in a deformation gradient. For an invertible orientation-preserving , the polar decomposition in continuum mechanics iswhere is a rotation matrix and are positive-definite symmetric matrices. The right stretch tensor acts on directions in the reference configuration; the left stretch tensor acts on directions in the current configuration.
The eigenvalues of either stretch tensor are the principal stretches. A value two means a material line in the corresponding principal direction doubles its length, while a value one means unchanged length. For a pure rotation, both stretch tensors are , even though need not equal . This makes them useful when an experiment includes substantial specimen rotation that should not be mistaken for strain.
Simple shear displaces one coordinate direction in proportion to a transverse coordinate while preserving volume. Its deformation gradient has a single off-diagonal component . Simple shear flow describes the corresponding velocity field rather than the finite deformation.
Pinned article: Introduction to the OurBigBook Project
Welcome to the OurBigBook Project! Our goal is to create the perfect publishing platform for STEM subjects, and get university-level students to write the best free STEM tutorials ever.
Everyone is welcome to create an account and play with the site: ourbigbook.com/go/register. We belive that students themselves can write amazing tutorials, but teachers are welcome too. You can write about anything you want, it doesn't have to be STEM or even educational. Silly test content is very welcome and you won't be penalized in any way. Just keep it legal!
Intro to OurBigBook
. Source. We have two killer features:
- topics: topics group articles by different users with the same title, e.g. here is the topic for the "Fundamental Theorem of Calculus" ourbigbook.com/go/topic/fundamental-theorem-of-calculusArticles of different users are sorted by upvote within each article page. This feature is a bit like:
- a Wikipedia where each user can have their own version of each article
- a Q&A website like Stack Overflow, where multiple people can give their views on a given topic, and the best ones are sorted by upvote. Except you don't need to wait for someone to ask first, and any topic goes, no matter how narrow or broad
This feature makes it possible for readers to find better explanations of any topic created by other writers. And it allows writers to create an explanation in a place that readers might actually find it.Figure 1. Screenshot of the "Derivative" topic page. View it live at: ourbigbook.com/go/topic/derivativeVideo 2. OurBigBook Web topics demo. Source. - local editing: you can store all your personal knowledge base content locally in a plaintext markup format that can be edited locally and published either:This way you can be sure that even if OurBigBook.com were to go down one day (which we have no plans to do as it is quite cheap to host!), your content will still be perfectly readable as a static site.
- to OurBigBook.com to get awesome multi-user features like topics and likes
- as HTML files to a static website, which you can host yourself for free on many external providers like GitHub Pages, and remain in full control
Figure 2. You can publish local OurBigBook lightweight markup files to either OurBigBook.com or as a static website.Figure 3. Visual Studio Code extension installation.Figure 5. . You can also edit articles on the Web editor without installing anything locally. Video 3. Edit locally and publish demo. Source. This shows editing OurBigBook Markup and publishing it using the Visual Studio Code extension. - Infinitely deep tables of contents:
All our software is open source and hosted at: github.com/ourbigbook/ourbigbook
Further documentation can be found at: docs.ourbigbook.com
Feel free to reach our to us for any help or suggestions: docs.ourbigbook.com/#contact








